使用 decic命令 错误分析

unfaithful2021

我的M文件：
y0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
yp0est = zeros(10,1);
  }1 E& ^( x6 ~; xopt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
. C6 [% x. l4 N* |[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)


5 _, l: p5 b! O( J- S  V运行后报错：
错误使用 decic>sls (line 170)
% s! w2 E* c2 v; S# {索引可能大于 1。

, \/ T& }6 k9 ~# {. n4 {4 C! ^出错 decic (line 77)
    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);

出错 DAEsL (line 39)
5 H  ]% _6 H1 B[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)

大家帮分析一下是什么问题？

I_believe

索引可能大于 1。
) M2 O( h8 w& [: p4 qy0est,yp0est  只能是1个 ？你的有10个。

unfaithful2021

I_believe 发表于 2022-9-15 13:58
6 R/ r* f7 p7 ?' x索引可能大于 1。
0 S2 f, P; O8 {. ny0est,yp0est  只能是1个 ？你的有10个。

7 X1 E6 k. y8 Z: ?7 \. @. Py0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个

peerless2021

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