| 查看: 233 | 回复: 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 ] |
回复此楼» 猜你喜欢
今年春晚有几个节目很不错，点赞！
已经有10人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有6人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有5人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
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