24小时热门版块排行榜    

查看: 653  |  回复: 2

huangshuang1

木虫 (著名写手)

[求助] Matlab ode45求解非线性方程 已有1人参与

function f=mmm(t,y)
f=zeros(4,1);
f(1)=y(2);
f(2)=(50*pi-y(2))/((50*pi-y(2))^2+595.8566)*1.858*10^6-1.52*10^5*(0.0625*y(1)-0.149*y(3))-1.43*10^7*(0.0625*y(1)-0.149*y(3))^2-2.86*10^6*(0.0625*y(1)-0.149*y(3))^3-32.36*(0.0625*y(2)-0.149*y(4));
f(3)=y(4);
f(4)=-86.994-2.437*10^4*(0.149*y(3)-0.0625*y(1))-2.29*10^6*(0.149*y(3)-0.0625*y(1))^2-4.58*10^5*(0.149*y(3)-0.0625*y(1))^3-5.185*(0.149*y(4)-0.0625*y(2));

y0=[0,50*pi,0,50*pi];
[t,y]=ode45(@mmm,[0:0.002:20],y0);
plot(t,y(:,2),t,y(:,4),':')


用ODE45一直处于busy状态,然后我换成ODE23S,可是结果振荡得很厉害,不收敛。
请问各位,这个该怎么解??
回复此楼

» 猜你喜欢

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

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

huangshuang1

木虫 (著名写手)

请问,有人帮我看看么??实在是找不到原因了
2楼2015-01-29 14:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
huangshuang1: 金币+10, ★★★很有帮助, 我试了,还是振荡的厉害,不收敛 2015-03-04 21:10:57
你震荡得厉害,把ode45改成ode15s
3楼2015-01-31 20:49:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huangshuang1 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程调剂 +3 大可digkids 2026-03-16 3/150 2026-03-16 09:09 by DDDddddmm
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
信息提示
请填处理意见