| 查看: 273 | 回复: 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 ] |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
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
20