[求助]那位大虾知道这段程序错在哪里了，能否帮忙看下，不胜感激

chenhua835 · 发表于 2006-5-23 10:24

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

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

x

function tODE clear all;clc,format long global Pmi Poi L E A d l R T Mi v Pmi=2.3e-014;%定义包装材料氧气的渗透系数 Poi=2.1232e+004;%定义包装袋外氧气的压强分压 A=176.7146;%定义包装材料的总面积 L=0.003;%定义薄膜的厚度 l=0.007;%定义孔的长度 T=293;%定义实验进行的温度 v=2.6507e+003;%定义实验用的包装自由体积 Mi=32; R=8.31e+006;%定义气体常数 D=30e-004;%定义孔的直径 n=1;%定义包装材料的孔个数 E=D^2*pi*n/(4*A);%定义包装材料的孔隙率 y0=1.01e+4; [t1,y1]=ode45(@fun,[0,20],y0); [t2,y2]=ode23(@fun,[0,20],y0); plot(t1,y1,'k-',t2,y2,'b--'),xlabel('t'),ylabel('Pin') %--------------------------------- function dydt=fun(t,y) global Pmi Poi L E A d l R T Mi v % simulataneous equations dydt=(Pmi*(Poi-y)*(1-E)*A/L+d*E*A*(Poi-y)*((8*R*T/pi/Mi)^0.5))*R*T/v

chenhua835 · 发表于 2006-5-23 11:29

帮帮忙哈，各位大虾

chenhua835 · 发表于 2006-5-23 16:23

[em01][em08]

flybaly · 发表于 2006-5-23 17:01

d*E*A*(Poi-x)*(sqrt(8*R*T/pi/Mi)) 这里好像有错

账号		自动登录	找回密码
密码			我要加入

[求助]那位大虾知道这段程序错在哪里了，能否帮忙看下，不胜感激

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

回复：（chenhua835）[求助]那位大虾知道这段程序错在...