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

使用matlab中的fsolve求解非线性方程时报错

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

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2023-2-7 15:47 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

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

    x
    本帖最后由 peerless2021 于 2023-2-7 15:48 编辑
    / a( l% X$ r2 L% W1 B
    & r. i: A; C" {# Q& d" T原程序为
    • function F = myeqn(alpha,L)
    • F = (0.3*((9.8)*L^1.83)/(2*cos(alpha)^2)+(1.96*L^3)/(3*cos(alpha)^3)-10)/(0.9*tan(alpha));
    • end
    • alpha0 = 0;
    • L0 = 0;
    • x0 = [alpha0,L0];
    • x = fsolve(@myeqn,x0);
      ) I0 q3 y: _+ B) E
    在调用时会报错截图如下:
    0 p) D7 F* J; m, ^# H6 ^/ E
    * u" w% e# q' c& Q4 ]" f6 Y; }
    原方程比较复杂,如图:5 F; q+ r, t9 i
    ' x' D* {6 P4 ~! ?) A+ k
    是调用的方式出现了问题吗?

      A2 v" i1 w1 u$ \- f0 S% R3 p& N5 ?, G" J

    ; R1 |: q( j. m+ J9 X* g

    该用户从未签到

    2#
    发表于 2023-2-7 16:37 | 只看该作者
    你这样写不符合 fsolve函数(超链接,自己点进去看) 的语法要求。
      x  S2 y" Z5 `" f+ X仔细看帮助文档里的示例,以 x = fsolve(@myfun,x0) 这种语法调用时,要求输入函数形为 function F = myfun(x) 这样的“单输入”函数,你有两个输入,当然不行了。要么改写为
    4 o/ q( C7 E7 ]1 }function F = myeqn( Input )
    9 J; Y7 c4 Q& Nalpha = Input( 1 );
      K6 Y/ a$ l" G  b+ s! nL = Input( 2 );
    ( n1 ]! L' P& ?3 ]; C8 Q5 aF = (0.3*((9.8)*L^1.83)/(2*cos(alpha)^2)+(1.96*L^3)/(3*cos(alpha)^3)-10)/(0.9*tan(alpha));
    # l5 t3 l5 d4 h! E, c& T, n1 Kend# D, m' b* y" f' w8 {; n! J, M
    这样的“单输入”,然后用 x0 = [ alpha0, L0 ]; x = fsolve( @myeqn, x0 ) 来求“满足方程的数值向量”(另外还得注意你这初值取得不合适,导致分母出现了零,新手要避免直接取全零或全1这种简单值作为初值,很容易导致各种无穷或零出现在分母之类位置上,导致迭代运算无法继续)
    " F3 @% n% P; Z/ i  |9 y! a要么不改动自定义函数,但需要另行指定一个输入的值,求另一个,需要改写为9 m/ U; q) n, ?
    x0 = [ alpha0 ] + eps; % 避免初值无效0 @7 G1 ^; Z/ }+ g& C5 z
    x = fsolve( @( t ) myeqn( t, L0 ), x0 )9 o' x4 r9 Z+ ]( e: p$ X' c
    disp( myeqn( x, L0 ) )# e$ K) H1 m: g6 m; {% f3 L# g
    最后,你得好好检查检查你的函数,你为什么要多此一举地把 0.9*tan(alpha) 扔到分母上去。原方程是1 d8 J0 J0 X% h7 o/ v! D7 s/ j
    0 == 0.9*tan(alpha) + Others
    8 t* \8 m# h  r7 v/ F' y怎么你改写就成了 0 == Others / ( 0.9*tan(alpha) ) 了?

    点评

    是的,我的函数搞错了,是我打的时候想多了。按照您的两个方式我都尝试了一下,程序可以正常运行了,但是我没明白是在哪一步修改初值,最后得到的结果数字都是没输入过的  详情 回复 发表于 2023-2-7 16:45
  • TA的每日心情
    开心
    2022-1-29 15:04
  • 签到天数: 1 天

    [LV.1]初来乍到

    3#
     楼主| 发表于 2023-2-7 16:45 | 只看该作者
    aid4her 发表于 2023-2-7 16:37. i; b" Z3 v' D5 \
    你这样写不符合 fsolve函数(超链接,自己点进去看) 的语法要求。
    $ D/ l* a8 {+ A6 k$ P6 ~8 c仔细看帮助文档里的示例,以 x = fsolv ...
    & M  d& ~+ b1 h
    是的,我的函数搞错了,是我打的时候想多了。按照您的两个方式我都尝试了一下,程序可以正常运行了,但是我没明白是在哪一步修改初值,最后得到的结果数字都是没输入过的
    # _* o: I& `2 |! a
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-10-28 00:10 , Processed in 0.140625 second(s), 27 queries , Gzip On.

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

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

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