| 查看: 1082 | 回复: 3 | |||
[交流]
【求助】如何用matlab求解这个微分方程？【已解决】 已有3人参与
|
请大家帮忙啊，谢谢！ [ Last edited by kuhailangyu on 2010-6-22 at 16:04 ] |
回复此楼» 猜你喜欢
材料工程085601，270求调剂
已经有25人回复
-
297求调剂
已经有18人回复
-
材料工程调剂
已经有11人回复
-
288求调剂
已经有15人回复
-
291求调剂
已经有4人回复
-
283电子信息求调剂
已经有4人回复
-
085801 总分275 本科新能源 求调剂
已经有7人回复
-
0860004 求调剂 309分
已经有6人回复
-
化学工程调剂289
已经有33人回复
-
08工科求调剂290分
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求解描述振动的二阶微分方程
已经有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
★ ★ ★ ★
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
