如何用SymPy求解含共轭复数的方程及方程组？

问题本质 #

直接用symbols('x')定义的复变量，SymPy的solve()无法自动识别共轭操作对应的连续解集（比如整个实轴/虚轴），只会尝试寻找离散解，导致输出不符合预期。

可靠解决方案：拆分实部与虚部 #

将复变量拆分为独立的实部、虚部变量，把共轭操作转化为实部不变、虚部取反的形式，再求解方程，就能得到清晰的约束条件。

代码示例 #

from sympy import symbols, solve, pprint

# 定义实部a和虚部b，明确约束为实数
a, b = symbols('a b', real=True)
x = a + b*1j  # 复变量x = a + bi

# 求解x + conjugate(x) = 0
eq1 = x + x.conjugate()
sol1 = solve(eq1, (a, b))
print("x + conjugate(x) = 0 的解：")
pprint(sol1)

# 求解x - conjugate(x) = 0
eq2 = x - x.conjugate()
sol2 = solve(eq2, (a, b))
print("\nx - conjugate(x) = 0 的解：")
pprint(sol2)

结果解释 #

• 第一个方程的解为{a: 0}：表示实部必须为0，虚部b可取任意实数，对应复平面的整个虚轴。
• 第二个方程的解为{b: 0}：表示虚部必须为0，实部a可取任意实数，对应复平面的整个实轴。

进阶：直接表示连续解集 #

如果需要用集合形式直接描述解，可以利用SymPy的实数域符号：

from sympy import S, I

# 实轴：所有虚部为0的复数
real_axis = S.Reals
# 虚轴：所有实部为0的复数
imaginary_axis = I * S.Reals

原方法失效原因 #

SymPy的solve()默认优先寻找离散解，而含共轭的方程常对应连续解集（如整条轴），此时直接调用solve()无法识别这类连续约束，必须通过拆分实虚部将问题转化为实变量方程组，才能得到准确的解描述。

内容的提问来源于stack exchange，提问作者Joshua Ruiter