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

flamesmnm

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

  r9 x" n4 b' F

x0= [1;0;0;1]
& b& N2 c2 D5 w使用ode15i进行求解，具体的命令如下：
$ f9 K/ d* O8 [& i' |( t! q' [
, b, v$ f9 ]8 y3 |% Y0 V" W
8 F7 D: S. z2 I2 S: yodefun=@(t,x,dx)[dx(1)-x(2)
$ M" P5 z! k; O3 ?3 y0 I1 idx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
1 F3 m- i. A; L' d( S' P* w& v) Rdx(3)-x(4)
# s1 T5 y3 W3 Y9 f: u& Xx(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];

, \6 J9 w) |. i: p. b  o5 `
3 A9 F5 R9 l7 f% et0=0
. C1 ^4 C2 G9 ~9 n& i7 o# K. O. Qx0=[1 0 0 1]';%本题初值 x0 的都给出了
fix_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);
. n0 h# h/ t* v% N) d  Z# n# N8 o
: A) R8 a# @4 N# r% A/ W  Z
% G7 z6 f9 |% Q! @. P$ s0 v, }但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
2 y! ?5 g3 W* ?1 Z6 B& o- E) u$ dError using decic (line 109)
Convergence failure in DECIC.
请问出错是为什么呢？
8 N1 f. `7 {, q# Y( W

qq666888qqw

你可以给一个稍微好一些的初值：
dx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
你可以给一个稍微好一些的初值：
dx0=[0 0 1 -0.5]';%

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

9 h# A7 G  Z% n' |* V7 K( ]8 a) o


2 t3 }) u  x) h, }6 G) o