ode45、ode15s、ode15i怎么选？

%主程序-求解器：ode45
clear all
close all
clc
y0 = 0;
tspan = [0:0.01:10];
[t,y] = ode45(@subfun,tspan,y0);
figure(1)
plot(t,y);% 画出每次积分结果与时间的图形

%主程序-求解器：ode15s
clear all
close all
clc
y0 = 0;
yp0 = 1;
tspan = [0:0.01:10];
options = odeset('RelTol',1e-12);
[t,y] = ode15s(@subfun,tspan,y0,yp0,options);
figure(1)
plot(t,y);% 画出每次积分结果与时间的图形

%主程序-求解器：ode15i
clear all
close all
clc
y0 = 0;
yp0 = 1;
tspan = [0:0.01:10];
options = odeset('RelTol',1e-12);
[t,y] = ode15i(@subfun,tspan,y0,yp0,options);
figure(1)
plot(t,y);% 画出每次积分结果与时间的图形

%以“subfun”命名的子函数
function f = ode(t,y,yp)
f = cos(t);

#ode45求解结果

#ode15s求解结果

#ode15i求解结果:积分公差要求无法满足：我暂时还没找到是什么原因产生的

最后

以上就是高高溪流为你收集整理的ode45、ode15s、ode15i怎么选？的全部内容，希望文章能够帮你解决ode45、ode15s、ode15i怎么选？所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。