使用 decic命令 错误分析

unfaithful2021

我的M文件：
y0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
6 L( {% x9 f& `5 I1 e' zyp0est = zeros(10,1);
opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));
$ r, E# f+ D3 _& b# X[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
2 U! v# Y! x. H3 [
( y. [5 K; A9 F2 Y3 g
运行后报错：
, U6 w9 P5 Q0 u" d2 ?8 L错误使用 decic>sls (line 170)
索引可能大于 1。

5 C8 H' P+ U& ^出错 decic (line 77)
    [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);
8 Q: P* g2 n( h9 I
# q  |7 _( H2 i; o0 `" s出错 DAEsL (line 39)
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)

- K3 ]* O2 L/ G6 |; n" h大家帮分析一下是什么问题？
, Z+ v& y6 R# q  I9 s% S0 }

I_believe

索引可能大于 1。
# s0 g. N( g' Ly0est,yp0est  只能是1个 ？你的有10个。

unfaithful2021

I_believe 发表于 2022-9-15 13:58
0 F) R4 _; z  }$ @; [6 ]- W' }索引可能大于 1。
, A$ d* T, ~6 J/ G' {y0est,yp0est  只能是1个 ？你的有10个。

y0est和yp0est的意义不是方程中参数的初值么？我这个方程有十个参数，是不是应该有十个
4 Z* l* Q$ v! ~3 a! U& t: n7 v

peerless2021

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