| 查看: 1259 | 回复: 3 | |||
[交流]
【求助】如何用matlab求解这个微分方程？【已解决】 已有3人参与
|
请大家帮忙啊，谢谢！ [ Last edited by kuhailangyu on 2010-6-22 at 16:04 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求解描述振动的二阶微分方程
已经有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 (幼儿园)
- 金币: 29359.8
- 散金: 160
- 红花: 15
- 帖子: 3406
- 在线: 200.8小时
- 虫号: 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
200