简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
代码如下：
% Define the given equations
syms pc qc N alpha po b pp qq q_ p_
* L# V5 X/ L4 D& a% ^' ?" Seq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
' Z. b6 r4 U% Z- v8 z5 Req2 = q_ == qc + b * (qq - qc);
. H/ \4 {8 U9 l1 Beq3 = p_ == pc + b * (pp - pc);
* v: ^/ X0 u6 P$ \! w
: N) D2 X% Y! X  ]% A4 a% Solve for b
eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
disp(eq4);
b = solve(eq4, b);
; u* B- d+ F8 k& j5 k1 u
运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
" ]# L1 D" D( g9 N改成
2 M0 p. l3 |. A$ Q5 a5 i/ R! Gsyms pc qc N alpha po b pp qq q_ p_
  D- E6 i9 W5 t5 oeq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0; % eq1是方程，所以才会有 eq1 =(赋值) 方程左侧 == 方程右侧 这样的写法
eq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”，只使用“=赋值”，不要混用==
6 W, \8 u/ |' T! v! y, @/ {4 j; D" seq3 = pc + b * (pp - pc); % 同理
eq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
disp( eq4 );
b = solve( eq4, b )