2、数学基础：自由曲面描述方法（多项式、Zernike、B样条、NURBS）、坐标变换与微分几何基础，以及曲率与法向量的计算
自由曲面光学设计的核心挑战在于：如何用数学语言精确、紧凑且可制造地描述一个非旋转对称的光学表面。与球面或标准非球面仅需曲率半径、圆锥系数等少数参数不同，自由曲面需要更强大的数学工具来表征其空间变化的面形。本章将系统梳理四类主流描述方法（多项式、Zernike、B样条、NURBS）的数学本质与适用场景，并深入讲解坐标变换、微分几何基础以及曲率与法向量的计算方法。这些内容构成了自由曲面设计、优化、公差分析与加工检测的数学基石。
2.1 自由曲面描述方法：从全局多项式到局部控制
2.1.1 核心原理：不同数学基底的面形表达能力
自由曲面的数学描述本质上是在二维参数域 (x, y) 或 (u, v) 上构造一个标量函数 z(x, y)，该函数需满足光学设计对连续性和可微性的要求。四种主流方法各有其数学基底与特性：多项式法（如XY多项式、Q-type）以幂级数或正交多项式为基底，具有全局支撑性，适合描述平滑、低阶的面形偏差；Zernike多项式在单位圆域上正交，其各项系数与Seidel像差（球差、彗差、像散等）有明确的对应关系，便于光学工程师理解面形误差的物理意义；B样条基于分段多项式，具有局部支撑性，即移动一个控制点仅影响附近区域，适合描述非对称、局部突变的面形；NURBS（非均匀有理B样条）在B样条基础上引入权重因子，可精确表示圆锥曲线和二次曲面，是CAD/CAM领域的工业标准。
| 描述方法 | 数学形式 | 支撑域 | 主要优势 | 典型应用场景 |
|---|---|---|---|---|
| XY多项式 | z = Σ cij xi yj | 全局 | 形式简单、求导方便 | 离轴非球面、自由曲面初值 |
| Q-type (Forbes) | z = f(ρ) + 正交多项式 | 全局 | 系数动态范围小、加工友好 | 高精度非球面、自由曲面 |
| Zernike | z = Σ anm Znm(ρ, θ) | 单位圆内全局 | 与像差对应、正交性 | 波前分析、面形误差表征 |
| B样条 | z = Σ Ni,p(u) Nj,q(v) Pij | 局部 | 局部控制、C2连续 | 复杂面形、自由曲面优化 |
| NURBS | z = Σ wij Ni,p Nj,q Pij / Σ wij Ni,p Nj,q | 局部 | 精确圆锥、权重调节 | CAD数据交换、注塑/模压 |
2.1.2 判断方法：如何选择描述方法
选择描述方法需综合考虑面形复杂度、优化自由度、制造约束和软件兼容性。判断流程如下：
- 面形复杂度评估：若面形偏离球面量小于数十微米且变化平缓，优先选择Q-type或XY多项式（低阶项即可）；若存在局部陡变或非对称特征，则需B样条或NURBS。
- 优化自由度需求：多项式方法系数少（通常10-30项），优化速度快但灵活性有限；B样条控制点数量可灵活调整（如10×10网格），适合复杂面形但需注意过拟合。
- 制造与检测兼容性：若后续采用单点金刚石车削（SPDT），Q-type多项式因系数动态范围小，对机床动态响应更友好；若采用磁流变抛光（MRF），Zernike系数可直接映射到驻留时间求解。
- 软件接口：Zemax/OpticStudio原生支持Zernike和XY多项式；Code V支持Q-type；CAD导入需NURBS格式（如STEP/IGES）。
2.1.3 案例拆解：Zernike与B样条在离轴三反系统中的应用对比
某离轴三反望远镜系统，主镜口径300mm，F数为8，视场角2°×2°，要求全视场波前误差RMS < 0.05λ（λ=632.8nm）。初始设计采用XY多项式，优化后RMS为0.12λ，主要残差为高阶彗差和像散。将主镜面形改用Zernike多项式（前20项）重新优化，RMS降至0.07λ，但系统在边缘视场仍存在局部面形误差。进一步将主镜和第三镜改用B样条（12×12控制点网格），经过200次迭代优化，RMS降至0.04λ，满足指标。关键参数对比：Zernike方法优化自由度为20，B样条为144，但B样条引入了局部修正能力，有效补偿了离轴系统固有的非对称像差。代价是优化时间从5分钟增加到40分钟，且需要额外约束控制点平滑性以避免面形突变。
2.1.4 常见误区与总结
- 误区一：Zernike系数越多越好。实际上，高阶Zernike项（如Z20以上）对制造误差极其敏感，且容易引入高频面形振荡，应优先保证低阶项精度。
- 误区二：B样条控制点越多越好。控制点过多会导致面形自由度冗余，优化时产生“波浪形”面形，且加工检测困难。建议控制点数量不超过面形特征波长的2倍。
- 误区三：忽略正交性。XY多项式在矩形域上不正交，系数间存在耦合，优化时需注意系数相关性；Zernike在圆域正交，但若光瞳为矩形（如某些离轴系统），需改用矩形正交多项式。
2.2 坐标变换与微分几何基础
2.2.1 核心原理：从全局坐标到局部坐标的桥梁
自由曲面光学系统通常包含多个光学元件，每个元件都有自己的局部坐标系。坐标变换用于描述光线在全局坐标系与元件局部坐标系之间的转换关系。刚体变换包含三个平移量（dx, dy, dz）和三个旋转量（α, β, γ），可用4×4齐次变换矩阵表示。微分几何则提供了分析曲面局部性质的数学工具：曲面上的点由参数方程 r(u, v) = (x(u,v), y(u,v), z(u,v)) 定义，其切平面由偏导数 ru 和 rv 张成，法向量 N 由两者的叉积归一化得到。第一基本形式（E, F, G）描述曲面上的度量性质，第二基本形式（L, M, N）描述曲面在空间中的弯曲程度。
2.2.2 操作流程：坐标变换与法向量计算步骤
- 建立坐标系：定义全局坐标系（通常以系统光轴为Z轴）和各元件的局部坐标系（原点在元件顶点，Z轴沿表面法线方向）。
- 构造变换矩阵：根据元件的位置和姿态，计算从全局到局部的齐次变换矩阵 T，包含旋转矩阵 R 和平移向量 t。
- 参数化面形：在局部坐标系中，将面形表示为 z = f(x, y) 或参数形式 r(u, v)。
- 计算偏导数：对 f(x, y) 求一阶偏导 fx, fy，二阶偏导 fxx, fxy, fyy。
- 求法向量：在显式形式下，法向量为 N = (-fx, -fy, 1) / √(fx² + fy² + 1)。
- 变换到全局坐标：将局部法向量通过旋转矩阵 R 变换到全局坐标系。
// 伪代码：计算自由曲面在全局坐标下的法向量
// 输入：局部坐标 (x, y)，面形函数 z = f(x, y)
// 输出：全局法向量 N_global
fx = df_dx(x, y); // 一阶偏导
fy = df_dy(x, y);
N_local = normalize([-fx, -fy, 1]); // 局部法向量
N_global = R * N_local; // 旋转到全局坐标
2.2.3 案例拆解：离轴抛物面反射镜的坐标变换与法向量计算
设计一个离轴抛物面反射镜，母镜焦距 f = 500mm，离轴量 h = 80mm。在全局坐标系中，反射镜中心位于 (80, 0, 0)，法线方向与全局Z轴夹角为 9.1°（即离轴角）。局部坐标系原点在反射镜中心，Z轴沿表面法线。面形在局部坐标系中表示为 z = (x² + y²) / (4f)。在局部坐标点 (10mm, 5mm) 处，计算得 fx = 0.01, fy = 0.005，局部法向量为 (-0.01, -0.005, 0.9999)。通过旋转矩阵（绕Y轴旋转9.1°）变换到全局坐标，得到全局法向量 (-0.0158, -0.005, 0.9999)。该法向量用于后续光线追迹中的反射定律计算，确保反射光线方向正确。
2.2.4 常见误区与总结
- 误区一：忽略坐标变换的累积误差。在多镜系统中，每个元件的坐标变换误差会累积，导致光线追迹结果偏差。建议使用双精度浮点数，并定期校验变换矩阵的正交性。
- 误区二：法向量方向不一致。光学设计中通常约定法向量指向入射光一侧，若计算时方向相反，会导致反射定律应用错误。建议在代码中统一法向量方向约定。
- 误区三：混淆局部与全局曲率。曲率是局部几何属性，但计算时必须明确是在哪个坐标系下进行。对于离轴元件，局部曲率与全局曲率存在差异，需通过坐标变换正确转换。
2.3 曲率与法向量的计算
2.3.1 核心原理：曲率描述面形的局部弯曲程度
曲率是自由曲面设计中最关键的几何量之一，直接影响光学系统的像差特性。对于显式曲面 z = f(x, y)，平均曲率 H 和高斯曲率 K 的计算公式为：
H = [(1+fy²)fxx - 2fxfyfxy + (1+fx²)fyy] / [2(1+fx²+fy²)3/2]
K = (fxxfyy - fxy²) / (1+fx²+fy²)²
主曲率 k1, k2 可通过 H 和 K 求得：k1,2 = H ± √(H² - K)。在光学设计中，我们更关心局部曲率半径 R = 1/k，因为它直接决定了光线在表面的偏折能力。法向量的计算已在2.2节详述，此处强调其在曲率计算中的基础作用——曲率本质上描述了法向量沿曲面的变化率。
2.3.2 判断方法：曲率分布分析流程
- 采样面形：在口径内均匀采样（如100×100网格），计算每个采样点的 fx, fy, fxx, fxy, fyy。
- 计算曲率：根据上述公式计算 H, K, k1, k2。
- 可视化曲率图：绘制曲率分布云图，检查是否存在曲率突变区域（通常对应面形不连续或加工风险点）。
- 与像差关联：将曲率分布与系统波前像差对比，识别曲率异常区域对成像质量的影响。
- 制造可行性评估：检查曲率半径是否在加工刀具可达范围内，曲率变化率是否超过机床动态响应能力。
2.3.3 案例拆解：自由曲面透镜曲率分析与工艺优化
某自由曲面准直透镜，口径25mm，采用XY多项式描述（10项）。设计优化后，在Zemax中模拟性能良好（RMS波前误差0.03λ），但试制样品出现局部面形误差，干涉仪检测显示在 (x=8mm, y=3mm) 附近存在一个直径约2mm的凹陷区域，深度约0.5μm。通过计算该区域的曲率分布，发现最小曲率半径仅为3.2mm，远小于SPDT刀具半径（通常为0.5-1mm），导致刀具无法完全进入凹陷区域，产生残留。解决方案：在优化函数中加入曲率约束（最小曲率半径≥5mm），重新优化后，该区域曲率半径提升至6.8mm，面形误差降至0.08μm，满足加工要求。此案例说明，曲率分析不仅是理论计算，更是连接设计与制造的桥梁。
2.3.4 常见误区与总结
- 误区一：只关注曲率大小，忽略曲率变化率。曲率变化率（即三阶导数）决定了加工时刀具加速度需求，过大的变化率会导致机床振动和表面波纹。
- 误区二：法向量计算精度不足。在光线追迹中，法向量误差会直接导致反射/折射方向偏差，建议使用解析求导而非数值差分（如中心差分），以保证精度。
- 误区三：忽略曲率的符号约定。光学设计中通常约定凸面朝入射光方向时曲率为正，但不同软件可能有不同约定，跨软件数据交换时需特别注意。