3、自由曲面数学描述：多项式（Zernike、XY多项式、Q-type）、NURBS、B样条及网格面形的数学表达与适用场景
自由曲面的“自由”二字，本质上是对传统球面与偶次非球面数学描述能力的突破。在光学设计中，面形数学描述不仅决定了设计变量空间的大小，更直接制约了后续的加工、检测与装配环节。选择错误的数学表达，轻则导致优化收敛缓慢，重则使设计结果在加工阶段完全不可实现。本章将系统梳理光学自由曲面最常用的四类数学工具：以Zernike多项式、XY多项式、Q-type为代表的正交/非正交多项式，以NURBS与B样条为代表的参数化样条曲面，以及直接由离散点云构成的网格面形。我们将从数学定义、适用场景、工程约束三个维度展开，帮助读者建立“设计意图—数学表达—制造可行性”的完整映射关系。
3.1 多项式类描述：Zernike、XY多项式与Q-type
核心原理
多项式类描述将自由曲面表示为基底函数的线性叠加。Zernike多项式在单位圆域内正交，其各项分别对应不同的像差类型（如离焦、像散、彗差、球差），因此在光学设计中具有天然的物理对应性。XY多项式则是x和y的幂次乘积组合，形式简单，适合描述具有明显矩形孔径或柱面特征的面形。Q-type多项式由Forbes提出，其核心思想是将面形分解为基底项与扰动项，基底项描述最佳拟合球面或二次曲面，扰动项则用正交多项式描述偏离量，从而显著改善系数动态范围与制造可行性。
判断方法或操作流程
- 选择依据：若系统孔径为圆形且关注像差平衡，优先Zernike；若孔径为矩形或存在非对称结构，XY多项式更直接；若面形偏离球面较大且需控制加工难度，Q-type是首选。
- 系数归一化：Zernike需将坐标归一化到单位圆，XY多项式需注意幂次项数值稳定性，Q-type需设定基底二次曲率。
- 梯度与雅可比：优化时需提供面形对系数的偏导数，多项式类均可解析求导，利于梯度优化器。
案例拆解：离轴三反系统主镜设计
某离轴三反系统主镜口径200mm，离轴量80mm，初始设计采用Zernike多项式前16项。优化后发现第7项（彗差）系数达到12.3μm，而第4项（离焦）仅0.2μm，系数动态范围过大导致加工检测困难。改用Q-type描述后，基底项吸收了大部分曲率，扰动项最大系数降至1.8μm，面形RMS从0.035λ降至0.021λ（λ=632.8nm），且加工验证一次通过。关键操作：将Zernike系数转换为Q-type基底+扰动，需重新拟合面形采样点，采样密度取128×128。
常见误区
- 误区一：认为Zernike正交性在离散采样下依然严格成立——实际需保证采样点均匀覆盖单位圆。
- 误区二：忽略Q-type基底项与扰动项的权重平衡，导致基底项过拟合。
- 误区三：XY多项式系数过多时（如超过20项）易产生龙格现象，需引入正则化。
3.2 样条参数类：NURBS与B样条
核心原理
B样条通过控制点网格与基函数定义曲面，具有局部支撑性——移动一个控制点仅影响邻近区域。NURBS在B样条基础上引入权重因子，可精确表示圆锥曲线与二次曲面。在光学设计中，NURBS常用于照明系统自由曲面透镜的构建，以及光机结构一体化设计中的光学面与机械基准面衔接。其数学形式为：S(u,v) = ΣΣ N_i,p(u) N_j,q(v) W_ij P_ij / ΣΣ N_i,p(u) N_j,q(v) W_ij。
判断方法或操作流程
- 确定控制点网格分辨率：光学面通常取8×8至16×16控制点。
- 设定节点向量：均匀节点适合简单边界，非均匀节点可增强局部控制。
- 权重初值设为1，优化时先固定权重仅调整控制点位置，收敛后再放开权重。
- 连续性检查：至少保证G2连续（曲率连续），否则影响光线追迹稳定性。
实战场景：LED自由曲面透镜设计
设计一款用于道路照明的LED透镜，目标配光为矩形光斑。采用NURBS描述透镜外表面，控制点网格12×12，初始权重全为1。优化目标为照度均匀度>0.8，光效>85%。经过50次迭代，控制点Z坐标最大变化4.2mm，权重变化范围0.7~1.3。最终均匀度0.83，光效88%。但发现边缘区域出现波纹，原因是节点向量在边界处分布过密。调整节点向量使边界处节点间隔增大后，波纹消失。
常见误区
- 误区一：控制点越多越好——过多控制点导致优化变量爆炸且曲面易产生局部抖动。
- 误区二：忽略权重对曲率的影响——权重突变会造成曲率不连续，影响注塑脱模。
- 误区三：将NURBS直接用于成像系统优化——其梯度计算复杂，且与像差理论脱节。
3.3 网格面形：点云、STL与OBJ格式
核心原理
网格面形由离散的三维坐标点（点云）及连接关系（三角面片或四边形面片）构成。其来源包括：逆向工程扫描数据、有限元分析变形结果、以及光学设计软件导出的离散采样。网格面形没有解析表达式，但可通过插值或拟合转换为多项式或样条形式。在光学仿真中，网格面形可直接用于光线追迹，但需注意法向量计算的准确性。
判断方法或操作流程
| 数据来源 | 网格密度建议 | 预处理操作 | 后续处理 |
|---|---|---|---|
| 干涉仪测量 | 256×256 | 去噪、去倾斜、去离焦 | 拟合Zernike或Q-type |
| 有限元变形 | 按应力梯度自适应 | 去除刚体位移 | 拟合B样条 |
| 设计软件导出 | 按需加密 | 检查法向一致性 | 直接追迹或拟合 |
案例拆解：空间相机反射镜重力变形分析
某空间相机主镜口径1.2m，在地面装调时受重力影响产生变形。有限元分析输出网格面形（节点数约50万），导入光学软件后追迹发现波前RMS从设计值0.02λ恶化至0.15λ。分析步骤：①提取镜面节点位移；②去除刚体平移与旋转（最小二乘拟合）；③将残余变形拟合为Zernike前36项，发现主要贡献为像散（0.09λ）与彗差（0.04λ）；④在光学模型中反向补偿该像差，验证了重力卸载方案的可行性。最终采用三点支撑优化，将RMS控制在0.03λ以内。
常见误区
- 误区一：直接使用原始点云追迹——未去除测量噪声导致追迹结果失真。
- 误区二：网格过度细化——追迹速度下降且内存占用剧增，建议在曲率变化大的区域局部加密。
- 误区三：忽略网格拓扑错误——如法向翻转、重复顶点、非流形边，需用专业工具修复。
3.4 四种描述方式的对比与选型总则
| 维度 | Zernike | XY多项式 | Q-type | NURBS/B样条 | 网格面形 |
|---|---|---|---|---|---|
| 连续性 | 无限阶可导 | 无限阶可导 | 无限阶可导 | 可达C^2以上 | 离散，需插值 |
| 局部控制 | 全局 | 全局 | 全局 | 局部 | 局部 |
| 像差对应性 | 强 | 中 | 中 | 弱 | 无 |
| 制造友好性 | 中 | 中 | 高 | 高 | 低 |
| 优化效率 | 高 | 高 | 高 | 中 | 低 |
| 典型应用 | 成像系统像差平衡 | 矩形孔径照明 | 成像/照明最终面形 | 照明透镜、光机结构 | 测量数据、FEA结果 |
选型总则：设计初期用Zernike或XY多项式快速探索像差空间；设计中期若面形偏离球面较大，切换Q-type以改善系数条件数；涉及复杂拓扑或机械接口时引入NURBS；最终验证阶段用网格面形承载测量或仿真数据。实际工程中，多种描述方式常共存——例如设计用Q-type，CAD模型用NURBS，检测用网格面形，三者之间需建立严格的转换与误差评估流程。