找回密码
 注册
关于网站域名变更的通知
查看: 1334|回复: 3
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-3-19 14:48 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
程序如下:3 r. T% R% E7 p. E' H
function out=mult_syg7(z,x)$ _- d/ Q4 [! X$ ]8 i
Pi=3.14159;) `( m  K1 A3 B" |
r0= 1.24*10^-3; % Outer diameter  of the membrane, [m]
$ S3 H6 q- g, }3 j% _! iri= 9.4*10^-4;% Inner diameter  of the membrane, [m]3 |8 X9 ~5 O$ j& `5 K
r=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]7 s3 L( E. e7 }8 j
L=0.05;% membrane length,[m]- p/ }8 G( ^* L9 b  q7 D' A
Ac=Pi*r*L;% membrane area [m^2]7 ], |4 e7 p  w; ^  S5 I; ]5 e
R=8.314;
3 v- s! F9 U) F& B" z0 D/ F; x1 zT0=1000;
  V. G  ~- C5 `, q. uT=T0+273;
8 G/ m# n; a; J$ h/ CP=101.325;%kpa
. e+ o/ L. S0 F7 ]: \F=2.678*10^-2*2;0 s  s3 P: Y) N* A
m=3;
* W+ x$ F* C/ C2 |8 k% On=1;
/ `/ T% A9 S- ]& C8 ]8 Y5 |2 iS=(1+m+n+2*x(1))/P;
7 D; ^( Y/ U8 {2 yPCH4=(1-x(1))/S;
, w. a8 Z/ ~, y" ~7 KPH2O=(m-x(1)-x(2))/S;
# D$ ^7 V9 }0 ]9 h4 Yout= ...* O4 B& N+ y  h: ~, F7 L; X, a  H
    [Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08# Y# u! I! T/ v$ X
    (3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];) [: p& {9 x8 x/ P- b. S+ o
" k" n9 L) H5 P
调用如下:' w: y& e5 @' x# I7 m3 U, f" ]
>> x0 = [0; 0];
; }2 v) M9 i6 k0 i0 G  Rtspan = [0:0.1:1];
. e" K& G$ ^$ ^: s. yM = [1 0 ; 0 0];( E* E- o" R: c# Z4 D
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
- y' d( N' {! u/ e3 o3 D: F7 {[t,x] = ode15s(@mult_syg7,tspan,x0,options)3 Y  X" W% H* u! t! Z% w1 a

; X5 |( {6 j8 t报错如下:
3 {/ g. y6 K& k/ p+ }7 g错误使用 daeic12 (line 76)6 k2 M0 o- _& G$ V) a
此 DAE 的索引大于 1。
# V7 c- @/ _, C0 \- s1 q) P- z2 o9 p
出错 ode15s (line 310)0 J) k8 E, ^( b. a" ~2 ?# ?
    [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
& e5 b& H. `: E) y: P

该用户从未签到

2#
发表于 2021-3-19 15:13 | 只看该作者
帮你顶一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-31 07:17 , Processed in 0.125000 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表