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

flamesmnm

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

# `% m$ o9 r& q4 [5 ^6 {x0= [1;0;0;1]
使用ode15i进行求解，具体的命令如下：


odefun=@(t,x,dx)[dx(1)-x(2)
* _) _0 m0 V  N# R8 O8 Jdx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
, E, j6 w7 K6 |0 j3 wdx(3)-x(4)
( S6 i* }/ Q; p* j5 V5 {% Sx(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];

! l; D6 L: m1 V! T6 w
) Y( A" w% R) k, J+ p( j4 s# Jt0=0
x0=[1 0 0 1]';%本题初值 x0 的都给出了
. q1 Z/ f" N) j7 l( Y8 cfix_x0=ones(4,1);%本题中 x0 都给出了
# `2 z, a' E- i
dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
fix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);
* u. `( ?* a7 D1 U9 J

但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
& }8 _: H, j. s3 ?8 KError using decic (line 109)
Convergence failure in DECIC.
请问出错是为什么呢？

qq666888qqw

你可以给一个稍微好一些的初值：
1 y) v7 G2 c2 M- Z2 rdx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
你可以给一个稍微好一些的初值：
$ D" U+ a' B- |1 a& L* H; g9 Jdx0=[0 0 1 -0.5]';%

$ l9 ^  t* z% ?' q! b' q, n1 A6 v
; q; w5 M0 H( J2 I

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

$ M% Y2 ?" y8 |