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

flamesmnm

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

' B' I' u3 n: h* ]7 b  nx0= [1;0;0;1]
2 b& G: ?* p. ^% u/ n  t. {1 B2 Z使用ode15i进行求解，具体的命令如下：


odefun=@(t,x,dx)[dx(1)-x(2)
1 Z( ]' ^5 Y7 ydx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
: w8 c( o2 z3 ddx(3)-x(4)
x(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];


t0=0
x0=[1 0 0 1]';%本题初值 x0 的都给出了
fix_x0=ones(4,1);%本题中 x0 都给出了

* U9 e) q+ @! l  f1 k! c4 s( X6 \. Q" |dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
0 ]5 R* D  G4 Nfix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);


但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
Error using decic (line 109)
7 a0 {2 r( W% _' ^( T- h; z- gConvergence failure in DECIC.
# W( t3 G* J* ?/ j9 A$ v( Y* ~请问出错是为什么呢？

qq666888qqw

你可以给一个稍微好一些的初值：
( R  C, z4 ^- Z9 c+ ]% jdx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
' a# |$ v0 Z% h- Y" g( k你可以给一个稍微好一些的初值：
! L( s* h) L8 F: U& y' idx0=[0 0 1 -0.5]';%

: N% Z: n) l1 u# r6 f* t

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

$ Z- E$ A, k/ u# f+ B5 W( B& |