| 查看: 223 | 回复: 1 | |||
| 当前主题已经存档。 | |||
fanbao_keke铁虫 (正式写手)
|
[交流]
【求助】求高人指导 Matlab 程序 已有1人参与
|
||
|
您好：我现在遇到一个二阶方程，想用龙格-库塔法求解，但是在第一个微分方程中 的右边也存在时间的微分项，不知如何处理，求高人指导一下，万分感谢！ y=[y(1);y(2)]; y0=[1,0]; [dy(1)/dt= - 4*sin(3*t)* sqrt(1-y(1)^2) * sin(y(2)) - 0.1*dy(2)/dt; dy(2)/dt = 2*y(1)+4*sin(3*t)*y(1)*cos(y(2)) / sqrt(1-y(1)^2) + 0.5*sin(3*t)]; 求y(1)随时间的演化。 [ Last edited by fanbao_keke on 2010-4-12 at 09:03 ] |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★
fanbao_keke(金币+10):非常谢谢！ 2010-04-12 10:54
adu886886(金币+2):谢谢交流 2010-04-13 08:30
fanbao_keke(金币+10):非常谢谢！ 2010-04-12 10:54
adu886886(金币+2):谢谢交流 2010-04-13 08:30
|
function test y0=[1-eps,0]; tspan=0:0.01:2; [t,y]=ode45(@ODEfun,tspan,y0); plot(t,y(:,1)) function dydt=ODEfun(t,y) dy2dt = 2*y(1)+4*sin(3*t).*y(1).*cos(y(2))./ sqrt(1-y(1).^2) + 0.5*sin(3*t); dy1dt= - 4*sin(3*t).* sqrt(1-y(1).^2) .* sin(y(2)) - 0.1*dy2dt; dydt=[dy1dt;dy2dt]; |
2楼2010-04-12 10:24:58