24小时热门版块排行榜    

查看: 658  |  回复: 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的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
huangshuang1: 金币+10, ★★★很有帮助, 我试了,还是振荡的厉害,不收敛 2015-03-04 21:10:57
你震荡得厉害,把ode45改成ode15s
3楼2015-01-31 20:49:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

huangshuang1

木虫 (著名写手)

请问,有人帮我看看么??实在是找不到原因了
2楼2015-01-29 14:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 材料080500调剂求收留 +4 一颗meteor 2026-03-13 4/200 2026-03-19 10:32 by 30660438
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见