|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
程序如下：; L4 |2 J3 q' U4 t
function out=mult_syg7(z,x)/ ]% j/ V$ p' J S4 B" I& \0 b
Pi=3.14159;7 ~/ z: ]2 `% A G. E, y1 b- H- U
r0= 1.24*10^-3; % Outer diameter of the membrane, [m]# m9 g4 [+ j& d
ri= 9.4*10^-4;% Inner diameter of the membrane, [m]! [; f/ e5 `( y# w/ \/ ^' ^4 b
r=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]( Y8 D* k1 S6 k3 c4 }& g
L=0.05;% membrane length,[m]
7 F; V: y4 ? k% V) DAc=Pi*r*L;% membrane area [m^2]
7 W9 [0 b; b' t, L' wR=8.314;
5 R6 l( p g4 }1 y R% ^T0=1000;
, v* x0 s( C1 b# r$ k/ rT=T0+273;* o5 O, t9 F4 ~$ b
P=101.325;%kpa
1 J C1 l5 }- n" e& q/ vF=2.678*10^-2*2;, o. b) \6 |8 V Y, }
m=3;
1 n$ S ^! |1 K' pn=1;
. w' c4 E5 G* C4 L+ e$ ~S=(1+m+n+2*x(1))/P;
/ `% m+ c* i/ t1 uPCH4=(1-x(1))/S;
1 c% K, q# i8 \2 y- A( |& X0 ]PH2O=(m-x(1)-x(2))/S;
4 l3 [" N( D9 \" |9 z* }out= ..., w6 r+ _/ L ~2 ^2 M6 t' O/ s6 Z
[Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08+ ]' L) q. z. a1 G! ~# M& u$ z
(3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];) L# k/ ^4 D: o
! J$ M J% m; x/ u W5 c: ]
调用如下:
: R$ ?) v$ V( s6 O>> x0 = [0; 0];0 @, {4 y4 N* J$ D' d: ]- m
tspan = [0:0.1:1];
v! B7 q1 V2 \# ^M = [1 0 ; 0 0];
7 G$ ^ |( N z/ u& O( |options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
: _- ]* `& R5 Z; Y; O6 |* N[t,x] = ode15s(@mult_syg7,tspan,x0,options)" f! X, u3 U7 U$ Y2 \8 X6 k! p
6 f9 V0 Q G; k8 U" [$ g
报错如下：8 P- _, u3 X- _ \/ E* y
错误使用 daeic12 (line 76)5 h' j4 |6 Z- a
此 DAE 的索引大于 1。
z/ {# d d9 I- o; _, ?' b; h4 ~6 m* f8 L, w
出错 ode15s (line 310)
4 H0 w* o* D' U- f" M6 b$ ? [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,.../ R- o$ E* T5 i% Z/ s
|
|
/1 
发表于 2021-3-19 14:48
收藏
淘帖
支持！
反对！