微分代数方程求解问题，运行报错：‘索引大于1’

勇往直前11

程序如下：
function out=mult_syg7(z,x)
8 N; z+ M7 c9 G1 Y* \6 mPi=3.14159;
r0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
ri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
/ H) g0 G: ~6 E7 I& Y/ L6 Fr=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
L=0.05;% membrane length,[m]
8 O5 f9 k. r" S0 ?; z2 XAc=Pi*r*L;% membrane area [m^2]
R=8.314;
1 }* ^& D: U  W9 w, w* D6 [T0=1000;
T=T0+273;
2 }  D  `) V0 x" h8 s& `P=101.325;%kpa
- N2 n8 z/ ?% r! S0 PF=2.678*10^-2*2;
m=3;
n=1;
0 s/ J8 J/ B& L/ Y: ^+ d0 B7 mS=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
PH2O=(m-x(1)-x(2))/S;
out= ...
    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];
3 x! u2 _, r* [1 w
% o- w& P! C5 O4 }调用如下:
$ u' e' Y$ W" t6 ?>> x0 = [0; 0];
5 m& f  p& k# D' `& Etspan = [0:0.1:1];
M = [1 0 ; 0 0];
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
[t,x] = ode15s(@mult_syg7,tspan,x0,options)

, Q. M7 x+ Z' V; h报错如下：
0 `, ?+ d7 I! i7 e8 o错误使用 daeic12 (line 76)
; o1 O) [& q+ j' D" F+ `此 DAE 的索引大于 1。

出错 ode15s (line 310)
    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
, _1 T! M4 R8 E

cichishia

帮你顶一下

kekek

shuddkk