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

勇往直前11

程序如下：
% F- K7 j1 I' Y& J% \function out=mult_syg7(z,x)
: A0 K+ ?6 y4 n, ?; U. KPi=3.14159;
+ b9 |! p0 {/ x0 ir0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
9 P9 ?$ {6 u% Z& T2 Pri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
) |1 r0 Y; N& w: gr=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
L=0.05;% membrane length,[m]
  q4 A) M, d. I6 l% d1 d- HAc=Pi*r*L;% membrane area [m^2]
R=8.314;
T0=1000;
T=T0+273;
" i3 o- F+ d+ ~0 v5 S% I: I: S8 B0 rP=101.325;%kpa
F=2.678*10^-2*2;
m=3;
/ S% b3 v" Y0 ]+ O" Q" Gn=1;
S=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
+ ?$ \6 X6 A0 g1 {: f) qPH2O=(m-x(1)-x(2))/S;
3 l4 D% c1 }" n* b" eout= ...
" Z6 ?, x6 e- C, B    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
4 l: B2 t8 B  |. N7 M) x+ o    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];

8 F7 O0 h8 R& d1 Y2 n9 y调用如下:
" j5 ^9 Y$ ~# G' y>> x0 = [0; 0];
tspan = [0:0.1:1];
M = [1 0 ; 0 0];
7 _' T- x4 [9 o; q( Y$ k% Ioptions = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
" f$ P4 W' S, m/ _) l. _: i[t,x] = ode15s(@mult_syg7,tspan,x0,options)

% l' j0 s; K3 {" |+ W报错如下：
0 W) }3 d) V% t8 C$ o8 k4 G错误使用 daeic12 (line 76)
7 {! X# R6 h& W9 B$ [, {( V% E此 DAE 的索引大于 1。

  Y/ @( b- ^' ^. X4 p出错 ode15s (line 310)
    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...

cichishia

帮你顶一下

kekek

shuddkk