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

flamesmnm

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

, u' @+ l, h9 a" Q$ L

5 ]  i, \3 L& |3 _* E+ W- Dx0= [1;0;0;1]
使用ode15i进行求解，具体的命令如下：


odefun=@(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)
dx(3)-x(4)
x(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];
+ {% p7 H- `' A* `0 g/ _& h# ~
! k/ P, e4 s' b2 M
t0=0
x0=[1 0 0 1]';%本题初值 x0 的都给出了
fix_x0=ones(4,1);%本题中 x0 都给出了

dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出，那么全部任意写
- u5 B( y  F( v  ?: ?fix_dx0=zeros(4,1);%本题中 dx0 一个没有给出，故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);

' B2 \: s5 m: d* m& l) G/ d- x
- b% r& F8 |$ t7 G9 D但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
Error using decic (line 109)
Convergence failure in DECIC.
# y$ |- j; N/ G& m( q$ t请问出错是为什么呢？

qq666888qqw

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

flamesmnm

qq666888qqw 发表于 2021-6-8 15:32
) o8 P9 c; v6 S) l9 d3 Z你可以给一个稍微好一些的初值：
' b- }) l  q4 a& |# Q5 Hdx0=[0 0 1 -0.5]';%

8 w) H; F6 i: A* Q

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

/ ^+ U" X  F& g. F" H
+ R/ q0 M9 w* f& Y) t- H) k+ W# ?