Question

matlab解微分方程，在线等，谢谢

求matlab代码，谢谢

Answer 1

%主函数，调用四阶龙格库塔
y0=[2,1]';
tf=10;h=0.001;
y = rk4(@eom,h,y0,tf/h+1);
y(end,1)%输出t=10时，y1的结果
plot(t,y)

%微分方程
function ydot = eom(y)
ydot(1) = y(2);
ydot(2) = 2*(1-y(1)^2)*y(2)-y(1);

%龙格库塔
function x=rk4(f,t,h,x0,n)
x(:,1)=x0;
for i = 1:n-1
L1=feval(f,t(i),x(:,i));
L2=feval(f,t(i)+h/2,x(:,i)+(h/2)*L1);
L3=feval(f,t(i)+h/2,x(:,i)+(h/2)*L2);
L4=feval(f,t(i)+h,x(:,i)+h*L3);
x(:,i+1)=x(:,i)+(h/6)*(L1+2*L2+2*L3+L4);
end

分开写的好处是，以后可以用这三个文件计算任意微分方程，只要把主函数中初值、步长、仿真时间和微分方程文件中的模型变换即可。