lhy 发表于 2011-5-17 16:41

求助：方程中参数不确定的微分方程该怎么用龙格——库塔法解？

dx=; % 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));

lhy 发表于 2011-5-17 16:44

如上，写完了才发现一些参数并不是确定的值，不能直接用ode45求解，请高手帮忙，这样的微分方程该怎么列？

lhy 发表于 2011-5-18 16:01

有没有高手可以帮下忙啊~~

gyq274706322 发表于 2011-5-26 10:38

你的参数是在变化的吗？那就给个范围，然后在ODE45调用

机械毛驴 发表于 2012-11-2 18:03

我也遇到了相同的问题，请问你解决了吗？

hsfy919 发表于 2012-11-8 16:02

机械毛驴 发表于 2012-11-2 18:03 static/image/common/back.gif
我也遇到了相同的问题，请问你解决了吗？

对于具有随机参数问题的求解可以试试随机有限元的思路
页: [1]
查看完整版本: 求助：方程中参数不确定的微分方程该怎么用龙格——库塔法解？