简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
: C% c6 z. R: `/ H5 K3 K$ v8 ]2 H代码如下：
% Define the given equations
5 i4 o# p- j4 O( `9 V' bsyms pc qc N alpha po b pp qq q_ p_
, @) K7 v4 A" Y7 `9 neq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
eq2 = q_ == qc + b * (qq - qc);
: K/ Z1 O! ?' M: A# L& ]eq3 = p_ == pc + b * (pp - pc);

9 D5 r: a9 D) U. b" `6 f% Solve for b
# K3 `& Y3 t% X( j) _eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
8 U# y9 U6 _5 T. \) n. Pdisp(eq4);
b = solve(eq4, b);

0 @; M5 {9 S* V( j0 V6 M运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
+ h2 w9 \1 s& S7 q改成
syms pc qc N alpha po b pp qq q_ p_
# }: d5 d0 N; Z" zeq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0; % eq1是方程，所以才会有 eq1 =(赋值) 方程左侧 == 方程右侧 这样的写法
; \% Q0 V4 x7 L6 O  r  f# Jeq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”，只使用“=赋值”，不要混用==
) n8 l; p9 h$ e& Weq3 = pc + b * (pp - pc); % 同理
eq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
disp( eq4 );
b = solve( eq4, b )