简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
* {, t; J5 ~# p& c/ w$ P代码如下：
: D; d& Y4 H$ f% Define the given equations
/ E$ C( i, k6 _' Usyms pc qc N alpha po b pp qq q_ p_
eq1 = (q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) == 0;
8 |3 _  M1 d" V- @: m+ A  qeq2 = q_ == qc + b * (qq - qc);
eq3 = p_ == pc + b * (pp - pc);

$ V7 S' {8 z, h/ G, R* S% Solve for b
9 j6 b# i. n+ X4 k: Ceq4 = subs(eq1, [q_, p_], [eq2, eq3]);
disp(eq4);
+ M. Z- M+ m' Ab = solve(eq4, b);
5 `+ Y, n- ]  e6 p  z# H
  G/ w% E. r6 m+ }3 n/ [. c运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！
% \- K3 Q! {% J2 ~4 y9 y; l% Z

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
改成
( v/ R# r$ @' ?9 p6 N4 z6 h% [syms 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 =(赋值) 方程左侧 == 方程右侧 这样的写法
eq2 = qc + b * (qq - qc); % eq2纯粹只是为了随后给 q_ 做“变量代换”，只使用“=赋值”，不要混用==
) I8 ]/ |4 I% ?1 i6 r0 Xeq3 = pc + b * (pp - pc); % 同理
* U, _% H+ \( U" b6 E$ veq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
8 a+ g( g  q4 K5 L  P' Z: K1 \disp( eq4 );
b = solve( eq4, b )