|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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)
就没有结果,结果溢出了,谁给解释下。。。 |
|