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

勇往直前11

程序如下：
' B% t6 I4 e$ N) a/ Ofunction out=mult_syg7(z,x)
7 ]8 W% ?/ H' Z9 Y6 b' `Pi=3.14159;
" }! e1 `' \/ A3 ?9 \! t7 qr0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
ri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
+ g' Z& D* N4 [2 L) B) Ar=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
L=0.05;% membrane length,[m]
Ac=Pi*r*L;% membrane area [m^2]
R=8.314;
T0=1000;
T=T0+273;
* D; G6 x9 ~9 f$ i, iP=101.325;%kpa
8 H/ h: F+ E& M& X. vF=2.678*10^-2*2;
m=3;
. G  k& ~! C& {: c0 p8 Kn=1;
S=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
3 t- D7 k- ?* s3 h" xPH2O=(m-x(1)-x(2))/S;
4 Z8 d* o2 f: N2 k* y6 Q1 q0 Sout= ...
    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
" T; h7 D: |' R% [% A* y* S    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];
) o, Q: L' S4 t1 o4 _1 |
调用如下:
>> x0 = [0; 0];
tspan = [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)

报错如下：
错误使用 daeic12 (line 76)
此 DAE 的索引大于 1。
' V3 R7 L: k* U5 r3 w2 Z. S  v' G
出错 ode15s (line 310)
- r  y" }; I7 M" X6 q- T    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
! O( ^4 A( h; \5 W+ |

cichishia

帮你顶一下

kekek

shuddkk