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

勇往直前11

程序如下：
& S: w- }( M! o" G& A+ Dfunction out=mult_syg7(z,x)
+ d) ~1 z2 Y! N% r2 dPi=3.14159;
1 B0 n: Y6 X. Fr0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
9 k) d  d$ ~& u* `- C9 l+ xri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
r=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
- \* ]* f  o$ _L=0.05;% membrane length,[m]
Ac=Pi*r*L;% membrane area [m^2]
! @( Z5 o& N8 H/ qR=8.314;
T0=1000;
1 r: J& v8 w0 ZT=T0+273;
! W$ x2 Q7 Q5 }/ E: b* uP=101.325;%kpa
F=2.678*10^-2*2;
6 O) m0 J8 Z& t: u3 w' q- }5 Zm=3;
n=1;
S=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
+ r- A7 [8 B) O% w7 \2 [7 LPH2O=(m-x(1)-x(2))/S;
: @) F% e; q) ~, S  i0 [. lout= ...
# }( P4 Y# q/ t* U3 |5 e- J    [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)];

调用如下:
>> x0 = [0; 0];
tspan = [0:0.1:1];
+ M7 G0 }7 @5 UM = [1 0 ; 0 0];
: S" `) A+ X$ w0 \/ }options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
- N1 `8 Q- W" |+ l! k- r1 k[t,x] = ode15s(@mult_syg7,tspan,x0,options)
% C; L! `: J; c% S& M' M
报错如下：
( p! H  i. ^. H3 Z- A: Z错误使用 daeic12 (line 76)
% V" p7 o# _1 s* s1 U4 D此 DAE 的索引大于 1。

' }5 V$ N; t2 K' {. V5 {出错 ode15s (line 310)
  u( F$ _( h+ e" \. L2 J& X    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
$ {$ n% i* V4 I5 b1 q: x" A

shuddkk

kekek

cichishia

帮你顶一下