3. 标量衍射理论
标量衍射理论是信息光学中描述光波传播与衍射行为的数学基石。在工程实践中，无论是设计傅里叶变换透镜、分析空间滤波系统，还是评估全息记录的分辨率极限，都需要对衍射场的复振幅分布进行定量计算。本章从惠更斯-菲涅尔原理的物理图像出发，引出基尔霍夫衍射公式的数学表达，并重点剖析菲涅尔近似与夫琅禾费近似的物理含义、适用判据及典型工程场景。理解这些近似条件的边界，是避免在系统设计中产生重大误差的关键。
3.1 惠更斯-菲涅尔原理：衍射的物理图像
惠更斯-菲涅尔原理是理解衍射现象的起点。该原理指出：波前上的每一点都可以视为一个新的次级球面波源，而空间中任意观察点的光场是所有这些次级波在该点相干叠加的结果。这一思想将连续波前的传播问题转化为无数点源的叠加积分问题，为后续的数学表述提供了物理基础。
在信息光学中，我们更关心的是衍射屏（如孔径、光栅）后的光场分布。设衍射屏位于 z=0 平面，其复振幅透过率为 t(x₀, y₀)，入射光场为 U₀(x₀, y₀)，则屏后紧贴处的光场为 U(x₀, y₀) = U₀(x₀, y₀)·t(x₀, y₀)。根据惠更斯-菲涅尔原理，观察平面 (x, y) 处的光场可表示为：
U(x,y,z) = (1/jλ) ∬ U(x₀,y₀) · [exp(jkr)/r] · K(θ) dx₀dy₀其中 r = √[(x-x₀)² + (y-y₀)² + z²] 是源点到观察点的距离，k = 2π/λ 为波数，K(θ) 为倾斜因子，表示次级波在不同方向上的辐射强度差异。在傍轴条件下，K(θ) ≈ 1，这一简化是后续近似的基础。
3.2 基尔霍夫衍射公式：从物理原理到数学表达
基尔霍夫利用格林定理严格求解了亥姆霍兹方程，给出了衍射场的精确积分表达式。基尔霍夫衍射公式为：
U(x,y,z) = (1/jλ) ∬ U(x₀,y₀) · [exp(jkr)/r] · [(cos(n,r) - cos(n,r₀))/2] dx₀dy₀式中 (cos(n,r) - cos(n,r₀))/2 是精确的倾斜因子，n 为衍射屏法线方向。该公式在数学上严格，但实际应用中往往过于复杂。在傍轴近似下（即观察点与源点均靠近光轴），倾斜因子趋近于 1，基尔霍夫公式退化为：
U(x,y,z) = (1/jλz) · exp(jkz) · ∬ U(x₀,y₀) · exp{jk[(x-x₀)² + (y-y₀)²]/(2z)} dx₀dy₀这一形式被称为菲涅尔衍射积分，是信息光学中最常用的传播公式。它本质上是一个卷积运算，将输入光场与一个二次相位因子（即菲涅尔核）进行卷积。
| 公式形式 | 数学表达 | 物理含义 | 计算复杂度 |
|---|---|---|---|
| 基尔霍夫衍射公式 | 含精确倾斜因子 | 严格标量解 | 高（数值积分） |
| 菲涅尔衍射积分 | 傍轴近似，二次相位 | 近场衍射 | 中（卷积/FFT） |
| 夫琅禾费衍射积分 | 相位线性化 | 远场衍射 | 低（直接FFT） |
3.3 菲涅尔近似：适用条件与工程判断
菲涅尔近似的核心是对球面波因子 exp(jkr)/r 中的 r 进行二项式展开，并保留到二次项。展开条件为：
z³ ≫ (π/4λ) · [(x-x₀)² + (y-y₀)²]²_max这一条件确保展开式中被忽略的高阶项（三次及以上）对相位的影响远小于 1 弧度。在实际工程中，更常用的判据是：
z > (D²/λ) 的 1/4 次方量级，其中 D 为孔径最大尺寸例如，对于直径 D = 10 mm 的孔径，λ = 632.8 nm（He-Ne 激光），则菲涅尔区的最小距离约为 z_min ≈ (D⁴/λ)^(1/3) ≈ 0.5 m。这意味着在 0.5 m 以外的区域，菲涅尔近似均适用。
3.4 夫琅禾费近似：远场条件与傅里叶变换
当传播距离进一步增大，使得观察平面上的二次相位因子近似为线性相位时，菲涅尔积分退化为夫琅禾费衍射积分：
U(x,y,z) = [exp(jkz)/(jλz)] · exp[jk(x²+y²)/(2z)] · ∬ U(x₀,y₀) · exp[-j2π(xx₀+yy₀)/(λz)] dx₀dy₀观察上式，积分部分正是 U(x₀,y₀) 的二维傅里叶变换在频率 (u = x/(λz), v = y/(λz)) 处的取值。这一发现是信息光学的核心——夫琅禾费衍射实现了光场的傅里叶变换，为光学信息处理提供了物理基础。
夫琅禾费近似的适用条件为：
z ≫ k(x₀² + y₀²)_max / 2即观察距离远大于孔径尺寸的平方与波长的比值。对于 D = 10 mm 的孔径和 λ = 632.8 nm，夫琅禾费距离约为 z_FF ≈ D²/λ ≈ 158 m。这个距离在实验室中往往难以直接实现，因此通常借助透镜的傅里叶变换性质，在透镜后焦面上获得夫琅禾费衍射图样。
3.5 案例拆解：菲涅尔近似在数字全息再现中的应用
数字全息是信息光学的重要应用领域，其再现过程通常采用菲涅尔衍射积分。以下是一个典型的工程案例：
背景与约束：某实验室使用 CMOS 相机（像素尺寸 3.45 μm，像素数 2048×2048）记录离轴菲涅尔全息图。记录距离 z = 150 mm，照明波长 λ = 532 nm，物体尺寸约 5 mm。
初始现象与目标：直接对全息图做傅里叶变换（即采用夫琅禾费近似）得到的再现像模糊不清，无法分辨物体细节。目标是通过正确的衍射模型获得清晰的再现像。
分析判断过程：
- 计算菲涅尔数：N_F = D²/(λz) = (5×10⁻³)² / (532×10⁻⁹ × 0.15) ≈ 313。由于 N_F ≫ 1，属于菲涅尔衍射区，夫琅禾费近似不成立。
- 检查采样条件：根据奈奎斯特采样定理，菲涅尔积分中的二次相位因子在频域的最大频率应小于相机采样频率的一半。经计算，最大空间频率为 f_max = D/(λz) ≈ 62.7 lp/mm，而相机的采样频率为 1/(2×3.45 μm) ≈ 145 lp/mm，满足采样要求。
- 选择正确的再现算法：采用卷积法（基于菲涅尔核的卷积）或角谱法进行数值再现。
实际操作步骤：
1. 读取全息图 H(x,y)，去除零级像和共轭像（频域滤波）
2. 构造菲涅尔传播核 h(x,y) = exp(jkz)/(jλz) · exp[jk(x²+y²)/(2z)]
3. 计算再现像：U_rec = IFFT{ FFT(H) × FFT(h) }
4. 对再现像进行振幅提取和相位解包裹结果对比：采用菲涅尔卷积法后，再现像的分辨率从原来的约 200 μm 提升至 15 μm，接近衍射极限。而使用夫琅禾费近似（直接 FFT）时，由于二次相位因子未被补偿，再现像严重散焦，分辨率仅约 500 μm。
风险与代价：卷积法计算量约为直接 FFT 的 3 倍（需要两次 FFT 和一次 IFFT），但现代计算机可在毫秒级完成，工程上完全可接受。主要风险在于传播核的采样不足可能导致混叠，需确保核函数的带宽满足采样定理。
复盘经验：在数字全息系统中，记录距离 z 的选择需综合考虑菲涅尔数、采样分辨率和视场大小。当 z 过小时，菲涅尔数过大，二次相位因子变化剧烈，容易产生采样混叠；当 z 过大时，分辨率下降。通常选择 z 使得菲涅尔数在 10~1000 之间，并配合适当的放大光路。
3.6 常见误区与总结复盘
在标量衍射理论的应用中，以下误区值得特别注意：
- 误区一：将菲涅尔近似与夫琅禾费近似混用。两者的适用距离相差数个数量级，错误使用会导致严重的相位误差。判断方法：计算菲涅尔数 N_F = D²/(λz)，若 N_F ≪ 1 则可用夫琅禾费近似，若 N_F ≥ 1 则必须用菲涅尔近似。
- 误区二：忽略倾斜因子的影响。在大数值孔径（NA > 0.1）系统中，倾斜因子不能简单取 1，否则会引入振幅误差。此时应使用更精确的瑞利-索末菲公式。
- 误区三：在近场区域使用标量理论。当特征尺寸接近波长量级（如亚波长光栅）时，标量衍射理论失效，必须使用矢量衍射理论（如严格耦合波分析 RCWA）。
- 误区四：混淆透镜后焦面与自由空间夫琅禾费衍射。透镜后焦面是精确的傅里叶变换面，不受远场条件限制；而自由空间夫琅禾费衍射需要满足 z ≫ D²/λ。
总结而言，标量衍射理论的三级近似构成了信息光学计算的基础框架。惠更斯-菲涅尔原理提供了物理图像，基尔霍夫公式给出了严格数学表达，而菲涅尔与夫琅禾费近似则在不同距离尺度上简化了计算。工程应用中，应根据具体系统的孔径尺寸、波长和传播距离，选择合适的近似阶次，并在必要时通过数值仿真验证近似误差。掌握这些适用条件，是进行光学系统设计、全息再现计算和空间滤波分析的前提。