简单的解函数问题

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
% N; J( Y2 h, D7 j- B/ Fsyms pc qc N alpha po b pp qq q_ p_
' Q# V! f4 _% l: A# Beq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
eq2 = q_ == qc + b * (qq - qc);
. P% f, C& _7 W5 Keq3 = p_ == pc + b * (pp - pc);
$ x* {, U# Z9 L2 F
% Solve for b
eq4 = subs(eq1, [q_, p_], [eq2, eq3]);
disp(eq4);
b = solve(eq4, b);
2 c0 d, F. e1 v9 d. r: [2 A
运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！
" z3 s8 o: w" `* s- j( y

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
# g8 l: {1 @  T/ E, v% h! D改成
syms pc qc N alpha po b pp qq q_ p_
8 `' _6 o7 f) ~7 v9 Peq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0; % eq1是方程，所以才会有 eq1 =(赋值) 方程左侧 == 方程右侧 这样的写法
eq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”，只使用“=赋值”，不要混用==
+ f+ _: s% u. u# X9 Ueq3 = pc + b * (pp - pc); % 同理
eq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
disp( eq4 );
3 R. N4 ~' }+ W! Wb = solve( eq4, b )