急向大侠们求助一个问题
向大侠们急求“变步长”龙格库塔法解下面方程组。function dxdt=fun(t,x)
k1=2;k2=3;k3=0.05;k4=20;k5=1;k6=1;k7=3;
dxdt=zeros(4,1);
dxdt(1)=x(3);
dxdt(2)=x(4);
dxdt(3)=(1-k1)*sin(k2*t+k4)-2*k3*x(3)+2*k3*x(4)-x(1)+x(2);
dxdt(4)=(k1*sin(k2*t+k4)+2*k3*x(3)-2*k3*(1+k5)*x(4)+x(1)-(1+k6)*x(2))/k7;
end
是碰撞系统的方程！急求，谢谢各位大侠，小弟初学，各种困难，望帮助，不胜感激。 回复 1 # 煜宸0922 的帖子
你直接用ode45解不就可以了吗！
回复 2 # meiyongyuandeze 的帖子
呵呵，大侠，能用ode做我就做啦，关键是这是要求呀
页:
[1]