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

使用 decic命令 错误分析

[复制链接]
  • TA的每日心情
    开心
    2022-1-29 15:04
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2022-9-15 13:15 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

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

    x
    我的M文件:
    - Q, i- H1 f3 x) }8 v; by0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
    6 ^! d; o) v0 ^/ [4 x5 Eyp0est = zeros(10,1);& V) P* t9 f1 w$ n' ~
    opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));  ^) @9 x2 L; b2 @+ j
    [y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)- z! J; f3 J8 `- r3 Q- A1 L  E, P( m
    5 b; s# s' e' K  u; J
    & @3 j" ~! q" G# U# o4 [
    运行后报错:! f! O+ h, c( m3 u* _* z' ~
    错误使用 decic>sls (line 170)
      p5 b$ N' m( t! f2 ]/ S- H索引可能大于 1。
    0 ~$ m" r$ q' ~  O% t2 r
    # [; g' d6 Z& o! T. S出错 decic (line 77)6 g  H$ z9 `4 X3 @9 s3 S
        [dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);" O7 ]5 O5 c  [7 q& L

    5 l& Z; U& P6 {  T出错 DAEsL (line 39)
    , K$ b0 C* d3 t" m# h[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
    4 Q7 f; a9 |  z! A) z: {, A1 P0 j4 [! \0 M5 B7 `5 e
    大家帮分析一下是什么问题?- ^, s8 w0 r7 A- k7 H7 K+ `6 F
  • TA的每日心情
    开心
    2022-1-29 15:03
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-9-15 13:58 | 只看该作者
    索引可能大于 1。: p- s+ b+ r0 D- v! f0 {
    y0est,yp0est  只能是1个 ?你的有10个。

    点评

    y0est和yp0est的意义不是方程中参数的初值么?我这个方程有十个参数,是不是应该有十个  详情 回复 发表于 2022-9-15 14:03
  • TA的每日心情
    开心
    2022-1-29 15:04
  • 签到天数: 1 天

    [LV.1]初来乍到

    3#
     楼主| 发表于 2022-9-15 14:03 | 只看该作者
    I_believe 发表于 2022-9-15 13:581 L' }* x8 X  x! v0 Z
    索引可能大于 1。
    6 u' ?+ C# V& j9 O+ ly0est,yp0est  只能是1个 ?你的有10个。

    * _3 A6 W9 ]% x- l4 R  X. py0est和yp0est的意义不是方程中参数的初值么?我这个方程有十个参数,是不是应该有十个8 g8 [: h# e! ], A
  • TA的每日心情
    开心
    2022-1-29 15:04
  • 签到天数: 1 天

    [LV.1]初来乍到

    4#
    发表于 2022-9-15 14:05 | 只看该作者
    Decic调用ode15ipdinit(),它产生几个结果,包括ode函数的偏导数。然后将第二个偏导数传递给decic内部例程sls(),该例程计算其秩。如果rank小于固定变量的数量,那么sls()将发出您看到的错误消息。简单说就是您的自变量,超过等式的数量了。
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 14:54 , Processed in 0.156250 second(s), 25 queries , Gzip On.

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

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

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