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

勇往直前11

程序如下：
# U$ U: L9 ~* w) L( Y9 f* D% bfunction out=mult_syg7(z,x)
Pi=3.14159;
r0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
ri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
3 s9 U# F! V8 c+ d. j. Y7 er=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
L=0.05;% membrane length,[m]
1 k( ?" N" X9 ], BAc=Pi*r*L;% membrane area [m^2]
R=8.314;
T0=1000;
; ~% P' T! T, aT=T0+273;
P=101.325;%kpa
F=2.678*10^-2*2;
: p! z% s" B8 ]( |8 \m=3;
3 w7 \0 a$ K+ k6 e% p4 tn=1;
6 O0 M2 x+ E* C$ ^) z% n$ nS=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
' x3 h9 B8 B) g# x: ?1 q, PPH2O=(m-x(1)-x(2))/S;
out= ...
  C5 O' r; l) N" \( }7 |    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
. ^9 R# |# f8 R+ c7 _) S    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];

( S- O5 O2 S' B调用如下:
' Q" I3 B6 }$ I! S0 a. L* p>> x0 = [0; 0];
- s. ?" U' U8 s& Ctspan = [0:0.1:1];
8 Z0 P/ Q+ Z& H: TM = [1 0 ; 0 0];
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
[t,x] = ode15s(@mult_syg7,tspan,x0,options)
! T% X8 d3 [) I4 b
报错如下：
4 N+ V  b* y" g& E  R错误使用 daeic12 (line 76)
, F# m9 w4 n# t此 DAE 的索引大于 1。

出错 ode15s (line 310)
0 Q; L+ g- c! h2 [5 f- I    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
; U% R8 ]& f/ f7 Z

cichishia

帮你顶一下

kekek

shuddkk