24小时热门版块排行榜    

查看: 755  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-16 23:28 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-16 22:40 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-16 21:39 by MHf40yA0Uhyb
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-16 21:26 by MHf40yA0Uhyb
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-16 21:26 by MHf40yA0Uhyb
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 16:12 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 14:27 by i3KBRD19ERkG
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:24 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:27 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:35 by DgNGHc3h5tPl
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
信息提示
请填处理意见