| 查看: 1210 | 回复: 3 | |||
[交流]
【求助】如何用matlab求解这个微分方程？【已解决】 已有3人参与
|
请大家帮忙啊，谢谢！ [ Last edited by kuhailangyu on 2010-6-22 at 16:04 ] |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab 线性常微分方程求解
已经有5人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- 隐式欧拉法求解一阶常微分方程
已经有7人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】Python 中 求解微分方程 odeint 的用法
已经有11人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
gitan
木虫 (正式写手)
Visiting Scholar
- 应助: 8 (幼儿园)
- 金币: 2259.8
- 红花: 1
- 帖子: 369
- 在线: 195.1小时
- 虫号: 831616
- 注册: 2009-08-21
- 性别: GG
- 专业: 光学信息获取与处理
2楼2010-06-21 23:30:46
photons
至尊木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 28533.4
- 散金: 160
- 红花: 15
- 帖子: 3340
- 在线: 200.2小时
- 虫号: 437799
- 注册: 2007-09-07
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★ ★
nono2009(金币+2):鼓励应助。 2010-06-22 07:46:55
qhlu788(金币+8): 2010-06-22 09:03:16
zzuwangshilei(金币+2):帮助他人解决问题 2010-06-22 19:47:15
nono2009(金币+2):鼓励应助。 2010-06-22 07:46:55
qhlu788(金币+8): 2010-06-22 09:03:16
zzuwangshilei(金币+2):帮助他人解决问题 2010-06-22 19:47:15
|
function ydot=DyDt(t,y) ydot=[y(2);(y(4)^2*sin(y(3)-y(1))-y(2)^2*sin(y(3)-y(1))*cos(y(3)-y(1))... -2*sin(y(1))+cos(t)*2*cos(y(1))+cos(y(3)-y(1))*sin(y(3))-cos(t)*... cos(y(3))*cos(y(3)-y(1)))/(2-cos(y(3)-y(1))^2);y(4);(1/2*y(4)^2*... sin(y(3)-y(1))*cos(y(3)-y(1))-sin(y(1))*cos(y(3)-y(1))+cos(t)*cos(y(1))... *cos(y(3)-y(1))-y(2)^2*sin(y(3)-y(1))+sin(y(3))-cos(t)*cos(y(3)))... /(1/2*cos(y(3)-y(1))^2-1)]; clear all tspan=[0,100]; y0=[1;0;1;0]; [tt,yy]=ode45(@DyDt,tspan,y0); subplot(2,1,1) plot(tt,yy(:,1),tt,yy(:,3),'r:') legend('\theta_{1}','\theta_{2}') grid on ttt=tt(1:10); yyy=yy(1:10,1); yyyy=yy(1:10,3); subplot(2,1,2) plot(ttt,yyy,ttt,yyyy,'r:') grid on legend('\theta_{1}','\theta_{2}')  |
3楼2010-06-22 07:35:17
4楼2010-06-22 09:03:38
