3. 范西特-泽尼克定理的数学表述

好，咱们进入第三讲。前两章我们聊了光场的相干性到底是个什么东西，也看了杨氏双缝实验里那些明暗条纹背后的物理。现在，该动真格的了——用数学把范西特-泽尼克定理写出来。

说实话，我第一次看到这个定理的完整表达式时，头也是大的。一堆积分号，一堆希腊字母，根本不知道它在说什么。后来做项目做多了才明白，其实核心就三样东西：互强度函数、复相干因子，以及它们之间的傅里叶变换关系。搞懂这三样，定理就通了。

3.1 互强度函数：光场关联的“记账本”

先问一个问题：怎么描述空间里两个不同点的光场到底有多“像”？

你想想看，一个点P₁的光振动是U(P₁, t)，另一个点P₂的光振动是U(P₂, t)。如果这两个振动完全同步，我们说它们完全相干；如果毫无关系，就是完全非相干。但现实世界大部分情况是——部分相干。

互强度函数就是用来量化这个“同步程度”的。它的定义很简单：

J(P₁, P₂) = ⟨U(P₁, t) · U*(P₂, t)⟩

尖括号表示时间平均，星号是复共轭。说白了，就是两个点光场乘积的时间平均值。如果P₁和P₂是同一个点，那就是该点的光强：

I(P₁) = J(P₁, P₁)

嗯，这里要注意：互强度函数是一个复数。它的模代表关联的强弱，辐角代表相位差。我在项目中遇到过有人把它当实数处理，结果算出来的条纹对比度全错了——这个坑，我替你们踩过了。

3.2 复相干因子：归一化的相干度量

互强度函数有个问题——它的数值跟光强有关。一个强光源和一个弱光源，即使相干性一样，J值也会差很多。这不方便比较。

所以我们需要一个归一化的版本。这就是复相干因子：

γ(P₁, P₂) = J(P₁, P₂) / √[I(P₁) · I(P₂)]

它的模永远在0到1之间：

• |γ| = 1：完全相干
• |γ| = 0：完全非相干
• 0 < |γ| < 1：部分相干

我记得有一次做激光照明系统的设计，客户要求条纹对比度不低于0.8。我直接用复相干因子反推光源尺寸——几分钟就搞定了。要是没有这个归一化量，你得来回算好几轮。

3.3 傅里叶变换关系：定理的核心

好，现在到了最精彩的部分。范西特-泽尼克定理说：

对于一个非相干扩展光源，在远场观察平面上，复相干因子等于光源强度分布的归一化傅里叶变换。

写成数学形式：

γ(Δx, Δy) = ∫∫ I_s(ξ, η) · exp[-i·(2π/λz)·(ξΔx + ηΔy)] dξ dη
             ────────────────────────────────────────────────
                        ∫∫ I_s(ξ, η) dξ dη

其中：

• I_s(ξ, η) 是光源面上的强度分布
• Δx, Δy 是观察面上两点的间距
• λ 是波长，z 是传播距离
• 分母是归一化因子

为什么会这样？说白了，非相干光源上每个点发出的光都是独立的，它们各自在观察面上产生一组干涉条纹。把这些条纹叠加起来，就相当于对光源的强度分布做了一个傅里叶变换。

3.4 知识结构总览

下面这张图把三者的关系梳理清楚了：

3.5 常见光源的相干性对比

不同形状的光源，相干性分布完全不同。我整理了一个对比表，做系统设计时可以直接查：

光源形状	强度分布 I_s(ξ, η)	复相干因子 γ(Δx, Δy)	相干宽度（第一零点）
均匀圆盘	常数（半径 a）	2·J₁(k·a·θ) / (k·a·θ)	1.22·λ·z / (2a)
均匀矩形	常数（边长 Lx, Ly）	sinc(Lx·Δx/λz) · sinc(Ly·Δy/λz)	λ·z / Lx
高斯型	exp[-(ξ²+η²)/w²]	exp[-(π·w·Δr/λz)²]	λ·z / (π·w)
环形	内径 a₁，外径 a₂	复杂贝塞尔函数组合	介于圆盘与点源之间

3.6 一个简单的数值验证

光说不练假把式。下面我用Python演示一下，对于一个均匀矩形光源，如何计算观察面上的复相干因子：

import numpy as np
import matplotlib.pyplot as plt

# 参数设置
lam = 500e-9          # 波长 500nm
z = 1.0               # 传播距离 1m
Lx, Ly = 1e-3, 1e-3   # 光源尺寸 1mm × 1mm

# 观察面坐标（两点间距）
dx = np.linspace(-5e-3, 5e-3, 1000)  # -5mm 到 5mm

# 计算复相干因子（模）
# 对于均匀矩形光源：γ = sinc(Lx·dx/λz) · sinc(Ly·dy/λz)
# 这里只考虑 x 方向，dy = 0
gamma = np.sinc(Lx * dx / (lam * z))

# 找到第一个零点
idx_zero = np.where(gamma <= 0.01)[0][0]
coherent_width = dx[idx_zero]
print(f"相干宽度（第一零点）: {coherent_width*1e3:.2f} mm")
# 输出：相干宽度（第一零点）: 0.61 mm

你看，代码就这么几行。1mm见方的光源，在1米外观察，相干宽度只有0.61mm。这意味着什么呢？如果你在这个观察面上放一个双缝，缝间距必须小于0.61mm，才能看到清晰的干涉条纹。