|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
程序如下：
9 c- a1 }1 m' _8 S2 {/ a7 L9 `function out=mult_syg7(z,x)
% G- Y" R$ a- {# H/ Q* OPi=3.14159;
: R; I H7 n* Mr0= 1.24*10^-3; % Outer diameter of the membrane, [m]
! o: v1 e* D( P( a7 _9 Jri= 9.4*10^-4;% Inner diameter of the membrane, [m]
1 M/ b' G$ o+ {9 cr=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
R" u, t* U: R+ X+ HL=0.05;% membrane length,[m]
! |8 T, {4 m- m3 w$ F9 y4 `Ac=Pi*r*L;% membrane area [m^2]
; k( J7 }3 O* N9 y- g6 ]9 i) FR=8.314;% B# ~! F$ ]4 y5 `+ J0 ~
T0=1000;: ^; r* z4 K$ \
T=T0+273;
O% f# n/ I* R3 gP=101.325;%kpa
0 x$ a o* W# E0 A( A/ P; Q4 pF=2.678*10^-2*2;# P& o2 ~5 }, p0 H U6 ?
m=3;
- ~, `! M! @* b% {" h0 Xn=1;
4 T% ? b! Z s% {S=(1+m+n+2*x(1))/P;; C2 n3 y% D7 Q# [9 q; x
PCH4=(1-x(1))/S;
4 n2 e6 u: I7 o7 S0 b+ L$ z0 MPH2O=(m-x(1)-x(2))/S;+ D" C6 v$ J* R4 Y( K
out= ...
: |7 C9 p! z' E7 J {) d [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08- G. {# f4 u1 F; Y; \$ w6 ~6 h
(3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];3 r; [ H# u" Y& i
Q' u7 _% ^+ T8 Q1 j( l" k调用如下:; M1 e; Y; o# H8 g
>> x0 = [0; 0];
% P3 N" X; c6 p/ Jtspan = [0:0.1:1];( J+ ~8 R- K% @9 V* c. b$ e
M = [1 0 ; 0 0];
9 z( \; T1 v# K3 y1 \1 o1 ~options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
8 U! \3 y) D0 p( T# Z3 | [[t,x] = ode15s(@mult_syg7,tspan,x0,options): g7 U3 _8 |0 N. U, @' `" \
6 G h7 H8 v$ L$ Q报错如下：
& l7 E! f" a% {( r0 i错误使用 daeic12 (line 76)! t t! J! B$ Y$ C
此 DAE 的索引大于 1。
A7 F0 W4 S! e$ A, o
7 c M! g. D, m. n1 H& i出错 ode15s (line 310)
1 e- {& ~3 h0 f4 q" k% m7 O" Q [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,..." h8 @/ ?1 }5 e2 b8 a
|
|
/1 
发表于 2021-3-19 14:48
收藏
淘帖
支持！
反对！