NumPy复相位角处理疑似Bug：叠加复平面波相位计算结果不一致问题

你遇到的问题并不是NumPy的bug，而是理论推导时忽略了振幅项的符号导致的相位偏移差异，下面详细解释：

问题本质分析 #

先看叠加波的正确数学推导：两个复平面波相加可以化简为：
$$
e^{jk_1x} + e^{jk_2x} = e^{j\frac{k_1+k_2}{2}x} \left( e^{-j\frac{k_2-k_1}{2}x} + e^{j\frac{k_2-k_1}{2}x} \right) = 2\cos\left( \frac{(k_2-k_1)x}{2} \right) e^{j\theta}
$$
其中 $\theta = \frac{k_1+k_2}{2}x$ 是你定义的平均相位。

你的理论实现中，rho = 2 + 2cos((k2-k1)*x) 其实等于 $4\cos^2\left( \frac{(k2-k1)x}{2} \right)$，所以 sqrt(rho) = 2|cos\left( \frac{(k2-k1)x}{2} \right)|——这里的绝对值是关键！你构造的superwave_theory = sqrt(rho)*exp(1j*theta) 等价于 $2|\cos(\dots)|e^{j\theta}$，但实际叠加波是 $2\cos(\dots)e^{j\theta}$：

• 当 $\cos\left( \frac{(k2-k1)x}{2} \right) > 0$ 时，两者完全一致，相位也相同；
• 当 $\cos\left( \frac{(k2-k1)x}{2} \right) < 0$ 时，实际叠加波等价于 $2|\cos(\dots)|e^{j(\theta+\pi)}$，相位比你的理论实现多了一个 $\pi$，而NumPy的angle()函数会把这个相位包裹到 $(-\pi, \pi]$ 区间，最终表现为相位差了 $-\pi$（比如 $\theta+\pi$ 超过 $\pi$ 时，会被减去 $2\pi$ 得到 $\theta-\pi$）。

这就是两种方法相位结果不一致的核心原因——你的理论实现丢掉了余弦项的负号，导致相位少了一个关键的偏移。

修正后的理论实现 #

把理论推导的符号考虑进去，修正代码如下：

from pylab import *
x = linspace(0,2*pi)
k1 = 1
k2 = 2
planewave1 = exp(1j*k1*x)
planewave2 = exp(1j*k2*x)
superwave = planewave1 + planewave2

# 修正后的理论计算
dk = k2 - k1
theta = (k1 + k2)/2 * x
cos_term = cos(dk * x / 2)
superwave_theory_correct = 2 * cos_term * exp(1j * theta)

phase_Numpy = angle(superwave)
phase_theory_correct = angle(superwave_theory_correct)

# 验证一致性
print(allclose(phase_Numpy, phase_theory_correct))  # 输出 True

这样修正后，两种方法计算的相位就完全一致了。

补充说明 #

NumPy的angle()函数工作正常：它会根据复数的实部和虚部的符号，严格将相位映射到 $(-\pi, \pi]$ 区间。你看到的“偏差$\pi$”，本质是因为错误的理论实现没有考虑振幅项的负号对应的相位偏移，而不是NumPy的复算术处理有问题。

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