| 查看: 962 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】如何用matlab求解这个微分方程？【已解决】 已有3人参与
|
|||
请大家帮忙啊，谢谢！ [ Last edited by kuhailangyu on 2010-6-22 at 16:04 ] |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求解描述振动的二阶微分方程
已经有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人回复
★ ★ ★ ★
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