|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
程序如下：
# ]7 k2 q9 }& Yfunction out=mult_syg7(z,x)
- }! O! C: E+ L7 HPi=3.14159;
' o. R( C8 v7 ~r0= 1.24*10^-3; % Outer diameter of the membrane, [m], q; E+ P0 o( \* C
ri= 9.4*10^-4;% Inner diameter of the membrane, [m]
+ ?5 n) ^$ Z( c- \: z3 Ur=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
2 D5 K! K! z$ ^4 m, g2 NL=0.05;% membrane length,[m]
4 f7 k7 h% w4 ^ j2 y0 jAc=Pi*r*L;% membrane area [m^2]: b" {6 {" y: T' |
R=8.314;
7 ]5 x( a7 B! V+ a3 ^7 vT0=1000;
/ j' }" b+ O" R; [T=T0+273;1 U: U1 d3 M! m3 l5 C7 U1 |
P=101.325;%kpa8 m1 U/ R3 h6 T+ P) F3 R
F=2.678*10^-2*2;- Q# B- J- N( P- f) u7 P4 _
m=3;
2 ?( u1 o7 m* r# ~* K& Rn=1;. Y6 H( P6 K; W: t# |
S=(1+m+n+2*x(1))/P;8 O( i ^! ~1 W8 `0 i+ T# ^$ u
PCH4=(1-x(1))/S;
( P1 }1 h! o H) X3 SPH2O=(m-x(1)-x(2))/S;( ~" [% }5 u' [6 y( K! x# p/ v- e
out= ...
8 J: o8 R; [0 t: t4 U% t8 f [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
b7 s: A/ |4 d8 K) k, d (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];2 f4 @3 \& z9 N; l6 l3 I6 y7 w+ v
1 j9 L0 z3 p( F9 f4 _. |. y调用如下:
% | I4 g+ I& u8 i3 X5 w% G, @>> x0 = [0; 0];! q" D9 O& a5 {' a/ s
tspan = [0:0.1:1];. A# F* g5 k1 j: x
M = [1 0 ; 0 0];
& V! Y1 x- c& h l* I- X8 `! ioptions = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
: ~" B$ j* l5 n1 v7 j7 y3 g) \[t,x] = ode15s(@mult_syg7,tspan,x0,options)
8 z$ N2 |) S8 s- d6 S% ^& \- S
$ [- b0 `9 n. ?1 [7 S3 o0 l+ H报错如下：) o% r: J' P/ j8 w2 i4 @' @
错误使用 daeic12 (line 76)% {% s3 \: [+ f3 J. n$ Z( l
此 DAE 的索引大于 1。- E. S% g! U9 J1 k6 y( ]- O
4 N# Q6 j9 s4 |' g9 c% \出错 ode15s (line 310)% k2 u8 S) g# m* M
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...* y- Z- c) B( K
|
|
/1 
发表于 2021-3-19 14:48
收藏
淘帖
支持！
反对！