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

flamesmnm

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

) r" Q; {8 M# \5 |, Qx0= [1;0;0;1]
使用ode15i进行求解，具体的命令如下：
2 _3 v* \* G: [+ k' I; I0 D. P

9 O, x7 `6 S# `2 iodefun=@(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)
+ [5 ?$ b% U7 U: ?% ~dx(3)-x(4)
* X) d- I6 J) Ax(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];


t0=0
5 @+ Y4 J; _. O# G/ _9 x! l. Lx0=[1 0 0 1]';%本题初值 x0 的都给出了
. r9 c) g' D# @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);
0 }$ T( i! B: y8 n$ \/ F0 I

但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
+ |: [  S: q4 S- [: ?8 MError using decic (line 109)
7 E9 _/ o1 L3 _! V! dConvergence failure in DECIC.
请问出错是为什么呢？
: X( K4 D0 i% x, S8 p- `5 P+ ~3 G

qq666888qqw

你可以给一个稍微好一些的初值：
! j* B# a' I- }( f5 gdx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
* V0 N; K3 w% w8 a6 I# P, b你可以给一个稍微好一些的初值：
/ j9 J5 x8 w6 `% x8 _dx0=[0 0 1 -0.5]';%

/ P, C- d! o: l, g( ?7 @7 K

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

" x# e& H4 f0 G) {8 u4 P


! \1 c0 _* L  v& N1 W  R