简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
代码如下：
) |; z* A+ R* k" E% Define the given equations
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;
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
% Solve for b
eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
  \3 p* g( ]" z& N4 W" Ldisp(eq4);
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

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
# f- z5 {' d, e: K9 Z* W  {改成
5 i" k" O) _) y" Tsyms pc qc N alpha po b pp qq q_ p_
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_ 做“变量代换”，只使用“=赋值”，不要混用==
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 )