使用 decic命令 错误分析

unfaithful2021

我的M文件：
y0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
+ r: d/ }4 W+ k( t; L: Lyp0est = zeros(10,1);
opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
: H1 a7 G; O3 Q3 z' _- K
$ P( E0 V+ i/ t% O
/ K: R/ Y  S4 x: G/ a" R5 H$ r  k运行后报错：
& I4 M- x. R& \" ?4 Z" E错误使用 decic>sls (line 170)
索引可能大于 1。
# h) l' e/ i! F, x* [4 d! f
出错 decic (line 77)
    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);

) C3 H" ^, z) F7 l出错 DAEsL (line 39)
4 _# M7 r7 A( r6 l1 R' v2 U[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
! P9 Q, y3 S5 F$ Q1 |( m  B% x
大家帮分析一下是什么问题？

I_believe

索引可能大于 1。
1 N% z, c. n+ e5 ky0est,yp0est  只能是1个 ？你的有10个。

unfaithful2021

I_believe 发表于 2022-9-15 13:58
$ s$ w% u- [, P2 L- u: B索引可能大于 1。
# N9 u3 J! K9 B$ ?  ^9 Py0est,yp0est  只能是1个 ？你的有10个。

y0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个
! a$ B# j& N- ~$ M% N% \" i

peerless2021

Decic调用ode15ipdinit()，它产生几个结果，包括ode函数的偏导数。然后将第二个偏导数传递给decic内部例程sls()，该例程计算其秩。如果rank小于固定变量的数量，那么sls()将发出您看到的错误消息。简单说就是您的自变量，超过等式的数量了。