24小时热门版块排行榜    

CyRhmU.jpeg
查看: 921  |  回复: 3

qhlu788

金虫 (正式写手)

[交流] 【求助】如何用matlab求解这个微分方程?【已解决】已有3人参与

请大家帮忙啊,谢谢!



[ Last edited by kuhailangyu on 2010-6-22 at 16:04 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gitan

木虫 (正式写手)

Visiting Scholar

★ ★
nono2009(金币+1):鼓励应助。 2010-06-22 07:46:44
qhlu788(金币+2): 2010-06-22 09:03:07
zzuwangshilei(金币+1):积极参与 2010-06-22 19:46:58
先降阶,把这个包含2个方程的二阶常微分方程组变成包含4个方程的一阶常微分方程组。接着看这个问题是初值问题(ode)还是边值问题(bvp)。
对于初值问题,可以采用ode45()编程求解;对于边值问题,则采用bvp4c()编程求解。
2楼2010-06-21 23:30:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

photons

至尊木虫 (职业作家)

★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qhlu788

金虫 (正式写手)

谢谢兄弟们帮忙,非常感谢!
4楼2010-06-22 09:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qhlu788 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见