EDA365电子论坛网
标题:
微分代数方程求解问题，运行报错：‘索引大于1’
[打印本页]
作者:
勇往直前11
时间:
2021-3-19 14:48
标题:
微分代数方程求解问题，运行报错：‘索引大于1’
程序如下：
: g) w" x% E" G6 p
function out=mult_syg7(z,x)
$ q2 R4 R! X1 k
Pi=3.14159;
% [- u, f6 s* p8 C
r0= 1.24*10^-3; % Outer diameter of the membrane, [m]
y6 k& ^$ b% o! h; N4 Z
ri= 9.4*10^-4;% Inner diameter of the membrane, [m]
" C: v" L% ^0 F; t3 A( R
r=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
) l3 {: P/ \" y
L=0.05;% membrane length,[m]
# x1 Q- Y5 L+ w8 a N( [- k+ `
Ac=Pi*r*L;% membrane area [m^2]
- g+ e" D& n4 [% W/ P3 X3 `
R=8.314;
, O' _4 P/ e4 Q: x. q+ z
T0=1000;
7 y5 {) b5 ]' f* a+ Y( x$ ^$ ` y
T=T0+273;
$ z' Y" p* k& e. Y
P=101.325;%kpa
$ r7 @/ z/ F6 @% L4 o9 \- ], e
F=2.678*10^-2*2;
) A# [) \7 o( ]# H7 I6 k. M
m=3;
: q- x8 Z% f3 o+ d
n=1;
" j) X5 ^6 y" t4 u8 y
S=(1+m+n+2*x(1))/P;
8 G1 K7 w* j0 t. T
PCH4=(1-x(1))/S;
3 X. P. k, I5 x. |+ l$ M3 ]+ {
PH2O=(m-x(1)-x(2))/S;
) H0 [7 l1 w) [3 q
out= ...
+ t# i- [- X0 O. `
[Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08
- j% W% `- Z I% I$ B
(3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];
+ W/ o3 r9 W- X& ?9 @- B$ I
" |8 b# P% ?; t* i* i7 B
调用如下:
3 w3 q+ R& S% r. d, y
>> x0 = [0; 0];
$ Q8 R9 r/ L. w7 s$ C
tspan = [0:0.1:1];
6 A& ^2 G- K; @$ r6 m# R
M = [1 0 ; 0 0];
& m' t8 u0 q# n' ]# y% C
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
0 R7 y2 a& R/ Z1 k1 C3 S
[t,x] = ode15s(@mult_syg7,tspan,x0,options)
* i+ ?7 s J" X$ r" q o; M
- g8 b! e/ z: F* [5 F! Y
报错如下：
- c+ A6 ^' C# M N& u5 S1 z
错误使用 daeic12 (line 76)
6 i( G2 l, y1 Q* ~4 o2 Y! U; V
此 DAE 的索引大于 1。
% I* ]1 N3 b# H, W; S' p w
; K" R2 M' \; k# e7 U' V
出错 ode15s (line 310)
5 Q1 [0 G, R* `0 P8 m6 y% m
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
- Z) [9 n9 c/ V
作者:
cichishia
时间:
2021-3-19 15:13
帮你顶一下
作者:
kekek
时间:
2021-3-19 15:18
作者:
shuddkk
时间:
2021-3-22 13:27
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2
