关于ode15i的一阶导数值求解函数decic函数不能收敛的问题

flamesmnm

最近看关于求解隐式微分方程的求解的内容，遇到了这样的问题，对于求解问题

x0= [1;0;0;1]
使用ode15i进行求解，具体的命令如下：


odefun=@(t,x,dx)[dx(1)-x(2)
dx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
dx(3)-x(4)
. D: e9 T! i/ h- i6 Yx(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];
! F( c6 z  t& e# l5 k

, d, N; j2 A3 l; T8 }t0=0
3 p; O" k. Y6 T; ~5 K( j& l* [x0=[1 0 0 1]';%本题初值 x0 的都给出了
+ r& e! h2 Z5 V  _0 Tfix_x0=ones(4,1);%本题中 x0 都给出了

dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
fix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);
4 u1 r) f$ d) I2 H
) f/ ~! h+ q. i1 N; J
但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
" n' i$ W' o- O8 xError using decic (line 109)
7 D+ }; w2 i2 [+ x& U4 z" f7 H) V$ GConvergence failure in DECIC.
请问出错是为什么呢？

qq666888qqw

你可以给一个稍微好一些的初值：
1 O4 u: e+ o! l6 a+ K6 hdx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
' j+ t( N7 _) r; B. ~& N3 b6 ?/ u你可以给一个稍微好一些的初值：
dx0=[0 0 1 -0.5]';%

" e6 O1 s/ ~3 H( w1 N

为什么正好初值就是真正解的时候反而就解不出来呢？

  d0 E) \6 I! x+ ]( B; U: O1 O2 O
( N* M3 O- V( m
5 H- f# e% W% [% |5 w1 Y