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

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

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

帮帮忙哈，各位大虾

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

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

d*E*A*(Poi-x)*(sqrt(8*R*T/pi/Mi))<BR><BR>这里好像有错
页: [1]
查看完整版本: [求助]那位大虾知道这段程序错在哪里了，能否帮忙看下，不胜感激