资讯动态

信息光学复习题解:傅里叶变换、衍射与空间滤波实战指南

发布时间:2026/9/19 7:42:59 来源:尧图企业网站定制
简介面向中科大及同类高校信息光学课程研究生与备考者这份复习题资料系统梳理了课程核心考点。内容聚焦傅里叶变换透镜与普通透镜的区别、相干光与非相干光处理优缺点、菲涅耳与夫琅和费衍射异同、光学传递函数性质等典型问答题并附平面波计算、衍射光栅、4F系统滤波等计算题解题思路帮助突破重点难点。资料为单个doc文档体积1.38MB方便打印或电子阅读。已有311人学习适合期末复习、考研冲刺及自学自测使用能够有效提升对光学信息处理原理的理解与应用能力。1. 拿到这份复习题先别急着背公式「中科大研究生信息光学复习题.doc」这个标题对正在备考的人来说是复习资料但对我们这些已经工作几年的光机电从业者来说它更像一张考点地图。信息光学这门课本质上就是傅里叶变换在光学里的应用从标量衍射到透镜成像从空间滤波到全息术最后都归结到「光场怎么描述、经过系统后怎么变化」这句话上。复习题的价值不在于题目本身而在于它把几十页公式压缩成了几十个高频考点告诉你考试想考什么、老师默认你掌握了什么。如果你正坐在实验室里对着这份 .doc 文件发愁我的建议是不要按章节顺序从头看到尾先把题目按「衍射计算、透镜变换、滤波与成像、全息与相干性」四类分好再针对每一类建立自己的解题模板。后文我会把每类的核心逻辑、常用公式和最容易踩的单位换算坑拆开讲并且给到可以直接抄的 MATLAB 验证脚本让你从「背题」变成「能推」。2. 信息光学的主线傅里叶变换、衍射与透镜的三角关系2.1 标量衍射理论是信息光学的地基信息光学里所有结论都建立在惠更斯-菲涅耳原理的标量近似之上。所谓标量近似就是把光场当成一个复标量函数 $U(x,y)$ 来处理忽略偏振分量之间的耦合。这个近似在绝大多数光学系统里成立因为探测器口径和光波长相比大得多场的横向变化尺度也远大于波长。衍射的数学表达你可以记成一句话夫琅禾费衍射就是一次傅里叶变换菲涅耳衍射就是一次带二次相位因子的傅里叶变换。很多复习题卡住的人其实是没想清楚「为什么要近似」。如果观察屏距离衍射屏足够远远场球面波可以近似为平面波二次相位因子趋于 1积分就退化成标准的傅里叶变换如果距离不够远近场就必须保留二次相位因子。判断远场还是近场就看菲涅耳数 $N_F a^2 / (\lambda z)$ 是否远小于 1其中 $a$ 是孔径特征尺寸。复习的时候我建议把下面这张表抄在笔记本扉页上它把「物理近似条件」和「数学表达式」一一对应近似名称适用条件数学本质典型考题场景菲涅耳近似$z^3 \gg \frac{\pi}{4\lambda}[(x-\xi)^2(y-\eta)^2]^2_{\max}$球面波保留二次项近场衍射、距离较近的成像夫琅禾费近似$z \gg \frac{k(\xi^2\eta^2)_{\max}}{2}$二次相位因子约等于 1远场衍射、频谱面分析傍轴近似光线与光轴夹角很小$\cos\theta \approx 1$薄透镜成像、矩阵光学薄透镜近似透镜厚度远小于焦距位相变换集中在中心平面透镜的傅里叶变换性质2.2 透镜为什么能当傅里叶变换器这是信息光学区别于普通物理光学的分水岭。复习题里凡是出现「透镜」「频谱面」「4f 系统」字样的题目核心都在问一个事透镜前表面的光场和透镜后焦面的光场是什么关系。结论很干净在傍轴条件下薄透镜的复振幅透过率是一个纯相位因子 $t(x,y) \exp[-jk(x^2y^2)/(2f)]$它正好抵消掉菲涅耳衍射积分里的二次相位项。结果是——透镜后焦面上的光场分布正比于透镜前表面光场的二维傅里叶变换只不过空间频率和坐标之间有确定关系 $f_x x_f / (\lambda f)$。注意这里的 $f_x$ 是空间频率单位线对/毫米不是焦距。这个结论的工程意义极大。只要有透镜就能用光学方法做傅里叶变换速度是光速带宽是光学频率远快于任何数字信号处理器。复习题里常让你证明「物体放在透镜前焦面后焦面得到标准傅里叶变换物体放在透镜后频谱面在何处」——前者是标准变换后者是带二次相位的变换。我的记忆方法是物体在透镜前焦面时后焦面是精确傅里叶变换因为相位因子完全抵消物体紧贴透镜前表面时后焦面频谱外加一个球面相位弯曲。2.3 复习时必须记住的 3 组傅里叶变换对研究生考试不是本科期末不会让你查表但也不会考偏门变换对。常考的翻来覆去就三组而且都和物理场景绑定第一组矩形函数和 $\text{sinc}$ 函数。单缝衍射就是典型缝宽 $a$夫琅禾费衍射光强分布为 $I I_0 \cdot \text{sinc}^2(a \cdot x_f / (\lambda f))$。第二组$\delta$ 函数和常数。透射光栅的频谱是若干 $\delta$ 峰对应衍射级次这个结论直接用在空间滤波的级次分离上。第三组高斯函数和高斯函数。高斯光束通过透镜后仍是高斯光束复习激光腔与光纤耦合的题目会反复用到。这三组变换对的价值在于拿它们当「基准信号」所有复杂孔径的衍射都能看成基元信号的线性叠加。做题时先用矩形函数判断主极大宽度再用 $\delta$ 函数判断峰位最后用高斯函数估算能量集中度三步走基本能覆盖 80% 的衍射计算大题。3. 把复习题拆成四类从衍射到全息的解题模板3.1 第一类标量衍射计算题先判断区域再选公式拿到一道衍射计算题第一件事不是套公式而是算菲涅耳数。我给你一个可复用的步骤先看题目给出的孔径尺寸 $a$、波长 $\lambda$、传播距离 $z$如果 $z$ 远大于 $a^2/\lambda$直接用夫琅禾费近似如果不够远只能用菲涅耳积分这时常见做法是把积分里的二次相位因子展开看能不能凑成傅里叶变换或卷积的形式。一道典型的填空题长这样单缝宽 $2\text{ mm}$波长 $632.8\text{ nm}$观察屏距缝 $1\text{ m}$求中央主极大半角宽度。计算很简单$\theta \approx \lambda / a 3.164 \times 10^{-4}\text{ rad}$。但很多人在这一步栽跟头——题目给的是缝宽 2 mm而你在草稿纸上可能下意识写成 $a 2 \times 10^{-3}$ 还是 $\lambda / a$ 的分子分母方向搞反。我的习惯是最后做一次量纲自检角宽度的单位是弧度弧度无量纲所以长度单位必须上下一致$\lambda(\text{m}) / a(\text{m})$。至于多缝衍射和光栅记住「单缝调制 多缝干涉」的乘积关系。光强公式 $I I_0 \cdot \text{sinc}^2(\beta) \cdot [\sin(N\gamma)/\sin(\gamma)]^2$前一项是衍射因子包络后一项是干涉因子细峰。缺级条件就是衍射极小和干涉极大重合$d \cdot \sin\theta m\lambda$ 与 $a \cdot \sin\theta n\lambda$ 联立推出 $d/a m/n$。这个考点反复出现因为它把两个物理过程揉在一起考正好测试你对公式物理含义的理解。3.2 第二类4f 系统与空间滤波画频谱图是得分关键空间滤波是信息光学最「工程」的一章因为它直接对应光学信息处理硬件。4f 系统的结构必须记牢两个焦距为 $f$ 的透镜相距 $2f$物体放在第一个透镜前焦面频谱面在第一个透镜后焦面也是两个透镜的共同焦面输出面在第二个透镜后焦面。复习题里常出现这样的问题「物面上是一个一维光栅缝间距 $d$缝宽 $a$频谱面上放一个低通滤波器只让零级通过求输出面光强分布」。解题时画频谱面示意图是最快的路径频谱面上有三个峰位于 $f_x 0, \pm 1/d$ 处低通滤波只保留 $f_x0$ 的峰逆傅里叶变换后输出是常数对应均匀光场光栅条纹完全消失。关于滤波器参数三个决定性参数要单独记忆滤波器类型传递函数输出效果高频考题问法低通圆孔/方孔截止频率 $f_c$平滑边缘去掉周期细节求截止频率最小值高通挡零频边缘增强、轮廓提取为什么零频被挡后背景变暗带通/方向滤波只让特定频率范围或方向通过提取特定周期结构让某级次衍射光单独通过一个常见的误用是认为高通滤波后对比度一定提高。实际上零频被挡住后输出场的平均值变为 0背景会变成灰色如果你只看强度图视觉上变暗不等于对比度提升。正确的对比度定义是 $(I_{\max} - I_{\min})/(I_{\max} I_{\min})$高通滤波对边缘的锐化来源于强度梯度变大而不是背景变亮。做题时如果题目问「为什么」答这个点能拿高分。3.3 第三类全息术的三种类型分清记录与再现过程全息术在复习题里出现的频率很高而且题型相对固定。核心逻辑是全息记录就是把物光波和参考光波干涉把振幅和相位信息编码成干涉条纹全息再现就是用一个照明光照射全息图通过衍射重建原始波前。三种类型需要分清边界。菲涅耳全息物到记录面的距离在菲涅耳区记录的是菲涅耳衍射图样重建时需要调节照明光波长或曲率来消除像差。傅里叶全息物位于透镜前焦面记录面上得到的是物光场的傅里叶频谱与参考光的干涉重建时再做一次逆傅里叶变换就能得到原物。像全息Image Hologram记录的是物体的像而不是物光本身观察时用白光就能看——因为不同波长的衍射像错位很小可以用白光再现。判断一道题考的是哪种全息就看题目描述里物和记录面的相对位置以及透镜是否存在。关于全息图分辨率和记录介质有一个实用的数量级判断记录介质的空间分辨率至少要是物光最高空间频率的两倍。空间频率上限大约为 $f_{\max} \approx \sin\theta_{\max}/\lambda$其中 $\theta_{\max}$ 是物光与参考光的最大夹角。如果物光和参考光的夹角是 30 度用 $632.8\text{ nm}$ 的 He-Ne 激光$f_{\max} \sin 30° / (632.8\times 10^{-9}) \approx 790\text{ 线对/mm}$。这就是为什么全息干板的分辨率要做到 1000 线对/mm 以上——普通胶片根本达不到这个量级。复习题如果问「为什么全息记录介质必须用高分辨率干板」上面这段话就是标准答案。3.4 第四类相干性与传递函数注意时间相干和空间相干别混淆相干性这块的考点集中在两个地方一个是杨氏双缝干涉条纹的可见度与光源尺寸、光谱宽度的关系另一个是光学传递函数OTF与相干传递函数CTF的区别。时间相干性看相干长度 $L_c \lambda^2 / \Delta\lambda$。如果用白光带宽宽相干长度只有微米量级所以白光干涉仪只能用于厚度极薄的薄膜测量。空间相干性看光源的横向尺寸与干涉孔径的匹配关系范西特-泽尼克定理给出复相干度是光源强度分布的归一化傅里叶变换。复习题考得多的结论是光源直径越小、离干涉屏越远空间相干性越好——这和「用针孔滤波改善激光光束质量」是同一个物理规律。CTF 和 OTF 的区别一句话讲清楚相干成像系统对复振幅是线性的传递函数是瞳孔函数本身CTF取值 0 或 1非相干成像系统对光强是线性的传递函数是瞳孔函数的自相关OTF取值连续下降。前者没有截止频率前的衰减后者从零频就开始下降。做题时如果题目给的是相干照明就用瞳孔函数直接判断截止频率 $f_c D/(2\lambda f_i)$如果是非相干照明就用自相关求 OTF截止频率是共焦时的两倍。4. 高频考题的实操单位、卷积与 MATLAB 验证4.1 统一单位微米、毫弧度与毫米的换算陷阱信息光学的计算题本身不难难的是单位。我做过的复习题里至少 30% 的低级错误来自单位混用波长用了纳米孔径用了毫米焦距用了厘米三个数直接相乘得到完全离谱的结果。我的统一规则是先把所有长度换算成米角度统一用弧度最后在答题纸上标注转换结果。比如波长 $632.8\text{ nm} 6.328 \times 10^{-7}\text{ m}$缝宽 $0.1\text{ mm} 10^{-4}\text{ m}$焦距 $50\text{ cm} 0.5\text{ m}$。在计算空间频率时$f_x x_f / (\lambda f)$ 的单位是 $\text{m}^{-1}$如果你希望得到「线对/毫米」最后再除以 1000。这道题是我经常拿来检验学生单位意识的一个矩形孔径宽 $a 50\text{ }\mu\text{m}$高 $b 2\text{ mm}$求夫琅禾费衍射图样的主极大尺寸。答案是 $2\lambda f / a$ 和 $2\lambda f / b$前者单位是毫米后者也是毫米但因为宽和高差了 40 倍两个方向的主极大尺寸相差 40 倍。题目看起来简单能做对的人不到一半——不是不会衍射是不习惯微米量级。4.2 卷积与相关的符号陷阱移位方向不要搞反复习题的简答题里经常出现「写出两个函数卷积的表达式」或「画出矩形函数与自身卷积的结果」。看似送分实则有坑很多人搞不清卷积的翻转操作把 $f(x) * g(x) \int f(\xi)g(x-\xi)d\xi$ 中的 $g(x-\xi)$ 写成 $g(\xi-x)$导致移位方向反转。相关运算的定义是 $f(x) \star g(x) \int f(\xi)g^*(\xi - x)d\xi$注意相关不翻转但共轭。这两个运算在光学里的物理意义完全不同卷积描述的是成像系统对输入光场的响应叠加相关描述的是两个信号的相似度度量匹配滤波器的基础。做题时我用一个土办法区分卷积的两个函数都要翻一个再滑相关只有一个滑、另一个不翻。在 MATLAB 里验证卷积的方向问题最直接的是用一维数组测试% 验证卷积与相关的方向差异 f [1, 2, 3]; % 物函数假设是三个点的振幅 g [0.5, 1, 0.5]; % 系统脉冲响应 % 卷积MATLAB 的 conv 自动处理翻转 y_conv conv(f, g, full); % 结果为 [0.5 2 4 4 1.5] % 相关需要把 g 翻转后做卷积或者用 xcorr y_corr xcorr(f, g); % 注意 xcorr 默认是零中心的互相关 disp(y_conv); disp(y_corr);conv函数自带翻转操作xcorr给出的是相关结果。如果你把这两个结果画在同一张图上会发现峰值位置不同——这个差异在匹配滤波器的题目里决定你能不能拿满分。匹配滤波器的本质是相关器所以系统输出峰值对应的是信号与模板最相似的位置而不是卷积意义上的「叠加响应」。4.3 用 MATLAB 验证一个 4f 系统的滤波结果复习到后期我强烈建议把高频计算题在 MATLAB 里做数值验证。既能加深印象又能发现公式记忆中的隐性错误。下面是一个 4f 系统低通滤波的完整模拟脚本你可以直接复制运行体会「频谱面操作」的实际效果% 4f 系统空间滤波模拟 % 参数物面采样 512x512波长 632.8nm透镜焦距 0.3m clear; clc; N 512; % 采样点数 L 0.01; % 物面尺寸 10mm dx L / N; % 空间采样间隔 x (-N/2 : N/2-1) * dx; [X, Y] meshgrid(x, x); lambda 632.8e-9; % 波长 f 0.3; % 透镜焦距 k 2*pi / lambda; % 物函数一维光栅周期 d0.4mm缝宽 a0.2mm d 0.4e-3; a 0.2e-3; obj (mod(X, d) a) .* (abs(Y) 2e-3); % 竖直条纹光栅 % 透镜前焦面 - 后焦面傅里叶变换 Uf fftshift(fft2(ifftshift(obj))) * dx^2; % 数值傅里叶变换 fx (-N/2 : N/2-1) / (N*dx); % 空间频率坐标 [FX, FY] meshgrid(fx, fx); % 低通滤波器半径对应截止空间频率 1/d fc 1/d; % 截止频率 H double(sqrt(FX.^2 FY.^2) fc*0.6); % 只让零级和部分一级通过 % 滤波后逆变换回物面 Uout ifftshift(ifft2(ifftshift(Uf .* H))) * N^2 * dx^2; Iout abs(Uout).^2; % 显示对比 figure; subplot(1,3,1); imshow(obj, []); title(原始物面); subplot(1,3,2); imshow(log(1abs(Uf)), []); title(频谱幅度(对数)); subplot(1,3,3); imshow(Iout, []); title(滤波后输出强度);这段脚本里三个细节值得重点说明。第一fftshift和ifftshift配对使用确保坐标原点在数组中心否则频谱会错位半个像素。第二Uf fft2(obj) * dx^2的dx^2是把离散求和归一化成连续积分的近似不做这步输出幅值的量纲就是错误的。第三滤波器半径设成fc*0.6而不是fc是因为实际空间滤波时一级衍射峰的中心位于 $f_x 1/d$主瓣有一定宽度留出余量才不至于把一级峰切掉一半。如果题目要求只让零级通过就应该把半径进一步缩小到fc*0.3以下。运行这段代码你会发现滤波后的输出强度变成均匀背景光栅条纹全部消失——这就是 4f 系统低通滤波的直观结果。全程不需要任何光学平台一台电脑就能把课本上的公式变成看得见的现象。5. 考前 48 小时的验证清单与快速抓分技巧5.1 用「量纲自检」做每道题的最后一关计算题的答案如果通过了量纲检验大概率是对的通不过必然有错。我给自己定的规矩是每道题算出数值后花十秒钟做一个量纲复核。空间频率的单位是 $\text{m}^{-1}$ 或「线对/mm」角度的单位是弧度无量纲电场强度的单位是 $\text{V/m}$——但信息光学题目里常把光场写成复振幅这时要特别注意题目问的是振幅还是强度。振幅是复数的模单位与初始光场一致强度是振幅平方的时间平均单位是 $\text{W/m}^2$。光强公式里出现 $|U|^2$ 时两个复数相乘后的单位是原单位平方。如果题目给的物函数是归一化的透过率函数取值 0 到 1那么计算得到的强度是相对值直接除以入射光强得到衍射效率如果题目给的是绝对振幅比如入射光强 $I_0 1\text{ mW/cm}^2$就需要做归一化换算。这一步是最容易被扣分的地方因为很多复习题的答案给的是相对强度你算出绝对数值后反而对不上标准答案。5.2 背下这张表模板题直接套信息光学期末考试和考研复试的题目高度模式化最后两天不要再看推导过程直接背高频结论题目特征核心公式/结论必用参数单缝夫琅禾费衍射$I I_0 \cdot \text{sinc}^2(a\sin\theta/\lambda)$缝宽 $a$波长 $\lambda$矩形孔衍射主极大尺寸$\Delta x 2\lambda f / a$$\Delta y 2\lambda f / b$焦距 $f$孔径尺寸 $a,b$光栅缺级条件$d/a m/n$整数比缝距 $d$缝宽 $a$透镜后焦面频谱坐标$f_x x_f/(\lambda f)$$f_y y_f/(\lambda f)$波长、焦距、位置坐标相干成像截止频率$f_c D/(2\lambda f_i)$入瞳直径 $D$像方焦距 $f_i$非相干成像截止频率$f_c D/(\lambda f_i)$比相干时大一倍全息图分辨率要求$f_{\max} \sin\theta_{\max}/\lambda$物参光夹角 $\theta_{\max}$相干长度$L_c \lambda^2/\Delta\lambda$中心波长、谱宽背这张表的时候注意区分相干和非相干成像的截止频率相干是 $1/(2\lambda F#)$$F#$ 是焦比非相干是 $1/(\lambda F#)$正好两倍关系。这个考点几乎是每年必考而且容易混。我的联想记忆法是非相干系统对强度成像MTF 是 CTF 的自相关自相关的支撑域比原函数大一倍所以截止频率翻倍。5.3 交叉验证用数值结果反推物理过程最后一个技巧是考前练高频计算题时用「先算数值、再讲物理」的方式自查。比如一道光栅衍射题你算出零级、一级、二级的相对强度分别是 1、0.4、0.1那么你要能说出零级最强是因为缝宽不为零衍射包络在零频处有最大值一级比二级强是因为 $\text{sinc}^2$ 包络随角度增大而衰减。如果你算出来二级比一级还强基本可以断定包络函数写错了。还有一类常考的证明题「证明矩形函数与自身卷积的结果是三角函数」。考前可以不用纸笔推用 MATLAB 验证一遍就够了% 矩形函数 self-convolution - 三角函数验证 x linspace(-3, 3, 1000); rect (x) double(abs(x) 0.5); y conv(rect(x), rect(x), same) * (x(2)-x(1)); % 离散卷积近似 tri max(0, 1 - abs(x)); % 理论三角函数 plot(x, y, x, tri, --); legend(数值卷积, 理论三角);注意这里的same选项保证输出长度和输入一致* (x(2)-x(1))是为了补偿离散近似的采样间隔。对比两条曲线重合你就直观理解了一个重要事实成像系统的分辨率受限于孔径的卷积效应这正是艾里斑的物理来源。复习到这一步信息光学的知识已经形成闭环——从衍射到滤波从滤波到成像从成像到分辨率极限所有考题都能归位到一条逻辑链上此时再看那份 .doc 里的任何一道题都能判断出它考的是链条上的哪个环节。本文还有配套的精品资源点击获取

读完文章,也想定制专属网站?

尧图设计师 24 小时内与您沟通定制方案

免费获取报价