急向大侠们求助一个问题

煜宸0922 · 发表于 2011-5-9 11:22

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

向大侠们急求“变步长”龙格库塔法解下面方程组。
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
是碰撞系统的方程！急求，谢谢各位大侠，小弟初学，各种困难，望帮助，不胜感激。

meiyongyuandeze · 发表于 2011-5-9 12:07

回复 1 # 煜宸0922 的帖子

你直接用ode45解不就可以了吗！

煜宸0922 · 发表于 2011-5-9 14:40

回复 2 # meiyongyuandeze 的帖子

呵呵，大侠，能用ode做我就做啦，关键是这是要求呀

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 急向大侠们求助一个问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。