使用 decic命令 错误分析

unfaithful2021

我的M文件：
6 ^' ?( i+ |: c: m$ L4 @y0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
yp0est = zeros(10,1);
+ ~- i8 t- Z5 N! }; g& i$ z7 |0 Oopt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
- \* s  W  g& q5 m! T$ N

. d5 r! @, A" ?' _  M* J运行后报错：
错误使用 decic>sls (line 170)
索引可能大于 1。

* o7 z1 _- i  t6 p  B1 l! H. q出错 decic (line 77)
% x/ i' M2 q9 U! c8 x# @; I- i6 A& I& c    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);
2 D6 E1 M& P/ J, X
! o2 \5 S- V7 k出错 DAEsL (line 39)
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)

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

I_believe

索引可能大于 1。
8 w4 `7 S& z4 M$ P4 sy0est,yp0est  只能是1个 ？你的有10个。

unfaithful2021

I_believe 发表于 2022-9-15 13:58
索引可能大于 1。
. p1 E! H, W% e/ I# h+ n# Iy0est,yp0est  只能是1个 ？你的有10个。

y0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个
: v: E8 z; _/ D5 i" r* c7 |

peerless2021

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