第二章 互强度函数的数学表示:从相干性到范西特-泽尼克定理

这一章,我们得把互强度函数这玩意儿从物理直觉拉到数学台面上来。说实话,我早年做干涉测量时,总觉得相干性是个玄学概念——直到我亲手推导了一遍互相干函数,才明白它其实就是光场涨落之间的“默契程度”。

2.1 互相干函数的定义与推导

先看一个最简单的场景:空间两点 \( P_1 \) 和 \( P_2 \),时间差 \( \tau \)。这两点的光场分别是 \( u(P_1, t) \) 和 \( u(P_2, t+\tau) \)。

互相干函数定义为:

Γ(P₁, P₂, τ) = ⟨u(P₁, t) · u*(P₂, t+τ)⟩

尖括号表示时间平均。说白了,就是看两个点的光场在时间偏移 τ 后,相关性有多强。

我在项目中遇到过一个问题:用迈克尔逊干涉仪测光源相干长度,结果条纹对比度一直上不去。后来发现是互相干函数的实部计算时忘了归一化——嗯,细节决定成败。

归一化后得到复相干度:

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

它的模长在 0 到 1 之间。0 表示完全不相关,1 表示完全相干。你想想看,这其实就是光场之间的“相关系数”。

核心要点: 互相干函数同时包含了空间相干性和时间相干性的信息。空间相干性看 P₁ 和 P₂ 的间距,时间相干性看 τ 的大小。

2.2 空间相干性与时间相干性

这两个概念经常被混为一谈,但它们的物理根源完全不同。

特性 空间相干性 时间相干性
物理来源 光源的横向尺寸 光源的频谱宽度
关键参数 相干面积 \( A_c \) 相干长度 \( L_c \)
数学表达 \( \Gamma(P_1, P_2, 0) \) \( \Gamma(P, P, \tau) \)
典型实验 杨氏双缝 迈克尔逊干涉仪

为什么会这样?空间相干性描述的是同一时刻不同位置的光场相关性;时间相干性描述的是同一位置不同时刻的光场相关性。两者在互相干函数中通过不同的自变量体现。

我的习惯: 做系统设计时,先判断光源是扩展光源还是窄带光源。扩展光源优先考虑空间相干性,窄带光源优先考虑时间相干性。这个判断能省掉一半的调试时间。

2.3 范西特-泽尼克定理

这个定理是空间相干性分析的基石。它告诉我们:对于一个非相干扩展光源,其远场观察面上的互相干函数,等于光源强度分布的傅里叶变换。

数学形式:

Γ(P₁, P₂, 0) = ∫∫ I(ξ, η) · exp[-j·(2π/λ·R)·(Δx·ξ + Δy·η)] dξ dη

其中 \( (\xi, \eta) \) 是光源平面坐标,\( (\Delta x, \Delta y) \) 是观察点间距,\( R \) 是距离。

我曾经用这个定理设计一个天文干涉仪的基线长度。当时需要估算两颗星之间的相干性,直接用范西特-泽尼克定理反推光源角直径——结果和实测数据吻合得很好。那一刻,我觉得傅里叶变换真是上帝给光学工程师的礼物。

注意: 范西特-泽尼克定理成立的前提是:光源是非相干的,且观察面在远场(夫琅禾费衍射区)。近场情况下,这个定理不成立,需要用更一般的瑞利-索末菲衍射公式。

2.4 知识结构图

下面这张图梳理了本章的核心逻辑:从互相干函数出发,分支出空间相干性和时间相干性,再通过范西特-泽尼克定理将空间相干性与光源强度分布联系起来。

互强度函数与相干性知识结构图 互相干函数 Γ(P₁,P₂,τ) 空间相干性 Γ(P₁,P₂,0) 杨氏双缝实验 时间相干性 Γ(P,P,τ) 迈克尔逊干涉仪 范西特-泽尼克定理 Γ = ℱ{ 光源强度分布 } 核心逻辑 互相干函数 → 空间相干性(τ=0) + 时间相干性(P₁=P₂) → 范西特-泽尼克定理连接光源与观察面

2.5 实战中的避坑指南

最后,分享几个我在项目中踩过的坑:

  • 归一化顺序: 先计算互相干函数,再归一化。不要先归一化再计算,否则会丢失幅度信息。
  • 远场条件: 使用范西特-泽尼克定理前,务必确认观察距离 \( R \) 满足 \( R \gg a^2/\lambda \),其中 \( a \) 是光源尺寸。我见过有人用近场数据套远场公式,结果完全对不上。
  • 时间平均窗口: 实际测量中,时间平均窗口要远大于光源的相干时间,否则测出来的互相干函数会抖动。我曾经因为探测器积分时间设得太短,导致条纹对比度忽高忽低。
  • 偏振匹配: 互相干函数要求两个光场的偏振方向一致。如果偏振态不匹配,即使空间和时间完全相干,条纹对比度也会下降。
我的建议: 刚开始做相干性分析时,先用 MATLAB 或 Python 写一个简单的模拟程序,把互相干函数、复相干度、条纹对比度都算一遍。数值实验能帮你建立直觉,比死磕公式高效得多。