简单的解函数问题

nocturne

已知(q_-p_*alpha)^2-(N^2-alpha^2)*p_*(po-p_) = 0， 其中q_ = qc + b * (qq - qc); p_ = pc + b * (pp - pc).  求出b的表达式。
  A5 [- @1 Y. D0 b代码如下：
! J4 _! `: o# E+ B  D. ~  Z- i  i8 y% 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);
8 x% }' l1 D7 S0 J5 u6 R: Eeq3 = p_ == pc + b * (pp - pc);

. @) I4 A1 V) n# X7 H4 k4 F/ \! `% Solve for b
$ Y' W3 x, w% b1 e1 y& _6 x* ]* _; A4 deq4 = subs(eq1, [q_, p_], [eq2, eq3]);
disp(eq4);
" @8 T: ?0 f! a/ t( _  H, q2 |6 lb = solve(eq4, b);
! b8 S1 {4 _* T2 C" D9 h- |" l
/ r9 w8 f9 c" i4 {" t, L运行会报错，求各位大佬帮我看看是哪里出错了，感谢！！
! P& m$ W$ E$ N. p8 y

unfaithful2021

你混淆了 = 赋值和 == 连接方程两侧的功能。
" M; R; A9 M5 t, m改成
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_ 做“变量代换”，只使用“=赋值”，不要混用==
eq3 = pc + b * (pp - pc); % 同理
eq4 = subs( eq1, [ q_, p_ ], [ eq2, eq3 ] );
% t, X" v1 x- p" r7 ~- adisp( eq4 );
  d& @2 B' W( k* d% X8 g! G9 ~$ Gb = solve( eq4, b )