使用 decic命令 错误分析

unfaithful2021

我的M文件：
3 q9 u8 T: U" q+ fy0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
yp0est = zeros(10,1);
/ i: A+ _+ M( [+ t+ ~* ~% P8 ^1 Yopt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
1 w1 S) F: o# l5 v[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
* V4 p6 v0 G+ [  r! U' ]0 Z
" t" ~4 @7 `% Y
运行后报错：
错误使用 decic>sls (line 170)
索引可能大于 1。

出错 decic (line 77)
; G- J2 e( p! a8 f    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);

出错 DAEsL (line 39)
+ Q: ]  J1 O% u0 ?* H$ G8 S[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
+ {! r  S" b/ g0 X7 K) }0 n' i! t6 O
3 Q( ], [" P& `8 V1 j. W! r大家帮分析一下是什么问题？
: l8 N3 H0 w& ?, G$ G

I_believe

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

unfaithful2021

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

2 E& ]7 l; W: u3 ?/ w0 q; d9 r9 X! y8 |y0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个

peerless2021

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