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

flamesmnm

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

x0= [1;0;0;1]
使用ode15i进行求解，具体的命令如下：
7 g9 x% D$ Q* r7 {3 m: f7 {' M! l

odefun=@(t,x,dx)[dx(1)-x(2)
& z3 D+ _/ B& `: B2 [! M$ K1 B& udx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
& _$ {' _! f! ^; w0 p+ ^" Jdx(3)-x(4)
x(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];
3 G% M8 I4 ?( c. x! X  p; c
7 G5 o5 Y5 C# [* m2 B  J7 t1 Y0 ^' c
& I! Q; @  w" Q! N* |& f7 `3 j% _t0=0
% \4 E$ E3 @; C" hx0=[1 0 0 1]';%本题初值 x0 的都给出了
fix_x0=ones(4,1);%本题中 x0 都给出了
% Z* H" D0 T: r4 k8 [& C; A' ]( x
dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
- h: Z9 s; Q' q/ b9 R0 @/ f0 Dfix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);
0 S/ ?* W& L2 t- k0 @! I! ^7 t

但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
Error using decic (line 109)
Convergence failure in DECIC.
2 C- o/ A0 k. F* q+ E  V请问出错是为什么呢？
" j$ v# R* n8 ?: A- ]( r! v

qq666888qqw

你可以给一个稍微好一些的初值：
" A. n1 @$ y* Edx0=[0 0 1 -0.5]';%

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
你可以给一个稍微好一些的初值：
) H3 a- T  M9 R7 x: P3 Cdx0=[0 0 1 -0.5]';%

3 m( O. W, M' m! g

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

8 [* C- [9 s  c1 h
4 z/ n  c; \4 j- h5 [