很多人学傅里叶光学时,最痛苦的点在于:公式能背下来,积分也能算,但完全不知道它在描述什么物理图景。我当年也是这样——夫琅禾费衍射积分写了无数遍,知道矩孔的衍射花样是sinc函数,但这只是“会做题”。直到某天我把衍射积分和傅里叶变换在纸上对着看了一遍,才发现光学里最值钱的东西根本不是那套新物理,而是换了一种语言重新描述光场。这篇《傅里叶光学(1)》就是把第一课里最绕的几个概念,用尽量朴素的方式拆开:空间频率到底是什么,夫琅禾费衍射为什么就是傅里叶变换,透镜凭什么能做傅里叶变换,以及4F系统和空间滤波这个经典应用怎么从频域去理解。适正在学光学的学生、刚入行的光学工程师,以及做计算成像和机器视觉、想搞清楚镜头和成像本质的人。
1. 光场为什么要“拆开”看:平面波是光学的积木块
1.1 平面波为什么是基元函数
信号处理里有个很成熟的做法:任何一段信号,都可以看成不同频率正弦波的线性叠加。你分析一个音频,不会去逐点看波形,而是看它的频谱,因为频谱把“这个信号里有什么成分”一次性说清了。
光学里完全一样。任意一个复杂的单色光场,都可以写成无穷多个平面波的叠加。平面波是最简单的光场:振幅处处均匀,等相位面是一个平面,传播方向固定。它就像一个“积木块”,所有复杂光场都能用这种积木拼出来。
为什么必须是平面波?因为波动方程是线性的,叠加原理成立。一个平面波是波动方程最基本的解,而线性系统的任何一个解都可以分解成基本解的叠加。这个逻辑和傅里叶级数展开一模一样,只是把基函数从正弦波换成了平面波。你看一个复杂光场觉得它很乱,说到底只是因为你在用一个固定的视角看它;换到平面波视角,一切就清楚多了。
1.2 拆完之后的传播问题:角谱法
平面波当积木块的最大好处,是传播过程极其简单。一个平面波沿某个方向传播,从z=0平面走到z=d平面,只需要乘一个相位因子exp(jk_z·d),振幅不变,方向不变。这就是整个角谱理论的核心。
所谓角谱,就是“按传播方向展开的谱”。一个光场在z=0处是什么样,我先把它分解成一堆平面波;每个平面波各自传一段距离,分别乘上各自的传播相位;最后再把它们叠加起来,就得到任意距离处的光场。
这个过程可以类比信号处理里的线性时不变系统:平面波是特征函数,自由空间传播是线性系统,系统对每个特征函数的作用就是乘一个复增益。信号处理里你会把信号分解成不同频率的复指数,看系统对每个频率的响应;光学里你把光场分解成不同方向的平面波,看每个方向在空间里怎么传。两者本质上是同一个套路,只是从时间换成了空间。
1.3 衍射的本质是频谱演化,不是“边缘效应”
教材里讲衍射,经常画那种边缘绕过的示意图,给人感觉衍射是“光碰到边之后拐弯了”。这个图像不能说错,但很容易让人忽略更本质的东西。
一个无限大的理想平面波在自由空间里传多久都不会衍射,因为它的角谱只有一个方向分量,没有任何展宽。但一旦遇到孔径,光场被空间截断,角谱就展宽了——出现了大量偏离原传播方向的高空间频率分量。这些分量传播速度在z方向上的分速度不一样,互相干涉,光场花样就会随距离变化。这就是衍射。
所以衍射的实质是:空间截断导致频谱展宽,随后各频谱分量在传播中相位演化不一致,叠加结果随距离改变。孔径越小,频谱展宽越厉害,衍射越明显。这正是傅里叶变换尺度定理的光学版本:空域越窄,频域越宽。想清楚这一点,后面所有衍射现象都会变得非常顺理成章。
2. 空间频率才是傅里叶光学的第一语言
2.1 空间频率和时间频率是亲兄弟
时间信号x(t)=A·cos(2πf·t),f是每秒振动多少次,单位Hz。空间信号里,一束平面波斜着照到一个平面上,它在x方向的投影是一个周期性分布,E(x)=A·cos(2πu·x),u就是每毫米有多少个明暗周期,单位常用“线对每毫米”(lp/mm)。
时间频率描述的是信号随时间变化的快慢,空间频率描述的是光场在空间上变化的快慢。一个光场里,有大片均匀亮区也有密集细纹理,低空间频率对应那些变化缓慢的区域,高空间频率对应细密的结构和边缘。
这个类比特别有用:做信号处理的工程师对“频率”有本能直觉,但一到光学就开始晕,其实就是没把这两个频率对应上。时间上,频率越高,信号变化越快;空间上,频率越高,图案越密。脑子里有了这个对应关系,傅里叶光学的门槛直接降一半。
2.2 二维空间频率:一个方向对应一个频率
光场是二维分布,所以空间频率也是二维的,写成(u,v)。u对应x方向的周期变化快慢,v对应y方向的。一束平面波沿某个方向传播,在x方向和y方向上都有相位梯度,所以它对应一个二维频率点(u,v)。
这里最关键的一个关系是:空间频率和传播方向是直接对应的。波长为λ的平面波,如果传播方向与z轴夹角为θ,那么它在垂直于z轴的平面上的空间频率就是u=sinθ/λ。方向越斜,频率越高。换句话说,频谱面上的一个点,就代表了一束特定方向的平面波。
这个视角会颠覆你理解“成像”的方式:物体反射的光进入镜头后,并不是简单地被“映射”到像面上,而是先被分解成各个方向的平面波,再由光学系统决定哪些方向能通过、哪些方向被拦掉。每一个方向的平面波,就是一个空间频率分量。
2.3 系统通频带:为什么看不清就是丢高频
所有光学系统的口径都是有限的,所以能收集的平面波方向范围也有限。数值孔径NA=n·sinθ定义了系统能收集的最大角度。对应的,系统能通过的最大空间频率是:
[ v_c = \frac{NA}{λ} ]
这说明了什么?说明分辨率极限的本质不是“像素不够”,而是系统从物理上丢掉了高频空间频率。显微镜看不清楚细节,不是因为探测器分辨率差,而是物镜只收集了一个有限锥角内的光,那些对应更细结构的更大角度衍射光根本没进到系统里。
给做相机的读者一个对应:F数≈F的镜头,出射光锥角近似对应NA≈1/(2F)。550nm绿光下,F/8的镜头截止频率大约227lp/mm。一个像元尺寸3.45μm的传感器,奈奎斯特频率约145lp/mm。也就是说,镜头本身的能力可能还没被传感器完全用满。这种匹配计算,在日常光学系统设计里非常常用。
3. 夫琅禾费衍射就是一次傅里叶变换
3.1 从菲涅耳积分到傅里叶变换:到底近似掉了什么
标量衍射理论里,菲涅耳衍射积分长这样:
[ U(x,y) = \frac{e^{jkz}}{jλz} \iint U_0(ξ,η) \exp\left[j\frac{k}{2z}\left((x-ξ)^2+(y-η)^2\right)\right] dξ dη ]
这个公式猛一看很吓人,但你只需要盯住平方项展开后的样子:
[ (x-ξ)^2 = x^2 - 2xξ + ξ^2 ]
三项里,x²在远场条件下近似对观察面只贡献一个整体相位因子,不影响强度;ξ²是物面自身的二次相位因子,当距离足够远时可以忽略;剩下的-2xξ刚好构成傅里叶变换的核exp(-j2πu·ξ)。
夫琅禾费条件的判据是:
[ z \gg \frac{a^2}{λ} ]
其中a是孔径的特征尺寸。满足这个条件,物面上每个点的自相位近似一致,衍射积分就退化成透射函数的一个傅里叶变换,只不过自变量是空间频率坐标u=x/(λz)。
所以夫琅禾费衍射从来不是什么“额外近似”,它就是傅里叶变换的物理实现方式之一。反过来想:傅里叶变换是一个数学操作,而光的远场传播天然就是做这个操作,这就是傅里叶光学最底层的立足点。
3.2 双缝、光栅、圆孔的频谱怎么一眼看穿
一旦接受夫琅禾费衍射就是傅里叶变换,很多衍射花样都可以直接“看”出来:
双缝的透射函数是两个δ函数的卷积乘以单缝宽度,它的傅里叶变换是一个宽的sinc包络乘上细的余弦条纹。所以你在屏上看到的双缝花样,其实是“单缝衍射调制的双缝干涉”——这句话在频域里根本不需要背,卷积定理一说就通。
光栅是梳状函数,梳状函数的傅里叶变换还是梳状函数,所以光栅衍射是离散的级次,级次间距对应空间频率间隔1/d。光栅方程sinθ=nλ/d,本质就是空间频率分量u=n/d在自由空间传播时对应的方向角。
矩孔的衍射花样是二维sinc函数;圆孔因为圆对称性,傅里叶变换是贝塞尔函数,第一暗环半径1.22λz/D就是艾里斑边界。这些经典结论,全部可以统一到“衍射花样=透射函数的傅里叶变换”这一句话里。以后看到任何衍射花样,第一反应应该是:这个是哪个函数的傅里叶变换?而不是去翻积分公式。
3.3 实验里怎么“看到”傅里叶变换:用透镜拉近远场
夫琅禾费条件z≫a²/λ在实验室里经常是不现实的。1mm的孔径,532nm绿光下,条件要求z远大于1.9米。在桌面上做实验,这个距离很难摆得开。
透镜解决了这个问题。一束平行光经过透镜,在后焦面上会聚;如果把孔径样品放在透镜前方,后焦面上的复振幅分布正好是孔径透射函数的傅里叶变换。无穷远处的夫琅禾费衍射花样,被透镜“拉”到了它的后焦面上。这就是为什么几乎所有光学信息处理实验都用透镜:它把抽象的频谱面变成了物理上可触及的一个平面。
实际操作时,你可以直接用激光扩束后照一个光栅,透镜后焦面放白纸或CMOS,几秒钟就能看到清晰的衍射级次。那个面上的花样,就是光栅透射函数的二维傅里叶变换的幅度分布。这比在远处架屏幕靠谱得多,也清晰得多。
4. 透镜凭什么是傅里叶变换器:后焦面的秘密
4.1 薄透镜的相位变换函数
要理解透镜为什么能做傅里叶变换,得先知道透镜对光场做了什么。理想薄透镜可以建模为一个相位变换:
[ t(x,y) = \exp\left(-j\frac{k}{2f}(x^2+y^2)\right) ]
这个式子看着抽象,但物理图像很简单:凸透镜中间厚、边缘薄,对离轴位置的光线相位延迟更小,等效于给光场乘一个会聚的二次相位因子。平面波经过透镜后变成了向焦点会聚的球面波,正是这个二次相位在起作用。
平时几何光学说“透镜把平行光会聚到焦点”,你已经听得太多了;但在傅里叶光学里,你需要记住的是它的另一面:透镜对入射光场施加了一个二次相位变换。这个相位变换和自由空间传播中的二次相位因子刚好符号相反,所以传播过程中积累的相位弯曲被透镜抵消,最终在傅里叶面上得到一个干净的变换结果。
4.2 前焦面放物,后焦面就是严格频谱面
物体放在透镜前焦面时,物光先传播一个焦距距离到透镜,被透镜乘上二次相位,再传播一个焦距距离到后焦面。推导出来的结果是:后焦面上的复振幅分布就是物透射函数的严格傅里叶变换,连二次相位因子都没有残留,是真正的“纯”频谱。
物体不放前焦面时,后焦面仍然近似是频谱面,但频谱会多一个与物位置有关的二次相位因子。如果你只是探测光强,这个相位不影响结果;但如果你要做相干叠加、滤波后再成像,这个相位就会改变后续的衍射行为。所以严格的光学信息处理实验,都要求物放在前焦面。
还有一个经常被忽略的细节:频谱面的坐标换算关系是
[ u = \frac{x}{λf} ]
波长不同,同一个空间频率对应的位置就不同。宽光谱光源照明的时候,频谱面会看到明显的彩色错位,零频点甚至会呈现彩虹色。这也是很多频谱面实验要用单色光或激光的原因。
4.3 4F系统:光学信息处理的“标准间”
4F系统是傅里叶光学里最经典的架构:物面到第一个透镜,距离为f;第一个透镜后焦面就是频谱面;频谱面到第二个透镜也是f;第二个透镜的后焦面就是输出像面。整个过程“变换-滤波-逆变换”,输入图像经过一次光学傅里叶变换,在频谱面上被人为修改,再做一次逆变换,得到处理后的像。
在这个架构下,滤波变成了一个极其直观的操作:在频谱面插一根针,挡掉零频,像面的背景就被去掉;插一个针孔,只让低频通过,图像立刻变平滑;挡一束高频方向的频谱,对应的纹理方向就会消失。频谱面上的“点”就是“方向”,也就是“频率”,这就是4F系统最有价值的地方。
这个架构不只是教学实验里玩玩。空间相干度测量、光学相关识别、图像边缘增强、暗场/相衬显微的硬件实现,很多都能在4F系统的框架下理解。哪怕你后面做的是数字图像处理,这个“频域滤波”的图像也值得刻在脑子里。
5. 阿贝成像与空间滤波:从频域重新理解成像
5.1 像不是直接生成的,是频谱“干涉”出来的
阿贝1873年研究显微镜分辨率时提出了一个在今天看来非常“傅里叶光学”的观点:成像不是光线简单地从物点映射到像点,而是经历了两次衍射过程。第一次衍射:光通过物结构后形成各级频谱;第二次衍射:这些频谱分量作为有效波源继续传播,在像面干涉叠加,重建出物体的像。
换句话说,成像的本质是频谱的合成。像的质量,取决于系统到底让多少频谱分量参与了最终干涉。高频频谱被系统切得越多,像的边缘就越模糊;低频频谱缺失,像的整体轮廓和背景亮度就会异常。
这个观点对现代光学特别重要。现在做计算成像的人天天说的“编码孔径”“相位恢复”“叠层成像”,本质上都是在跟频谱打交道:要么想方设法恢复被系统丢掉的高频,要么故意对频谱做调制来获取更多信息。
5.2 挡频谱等于改图像:阿贝-波特实验
阿贝-波特实验是理解空间滤波的最佳方式,实验布置就是一个标准的4F系统,物面用光栅,频谱面放各种滤波器件。做完这个实验,几乎所有频域结论都能变成直观记忆:
- 只让零频通过:像面变成均匀亮场。因为零频表示全场平均值,它不能携带任何空间变化信息。
- 挡掉零频:像的背景变暗,边缘变成亮线。因为高频对应边缘和细节,滤掉零频相当于做了一次高通/微分。
- 只让水平方向的频谱通过:原本的二维图像,输出只保留竖直方向的结构。频谱方向对应图像纹理方向,这一点体现得淋漓尽致。
这个实验为什么值得亲手做一次?因为做完以后你看任何光学成像问题,脑子里都会自动浮现频谱面的样子:一个图像的内容落在低频中心,边缘和纹理分布在高频,纹理的方向对应频谱的方向。这个直觉比背一百条成像公式都有用。
5.3 从频域看暗场、相衬和全息
暗场显微的原理就是频谱面滤波:普通显微镜下透明样品几乎看不见,因为它不吸收光,只改变相位;暗场方法用挡光板挡住直接透射的零级光,只让被样品散射的高频光进入成像系统,结果背景变暗而样品细节变亮。本质上就是“挡零频”这个操作的显微镜版本。
相衬显微镜更进一步。透明样品不改变振幅只改变相位,而探测器只能测强度;泽尼克的方法是在频谱面的零频位置放一个移相片,让零级光相位移动约π/2。这样零级光和衍射光在像面干涉时,相位差就转化成了强度差。这也是频谱面滤波直接改善成像能力的经典案例。
全息术就更不用说了,记录的时候用参考光和物光干涉,本质上就是在频域或空域对复振幅做编码,用强度记录下来;重现时靠衍射把相位信息从编码中解出来。理解了“傅里叶光学=频域视角”,再去看全息、相干衍射成像甚至最新的傅里叶叠层显微(FPM),都不会觉得它们是独立的新技术,只是同一套思想的延伸。
6. 实验笔记:找频谱面、估极限分辨、绕开常见翻车点
6.1 桌面实验三步走
想亲眼验证“频谱面就是傅里叶变换面”,不需要复杂设备。一套最简配置:一台激光笔或小型激光器,一个扩束准直透镜组,一片透射光栅或针孔样品,一个焦距50~200mm的凸透镜,一块白卡纸或任意CMOS相机。
第一步,把激光扩束并准直成直径1~2cm的平行光。准直质量非常重要,入射光如果不是平面波,波前带曲率,就会给物光额外引入一个二次相位,频谱面位置和花样都会受影响。
第二步,样品放在透镜前焦面附近。样品可以是一维光栅、单缝、双缝或者打印在透明片上的任意图案。光线经过样品后,紧接着进入透镜。
第三步,在透镜后焦面放屏幕或相机。怎么确认那是后焦面?很简单,前后移动屏幕,找到光斑最锐利、衍射级次最清晰的那个位置。这个位置就是频谱面。别嫌这点时间省不掉,我见过太多人把屏幕放在“大约一个焦距”的位置就开拍,结果花样边缘发虚还找不到原因,其实只是焦面偏了两厘米。
6.2 频率估算公式与一组常用换算
傅里叶光学里最常用的三个估算公式,值得贴在实验本上:
- 系统截止频率:(v_c = NA/λ),数值孔径越大、波长越短,能通过的空间频率越高。
- 传感器奈奎斯特频率:像元尺寸p,奈奎斯特频率 (f_N = 1/(2p))。
- 频谱面坐标换算:(x = u·λ·f),反过来 (u = x/(λf))。
举个例子:532nm绿光,100mm焦距透镜,你在频谱面上把针孔移到离光轴5mm处,对应的空间频率是多少?算一下,(u = 0.005/(532×10^{-9}×0.1) ≈ 94lp/mm)。这就是你亲手滤掉的频率。如果你换200mm焦距透镜,同样是5mm位移,对应只有47lp/mm,频谱被摊得更开,更容易精确操作。这就是长焦透镜在4F系统里受青睐的原因。
下面这组常见数值可以临时查:
| 系统/参数 | 数值 | 截止或奈奎斯特频率 |
|---|---|---|
| 显微镜物镜 NA=0.65,绿光550nm | 截止频率 | 约1182 lp/mm |
| F/8镜头,550nm | 截止频率 | 约227 lp/mm |
| 像元3.45μm的传感器 | 奈奎斯特频率 | 约145 lp/mm |
| 像元1.4μm的传感器 | 奈奎斯特频率 | 约357 lp/mm |
把镜头截止频率和传感器奈奎斯特频率放在一张表里对比,很多“为什么这个镜头配这个相机合适/不合适”的问题会立刻有答案。
6.3 四个真实踩坑记录
第一个坑:照明光没准直。我早期做4F系统实验时,图省事直接用点光源照明,结果频谱面位置会随样品轴向移动而漂移,换上准直光之后所有现象都变得干净利落。原因是波前曲率会给物光加一个二次相位,等于额外引入了一个等效透镜。
第二个坑:透镜口径不够。透镜本身口径有限,相当于一个低通滤波器,高频分量还没到频谱面就被切掉了。这时候你在频谱面上看到的频谱和理论频谱不一样,滤波效果也混杂了系统本身的响应。做实验最好选口径明显大于样品尺寸的透镜。
第三个坑:频谱面光强反差极大。零级亮斑常常刺眼,而高频细节暗得看不见。这是真实频谱的样子,不是设备问题。处理办法是用中性密度衰减片把零级压暗,或者分多次曝光。强行一次拍全,不是过曝就是欠曝。
第四个坑:把频谱面强度当成完整频谱。探测器只能测强度,测不到相位。只凭频谱面的强度分布,无法唯一反推出物光场,这也是相位恢复领域存在的原因之一。做实验时千万别忽略相位信息,很多现象解释不了,就是因为只看了一半。
6.4 先养出“方向=频率”的直觉
如果你只从这篇文章里带走一个东西,我希望是这句:空间频率等于传播方向,频谱面里的一个点就是一束平面波。你什么时候看到衍射花样第一反应是“这是一次傅里叶变换”,而不是去翻衍射积分公式,这个学科就算正式入门了。
这一篇先把平面波拆解、空间频率、夫琅禾费衍射和4F系统这条主线立住。下一篇会深入菲涅耳衍射不做远场近似时怎么处理、PSF与OTF的关系,以及部分相干光场里的傅里叶光学。系列内容不追求堆公式,重点是把物理图像讲顺。毕竟公式可以随时翻书,图景一旦建立起来,就再也丢不掉了。