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

简单的解函数问题

[复制链接]
  • TA的每日心情
    奋斗
    2022-1-21 15:15
  • 签到天数: 1 天

    [LV.1]初来乍到

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

    EDA365欢迎您登录!

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

    x
    已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0, 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。6 W  @8 h+ |. N+ j, B  \9 S0 M
    代码如下:
    ) |; z* A+ R* k" E% Define the given equations  {3 o. P! m; d- w* ~. c
    syms pc qc N alpha po b pp qq q_ p_
    - U! j1 V6 g1 H+ ^5 Ieq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;- u1 p* B8 \! }
    eq2 = q_ == qc + b * (qq - qc);
    & l8 O) Z8 S, q, I, O- ]9 weq3 = p_ == pc + b * (pp - pc);
    3 P' j* K* Z# z/ ~4 |' Y! n/ Y" S: b: M8 L3 S+ o1 ]
    % Solve for b; T& G' I* H' C! q
    eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
      \3 p* g( ]" z& N4 W" Ldisp(eq4);7 _& I" F# H3 Y+ `
    b = solve(eq4, b);
    . e4 s! j9 K0 e0 U2 Z! {2 u+ C5 ?% d
    : U7 b  q# }0 {& a2 X运行会报错,求各位大佬帮我看看是哪里出错了,感谢!!
    2 q( `- t2 A5 F/ L
  • TA的每日心情
    开心
    2022-1-29 15:04
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2023-2-16 16:25 | 只看该作者
    你混淆了 = 赋值和 == 连接方程两侧的功能。
    # f- z5 {' d, e: K9 Z* W  {改成
    5 i" k" O) _) y" Tsyms pc qc N alpha po b pp qq q_ p_4 _4 y5 y# B+ S- R/ d* C
    eq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0; % eq1是方程,所以才会有 eq1 =(赋值) 方程左侧 == 方程右侧 这样的写法
      ^( j5 }8 {2 c2 seq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”,只使用“=赋值”,不要混用==1 `( r$ Z2 l% ^
    eq3 = pc + b * (pp - pc); % 同理
    8 [7 M2 h( D* jeq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
    " ?0 d! T  p" |( hdisp( eq4 );
    , H$ C( [" g5 Y/ T. Z2 Rb = solve( eq4, b )
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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