作业帮matlab中的ode45怎样设置定步长解微分方程?
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/23 20:18:26
matlab中的ode45怎样设置定步长解微分方程?
matlab中的ode45是变步长解法,定步长在simulink里有ode4和ode5,但是在matlab里面好像没有现成的函数,听说ode45可以实现定步长的输出,如何在matlab中（不是simulink）实现定步长解微分方程,
matlab中的ode45是变步长解法,定步长在simulink里有ode4和ode5,但是在matlab里面好像没有现成的函数,听说ode45可以实现定步长的输出,如何在matlab中（不是simulink）实现定步长解微分方程,
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
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