24小时热门版块排行榜    

查看: 614  |  回复: 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的回帖
信息提示
请填处理意见