4. 自由曲面建模：常用自由曲面方程（XY多项式、Q-type、径向基函数等）的构建与参数化方法，以及曲面连续性要求
自由曲面光学设计的核心在于“面型自由度”的释放。与传统球面或非球面相比，自由曲面不再受限于对称轴或固定曲率半径，而是通过数学方程直接定义面型在三维空间中的矢高分布。本章聚焦于三种最主流的自由曲面数学表达形式——XY多项式、Q-type多项式、径向基函数（RBF），并深入探讨曲面连续性的工程判定标准。理解这些方程的数学结构、参数化策略以及连续性约束，是后续进行优化设计和公差分析的前提。
4.1 核心原理：三种数学表达式的本质与适用边界
XY多项式是最直观的自由曲面表达方式，其数学形式为：
z(x,y) = c·ρ² / (1 + √(1 - (1+k)·c²·ρ²)) + Σ Cij·x^i·y^j其中第一项为二次曲面基底（c为顶点曲率，k为圆锥系数），第二项为多项式扰动项。XY多项式的优势在于与主流光学设计软件（Zemax OpticStudio、CodeV）原生兼容，设计者可以直接在软件中定义系数并进行优化。然而，XY多项式存在一个显著缺陷：当阶次升高时（如超过6阶），各系数之间的耦合效应增强，面型容易出现非物理的波浪形振荡，且对加工数据点的拟合稳定性下降。
Q-type多项式由Forbes于2007年提出，其核心思想是将面型分解为基底项与正交扰动项的叠加。Q-type分为两种：Qbfs（Q-polynomial for best-fit sphere）和Qcon（Q-polynomial for conic）。以Qbfs为例：
z(ρ) = c·ρ² / (1 + √(1 - c²·ρ²)) + (1/√(1 - c²·ρ²)) · Σ Qn(ρ²)·An其中Qn为Jacobi多项式，ρ为归一化径向坐标。Q-type的最大优势在于其正交性——各阶系数对面型的影响相互独立，优化时不会出现系数“打架”的现象。这使得Q-type在光学设计迭代中收敛速度更快，且对加工检测环节更友好。但Q-type仅适用于圆形或环形孔径，对于矩形或异形孔径无能为力。
径向基函数（RBF）采用完全不同的思路：面型由一系列基函数的加权叠加构成：
z(x,y) = Σ wi·φ(||r - ri||)其中ri为控制点位置，φ为径向对称的基函数（如薄板样条TPS、高斯函数、多二次函数等）。RBF的突出优势是局部控制能力——调整某个控制点的权重只会影响其邻近区域的面型，非常适合局部像差校正。但RBF的参数数量远多于前两种方法（每个控制点至少3个参数），且需要正则化处理以避免过拟合，计算开销较大。
| 对比维度 | XY多项式 | Q-type多项式 | 径向基函数(RBF) |
|---|---|---|---|
| 数学基础 | 幂级数展开 | 正交多项式（Jacobi） | 径向对称核函数 |
| 适用孔径 | 任意形状 | 圆形/环形 | 任意形状 |
| 参数独立性 | 较差（高阶耦合） | 优秀（正交） | 中等（需正则化） |
| 局部控制能力 | 弱 | 中等 | 强 |
| 软件兼容性 | Zemax/CodeV原生 | Zemax 19+支持 | 需自定义DLL |
| 加工检测适配 | 需额外拟合 | 直接适配干涉仪 | 需点云拟合 |
4.2 判断方法：如何选择适合当前设计场景的曲面方程
选择曲面方程类型时，应遵循以下决策流程：
- 确认孔径形状：若为圆形或环形孔径，优先考虑Q-type；若为矩形、自由曲面或异形孔径，则选择XY多项式或RBF。
- 评估优化自由度需求：若需要精细控制局部像差（如视场边缘的畸变），RBF更具优势；若只需整体面型调整，Q-type或低阶XY多项式即可满足。
- 检查加工与检测约束：若后续采用磁流变抛光或离子束修形，Q-type的正交特性可减少面形误差的频谱混叠；若采用单点金刚石车削，XY多项式的刀具路径生成更直接。
- 考虑软件与算法兼容性：若团队使用Zemax进行优化，建议优先使用内置的Q-type或XY多项式；若使用Python/Matlab自研优化算法，RBF的灵活性更高。
4.3 案例拆解：离轴三反系统自由曲面设计中的方程选择与参数化
背景与约束：某航天遥感项目需要设计一个离轴三反望远镜，视场角为5°×5°，F数为4，工作波段为可见光至近红外（0.4-1.0μm）。系统要求MTF在50lp/mm处大于0.3，且体积限制在300mm×200mm×200mm以内。由于离轴结构导致像差分布复杂，传统非球面无法满足全视场成像质量要求。
初始现象与目标：初步采用XY多项式（8阶）进行优化，但优化过程中发现高阶系数（x^6、y^6等）出现大幅振荡，面型在边缘区域产生非物理的波浪形，导致加工检测困难。目标是在保证成像质量的前提下，获得平滑、可加工的面型。
关键输入数据与参数：三片反射镜均采用Qbfs多项式，基底曲率半径分别为-500mm、-300mm、-200mm，圆锥系数k分别为-1.2、-0.8、-0.5。Qbfs阶次选择为6阶（即包含Q1至Q6共6个正交项），每片镜子共7个优化变量（1个曲率+6个Q系数）。
分析判断过程：对比XY多项式与Q-type的优化结果：
- XY多项式（6阶）优化后，RMS波前误差为0.035λ，但面型在边缘出现±2μm的波浪形起伏，超出加工允差（±0.5μm）。
- Q-type（6阶）优化后，RMS波前误差为0.028λ，面型平滑度显著改善，最大面形偏差仅为±0.3μm。
实际操作步骤：
1. 在Zemax中设置三片反射镜的面型为"Qbfs"类型；
2. 设定基底曲率半径和圆锥系数为初始值；
3. 将Q1-Q6系数设为变量，初始值设为0；
4. 使用Damped Least Squares算法进行优化；
5. 监控优化过程中的面型平滑度（通过Zemax的Surface Sag分析）；
6. 优化收敛后，导出面型数据点进行干涉仪仿真验证。结果对比：采用Q-type后，系统MTF在50lp/mm处达到0.35，满足指标要求；面型最大梯度变化率降低至0.8μm/mm，远优于XY多项式的2.5μm/mm。加工阶段采用磁流变抛光，实际加工周期缩短了30%。
风险与代价：Q-type仅适用于圆形孔径，若后续需要将反射镜切边为矩形以适应结构布局，则需重新拟合面型。此外，Q-type在Zemax中的优化速度略慢于XY多项式（约慢15%），因为正交多项式的计算涉及更多数学运算。
复盘经验：对于离轴反射系统，Q-type的正交性不仅提升了优化效率，更重要的是保证了面型的物理可实现性。建议在项目初期就根据孔径形状和加工条件确定方程类型，避免后期返工。
4.4 曲面连续性要求：从数学定义到工程实现
曲面连续性（Continuity）是自由曲面设计中不可忽视的约束。在光学设计中，连续性直接影响光线追迹的稳定性、加工刀具路径的平滑性以及检测数据的可解释性。
G0（位置连续）是最低要求，即相邻曲面片在边界处无缝隙或重叠。在光学设计中，G0不连续会导致光线在边界处发生不合理的偏折，通常出现在拼接面或分段定义的面型中。
G1（切线连续）要求相邻曲面在边界处的切线方向一致。G1连续的面型在视觉上无折痕，但曲率可能突变。对于光学面，G1连续通常足以保证光线追迹的稳定性，但在高精度干涉检测中，曲率突变会产生明显的干涉条纹畸变。
G2（曲率连续）是最高等级的连续性要求，即相邻曲面在边界处的曲率也保持一致。G2连续的面型在光学上表现为无像散突变，在加工上表现为刀具加速度无跳变。对于自由曲面光学元件，G2连续是理想目标，但实现难度较大。
工程实现方法：
- 多项式阶次控制：限制XY多项式的最高阶次（建议≤6阶），可自然避免高阶振荡导致的连续性下降。
- NURBS曲面拼接：在CAD软件中将自由曲面转换为NURBS表示，通过调整控制点权重实现G2连续。
- 局部光顺算法：对优化后的面型数据进行低通滤波或样条平滑，消除高频噪声。