数值解法与傅里叶变换的应用

在光学领域，描述光波如何在空间中演化是核心任务。传统的几何光学处理光线路径，而波动光学则关注光的干涉与衍射现象。当我们需要精确计算复杂孔径或任意形状障碍物后的光场分布时，解析解往往难以求得，此时数值解法便成为不可或缺的工具。本文将深入探讨如何利用计算机算法模拟光传播过程，并揭示傅里叶变换在其中扮演的关键角色。

菲涅尔与夫琅禾费衍射的数值逼近基础

要实施数值计算，首先必须明确所研究的物理场景属于哪种衍射 regime（区域）。在光学中，通常将远场近似为夫琅禾费衍射，而近场则对应菲涅尔衍射。这两种情形决定了我们选择何种算法进行迭代或变换。

对于菲涅尔衍射，光波前的相位变化与传播距离 $z$ 呈线性关系，其传递函数较为复杂。数值模拟通常采用角谱法（Angular Spectrum Method）或分步傅里叶变换法。在这些方法中，我们将复振幅分布表示为空间频率的叠加，通过引入系统的传输矩阵来推进光场演化。

菲涅尔衍射公式：在傍轴近似下，传播距离 $z$ 后的光场 $U(x, y)$ 可表达为输入光场与一个二次相位因子的卷积积分形式。
$$ U(x,y,z) = \frac{e^{ikz}}{i\lambda z} \iint U_0(\xi,\eta) e^{\frac{i\pi}{\lambda z}[(x-\xi)^2 + (y-\eta)^2]} d\xi d\eta $$
数值离散化：在实际编程中，我们将上述连续积分转化为离散的求和运算。利用快速傅里叶变换（FFT）可以将卷积操作高效地转换为频域的乘法操作，从而大幅降低计算复杂度。

角谱法的核心机制与实现步骤

角谱法是处理光波传播最经典且高效的数值方法之一，其核心思想是将空间域的光场分解为不同方向的平面波分量，分别计算每个分量的相位延迟后重新合成。这种方法在处理大孔径或复杂介质时尤为优势明显。

具体实施过程通常包含以下关键步骤：

光场采样与傅里叶变换：首先对初始时刻 $t=0$ 的光场复振幅分布进行离散化，并计算其二维快速傅里叶变换（FFT），得到空间频谱 $\tilde{U}(f_x, f_y)$。此时，每个频率分量对应一个特定传播方向的平面波。
频域相位调制：根据波动方程推导出的色散关系 $k_z = \sqrt{k^2 - k_x^2 - k_y^2}$（其中 $k$ 为真空中的总波数），计算每个空间频率分量在传播距离 $z$ 后的纵向波矢变化。这一步相当于给频谱乘以相位因子 $\exp(i k_z z)$，实现了光场向前的“时间”推进或“空间”演化。
逆傅里叶变换重构：将调制后的频谱通过逆变换（IFFT）转换回空间域，即得到距离 $z$ 处的光场分布。

这种方法的优势在于其并行计算能力强，且能够自然地处理非傍轴条件下的精确传播问题。例如，在模拟激光束经过透镜聚焦的过程时，角谱法可以准确描述球面波的会聚行为，而无需引入额外的近似条件。

傅里叶变换作为光学系统的传递函数解析器

除了直接的时间步进算法外，许多复杂的衍射系统可以通过傅里叶变换的性质进行简化分析。根据菲涅尔 - 基尔霍夫衍射理论，在夫琅禾费区域（即远场），光场的分布实际上等同于孔径函数的傅里叶变换。这一性质使得光学系统设计变得直观且强大。

当平行光照射到具有复杂透射函数 $T(x,y)$ 的透镜或掩模上时，其焦平面上的光强分布直接由该函数的频谱决定：

空间滤波：如果我们在频域放置一个特定的滤波器（如低通、高通带通），就可以选择性地保留或抑制某些空间频率成分。这在图像处理和光学信息处理中应用广泛。
透镜的傅里叶变换器角色：理想薄透镜可以在其焦平面上实现输入光场的傅里叶变换操作，前提是入射光是平面波且满足傍轴近似。这意味着我们可以通过调整透镜的位置和焦距来改变观察到的频谱结构。

例如，在光学全息术或相干成像系统中，利用这一原理可以设计出特定的衍射元件（如计算机生成全息图 CGH），通过编程计算其透射函数 $T(x,y)$ 的傅里叶逆变换，进而合成出期望的光场分布。这种“逆向工程”的设计思路极大地扩展了光学的功能边界。

总结与展望：数值模拟在光学设计中的价值

综上所述，数值解法结合傅里叶变换为我们提供了一种强大的工具箱来理解和预测光的传播行为。从菲涅尔区域的逐步推进到远场的频谱分析，这些方法不仅解决了传统解析公式无法覆盖的复杂场景问题，也为现代自适应光学、光刻技术和量子光学实验提供了理论支撑。随着计算机算力的提升和算法精度的优化，未来的数值模拟将更加精细地刻画非线性效应和多模态光的传播特性，推动光学技术向更高维度和更智能的方向发展。