3. 互强度传播定律：自由空间中的演化

这一章，我们来啃互强度理论里最核心的一块——自由空间传播。说白了，就是光从一个平面跑到另一个平面，它的互强度怎么变。当年我推导这个公式的时候，可是踩了不少坑，今天一并讲给你们听。

3.1 从衍射积分到互强度传播

先回忆一下标量衍射理论。一个光场 \( U(x,y) \) 从平面 \( z=0 \) 传播到 \( z=z_0 \)，可以用菲涅尔-基尔霍夫衍射积分描述：

U(x',y') = ∬ U(x,y) · h(x'-x, y'-y) dx dy

其中 \( h \) 是自由空间的脉冲响应。这个大家应该都熟。

那互强度呢？互强度定义是 \( J(x_1,x_2) = \langle U(x_1) U^*(x_2) \rangle \)。把衍射积分代进去，你会得到：

J(x1',x2') = ∬∬ J(x1,x2) · h(x1'-x1) · h*(x2'-x2) dx1 dx2

嗯，这就是互强度传播的基本公式。看起来就是衍射积分的“双份版”——一个场用 \( h \)，另一个场用 \( h^* \)。

3.2 我当年踩过的坑

第一次推导这个公式时，我犯了个低级错误——直接把衍射积分平方了。心想反正互强度是场的二阶矩，平方不就完了？

错！互强度是两个不同点的场乘积的平均，不是同一点。平方得到的是强度，不是互强度。

第二个坑：积分顺序。我当时图省事，先做系综平均再做积分。但系综平均和积分交换顺序是有条件的——场必须是广义平稳的。实际光源大多满足，但理论上你得先验证。

3.3 自由空间传播的简化形式

如果光源是准单色的，且传播距离满足菲涅尔近似，互强度传播可以大大简化。我直接给结果：

J(x1',x2') = (1/λ²z²) · exp[jk(x1'² - x2'²)/2z]
             · ∬∬ J(x1,x2) · exp[-jk(x1'x1 - x2'x2)/z] dx1 dx2

这个形式看着复杂，但有个漂亮的性质——它本质上是一个四维傅里叶变换。把互强度从空间域变换到空间频率域，传播就变成了简单的相位因子相乘。

传播条件	近似形式	适用场景
近场（菲涅尔）	四维傅里叶变换 + 二次相位因子	显微成像、近场光学
远场（夫琅禾费）	纯四维傅里叶变换	望远镜、远场衍射
完全非相干	退化为强度传播（照度定律）	普通照明系统

3.4 互强度传播的物理图像

你想想看，互强度传播到底在描述什么？

它描述的是：两个点之间的相干性，如何随着传播距离变化。近场时，两个点可能高度相干；跑远了，相干性逐渐丧失。这就是著名的范西特-泽尼克定理的前奏。

我做个项目时遇到过这种情况：设计一个激光照明系统，近场条纹清晰，远场却模糊了。当时以为是光学元件问题，查了半天。后来用互强度传播一算，发现是相干长度随传播距离增大而减小，导致干涉条纹对比度下降。嗯，这就是互强度理论的威力——它告诉你“为什么”会这样。

3.5 知识结构图

下面这张图总结了互强度传播的核心逻辑：

3.6 工程实用建议

最后，给几个实际干活时的建议：

• 先判断光源类型：激光还是LED？相干长度多少？这决定了你用哪个近似。
• 用数值积分时注意采样：互强度是四维函数，采样点稍微多一点，内存就爆了。我一般用稀疏采样加插值。
• 对称性别浪费：如果系统是旋转对称的，互强度可以降维成径向坐标的函数，计算量直接少一个数量级。
• 验证用边缘情况：把互强度退化成强度（令 x₁=x₂），看传播结果是否和强度传播一致。这是最直接的 sanity check。