找回密码
 注册
关于网站域名变更的通知
查看: 705|回复: 2
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-6-8 14:53 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

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

x
最近看关于求解隐式微分方程的求解的内容,遇到了这样的问题,对于求解问题
/ U. L& n: ?  z, j4 o

' B' I' u3 n: h* ]7 b  nx0= [1;0;0;1]
2 b& G: ?* p. ^% u/ n  t. {1 B2 Z使用ode15i进行求解,具体的命令如下:) r3 W3 U, O6 O7 s0 b1 [
6 g+ Z9 Z) I/ S  W. L0 i
+ O; C0 c/ `& Q* L  `
odefun=@(t,x,dx)[dx(1)-x(2)
1 Z( ]' ^5 Y7 ydx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
: w8 c( o2 z3 ddx(3)-x(4), c1 n6 a  d8 d/ z1 ^& @
x(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];
9 l2 S/ p" ]" _
+ E1 Y$ Q4 G: H* j- O4 z
$ [% I  M2 B) I: u7 E
t0=0* i* |2 v7 c# b- Y
x0=[1 0 0 1]';%本题初值 x0 的都给出了0 ?' K  D! I, J: Q
fix_x0=ones(4,1);%本题中 x0 都给出了" a% H. d* T! s( r6 e+ u' l

* U9 e) q+ @! l  f1 k! c4 s( X6 \. Q" |dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出,那么全部任意写
0 ]5 R* D  G4 Nfix_dx0=zeros(4,1);%本题中 dx0 一个没有给出,故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);/ B: m  A) Y  Q$ x" J
+ j, V# U6 p! Q5 Y
2 k  L* |6 f- n' Q9 d8 F: z  @
但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示* a8 ~2 b8 O$ U+ P( }5 ~; ?' j! ?
Error using decic (line 109)
7 a0 {2 r( W% _' ^( T- h; z- gConvergence failure in DECIC.
# W( t3 G* J* ?/ j9 A$ v( Y* ~请问出错是为什么呢?
  P5 L, G. w9 k9 D+ f1 ]
  • TA的每日心情
    开心
    2022-12-27 15:46
  • 签到天数: 4 天

    [LV.2]偶尔看看I

    2#
    发表于 2021-6-8 15:32 | 只看该作者
    你可以给一个稍微好一些的初值:
    ( R  C, z4 ^- Z9 c+ ]% jdx0=[0 0 1 -0.5]';%

    该用户从未签到

    3#
     楼主| 发表于 2021-6-8 16:31 | 只看该作者
    qq666888qqw 发表于 2021-6-8 15:32
    ' a# |$ v0 Z% h- Y" g( k你可以给一个稍微好一些的初值:
    ! L( s* h) L8 F: U& y' idx0=[0 0 1 -0.5]';%

    : N% Z: n) l1 u# r6 f* t; l5 e4 f# b$ v0 `3 n) n
    为什么正好初值就是真正解的时候反而就解不出来呢?
    : s, y4 `% S) }" C8 @9 O+ q
    # A0 }8 m3 B, K$ f
    , h, s1 w3 K5 K

    $ Z- E$ A, k/ u# f+ B5 W( B& |$ z1 ]2 Y# H  R& T/ A6 x& x- u5 m4 Z
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-29 10:21 , Processed in 0.156250 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表