自由曲面面型数学表示与数值建模 - 第4节
4.1 XY多项式的数学形式
XY多项式，说白了就是用x和y的幂次组合来拟合自由曲面。它的数学形式长这样：
z(x,y) = Σ Cij * x^i * y^j其中i和j是非负整数，Cij是对应的系数。i+j就是多项式的阶数。举个例子，i=2,j=0就是x²项，i=0,j=2就是y²项。
实际工程中，我们通常用到5阶或7阶。阶数太低拟合不准，阶数太高又容易过拟合。我个人习惯先试5阶，不够再加。
4.1.1 完整展开形式
拿5阶XY多项式来说，展开后是这样的：
z(x,y) = C00
+ C10*x + C01*y
+ C20*x² + C11*x*y + C02*y²
+ C30*x³ + C21*x²y + C12*x*y² + C03*y³
+ C40*x⁴ + C31*x³y + C22*x²y² + C13*x*y³ + C04*y⁴
+ C50*x⁵ + C41*x⁴y + C32*x³y² + C23*x²y³ + C14*x*y⁴ + C05*y⁵一共21项。你想想看，每个系数都是一个自由度。自由度越多，拟合能力越强，但优化难度也越大。
4.2 各项的物理意义
很多新手拿到XY多项式，觉得就是一堆数学符号。其实每一项都有明确的物理含义。搞懂这个，你才能有针对性地调整系数。
| 项 | 物理意义 | 对应像差 |
|---|---|---|
| C00 | 常数项（平移） | 轴向偏移 |
| C10*x, C01*y | 线性项（倾斜） | 像面倾斜 |
| C20*x², C02*y² | 二次项（离焦） | 离焦、场曲 |
| C11*x*y | 交叉项 | 像散 |
| C30*x³, C03*y³ | 三次项 | 彗差 |
| C40*x⁴, C04*y⁴ | 四次项 | 球差 |
4.2.1 二次项——离焦与场曲
x²项和y²项，说白了就是离焦。如果系统有场曲，你会发现边缘视场的离焦量和中心不一样。这时候调整C20和C02，就能补偿场曲。
我记得有一次做投影镜头，边缘画质总是模糊。折腾了半天，最后发现就是C20和C02没调好。两个系数差了0.002，边缘MTF直接掉了15%。
4.2.2 三次项——彗差
x³项和y³项对应彗差。彗差的特点是：视场越大越明显，而且呈不对称分布。如果你发现点列图有彗尾，优先调C30和C03。
这里要注意：x²y项和xy²项也会产生类似彗差的变形，但方向不同。实际优化时，我建议把同阶的所有项一起放开，让优化器自己平衡。
4.2.3 四次项——球差
x⁴项和y⁴项主要影响球差。球差是轴上像差，跟视场关系不大。如果你发现中心视场的点列图有弥散，先看C40和C04。
嗯，这里有个坑：高次项之间会互相耦合。比如C40和C22都会影响球差，但C22还会影响像散。我曾经吃过这个亏，单独优化C40，结果像散跑飞了。
4.3 拟合与优化技巧
搞清楚了各项的物理意义，接下来就是怎么用。这部分我踩过不少坑，分享几个实战经验。
4.3.1 拟合流程
拟合XY多项式，我一般按这个步骤来：
- 数据准备：从设计目标提取离散点云，通常1000-5000个点
- 归一化：把x和y归一化到[-1,1]区间，避免数值溢出
- 最小二乘拟合：用SVD分解求解系数
- 残差分析：检查拟合误差，如果RMS大于λ/20，考虑升阶
- 边界检查：确保拟合曲面在口径边缘没有翘曲
4.3.2 优化策略
在光学设计软件里优化XY多项式，有几点要注意：
- 从低阶开始：先固定高阶项为0，只优化低阶项。等低阶收敛了，再逐步放开高阶
- 加权重：边缘视场的权重可以设高一些，因为边缘更容易出问题
- 约束边界：给系数加边界约束，防止优化器把曲面拉变形
- 监控斜率：曲面斜率不能太大，否则加工不了。一般控制在±15°以内
我有个习惯：每轮优化完，都会把曲面斜率分布画出来看一眼。有一次优化出来的曲面RMS很好，但斜率局部超过20°，加工师傅直接拒了。
4.3.3 常见问题与对策
| 问题 | 原因 | 对策 |
|---|---|---|
| 拟合误差大 | 阶数不够 | 升阶到7阶或9阶 |
| 曲面边缘翘曲 | 高阶项过强 | 降低高阶项权重或加边界约束 |
| 优化不收敛 | 系数数量太多 | 减少自由变量，先优化主要项 |
| 加工困难 | 斜率过大 | 加斜率约束，或改用Zernike多项式 |
4.4 XY多项式 vs Zernike多项式
说到这儿，你可能想问：为什么不用Zernike多项式？两者确实很像，但有个关键区别。
Zernike多项式在单位圆内正交，适合圆形口径。XY多项式没有正交性，但形式简单，适合矩形口径或自由形状。
我个人习惯：圆形口径用Zernike，矩形口径或异形口径用XY多项式。如果你做的是矩形探测器配光系统，XY多项式是更好的选择。
4.5 实战案例：投影镜头自由曲面
最后分享一个实际案例。去年我做一款超短焦投影镜头，反射镜需要自由曲面。初始设计用的是偶次非球面，但边缘视场畸变太大。
我换成5阶XY多项式，只优化了C20、C02、C40、C04这四个主要项。结果畸变从8%降到了2%以内。你想想看，就四个系数，效果这么明显。
但要注意：不要贪心。我当时试着把全部21个系数都放开，结果优化器跑了三天没收敛。后来老老实实只放开8个关键系数，两天就搞定了。