课程导入：学习目标、课程结构、像差与多项式的基本关系概述

光学系统设计与像差分析是精密仪器、半导体制造、天文观测及消费电子镜头开发中的核心环节。无论是设计一枚手机摄像头，还是搭建一套深紫外光刻投影物镜，工程师最终面对的都是同一个数学问题：如何用一组正交多项式去描述波前偏离理想球面的程度，并从中分离出不同几何像差的贡献量。本章作为整个课程的起点，将明确学习目标、梳理课程结构，并建立“像差—多项式—分解方法”三者之间的基本逻辑框架。我们将从物理直觉出发，理解为什么Zernike多项式在圆形孔径系统中占据主导地位，以及为何在离轴或矩形孔径系统中需要引入其他多项式基底。

1.1 学习目标与课程定位

本课程面向具备基础光学和高等数学背景的工程师与研究人员，旨在建立从物理像差到数学分解的完整知识链条。完成本章及后续章节学习后，学习者应能够：第一，准确识别光学系统中各类初级像差（球差、彗差、像散、场曲、畸变）的波前特征；第二，理解Zernike多项式、Seidel系数、Legendre多项式等不同数学基底的定义域、正交性及适用场景；第三，掌握从离散波前采样数据中提取多项式系数的数值方法，包括最小二乘拟合、Gram-Schmidt正交化及奇异值分解；第四，能够将分解得到的系数反向映射为光学设计中的物理像差量，指导系统优化与公差分析。

课程结构按照“基础理论→方法实现→工程应用”三个层次递进。前五章聚焦于多项式数学基础与像差物理对应关系；第六至十五章深入各类分解算法的原理与实现细节；第十六至二十五章讨论实际光学系统（如光刻物镜、天文望远镜、车载镜头）中的像差分解案例；最后五章涉及高阶像差、非圆形孔径、自由曲面等前沿话题。每章均配有数值实验与代码示例，帮助学习者在实践中巩固理论认知。

1.2 像差与多项式的基本关系

光学像差本质上是对理想球面波前的偏离。在旋转对称光学系统中，波前像差函数W(ρ,θ)可以表示为极坐标下的二维函数，其中ρ为归一化径向坐标（0≤ρ≤1），θ为方位角。数学上，任何连续波前都可以在完备正交基底下展开为无穷级数。选择不同的基底函数，就对应不同的分解策略。Zernike多项式之所以在光学领域占据统治地位，是因为它在单位圆内满足正交性，且其低阶项与Seidel像差存在一一对应关系。例如，Zernike第4项（Z₄）对应离焦，第5、6项对应像散，第7、8项对应彗差，第11项对应球差。

然而，Zernike多项式并非万能。对于矩形孔径（如CMOS图像传感器感光面）或环形孔径（如反射式望远镜），Zernike基底不再正交，此时需要引入Legendre多项式或环形Zernike多项式。理解不同基底的适用边界，是正确选择分解方法的前提。下表总结了常见多项式基底及其适用场景：

从数学角度看，像差分解的本质是求解线性方程组：给定M个采样点上的波前值W(ρᵢ,θᵢ)，寻找N个系数Cⱼ使得ΣCⱼZⱼ(ρᵢ,θᵢ)在最小二乘意义下逼近W(ρᵢ,θᵢ)。当采样点数M远大于系数个数N时，问题为超定方程组，可通过正规方程或SVD求解。当M=N时，问题退化为插值问题。实际工程中，采样点通常包含噪声，因此最小二乘估计是最常用的方法。

1.3 案例拆解：从干涉图到Zernike系数

为建立直观认识，我们以一个典型的Fizeau干涉仪测量球面镜案例来说明完整流程。被测球面镜口径为50mm，曲率半径R=200mm，设计波长为632.8nm（He-Ne激光）。干涉仪输出一幅包含256×256像素的干涉图，通过相位解包裹算法得到波前数据W(x,y)。

初始现象与目标：干涉图显示中心区域存在明显的环形条纹畸变，边缘出现“蝴蝶形”图案。目标是通过Zernike分解确定主要像差成分及其大小，为后续抛光提供依据。

分析过程：首先将直角坐标数据转换为极坐标，并归一化到单位圆内。选择前15项Zernike多项式（对应到第5阶径向）作为基底。构建设计矩阵A（尺寸为65536×15），其中Aᵢⱼ = Zⱼ(ρᵢ,θᵢ)。求解最小二乘问题min||Ac - W||²，得到系数向量c。

# 伪代码示例：Zernike拟合核心步骤
import numpy as np
from scipy.special import zernike

# 假设已获得波前数据 W (256x256) 和坐标网格 rho, theta
A = np.zeros((W.size, 15))
for j in range(15):
    n, m = zernike_index(j)  # 返回径向阶数和角向频率
    A[:, j] = zernike(n, m, rho.flatten())(theta.flatten())

# 最小二乘求解
c, residuals, rank, s = np.linalg.lstsq(A, W.flatten(), rcond=None)
print("Zernike系数:", c)

结果与解读：分解结果显示，Z₄（离焦）系数为0.32λ，Z₅（像散0°）系数为0.18λ，Z₆（像散45°）系数为-0.05λ，Z₈（彗差）系数为0.07λ。这表明主要像差为离焦和0°方向像散，与干涉图观察到的环形条纹和蝴蝶形图案吻合。离焦可能源于测试光路调整不当，而像散则暗示镜面存在非旋转对称的加工误差。

风险与代价：仅使用前15项Zernike多项式可能遗漏高阶像差信息。若镜面存在局部高频误差（如抛光痕迹），这些误差会泄漏到低阶系数中，导致误判。因此，实际工程中应检查拟合残差：若残差RMS大于λ/20，则需要增加多项式项数或考虑其他基底。

1.4 常见误区与总结复盘

初学者在像差分解中常犯以下错误：第一，盲目使用Zernike多项式而不检查孔径形状。对于矩形孔径，Zernike基底不再正交，拟合结果会产生伪像差。第二，忽略多项式项数的选择。项数过少会丢失高频信息，项数过多则可能过拟合噪声。第三，混淆归一化半径的定义。Zernike多项式定义在单位圆内，若实际数据超出该范围，必须进行坐标缩放。第四，将数学分解系数直接等同于物理像差量。例如，Zernike离焦系数包含实际离焦与最佳像面位移的混合贡献，需要结合光学系统参数进行换算。

总结本章核心要点：像差分解是连接物理光学与数值计算的桥梁。正确的分解流程应遵循“孔径形状→基底选择→系数求解→残差验证→物理映射”的完整链路。本章建立的框架将贯穿后续所有章节，后续内容将逐一深入各类多项式的数学性质、高效算法及工程应用细节。