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

勇往直前11

程序如下：
1 \3 I9 Z- O. Y& _function out=mult_syg7(z,x)
Pi=3.14159;
% M/ n- m4 b, C# vr0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
% b  v5 x+ X" s6 r. c: Jri= 9.4*10^-4;% Inner diameter  of the membrane, [m]
r=(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;
P=101.325;%kpa
F=2.678*10^-2*2;
+ l; g) v/ ?9 v$ K& J  s$ ?5 Fm=3;
& T% N' M2 U+ z6 z+ ?n=1;
S=(1+m+n+2*x(1))/P;
PCH4=(1-x(1))/S;
PH2O=(m-x(1)-x(2))/S;
out= ...
. \* d" L; X% M6 B# V$ }    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
7 n- M- c" q. s2 A    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];
* n% U5 M9 `. ^, C
+ q$ u3 R# y! U$ d! P7 E' S调用如下:
>> x0 = [0; 0];
tspan = [0:0.1:1];
M = [1 0 ; 0 0];
5 L# x) c; `3 [' E; |. Loptions = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
[t,x] = ode15s(@mult_syg7,tspan,x0,options)
# P: x6 H% h# P6 s2 q! R( o
报错如下：
错误使用 daeic12 (line 76)
1 r5 Y2 {1 @( G' G此 DAE 的索引大于 1。

7 L$ W3 T$ k* {: p5 E/ {出错 ode15s (line 310)
    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
. @* z9 E* l+ v5 D/ s2 q/ M8 |7 x+ g6 _

cichishia

帮你顶一下

kekek

shuddkk