使用 decic命令 错误分析

unfaithful2021

我的M文件：
0 |- f, |& g3 \. i# a0 L6 P! E: u% Wy0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
& `3 L1 x& _" i/ i/ |9 L" Z2 ^yp0est = zeros(10,1);
opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)

- h4 E8 x* P" G6 D$ l- ^
运行后报错：
) H& ]6 s1 D# \) I8 {4 B错误使用 decic>sls (line 170)
4 L# @. }& z( y, }' R0 @索引可能大于 1。

出错 decic (line 77)
* }8 q: }  J* x6 v6 p    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);

' \* a2 {3 U( r/ p6 C出错 DAEsL (line 39)
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
2 J( c9 f$ k- h& i( ^8 z) F) f
大家帮分析一下是什么问题？

peerless2021

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

unfaithful2021

I_believe 发表于 2022-9-15 13:58
索引可能大于 1。
4 v- F7 q, Z1 S- I1 j( sy0est,yp0est  只能是1个 ？你的有10个。

7 ?, l3 [2 p1 s3 F- [2 Uy0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个

I_believe

索引可能大于 1。
y0est,yp0est  只能是1个 ？你的有10个。