使用 decic命令 错误分析

unfaithful2021

我的M文件：
y0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
yp0est = zeros(10,1);
opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
8 v# ]& q& X' w) ^. p" Q% A

, L  I2 w5 X3 K$ A+ j1 E4 N& G运行后报错：
" L" u6 u: q* L. J& c! N) @错误使用 decic>sls (line 170)
索引可能大于 1。

出错 decic (line 77)
    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);
! G4 u1 N. |  W1 _$ P
* m( H* R/ }! A+ _3 R/ m出错 DAEsL (line 39)
$ P/ {, i0 j1 X[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)

- P3 C2 c/ N; a. D( H大家帮分析一下是什么问题？
2 j1 c. \! c. \, U

I_believe

索引可能大于 1。
+ ~+ H/ e9 Q2 ]1 C3 f0 Ay0est,yp0est  只能是1个 ？你的有10个。

unfaithful2021

I_believe 发表于 2022-9-15 13:58
索引可能大于 1。
y0est,yp0est  只能是1个 ？你的有10个。

y0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个
) [* }; {5 h: {/ @

peerless2021

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