MATLAB中用ode45求解微分方程,如何设置最大步长?
答案:2 悬赏:60
解决时间 2021-04-18 03:23
- 提问者网友:✐ۖ﹏ℳ๓北风
- 2021-04-17 10:17
MATLAB中用ode45求解微分方程,如何设置最大步长?
最佳答案
- 二级知识专家网友:而你却相形见绌
- 2021-04-17 11:23
[t,y]=ode45('odefun',[t0:dt:tend],y0);
这里的dt就是最大步长,但是因为ode45是变步长的算法,如果需要计算的时间范围 [t0 tend] 不是很大的话,那么返回的时间矩阵 t 里面的时间间隔不一定全是 dt,会包含更小的时间步。
这里的dt就是最大步长,但是因为ode45是变步长的算法,如果需要计算的时间范围 [t0 tend] 不是很大的话,那么返回的时间矩阵 t 里面的时间间隔不一定全是 dt,会包含更小的时间步。
全部回答
- 1楼网友:野性且迷人
- 2021-04-17 11:36
function [x,y]=runge_kutta1(ufunc,y0,h,a,b) %参数表顺序依次是微分方程组的函数名称,初始值向量,步长,时间起点,时间终点,
n=floor((b-a)/h); %求步数
x(1)=a; %时间起点
y(:,1)=y0; %赋初值,可以是向量,但是要注意维数
for ii=1:n
x(ii+1)=x(ii)+h;
k1=ufunc(x(ii),y(:,ii));
k2=ufunc(x(ii)+h/2,y(:,ii)+h*k1/2);
k3=ufunc(x(ii)+h/2,y(:,ii)+h*k2/2);
k4=ufunc(x(ii)+h,y(:,ii)+h*k3);
y(:,ii+1)=y(:,ii)+h*(k1+2*k2+2*k3+k4)/6; %按照龙格库塔方法进行数值求解
end
这是给你编的定步长龙哥库塔,即ode4
我要举报
如以上问答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯