4. 部分相干光的传输理论：广义惠更斯-菲涅尔原理、交叉谱密度函数的传输方程、ABCD光学系统的传输矩阵
好，我们进入第四章。这一章是部分相干光理论从“描述”走向“应用”的关键桥梁。说白了，前面我们学会了怎么描述一束部分相干光——用交叉谱密度函数（CSD）。但光是要传播的，它经过透镜、自由空间、光阑之后，CSD 会变成什么样？这就是本章要解决的问题。
我个人习惯把这一章看作“工具箱”。你手里有 CSD 这个“原材料”，本章教你用三种工具去加工它：广义惠更斯-菲涅尔原理（最根本的物理图像）、CSD 传输方程（最直接的数学工具）、ABCD 矩阵法（最实用的工程手段）。三者本质相通，但适用场景不同。
4.1 从经典惠更斯原理说起
还记得大学物理里的惠更斯-菲涅尔原理吗？波前上的每一点都可以看作一个新的球面波源，这些次级波相干叠加，就形成了后面的波场。对于完全相干光，这个原理用起来很漂亮。
但部分相干光呢？它的波前上各点不是完全相干的。你不能简单地把所有次级波振幅加起来——你得考虑它们之间的统计相关性。
嗯，这里就是关键。广义惠更斯-菲涅尔原理的核心思想是：把 CSD 函数看作“统计相关的次级波源对”的传输。你想想看，在输入平面上，任意两点 (ρ₁, ρ₂) 的 CSD 是 W(ρ₁, ρ₂)。这两点各自发出球面波，传播到输出平面上的两点 (r₁, r₂)。由于 ρ₁ 和 ρ₂ 是部分相干的，它们的贡献不能直接加振幅，而是加互强度。
数学上，这个原理可以写成：
W(r₁, r₂) = ∫∫ W(ρ₁, ρ₂) · K*(r₁, ρ₁) · K(r₂, ρ₂) d²ρ₁ d²ρ₂其中 K(r, ρ) 是自由空间的传播核函数。对于菲涅尔衍射区，它长这样：
K(r, ρ) = (1 / iλz) · exp[ iπ/(λz) · |r - ρ|² ]我在项目中遇到过一个问题：直接用这个四重积分去算，计算量非常大。对于 1024×1024 的网格，那就是 1024⁴ 次运算，任何计算机都扛不住。所以实际工程中，我们很少直接算这个积分——后面会讲怎么简化。
核心要点：广义惠更斯-菲涅尔原理告诉我们，部分相干光的传输不是“场”的传输，而是“相关函数”的传输。输入平面上每一对点，对应输出平面上每一对点。这种“点对点”的映射关系，是理解一切传输问题的基础。
4.2 交叉谱密度函数的传输方程
有了广义惠更斯-菲涅尔原理，我们可以推导出 CSD 在自由空间中的传输方程。这个方程更紧凑，也更适合数值计算。
对于部分相干光在近轴近似下的传输，CSD 满足以下方程：
∂W/∂z = (iλ / 4π) · (∇₁² - ∇₂²) W其中 ∇₁² 和 ∇₂² 分别是对 r₁ 和 r₂ 的横向拉普拉斯算子。这个方程看起来像薛定谔方程，但注意符号——它是两个拉普拉斯算子的差，而不是和。
为什么会这样？我个人的理解是：CSD 描述的是两个点的相关性，每个点各自独立地衍射。r₁ 点的衍射用 ∇₁² 描述，r₂ 点的衍射用 ∇₂² 描述，但因为是互强度，所以是差而不是和。
这个方程有几个重要的特解：
- 高斯-谢尔模型光束：传输后仍然是高斯-谢尔模型，只是束腰和相干长度会变化
- 均匀相干光：如果初始 CSD 是空间不变的，传输后仍然保持空间不变性
- 完全相干极限：当相干长度趋于无穷时，方程退化为标准的波动方程
实用技巧：在实际仿真中，我建议用“角谱法”来解这个方程。把 CSD 变换到空间频率域，传输就变成了简单的相位相乘。这样计算复杂度从 O(N⁴) 降到 O(N² log N)，快了好几个数量级。
4.3 ABCD 光学系统的传输矩阵
好了，前面讲的是自由空间。但实际光学系统有透镜、有反射镜、有光阑。怎么处理？
这里我要引入一个非常强大的工具——ABCD 矩阵法。如果你做过激光谐振腔设计，对这个应该不陌生。对于完全相干光，ABCD 矩阵描述的是光线从输入面到输出面的几何变换。对于部分相干光，它描述的是 CSD 的变换。
具体来说，对于一个由 ABCD 矩阵描述的光学系统：
[ r' ] [ A B ] [ r ]
[ θ' ] = [ C D ] [ θ ]CSD 的传输公式可以写成：
W(r₁', r₂') = (1 / |B|²) · ∫∫ W(ρ₁, ρ₂) · exp{ iπ/(λB) · [ A|ρ₁|² - 2ρ₁·r₁' + D|r₁'|² - (A|ρ₂|² - 2ρ₂·r₂' + D|r₂'|²) ] } d²ρ₁ d²ρ₂这个公式看起来很吓人，但别怕。我来拆解一下它的物理含义：
- A 项：描述输入平面的二次相位因子（相当于透镜效应）
- B 项：描述自由空间的菲涅尔衍射
- D 项：描述输出平面的二次相位因子
- 1/|B|²：能量守恒的归一化因子
我曾经在一个激光照明系统项目中踩过坑。当时我用 ABCD 矩阵算出来的光强分布和实验对不上，查了两天才发现——我忘了考虑光阑的截断效应。ABCD 矩阵法只适用于无光阑的“理想”系统，一旦有硬边光阑，就必须用广义惠更斯-菲涅尔原理或者数值仿真。
避坑指南：ABCD 矩阵法有两个前提条件——(1) 系统是近轴、无像差的；(2) 没有硬边光阑或孔径。如果系统中有光阑，CSD 会被截断，产生衍射环，这时候 ABCD 公式就不准了。我曾经吃过这个亏，后来改用“扩展的 ABCD 矩阵法”加上数值积分才解决。
4.4 三种方法的对比与选择
讲了三种方法，你可能会问：实际工作中到底用哪个？我根据经验整理了一个对比表：
| 方法 | 适用场景 | 计算复杂度 | 精度 | 我的推荐 |
|---|---|---|---|---|
| 广义惠更斯-菲涅尔原理 | 任意系统（含光阑、像差） | O(N⁴) | 最高（严格解） | 理论分析、验证基准 |
| CSD 传输方程 | 自由空间、渐变折射率介质 | O(N² log N) | 高（近轴近似下精确） | 数值仿真首选 |
| ABCD 矩阵法 | 无光阑的理想系统 | O(N²) | 中等（忽略衍射细节） | 快速估算、系统设计 |
我个人习惯是：先用 ABCD 矩阵法做快速估算，确定系统参数的大致范围；再用 CSD 传输方程做精确仿真，验证设计是否满足要求；如果遇到光阑或复杂像差，最后用广义惠更斯-菲涅尔原理做严格验证。三步走，既快又稳。
4.5 知识结构总览
下面这张图总结了本章的核心逻辑。你可以把它当作一个“导航图”，以后遇到部分相干光传输问题时，先看这张图，就知道该走哪条路。
这张图把本章的三个核心方法放在一起对比。你可以看到，它们共享同一个输入和输出——都是 CSD 函数。区别在于中间的“传输算子”不同。选择哪个方法，取决于你的系统复杂度和精度要求。