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

flamesmnm

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

x0= [1;0;0;1]
& q5 U1 g% d  s使用ode15i进行求解，具体的命令如下：

" N. N8 ?8 P1 h7 V$ P
odefun=@(t,x,dx)[dx(1)-x(2)
: Q/ u3 b% E7 r9 Y, J$ z6 Tdx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
dx(3)-x(4)
2 n8 g; [5 d! `! kx(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];
) B$ N! {) D- i" Q& o

8 p( I9 q) b) p" H! it0=0
* c9 n1 E- \% C4 D( Px0=[1 0 0 1]';%本题初值 x0 的都给出了
fix_x0=ones(4,1);%本题中 x0 都给出了
) w8 W1 w5 `9 _/ [- `
/ `- i6 N: D# a6 F& s, xdx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
fix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);
. o0 s4 D2 {. X. @# d

/ i  G- E+ v6 b' M. y但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
. e2 Q1 _/ J( S! r( K" QError using decic (line 109)
5 B: K5 p) s7 V+ ?- B$ u# oConvergence failure in DECIC.
请问出错是为什么呢？
2 P' A* Z, ~: o, V

qq666888qqw

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

flamesmnm

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

, k3 S" ]+ [4 M5 h( C  o2 q& A

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

0 z4 F9 \( [- m
5 n( K5 r8 P# B: L