复数表示法在光波描述中的应用

在现代光学研究与工程设计中，精确描述光波的传播、干涉、衍射等行为是不可或缺的基础。光本质上是一种高频电磁波，其电场强度随空间和时间做快速振动。如果直接使用三角函数（如正弦和余弦）来描述光场的振动，在进行复杂的数学运算（如多光束干涉叠加、傅里叶光学衍射计算）时，计算过程会变得异常繁琐。

为了化繁为简，复数表示法成为了波动光学中最核心的数学工具之一。通过引入复振幅，原本复杂的三角函数运算被转化为了更为简便的代数运算。
一束单色平面光波的电场强度 $E(r, t)$ 通常可以表示为空间位置 $r$ 和时间 $t$ 的函数：

$$E(r, t) = A(r) \cos(\omega t - k \cdot r + \varphi_0)$$

其中，$A(r)$ 是振幅，$\omega$ 是角频率，$k$ 是波矢，$\varphi_0$ 是初相位。在处理光学问题时，时间因子 $\cos(\omega t)$ 在大多数稳态干涉和衍射实验中可以通过略去或分离来简化。因此，我们更关注光波的空间分布特征。

利用欧拉公式：

$$e^{i\theta} = \cos\theta + i\sin\theta$$

我们可以将实数域的光波表达式嵌入到复数域中。定义光波的复振幅（Complex Amplitude）为：

$$\tilde{U}(r) = A(r) e^{i[\varphi(r)]}$$

此时，瞬时电场可以表示为复振幅与其时间因子的乘积的实部：

$$E(r, t) = \text{Re} \left{ \tilde{U}(r) e^{-i\omega t} \right}$$

通过这种转换，我们成功地将包含时间和空间的动态波动方程，转化为仅处理空间复振幅的稳态数学模型。

复数表示法的核心优势

在波动光学分析中，采用复数表示法主要具备以下几大优势：

• 化乘积为加和（简化相位叠加）： 在处理多列光波叠加（如干涉现象）时，若使用三角函数，需要运用繁杂的和差化积公式。而采用复数表示，光波的叠加直接表现为复振幅的线性相加：
  $$\tilde{U}_{\text{total}} = \tilde{U}_1 + \tilde{U}_2 + \dots + \tilde{U}_N$$
  这极大地降低了多光束干涉和复杂衍射积分的计算难度。
• 方便处理微分方程： 波动方程本质上是二阶偏微分方程。指数形式的复函数求导极其简便（例如对空间坐标求导只需乘以波数常数），这使得研究光在介质中的传播（如亥姆霍兹方程）变得直观高效。
• 直观体现光强与相位的关系： 在物理实验中，光电探测器响应的是光强 $I$，它正比于电场振幅平方的时间平均值。在复数表示法中，光强可以非常简洁地表示为复振幅与其共轭复数的乘积：
  $$I(r) \propto |\tilde{U}(r)|^2 = \tilde{U}(r) \cdot \tilde{U}^*(r)$$
  而相位相对于参考面的滞后或超前，则直接体现在复数的辐角中。

波动光学中的通用应用全景

复数表示法贯穿了整个波动光学体系，是连接基础物理原理与现代光学工程的桥梁。

• 干涉与衍射的统一处理： 无论是在杨氏双缝干涉、法布里-珀罗干涉中，还是在夫琅禾费衍射、菲涅耳衍射的积分计算中，复振幅叠加原理都是核心。通过将孔径函数或透过率函数定义为复数（如复振幅光栅、全息图），可以同时对光波的振幅调制和相位调制进行建模。
• 傅里叶光学与系统建模： 在现代光学中，透镜对光波的变换作用可以用傅里叶变换来描述。空间频域分析离不开复指数基底，复数表示法使我们能够将光学系统视为线性时不变系统（LSI系统），利用传递函数全面评估成像质量。
• 光波与物质相互作用（色散与偏振）： 当光穿过介质时，介质的折射率常被处理为复折射率 $\tilde{n} = n + i\kappa$，其中虚部直接对应光的吸收或增益（与光的色散和衰减相关）。同样，在描述光的偏振态时，正交的两个电场分量也可以通过复数形式组合成琼斯矢量（Jones Vectors），用于优雅地计算各种偏振器件（如波片、起偏器）对光束的作用。

编程实践示例：复振幅叠加计算

在计算光学仿真中，复数表示法是编写代码的基础。以下是一个使用 Python 模拟两束相干光干涉的简单示例，展示了复振幅叠加的实际应用：

import numpy as np
import matplotlib.pyplot as plt

# 定义空间网格
x = np.linspace(-0.01, 0.01, 500)
y = np.linspace(-0.01, 0.01, 500)
X, Y = np.meshgrid(x, y)

# 光波参数
wavelength = 632.8e-9  # 氦氖激光波长 (m)
k = 2 * np.pi / wavelength
z = 0.5                # 传播距离 (m)

# 模拟两个点光源产生的球面波复振幅
# 源1位置 (0, -d/2, 0)，源2位置 (0, d/2, 0)
d = 1e-3
r1 = np.sqrt(X**2 + (Y - d/2)**2 + z**2)
r2 = np.sqrt(X**2 + (Y + d/2)**2 + z**2)

# 计算复振幅 (忽略常数项，仅看相位与距离衰减)
U1 = (1 / r1) * np.exp(1j * k * r1)
U2 = (1 / r2) * np.exp(1j * k * r2)

# 复振幅叠加
U_total = U1 + U2

# 计算光强 (振幅平方)
Intensity = np.abs(U_total)**2

# 绘图展示干涉条纹
plt.figure(figsize=(6, 5))
plt.imshow(Intensity, extent=[-10, 10, -10, 10], cmap='gray')
plt.title('Interference Pattern via Complex Amplitude')
plt.xlabel('X (mm)')
plt.ylabel('Y (mm)')
plt.colorbar(label='Intensity')
plt.show()

通过上述代码可以看出，利用 Python 的复数计算能力，原本复杂的三角函数叠加被简化为了纯粹的复数加法，这正是复数表示法在波动光学中强大生命力的体现。