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

flamesmnm · 发表于 2021-6-8 14:53

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

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

x0= [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)];

t0=0
x0=[1 0 0 1]';%本题初值 x0 的都给出了
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);

但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示
Error using decic (line 109)
Convergence failure in DECIC.
请问出错是为什么呢？

qq666888qqw · 发表于 2021-6-8 15:32

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

flamesmnm · 发表于 2021-6-8 16:31

qq666888qqw 发表于 2021-6-8 15:32/ l, a6 ]" \" T E! l: s& J- k
你可以给一个稍微好一些的初值：( F: |. W8 _( a1 r7 D. k
dx0=[0 0 1 -0.5]';%

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

帐号		自动登录	找回密码
密码			注册

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

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1