声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2053|回复: 5

[稳定性与分岔] 求助:方程中参数不确定的微分方程该怎么用龙格——库塔法解?

[复制链接]
发表于 2011-5-17 16:41 | 显示全部楼层 |阅读模式

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

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

x
dx=[x(2);fx/M+md*cos(tao+fai);x(4);fy/M+md*sin(tao+fai)-G]; % x=x(1) x'=x(2) y=x(3) y'=x(4)
fx=(sqrt((x(1)-2*x(4))^2+(x(1)+2*x(2))^2)/(1-x(1)^2-x(3)^2))*(3*x(1)*V-G1*sin(a)-2*S*cos(a))+(RC/L^2)*(8*pi/sqrt(x(1)^2+x(3)^2))*(1-1/sqrt(1-x(1)^2-x(3)^2))*cos(a);
fy=(sqrt((x(1)-2*x(4))^2+(x(1)+2*x(2))^2)/(1-x(1)^2-x(3)^2))*(3*x(3)*V-G1*cos(a)-2*S*sin(a))+(RC/L^2)*(8*pi/sqrt(x(1)^2+x(3)^2))*(1-1/sqrt(1-x(1)^2-x(3)^2))*sin(a);
V=(2+(X(3)*cos(a)-x(1)*sin(a))*G1)/(1-x(1)^2-x(3)^2);
S=(x(1)*cos(a)+x(3)*sin(a))/(1-(x(1)*cosa+x(3)*sin(a))^2);
G1=(2/sqrt(1-x(1)^2-x(3)^2))*(pai/2+atan((x(3)*cos(a)-x(1)*sin(a))/(1-x(1)^2-x(3)^2)));
a=atan((x(3)+2*x(2))/(x(1)-2*x(4)))-(pi/2)*sin((x(3)+2*x(2))/(x(1)-2*x(4)))-(pai/2)*sin(x(3)+2*x(2));
回复
分享到:

使用道具 举报

 楼主| 发表于 2011-5-17 16:44 | 显示全部楼层
如上,写完了才发现一些参数并不是确定的值,不能直接用ode45求解,请高手帮忙,这样的微分方程该怎么列?
 楼主| 发表于 2011-5-18 16:01 | 显示全部楼层
有没有高手可以帮下忙啊~~
发表于 2011-5-26 10:38 | 显示全部楼层
你的参数是在变化的吗?那就给个范围,然后在ODE45调用
发表于 2012-11-2 18:03 | 显示全部楼层
我也遇到了相同的问题,请问你解决了吗?
发表于 2012-11-8 16:02 | 显示全部楼层
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 19:23 , Processed in 0.059114 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表