Learning objectives / 学习目标
- 在 linear shift-invariant 近似下由 PSF 写出卷积成像，并由 PSF 定义 OTF/MTF。
- 用 modulation 与 spatial frequency 解释“锐度”不是单一数字。
- 从 ray/wavefront 图景区分 spherical、chromatic、coma、astigmatism、field curvature 与 distortion。
- 判断系统是 diffraction-limited 还是 aberration-limited，并用 Strehl approximation 做定量估算。
- 识别 PSF/MTF 测量的视场、波长、焦面与相干性依赖。
Prerequisites / 先修知识
- Paraxial imaging and pupils。
- Single-aperture diffraction 的 Airy/sinc 图景。
- Fourier transform 与 convolution 的基本定义；复振幅与强度 \(I\propto\lvert U\rvert^2\)。
Scope and conventions / 范围与约定
- 先讨论一个小视场内的 linear shift-invariant (LSI) imaging；强烈 field-dependent aberration 时必须分区使用 PSF。
- \(h_I(x,y)\) 表示 incoherent intensity PSF，归一化使 \(\iint h_I\,dx\,dy=1\)。
- \(\mathrm{OTF}=\mathcal F\{h_I\}\)，\(\mathrm{MTF}=\lvert\mathrm{OTF}\rvert/\mathrm{OTF}(0)\)；coherent imaging 的 field transfer function 不是同一对象。
- wavefront error \(W\) 是实际 wavefront 相对 reference sphere 的 optical path difference，单位 m；相位误差 \(\phi=2\pi W/\lambda\)。
- Seidel aberrations 在 rotationally symmetric、small-field、third-order 近似下分类；本课不推导所有系数。
- 图中 spot sizes 多为 schematic；声称比例时才按比例。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit |
|---|---|---|
| \(h_I(x,y)\) | normalized incoherent intensity PSF | m\(^{-2}\) |
| \(I_o,I_i\) | object and image intensity distributions | W m\(^{-2}\) |
| \(\nu_x,\nu_y\) | spatial frequencies | cycles m\(^{-1}\) |
| OTF, MTF | optical and modulation transfer functions | 1 |
| \(C\) | Michelson contrast | 1 |
| \(W,\sigma_W\) | wavefront error and its pupil RMS | m |
| \(S\) | Strehl ratio | 1 |
| \(\lambda,D_a\) | wavelength and pupil diameter | m |
Knowledge 1 — Point-spread function and modulation transfer / 点扩散函数与调制传递
Knowledge ID: P05-U03-L03-K01
PSF as the response to a point
有限 aperture、defocus 与 aberration 使 point object 形成有限斑。对小视场内的 incoherent LSI system：
\[I_i(x,y)=\iint I_o(\xi,\eta)\, h_I(x-M\xi,y-M\eta)\,d\xi\,d\eta,\]
即（忽略 magnification coordinate change 后）\(I_i=I_o*h_I\)。energy-normalized PSF 满足
\[\iint h_I(x,y)\,dx\,dy=1.\]
图 1 把“点变斑”定义为系统冲激响应；若 PSF 随物点位置变化，简单全图卷积失效。
OTF and MTF
取 Fourier transform：
\[\mathrm{OTF}(\nu_x,\nu_y)=\mathcal F\{h_I\}, \qquad \mathrm{MTF}=\frac{\lvert\mathrm{OTF}\rvert}{\mathrm{OTF}(0)}.\]
对正弦目标 \(I_o=\bar I[1+C_o\cos(2\pi\nu x)]\)，若系统在该方向 LSI，则
\[C_i(\nu)=\mathrm{MTF}(\nu)C_o(\nu).\]
Michelson contrast 定义为
\[C=\frac{I_{\max}-I_{\min}}{I_{\max}+I_{\min}}.\]
图 2 表示 MTF 传递对比度，不表示输出“频率变慢”。phase transfer 若非零，还会移动/扭曲结构，需看完整 OTF。
Knowledge 2 — Spherical and chromatic aberration / 球差与色差
Knowledge ID: P05-U03-L03-K02
Spherical aberration
球面折射面只在 paraxial limit 近似把同轴平行光聚到同一点。marginal rays 与 paraxial rays 的折射量不同，焦点沿轴分散，形成 caustic 和 circle of least confusion。缩小 aperture 可减小高阶 ray contribution，但会增强 diffraction。
图 3 画出典型正球差；具体焦点先后可随 optical design 改变，不能把图示方向当作普遍符号定律。
Chromatic aberration
折射率 \(n(\lambda)\) 使 lens power 随波长变化。正常色散下短波 \(n\) 较大，正透镜通常 blue focus 更近，red focus 更远。轴向色差改变焦位；lateral color 使 magnification 也随波长变。
\[\Phi(\lambda)=\left[n(\lambda)-1\right]K\]
是空气中 fixed-shape thin lens 的简化表达，\(K=1/R_1-1/R_2\)。
图 4 明确限定 normal dispersion 与 positive lens；mirror focusing 本身不产生这种 refractive chromatic power。
achromatic doublet 通过不同 glass dispersions 与相反 powers 使两个选定波长共焦；不等于所有波长、所有像差完全消失。
Knowledge 3 — Coma, astigmatism, field curvature, and distortion / 彗差、像散、场曲与畸变
Knowledge ID: P05-U03-L03-K03
这些主要是 off-axis field 的 third-order Seidel aberrations：
- coma： 同一 off-axis point 的不同 pupil zones 具有不同 magnification，PSF 呈 comet-like tail。
- astigmatism： tangential 与 sagittal ray fans 焦在不同 axial positions，中间有 circle of least confusion。
- field curvature： 最佳焦面为 curved Petzval surface，flat sensor 边缘失焦。
- distortion： point image 可仍小，但 magnification 随 field height 改变；直线弯成 barrel 或 pincushion。
图 5 是 coma spot 的结构示意，不是 ray trace 的定量 spot diagram。
图 6 合并展示 tangential/sagittal 分离与 curved best-focus surface；二者相关但不是同一个定义。
图 7 强调 distortion 与 blur 的区别：几何坐标错误不必伴随 PSF 变宽。
Knowledge 4 — Diffraction-limited versus aberration-limited imaging / 衍射极限与像差极限成像
Knowledge ID: P05-U03-L03-K04
理想 finite pupil 即使 \(W=0\) 也有 diffraction PSF；aberration 添加 pupil phase \(\exp[i2\pi W/\lambda]\)，把能量从中央峰移向旁瓣或不对称结构。系统称 diffraction-limited，通常表示 residual aberration 足够小，使性能接近同一 pupil 的理想 diffraction baseline，而不是“没有衍射”。
Strehl ratio 定义为相同总通量和 pupil 下，实际 on-axis PSF peak 与 ideal peak 之比。小 RMS phase error 的 Maréchal approximation：
\[\boxed{S\approx \exp\left[-\left(\frac{2\pi\sigma_W}{\lambda}\right)^2\right]}.\]
\(S\gtrsim0.8\) 常作工程判据，但不是所有任务的普遍质量阈值。
图 8 示意 aberration 降低峰值并扩宽/偏移 PSF；Gaussian 曲线不是 Airy 的精确公式。
图 9 比较同类系统的 MTF；必须在同一 wavelength、field、focus、orientation 和 normalization 下比较。
图 10 是审计流程：先定义理想 pupil baseline，再加入 wavefront error，最后才叠加 detector 与 noise。
Worked Examples / 完整例题
Worked Example 1 — Gaussian PSF width and MTF
已知： 一维归一化 intensity PSF 近似 \(h(x)\propto e^{-x^2/(2\sigma^2)}\)，\(\sigma=8.00\ \mu\mathrm m\)；求 FWHM 与 \(\nu=20.0\ \mathrm{cycles\,mm^{-1}}\) 的 MTF。
研究对象/系统： 局部 incoherent LSI imaging。
参考系/正向： sensor \(x\) 轴；MTF 非负并归一化。
模型： Gaussian approximation；真实 Airy/aberrated PSF 仅作局部拟合。
\[\mathrm{FWHM}=2\sqrt{2\ln2}\,\sigma =2.35482(8.00)=18.8\ \mu\mathrm m.\]
Gaussian transform 给
\[\mathrm{MTF}(\nu)=\exp(-2\pi^2\sigma^2\nu^2).\]
用 \(\sigma=0.00800\ \mathrm{mm}\)：
\[\mathrm{MTF}(20.0) =\exp[-2\pi^2(0.00800)^2(20.0)^2] =0.603.\]
结论： FWHM \(18.8\ \mu\mathrm m\)，20 cycles/mm 处传递 60.3% modulation。
检查： \(\sigma\nu\) 无量纲；\(\nu\to0\) 时 MTF \(\to1\)，高频趋 0。
Worked Example 2 — Measured contrast transfer
已知： object sinusoid \(I_{\max}=100\)、\(I_{\min}=20\)；image 为 70 与 30，同一 spatial frequency、线性响应。
研究对象/系统： 单频正弦 target 与局部 imaging system。
参考系/正向： 强度非负；contrast 取正。
模型： 无 additive offset drift，phase shift 不影响 extrema contrast。
\[C_o=\frac{100-20}{100+20}=0.6667, \qquad C_i=\frac{70-30}{70+30}=0.4000.\]
\[\mathrm{MTF}=\frac{C_i}{C_o}=0.600.\]
结论： 该频率 MTF 为 \(0.600\)。
检查： \(0\le C,\mathrm{MTF}\le1\)；若输出 offset 未校正，直接使用 extrema 会偏置结果。
Worked Example 3 — Longitudinal chromatic focal shift
已知： 空气中 thin positive lens 的 curvature factor \(K=4.00\ \mathrm{m^{-1}}\)；\(n_B=1.530\)，\(n_R=1.510\)。
研究对象/系统： 固定形状的单材料薄透镜。
参考系/正向： 光向右；positive focal length。
模型： \(\Phi=(n-1)K\)，忽略 thickness 与 secondary spectrum。
\[\Phi_B=(1.530-1)(4.00)=2.120\ \mathrm D, \qquad f_B=0.4717\ \mathrm m,\]
\[\Phi_R=(1.510-1)(4.00)=2.040\ \mathrm D, \qquad f_R=0.4902\ \mathrm m.\]
\[\Delta f=f_R-f_B=0.0185\ \mathrm m=18.5\ \mathrm{mm}.\]
结论： blue focus 比 red focus 靠近 lens，longitudinal separation \(18.5\ \mathrm{mm}\)。
检查： normal dispersion \(n_B>n_R\) 给 \(\Phi_B>\Phi_R\)；若无色散，两焦点重合。
Worked Example 4 — Strehl estimate from RMS wavefront error
已知： \(\sigma_W=\lambda/20\)；求 Maréchal Strehl estimate 并判断是否超过 0.8。
研究对象/系统： 与 ideal system 同 pupil、同总通量的 on-axis point image。
参考系/正向： \(W\) 的 piston 已去除，\(\sigma_W\ge0\)。
模型： small RMS phase error, approximately random/smooth residual aberration。
\[S\approx\exp\left[-\left(2\pi\frac1{20}\right)^2\right] =\exp(-0.09870)=0.906.\]
结论： \(S\approx0.906\)，超过常用 0.8 engineering criterion。
检查： \(\sigma_W\to0\) 时 \(S\to1\)；结果不代表所有 field/wavelength 都相同，也不替代完整 MTF。
Lesson synthesis / 本课综合
质量分析从 point object 开始：PSF 描述空间扩散，OTF/MTF 描述各 spatial frequency 的幅度和相位传递。aberration 不是一个统一“模糊量”：spherical/chromatic 主要含 on-axis/aperture/wavelength 依赖，coma/astigmatism/field curvature/distortion 强烈依赖 field。最终必须与相同 pupil 的 diffraction baseline 比较，再加入 detector、sampling 与 noise。
Formula-and-decision table / 公式与决策表
| Need | Formula | Conditions | Caution |
|---|---|---|---|
| incoherent image | \(I_i=I_o*h_I\) | local LSI system | PSF may vary with field |
| OTF/MTF | \(\mathrm{OTF}=\mathcal F\{h_I\}\); \(\mathrm{MTF}=\lvert\mathrm{OTF}\rvert/\mathrm{OTF}(0)\) | incoherent intensity imaging | coherent transfer differs |
| contrast | \(C=(I_{\max}-I_{\min})/(I_{\max}+I_{\min})\) | sinusoidal/periodic target | correct offsets |
| fixed-shape lens power | \(\Phi(\lambda)=[n(\lambda)-1]K\) | thin lens in air | simplified dispersion model |
| phase error | \(\phi=2\pi W/\lambda\) | scalar wave | piston does not affect image |
| Strehl estimate | \(S\approx e^{-(2\pi\sigma_W/\lambda)^2}\) | small RMS aberration | engineering approximation |
Model boundaries / 模型边界
- Global convolution requires shift invariance; off-axis aberrations often force field-dependent PSFs.
- incoherent intensity PSF/OTF and coherent field transfer must not be interchanged.
- third-order Seidel taxonomy is an approximation; real design uses full ray trace/wavefront maps and often Zernike bases.
- distortion moves coordinates without necessarily reducing point sharpness; geometric correction can remap it but changes sampling.
- \(S\ge0.8\) is a conventional criterion, not universal task performance.
- diffraction-limited means near ideal for the stated pupil/wavelength/field/focus, not “no blur.”
Misconception clinic / 误区诊所
- “PSF 就是 lens 的固定属性。” 它依赖 wavelength、field, focus、aperture、polarization 和 coherence。
- “MTF 只要一个数字。” MTF 是 spatial-frequency（通常还含方向）的函数。
- “缩光圈总能修好球差。” 它减少 marginal rays，也扩大 diffraction PSF。
- “reflecting telescope 没有像差。” 无 primary chromatic power，但仍有 monochromatic aberrations。
- “distortion 就是 blur。” distortion 可保持每个点很尖，却把坐标映射错。
- “Strehl 0.9 说明所有细节传得很好。” 它只比较 central peak；完整任务仍需 MTF、field 与 noise。
Practice and retrieval / 练习与提取
- K01： PSF 的 integral 为什么归一为 1？
Answer
对无增益归一化系统，它表示一个点源的总能量在 image plane 重分布而非凭空增加；若有 throughput loss，可先除以总通量再讨论形状。
- K01： object contrast 0.80，某频率 MTF 0.25，image contrast？
Answer
\(C_i=0.25(0.80)=0.20\)。
- K01： \(I_{\max}=90\)、\(I_{\min}=30\)，contrast？
Answer
\(C=(90-30)/(90+30)=0.500\)。
- K02： positive lens 在 normal dispersion 下，blue 与 red 哪个通常焦得更近？
Answer
blue；通常 \(n_B>n_R\)，故 blue power 更大、focal length 更短。
- K02： achromatic doublet 是否消除全部色差？
Answer
否；它通常令两个选定 wavelengths 共焦，仍可能有 secondary spectrum、lateral color 和其他 aberrations。
- K03： 哪种 aberration 可使 straight grid 弯曲但 point images 仍小？
Answer
distortion；magnification 随 field height 改变。
- K03： tangential 与 sagittal foci 分离对应什么？
Answer
off-axis astigmatism。
- K04： 若 \(\sigma_W=\lambda/14\)，估算 Strehl。
Answer
\[S\approx\exp[-(2\pi/14)^2]=0.818.\]
- K04： 为何理想无像差 lens 仍有有限 PSF？
Answer
pupil 有限，截断波前导致 diffraction；point image 是 Airy-like pattern 而非 delta function。
- 综合： 中心 field MTF 高、边缘呈 comet-tail PSF，最可能是哪类问题，为什么不能用一个全图 convolution？
Answer
典型为 field-dependent coma；PSF 随 field position 改变，系统不再在全图范围 shift invariant。
Exit ticket / 离堂检核
- 从 point object 到 PSF，再到 sinusoidal target 的 MTF，写出完整逻辑链。
- 给一幅边缘模糊且线条弯曲的图，分别设计实验区分 field curvature、distortion 与 sampling alias。
References and source notes / 参考资料与来源说明
Source-to-content map
| Source | Used for | Knowledge ID / section | Model boundary or convention |
|---|---|---|---|
| MIT OCW 2.71 Lecture 8 | spherical, chromatic, coma and Seidel picture | K02–K03 | rotationally symmetric third-order view |
| MIT OCW 2.71 Lecture 9 transcript | astigmatism and field curvature | K03 | off-axis focal surfaces |
| MIT OCW 2.71 Lecture 22 | PSF, OTF/MTF, coherent/incoherent distinction | K01, K04 | local LSI approximation |
| OpenStax `4.5 Circular Apertures | diffraction baseline and Rayleigh | K04 | ideal circular aperture |
| Optica/OSA Aberrations tutorial | wavefront reference and aberration taxonomy | K02–K03 | pedagogical ray/wavefront descriptions |
| University of Arizona OPTI 415R | OPD, PSF, Seidel/Zernike and MTF metrics | K01–K04 | declared normalization required |
资料和链接实际打开核对日期：2026-08-28。图示原创；例题与练习独立复算。
复习与拓展
- Optical Instruments：aperture 和 magnification。
- Two-Source Interference：phase/OPD。
- Single-Aperture Diffraction：ideal PSF。
- Resolution and Fourier Optics：OTF、convolution 与 filtering。