第三章 范西特-泽尼克定理：从推导到应用
这一章，我们终于要啃下光学相干性理论里最硬的一块骨头——范西特-泽尼克定理。说实话，我第一次看到这个定理的完整形式时，心里也犯过嘀咕：这玩意儿到底在说什么？后来在项目中反复用了几年，才真正体会到它的精妙。
简单来说，这个定理回答了一个核心问题：一个扩展的非相干光源，在远场产生的光场，其互强度函数与光源本身的强度分布之间，存在什么关系？
答案很漂亮——它们之间是傅里叶变换关系。没错，就是那个你在信号处理里天天见的傅里叶变换。
3.1 定理的数学推导
我们先从最基础的物理图像说起。假设有一个扩展的非相干光源，分布在某个平面上。光源上每个点都在独立地发光，彼此之间没有固定的相位关系。这就是“非相干”的含义。
现在，我们在距离光源足够远的观察平面上，放两个点 P₁ 和 P₂。我想知道这两个点上的光场之间的相关性——也就是互强度 J(P₁, P₂)。
推导的起点是惠更斯-菲涅尔原理。每个光源点发出的球面波，传播到观察点，然后叠加起来。数学上可以写成：
U(P₁) = ∫∫ Σ(ξ,η) · [exp(ikr₁) / r₁] · dξ dη
U(P₂) = ∫∫ Σ(ξ,η) · [exp(ikr₂) / r₂] · dξ dη其中 Σ(ξ,η) 是光源面上的复振幅分布，r₁ 和 r₂ 分别是光源点到 P₁ 和 P₂ 的距离。
互强度的定义是：
J(P₁, P₂) = ⟨U(P₁) · U*(P₂)⟩这里的尖括号表示时间平均。把上面的表达式代进去，再交换积分和平均的顺序，你会得到：
J(P₁, P₂) = ∫∫ I(ξ,η) · [exp(ik(r₁ - r₂)) / (r₁ r₂)] · dξ dη注意，这里 I(ξ,η) 是光源的强度分布。为什么交叉项消失了？因为光源是非相干的，不同点发出的光之间没有固定的相位关系，时间平均后交叉项为零。嗯，这一步是关键，也是很多初学者容易忽略的地方。
核心结论：非相干光源的互强度，等于光源强度分布与传播相位因子的卷积积分。这个积分形式，就是范西特-泽尼克定理的雏形。
3.2 远场近似：让公式变简单
上面的公式虽然精确，但用起来太麻烦。我们需要做近似。我个人习惯把远场近似拆成两步来理解：
- 振幅近似：在分母中，r₁ ≈ r₂ ≈ z（观察距离）。这个近似对振幅影响不大，因为分母的变化相对缓慢。
- 相位近似：在指数项中，r₁ - r₂ 必须精确到一阶项。因为相位对距离变化极其敏感，差半个波长就会导致完全不同的结果。
具体来说，在远场条件下，我们可以把 r₁ 和 r₂ 展开为：
r₁ ≈ z + [(ξ - x₁)² + (η - y₁)²] / (2z)
r₂ ≈ z + [(ξ - x₂)² + (η - y₂)²] / (2z)相减后，常数项 z 消掉了，得到：
r₁ - r₂ ≈ [ (x₂² + y₂²) - (x₁² + y₁²) ] / (2z) - [ (x₂ - x₁)ξ + (y₂ - y₁)η ] / z把这一项代回互强度公式，整理后得到：
J(P₁, P₂) = [exp(iψ) / (λ² z²)] · ∫∫ I(ξ,η) · exp[-i(2π/λz)(Δx·ξ + Δy·η)] · dξ dη其中 Δx = x₂ - x₁，Δy = y₂ - y₁，ψ 是一个与观察点位置有关的二次相位因子。
我的经验：在实际工程中，如果你只关心互强度的模（即相干度的绝对值），那个二次相位因子 ψ 可以忽略。但如果你需要完整的复相干度，ψ 必须保留。我曾经在一个干涉测量项目中忽略了它，结果相位反演全错了，排查了两天才找到原因。
3.3 空间相干性的传播
范西特-泽尼克定理告诉我们一个非常重要的物理图像：非相干光源在远场产生的光场，其空间相干性由光源的强度分布通过傅里叶变换决定。
换句话说，光源的尺寸越小，远场相干区域越大；光源的尺寸越大，远场相干区域越小。这就像是一个“光学不确定性原理”——你无法同时让光源又小又亮，同时又让远场光斑完全相干。
我们定义一个重要的量——复相干度 μ(P₁, P₂)：
μ(P₁, P₂) = J(P₁, P₂) / √[ I(P₁) · I(P₂) ]在远场近似下，复相干度就是光源强度分布的归一化傅里叶变换：
μ(Δx, Δy) = [∫∫ I(ξ,η) · exp(-i(2π/λz)(Δx·ξ + Δy·η)) dξ dη] / ∫∫ I(ξ,η) dξ dη这个形式是不是很眼熟？没错，它和夫琅禾费衍射的公式几乎一模一样。区别在于，衍射描述的是振幅的传播，而范西特-泽尼克定理描述的是相干性的传播。
| 物理量 | 夫琅禾费衍射 | 范西特-泽尼克定理 |
|---|---|---|
| 输入 | 孔径的复振幅透过率 | 光源的强度分布 |
| 输出 | 远场衍射光强 | 远场复相干度 |
| 数学关系 | 傅里叶变换 | 归一化傅里叶变换 |
| 物理意义 | 光能量的重新分布 | 光场相关性的空间分布 |
避坑指南：我曾经在分析一个LED阵列的相干性时，直接套用了夫琅禾费衍射公式来计算相干度，结果发现实验数据对不上。后来才意识到，夫琅禾费衍射用的是复振幅透过率，而范西特-泽尼克定理用的是强度分布。两者虽然数学形式相似，但物理输入完全不同。
3.4 远场近似条件
说了这么多，到底什么条件下才能用远场近似？我给大家一个实用的判断标准：
- 观察距离 z 必须远大于光源的最大线度 D：通常要求 z > 10D，这是最基本的条件。
- 观察距离 z 必须远大于 λ：这个一般都能满足，可见光波长才几百纳米。
- 更严格的条件：z > D²/λ。这个条件来自菲涅尔数，确保二次相位因子可以线性化。
举个例子，一个直径 1mm 的光源，用 500nm 的光照明，远场条件要求 z > (1mm)² / 500nm = 2m。也就是说，观察距离至少要 2 米以上，才能用范西特-泽尼克定理的简单形式。
如果观察距离不够远怎么办？那就得用更一般的菲涅尔衍射形式的互强度传播公式，计算会复杂很多。我在实际项目中通常先用远场近似做快速估算，再用数值方法验证关键参数。
3.5 一个简单的例子
假设光源是一个均匀发光的圆形孔径，半径为 a，强度分布为：
I(ρ) = I₀, ρ ≤ a
I(ρ) = 0, ρ > a其中 ρ = √(ξ² + η²)。那么远场复相干度就是圆域函数的傅里叶变换，结果是一阶贝塞尔函数：
μ(Δr) = 2 J₁(k a Δr / z) / (k a Δr / z)其中 Δr = √(Δx² + Δy²)，J₁ 是一阶贝塞尔函数。
这个结果和圆孔夫琅禾费衍射的振幅分布形式完全一样。但注意，这里描述的是相干度，不是光强。相干度第一次降到零的位置对应：
Δr₀ = 1.22 λ z / (2a)这个 Δr₀ 就是相干长度，它决定了在观察平面上多大范围内光场是部分相干的。
实用技巧：如果你需要快速估算一个扩展光源的相干面积，直接用这个公式。我在设计空间光通信系统时，经常用它来估算接收孔径的尺寸上限——超过相干长度的孔径，接收到的光场相干性会显著下降，影响干涉效率。
好了，范西特-泽尼克定理的核心内容就这些。从推导到近似条件，再到物理图像和实际应用，我希望你能感受到这个定理的美妙之处。它把看似复杂的相干性问题，简化成了一个傅里叶变换——而傅里叶变换，恰恰是我们光学工程师最熟悉的工具之一。