声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 995|回复: 0

[综合讨论] 常微分方程方程组的解为什么延长积分区间结果会溢出。。

[复制链接]
发表于 2010-10-11 17:46 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
函数
function dy=ODE1(z,y)
dy=zeros(3,1);
r=0.05;H1=206000;H2=165000;Cp=35.77;F=0.096;Fco=0.00672;Fco2=0.00192;pb=1175000;
dy(1)=(pi*r*r/Fco)*22656346.5*exp(-84816.7/8.314/y(3))*(70*(1-y(1))/(100-14*y(1)-4*y(2)))^0.4734;
dy(2)=(pi*r*r/Fco2)*1070043499*exp(-105802.4/8.314/y(3))*(20*(1-y(2))/(100-14*y(1)-4*y(2)))^0.000945;
dy(3)=(pi*r*r/F/Cp)*(22656346.5*exp(-84816.7/8.314/y(3))*(70*(1-y(1))/(100-14*y(1)-4*y(2)))^0.4734*H1+1070043499*exp(-105802.4/8.314/y(3))*(20*(1-y(2))/(100-14*y(1)-4*y(2)))^0.000945*H2);
当命令为
X0=[0;0;500]; tspan=[0 1.6]; [x,y]=ode45(@ODE1,tspan,X0)
还有结果,当命令为
X0=[0;0;500]; tspan=[0 1.7]; [x,y]=ode45(@ODE1,tspan,X0)
就没有结果,结果溢出了,谁给解释下。。。

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 23:23 , Processed in 0.056263 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表