24小时热门版块排行榜    

查看: 1058  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qhlu788

金虫 (正式写手)

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

请大家帮忙啊,谢谢!



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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 4 个回答

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的回帖

qhlu788

金虫 (正式写手)

谢谢兄弟们帮忙,非常感谢!
4楼2010-06-22 09:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 085600材料与化工 求调剂 +14 enenenhui 2026-03-13 15/750 2026-03-19 21:53 by allen-yin
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 一志愿南昌大学,327分,材料与化工085600 +3 Ncdx123456 2026-03-19 3/150 2026-03-19 13:18 by houyaoxu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
请填处理意见