简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
8 c# j6 }8 E* X$ n- t3 j' Z$ a, V  C代码如下：
% Define the given equations
syms pc qc N alpha po b pp qq q_ p_
eq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
eq2 = q_ == qc + b * (qq - qc);
5 B2 a8 X7 {& h" i/ r4 Req3 = p_ == pc + b * (pp - pc);
0 |3 k" o  J+ z2 \0 p- ?2 A3 ]
% Solve for b
& X6 _5 P) g4 A/ V- B, Yeq4 = subs(eq1, [q_, p_], [eq2, eq3]);
disp(eq4);
* W5 n8 U& I1 f4 P' Q% H4 s+ rb = solve(eq4, b);
7 @: o' k# h9 L+ H0 K
- \0 Y) G7 ~- n' s/ F运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！
8 X0 E: v/ t9 d5 c

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
改成
* N5 m$ Q, ~4 j, f& t2 Hsyms pc qc N alpha po b pp qq q_ p_
+ l4 d# U+ U  V* X- Req1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0; % eq1是方程，所以才会有 eq1 =(赋值) 方程左侧 == 方程右侧 这样的写法
) {3 u/ a, u# d* ]2 m2 ueq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”，只使用“=赋值”，不要混用==
( `% _8 {. j1 u% @' Z6 ~/ l# beq3 = pc + b * (pp - pc); % 同理
eq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
disp( eq4 );
( H8 k4 o5 p! ^. Y) i9 B" eb = solve( eq4, b )