二、衍射光学基础：惠更斯-菲涅尔原理、标量衍射理论、基尔霍夫衍射公式、菲涅尔与夫琅禾费近似

衍射光学元件（DOE）的设计本质上是求解一个逆衍射问题：给定期望的输出光场分布，反推元件表面的微结构相位分布。要完成这一过程，必须首先建立正向衍射计算的数学框架。本章将沿着光学衍射理论的历史发展脉络，从惠更斯-菲涅尔原理的物理图像出发，逐步过渡到标量衍射理论的数学表述，重点解析基尔霍夫衍射公式的物理意义与适用边界，并系统梳理菲涅尔近似与夫琅禾费近似的判据、适用范围及在DOE设计中的典型应用场景。这一基础是否牢固，直接决定了后续相位恢复算法（如GS算法）的迭代精度与设计结果的可靠性。

上图给出了衍射光学理论体系的整体框架与工程判断路径。从惠更斯-菲涅尔原理的物理图像出发，经由基尔霍夫公式的严格数学化，再到标量衍射理论的简化处理，最终根据传播距离与元件特征尺寸的相对关系，选择菲涅尔近似或夫琅禾费近似作为实际计算工具。这张图是本章内容的地图，后续所有公式推导与案例讨论都将围绕这一框架展开。

2.1 核心原理：从次波叠加到衍射积分

惠更斯-菲涅尔原理是衍射光学最根本的物理基础。该原理指出：波前上的每一点都可以看作一个新的次级球面波源，而空间中任意一点的振动是所有这些次波在该点相干叠加的结果。这一思想将光的传播问题转化为积分问题，为衍射计算提供了物理图像。然而，原始的惠更斯原理只解决了波的传播方向问题，并未涉及振幅与相位分布；菲涅尔补充了次波相干叠加的思想，引入了倾斜因子和相位延迟，使得衍射计算成为可能。

基尔霍夫在此基础上给出了严格的数学表述。他利用格林定理求解亥姆霍兹方程，推导出衍射场中任意一点P的复振幅表达式：

U(P) = (1/jλ) ∬Σ U(Q) · [exp(jkr)/r] · K(θ) dΣ

其中，U(Q)是衍射屏上Q点的复振幅分布，r是Q到P的距离，k=2π/λ是波数，K(θ)是倾斜因子，θ是衍射角。基尔霍夫公式的贡献在于：它明确给出了倾斜因子的具体形式K(θ) = (cosθ₀ + cosθ)/2，并规定了边界条件。但基尔霍夫公式在数学上存在自洽性问题——它在边界上的边界条件假设与物理实际不完全一致，后来由瑞利-索末菲理论进行了修正。对于DOE设计而言，基尔霍夫公式的精度已经足够，因为DOE的特征尺寸通常远大于波长，边界效应的影响可以忽略。

标量衍射理论的核心假设是：将光场视为标量场，忽略其矢量偏振特性。这一假设在以下条件下成立：衍射结构的特征尺寸远大于波长（通常要求大于10λ），且观察点离衍射屏的距离远大于波长。对于大多数DOE应用——如分束器、光束整形器、扩散器——这些条件均能满足。但当DOE的特征尺寸接近波长（如亚波长光栅）或涉及偏振敏感应用时，必须采用严格的矢量衍射理论（如RCWA、FDTD），标量理论将不再适用。

2.2 判断方法：菲涅尔数与近似选择流程

在实际DOE设计中，选择菲涅尔近似还是夫琅禾费近似，直接决定了计算复杂度和算法实现路径。判断的核心指标是菲涅尔数（Fresnel Number），其定义为：

F = a² / (λz)

其中，a是衍射孔径或DOE有效区域的半宽度，z是衍射屏到观察面的距离，λ是波长。菲涅尔数的物理意义是：孔径边缘与中心到观察点的光程差相对于波长的比值。当F远大于1时，光程差显著，必须使用菲涅尔近似；当F远小于1时，光程差可忽略，夫琅禾费近似成立。

具体判断流程如下：

1. 确定参数：明确DOE的口径D（或半宽a）、工作波长λ、目标距离z。
2. 计算菲涅尔数：F = a²/(λz)。
3. 判断近似类型：若F < 0.1，采用夫琅禾费近似；若0.1 ≤ F ≤ 10，采用菲涅尔近似；若F > 10，需考虑更严格的衍射积分或矢量方法。
4. 验证傍轴条件：无论哪种近似，都要求观察点靠近光轴，即横向偏移远小于z。

夫琅禾费近似的数学本质是：将衍射积分中的二次相位因子近似为1，使得积分变为标准的傅里叶变换形式。这意味着观察面上的光场分布是衍射屏透射函数的傅里叶变换。这一性质使得夫琅禾费近似下的DOE设计可以借助快速傅里叶变换（FFT）高效实现，是远场分束器、结构光投影器等元件设计的标准工具。

菲涅尔近似则保留了二次相位因子，其积分形式为：

U(x,y,z) = [exp(jkz)/(jλz)] · exp[jk(x²+y²)/(2z)] · ∬ U(ξ,η) · exp[jk(ξ²+η²)/(2z)] · exp[-j2π(xξ+yη)/(λz)] dξdη

该式可以看作一个二次相位因子与傅里叶变换的组合，在数值实现上通常采用单次FFT加相位相乘的方式完成。菲涅尔近似适用于近场DOE设计，如光束整形、焦深扩展、激光加工中的平顶光束生成等场景。

对比维度	菲涅尔近似	夫琅禾费近似
适用距离	近场（z ~ a²/λ）	远场（z >> a²/λ）
菲涅尔数范围	0.1 ~ 10	< 0.1
相位因子	保留二次相位	忽略二次相位
数学形式	菲涅尔积分（FFT+相位相乘）	傅里叶变换（纯FFT）
典型应用	光束整形、近场成像、焦深扩展	分束器、结构光、远场扩散
数值实现	单次FFT + 相位因子	单次FFT
误差来源	傍轴近似、采样率不足	远场条件不满足、孔径衍射

2.3 案例拆解：激光分束器的近似选择与仿真验证

以一个典型的1×5激光分束器设计为例，展示近似选择对设计结果的影响。设计目标：将波长为532nm的绿光激光束分束为5×5的均匀点阵，分束角为±5°，DOE口径为10mm×10mm。

背景与约束：该分束器用于激光雷达中的扫描照明模块，要求分束效率大于85%，均匀性优于±5%。DOE基底材料为熔融石英（折射率n=1.46），采用二元光学工艺刻蚀。

关键参数计算：分束角θ=5°对应空间频率f = sinθ/λ ≈ 0.0872/532nm ≈ 164 lp/mm。DOE口径10mm，则半宽a=5mm。根据夫琅禾费衍射条件，观察距离z需满足：z >> a²/λ = (5×10⁻³)² / 532×10⁻⁹ ≈ 47m。这意味着若要在1m处观察分束效果，菲涅尔数F = a²/(λz) = 25×10⁻⁶ / (532×10⁻⁹ × 1) ≈ 47，远大于0.1，夫琅禾费近似完全不成立。

分析判断：然而，分束器的设计通常是在傅里叶平面上进行的——即在透镜的后焦面上观察分束图案。此时，透镜将无穷远处的夫琅禾费衍射转换到焦平面上，等效距离为透镜焦距f。若选用f=100mm的傅里叶透镜，则等效菲涅尔数F = a²/(λf) = 25×10⁻⁶ / (532×10⁻⁹ × 0.1) ≈ 470，仍然很大。但这里的关键在于：透镜的傅里叶变换性质是精确的，不依赖于菲涅尔数的大小。因此，在透镜后焦面上，夫琅禾费近似严格成立，DOE的相位分布可以直接通过目标光场的逆傅里叶变换获得。

实际操作步骤：

1. 设定目标光场为5×5均匀点阵，每个点的振幅相等，相位随机初始化。
2. 对目标光场做逆傅里叶变换，得到DOE平面的复振幅分布。
3. 提取相位分布，量化到8个台阶（二元光学工艺）。
4. 使用菲涅尔衍射积分（而非夫琅禾费近似）对量化后的DOE进行正向仿真，验证在透镜焦平面上的实际分束效果。
5. 若均匀性不满足要求，迭代优化相位分布（如采用GS算法或梯度下降法）。

结果对比：采用夫琅禾费近似设计得到的相位分布，在透镜焦平面上的仿真结果与目标光场的偏差小于2%；而若直接在1m距离处用菲涅尔积分仿真，则由于菲涅尔数过大，分束图案会出现明显的散焦和畸变。这说明：夫琅禾费近似在透镜傅里叶变换系统中是精确的，但若试图在自由空间远场直接应用，则必须严格满足远场条件。

风险与代价：该设计方案的代价是需要额外的傅里叶透镜，增加了系统体积和成本。若希望取消透镜，直接在自由空间实现分束，则必须采用菲涅尔近似设计，此时DOE的相位分布不再是简单的傅里叶变换关系，而需要通过迭代优化算法求解，计算复杂度显著增加。

2.4 常见误区与总结复盘

在衍射光学基础理论的学习与应用中，以下误区最为常见，需要特别警惕：

• 误区一：混淆“远场”与“夫琅禾费近似”。远场是一个物理概念，指距离足够远使得光场成为以观察点为中心的球面波；夫琅禾费近似是数学上的简化，要求二次相位因子可以忽略。在透镜焦平面上，即使物理距离很近，夫琅禾费近似也严格成立。
• 误区二：忽视标量理论的适用边界。当DOE的最小特征尺寸接近波长时，标量理论预测的衍射效率会与实际值产生显著偏差。例如，对于特征尺寸为2λ的二元光栅，标量理论预测的+1级衍射效率为40.5%，而严格矢量计算可能只有30%左右。
• 误区三：采样率不足导致混叠。在数值实现菲涅尔衍射积分时，采样间隔必须满足奈奎斯特条件。对于DOE口径D和距离z，采样间隔Δx应满足Δx ≤ λz/D，否则高频成分会发生混叠，导致仿真结果失真。
• 误区四：忽略傍轴条件的限制。菲涅尔近似和夫琅禾费近似都要求衍射角较小（通常小于10°~15°）。对于大角度分束（如分束角超过20°），傍轴近似不再成立，需要采用非傍轴衍射模型或矢量理论。

总结本章内容，衍射光学设计的理论基础可以概括为一条主线：从惠更斯-菲涅尔原理的物理图像出发，经由基尔霍夫公式的严格数学化，在标量衍射理论的框架下，根据菲涅尔数选择合适的近似模型。夫琅禾费近似对应傅里叶变换关系，适合远场分束与结构光投影；菲涅尔近似保留二次相位因子，适合近场光束整形与成像应用。在实际工程中，设计者应始终带着两个问题审视自己的计算模型：一是标量近似是否成立？二是菲涅尔数是否支持所选的近似类型？只有这两个问题都得到肯定回答，后续的相位恢复算法才能建立在可靠的基础上。