积分工具CVODE 使用简介 sundials工具使用说明书

Step 1：网上下载sundialsTB（goole 搜索）包含在在文件sundials-2.6.1.tar中

Step2: 根据说明书进行安装

◆running the matlab script install STB.m which is present in the sundialsTB top

directory

◆add the sundialsTB startup to the system-wide startup file (by linking or copying):

◆说明：编译成mex文件可能会把报错compiler 不存在，最好的解决办法是安装visual

studio 2010版本，2013版本过高

Step 3: 算例测试

与ode15有些微的区别

◆微分函数返回变量，CVODE有三个变量：[y, flag, new_data]

◆对于传入参数的处理，CVODE主要是利用options中userdata进行传递：

主要代码如下：

参数传入：

data.f_fcn = f_fcn;

data.gh_fcn = gh_fcn;

t0 = 0.0;

tf = 13.0;

y0 = x0;

options = CVodeSetOptions('UserData',data,...

'RelTol',1.e-8,...

'AbsTol',1e-6);

子函数，参数提取：

function [y, flag, new_data] = QQGGSFunc_CVODE(t, xk,

data)

% function of ode -QGS or CDF

f_fcn = data.f_fcn ;

gh_fcn=data.gh_fcn;

…

其余的就参考附件的CVODE使用说明书