从理论到实践：Zemax驱动下的离轴三反光学系统全流程构建

1. 从零开始：为什么选择离轴三反系统？

如果你是一位光学工程师，接到一个高分辨率空间遥感相机的设计任务，你的第一反应是什么？是立刻打开Zemax开始画镜片吗？我刚开始做设计的时候也这么干过，结果往往是折腾半天，系统性能一塌糊涂，还得从头再来。后来踩过几次坑才明白，一个优秀的光学设计，尤其是像离轴三反这样的复杂系统，七分靠理论，三分靠软件。软件只是工具，真正决定系统上限的，是你对初始结构参数的计算和推导。

那么，为什么空间相机偏爱离轴三反系统呢？这得从它的“天赋”说起。首先，它是个纯反射系统，这意味着它天生没有色差。无论你的工作波段是从深紫外到远红外，它都能一视同仁，保持稳定的成像质量，这对于需要宽谱段工作的遥感卫星来说简直是福音。其次，反射镜可以用金属或特殊陶瓷制作，热膨胀系数低，在太空剧烈的温度变化下依然能保持形状稳定，不会像折射镜片那样因为热胀冷缩导致像质崩坏。再者，光路可以折叠，能用更短的总长实现长焦距，这对于寸土寸金的卫星载荷舱来说太重要了。最后，镜片可以做得很大、很轻，透过率接近100%，还没有辐射导致玻璃变色的风险。这些优点，让离轴三反系统成为了高分辨率对地观测、天文探测等领域的“明星选手”。

但是，离轴三反的设计门槛也相当高。它不像简单的双胶合透镜，有现成的库克三片式初始结构可以套用。你需要从一个共轴的三反系统出发，通过严谨的数学计算，解出一组合理的初始曲率、间隔和圆锥系数，然后再在软件里进行“离轴化”操作，避开中心遮拦，最终优化到接近衍射极限。这个过程，就像盖房子，理论计算是打地基，Zemax优化是精装修。地基打歪了，后面装修再漂亮，房子也是危房。接下来，我就带你完整走一遍这个从理论推导到Zemax实现的全流程，分享我这些年总结下来的实战经验和避坑指南。

2. 理论基石：手把手推导共轴三反初始结构

拿到设计指标（比如焦距f‘=137mm，入瞳直径D=90mm，视场角等）后，千万别急着打开Zemax。我们的第一步，是在纸面上（或者Matlab、Python里）把共轴系统的初始参数算出来。原始文章里给出了一组距离参数（d1=-130mm, d2=130mm, l3'=-130mm），但我要告诉你，这些参数不是拍脑袋定的，它们与系统的总长、后工作距以及像差校正潜力密切相关。在实际项目中，我们需要根据系统的总长约束、后截距要求来反复调整试探。

计算的核心是解一组基于光焦度分配和像差理论的方程。系统总光焦度Φ（就是焦距的倒数）是已知的。我们设三个反射镜的曲率分别为c1, c2, c3（曲率=1/半径）。根据近轴光学的高斯光学公式，可以建立关于c1, c2, c3的方程组。原始文章给出的方程组就是干这个的。解这个方程组，通常会得到多组解，我们需要根据工程常识进行筛选：比如，曲率不能太大（否则镜片太弯难加工），镜片间的间隔要合理（避免机械干涉），像面位置要合适等。

我个人的习惯是用Python的Sympy库或者Matlab来解这个方程，比手算或者计算器快得多，也便于批量尝试不同的间距组合。比如，你可以写个循环，让d1, d2在某个范围内变化，自动求解并过滤出合理的解组。这里有个小技巧：解出来的c1, c2, c3，最好能保证三个反射镜的遮拦比在一个可接受的初始范围内（例如次镜对主镜的遮拦不要一开始就超过0.4），不然后续离轴化会非常困难。

解出曲率后，就要算更关键的参数：圆锥系数（Conic Constant）。对于三反系统，我们通常使用二次圆锥曲面（包括抛物面、双曲面、椭球面）。这些系数（记为e1^2, e2^2, e3^2）决定了镜面的非球面程度，是校正球差、彗差等像差的关键。求解它们需要引入像差平衡条件，通常是令系统的三级球差（Spherical Aberration）和彗差（Coma）同时为零。这会引出一组更为复杂的非线性方程，如原文中所示。

注意：解这组方程时，数值稳定性很重要。不同的数学工具（如Matlab的solve和vpasolve）可能给出精度略有差异的解。建议将解代回原方程验证，确保误差足够小。

下面我给出一个比原文更清晰的、可直接运行的Matlab代码片段示例，用于求解其中一组中间变量（a1, a2, b1, b2）已知后的圆锥系数：

% 假设已通过前面计算得到以下中间参数 a1 = 0.4873380445; a2 = 1.947118885; b1 = 1.110698735; b2 = 0.4864206763; syms e12 e22 e32 % e1^2, e2^2, e3^2 % 根据像差理论（球差、彗差为零）建立三个方程 eq1 = e12 == 1 + 1/(b1^3*b2^3) * (e22*a1*b2^3*(1+b1)^3 - e32*a1*a2*(1+b2)^3 - a1*b2^3*(1+b1)*(1-b1)^2 + a1*a2*(1+b2)*(1-b2)^2); eq2 = e22*(a1-1)*b2^3*(1+b1)^3 - e32*(a2*(a1-1)+b1*(1-a2))*(1+b2)^3 == (a1-1)*b2^3*(1+b1)*(1-b1)^2 - (a2*(a1-1)+b1*(1-a2))*(1+b2)*(1-b2)^2 - 2*b1*b2; eq3 = e22*b2*(a1-1)^2*(1+b1)^3/(4*a1*b1^2) - e32*(a2*(a1-1)+b1*(1-a2))^2*(1+b2)^3/(4*a1*a2*b1^2*b2^2) == b2*(a1-1)^2*(1+b1)*(1-b1)^2/(4*a1*b1^2) - (a2*(a1-1)+b1*(1-a2))^2*(1+b2)*(1-b2)^2/(4*a1*a2*b1^2*b2^2) - b2*(a1-1)*(1+b1)*(1-b1)/(a1*b1) - (a2*(a1-1)+b1*(1-a2))*(1+b2)*(1-b2)/(a1*a2*b1*b2) - b1*b2 + b2*(1+b1)/a1 - (1+b2)/(a1*a2); % 求解方程组 [e12_sol, e22_sol, e32_sol] = solve([eq1, eq2, eq3], [e12, e22, e32]); % 转换为数值 e12_val = double(e12_sol); e22_val = double(e22_sol); e32_val = double(e32_sol); disp('求解得到的圆锥系数平方值(e^2):'); disp(['e1^2 = ', num2str(e12_val)]); disp(['e2^2 = ', num2str(e22_val)]); disp(['e3^2 = ', num2str(e32_val)]);

运行后，我们会得到一组e^2的值。这里要特别注意：圆锥系数k = -e^2。所以如果e1^2 = 2.5，那么第一个反射镜的圆锥系数（在Zemax中输入的“Conic”值）就是 -2.5。算完这些，恭喜你，共轴三反系统的“地基”参数就全部到手了。我们可以把这些参数整理成一张表，方便输入Zemax。

3. Zemax实战：从共轴到离轴的建模与优化

理论参数齐备，现在可以打开Zemax了。新建一个序列模式文件，按照上表的参数，在透镜数据编辑器（Lens Data Editor）里依次插入3个面（标准面），并将它们的“Glass”类型改为“MIRROR”。准确输入曲率半径、厚度和圆锥系数。这里有个巨坑需要注意：反射镜的厚度。在序列模式下，光线遇到反射镜会折返，因此下一个面的厚度（即到下一个光学面的距离）应该输入负值，或者使用“折叠”镜筒配合坐标断点面。更稳妥的方法是，在输入厚度时，直接按照光线实际传播方向的距离和正负号来输入。对于简单的三反，常见设置是：面1（物面）到面2（M1）为无限远（或入瞳距离），厚度为无穷大；面2到面3（M2）的厚度为-d1（因为光线从M1反射后是往回走的）；面3到面4（M3）的厚度为+d2（光线从M2反射后再次向前）；面4到面5（像面）的厚度为-l3‘（光线从M3反射后折向像面）。正负号取决于你定义的光线方向，务必画个光路草图核对一下！

输入完毕后，设置系统孔径为“入瞳直径”，输入90mm。波长根据你的设计指标来，比如设置8μm, 10μm, 12μm三个波长，10μm为主波长。此时，在0视场下，你的系统应该是一个完美的共轴三反系统，点列图可能很大，因为还没优化，但光路应该是通的。

接下来是关键一步：视场离轴。我们的目的是让有用的视场偏离光轴，从而避开次镜及其支撑结构造成的中心遮拦。在Zemax的视场设置里，选择“矩形视场”（Rectangular）。假设我们最终想要一个在Y方向偏转的离轴系统。X方向视场对称，设为0到2度；Y方向视场，我们想要一个宽度为3度的条带，但这个条带不是以0度为中心，而是以15度为中心。那么，Y视场就应该设置为13.5度到16.5度。这样设置后，Zemax会计算这个矩形网格内所有视场点的成像情况。

提示：离轴量（这里是15度）不是随便设的。它需要足够大，以确保在整个视场内，主光线及其边缘光线都能完全避开次镜的投影区域（即无遮拦）。但离轴量太大又会导致像差（特别是像散和场曲）急剧增大，增加优化难度。通常需要根据初始结构的遮拦情况，通过光线追迹图反复试探找到一个平衡点。

视场设置好，光路看上去应该已经“歪”了。但此时系统性能肯定很差，我们需要启动优化。优化的核心思想是：在实现无遮拦的同时，把像质拉回到衍射极限附近。原始文章提到了两种设置优化操作数的方法，我都用过，这里详细说说我的体会。

方法一（2x3视场点控制边缘像高）：这是原文的“法1”。它的思路是控制矩形视场在X和Y方向的最大、最小视场点在像面上的位置。你需要先用公式计算出这些“理想像高”。例如，Y方向视场角为13.5°和16.5°，焦距137mm，那么理想像高就是f * tan(θ)，计算出来大约是33.0mm和40.7mm。在Zemax的优化函数编辑器（Merit Function Editor）里，你需要为这些视场点配置REAY（实际Y方向像高）和REAX（实际X方向像高）操作数，然后用DIFF操作数计算实际值与理想值的差，最后用OPGT和OPLT操作数把这个差值约束在一个很小的范围内（比如±10μm）。这种方法直接控制像面尺寸，对于约束畸变很有效，但控制点较少，对于视场内部的像质控制可能不够细腻。

方法二（3x3视场点网格控制）：这是原文的“法2”，也是我更推荐的方法。我们在视场设置里就定义3x3的网格点（共9个视场）。优化时，不仅控制边缘视场的像高，还控制相邻视场点之间的像高差，使得像面上的像点分布更均匀。具体操作上，你可以用REAY和REAX分别读取每个视场点的实际像高，然后用DIFF操作数让相邻点的像高差等于根据理想光学公式计算出的差值。这种方法相当于用多个“锚点”把整个像面“绷”住了，优化出来的系统像面更平坦，畸变更均匀。但是，正如原文提到的，这种方法有时会导致优化过程中出现新的、意想不到的遮拦。我就遇到过优化后，边缘视场的光线打到了次镜背面的情况。

当出现遮拦时怎么办？原文给出了一个很好的思路：释放更多的变量，特别是倾斜变量。我的操作步骤通常是：1）先固定主镜厚度，让次镜随之固定，优化一轮；2）如果遮拦消失但像质不佳，在像面前（或三镜后）插入一个“坐标断点”面，将这个面的“倾斜X”和“倾斜Y”设为变量，这相当于允许像面有一个倾斜来更好地接收光线；3）同时，重新将主镜厚度设为变量。这样，优化器就有了更大的自由度去同时调整光路布局和像面位置，往往能跳出局部最优解，找到一个既无遮拦、像质又好的解。

在优化时，评价函数（Merit Function）我习惯用“波前-RMS-质心”作为主要目标，因为它直接优化波像差，与MTF关联紧密。变量方面，除了曲率半径、圆锥系数，我强烈建议把偶次非球面的高阶项（比如6阶、8阶系数）也作为变量。对于离轴三反这种大视场、大孔径系统，仅靠二次圆锥曲面很难校正所有高级像差，引入少量高阶项可以极大地提升优化潜力。但要注意，一开始优化时先别放开高阶项，等系统结构基本稳定、低阶像差校正得差不多了，再逐步加入高阶项进行微调，否则优化器容易跑飞。

4. 性能评估与设计结果分析

经过一番“折腾”，优化评价函数值不再显著下降，光线追迹图清晰无遮拦，这时候就可以停下来，好好看看我们的“作品”到底怎么样了。评估一个离轴三反系统，我主要看四样东西：点列图、MTF曲线、场曲/畸变图、以及光线足迹图。

首先看点列图（Spot Diagram）。在Zemax的分析菜单里打开它，选择显示艾里斑（Airy Disk）。一个好的设计，所有视场点的弥散斑RMS半径应该接近甚至小于艾里斑半径。原始文章说“优化至半径最小状态”，这是一个定性的说法。更定量地，我们可以查看点列图数据表格，看看RMS半径的几何平均值，以及最大视场的RMS值。对于我们的例子（工作波长10μm，F数约为1.5），艾里斑半径大约在16.5μm左右。如果全视场RMS点半径能控制在20μm以内，那就非常优秀了。

MTF曲线是衡量系统分辨率的金标准。打开“MTF”分析，选择“衍射”MTF，空间频率可以设到系统奈奎斯特频率附近。关键是要和衍射极限（Diffraction Limit）曲线对比。理想情况下，我们的MTF曲线应该紧紧贴着那条黑色的衍射极限曲线。原文提到“在空间频率25 lp/mm时接近衍射极限”，这是一个具体指标。在实际项目中，我们更关注在特定空间频率下（比如探测器像元尺寸对应的奈奎斯特频率），全视场的MTF值是否都高于一个阈值（例如0.3）。MTF曲线还能帮你诊断像差类型：如果弧矢和子午曲线分开严重，说明像散大；如果曲线整体下垂，说明离焦或球差大。

场曲和畸变（Field Curvature / Distortion）图对于遥感相机至关重要。场曲过大会导致探测器靶面无法同时清晰成像，通常需要后续的电子校正面阵探测器或采用曲面探测器。畸变会影响几何定位精度。原文通过操作数控制全视场畸变在40μm以内，换算成相对畸变可能不到0.1%，这对于高精度测绘相机是必须的。在Zemax中，你可以直接查看“畸变（Distortion）”分析图，确认其是否符合要求。

最后，别忘了再看一眼光线足迹图（Ray Fan）和系统三维布局图。光线足迹图可以直观看到每个视场、每个孔径的光线像差贡献。三维布局图则是检查机械可行性的第一步：镜片之间有没有空间留给支撑结构？镜片会不会太大导致加工成本飙升？光线在镜筒内有没有发生意外的杂散光反射？这些都是在软件里就要提前考虑的问题。

当所有这些分析结果都满足设计指标时，你的离轴三反光学系统设计才算真正完成。但这还不是终点，接下来还要进行公差分析、热分析、鬼像分析等一系列工程化工作。不过，那又是另一个漫长而有趣的故事了。回顾这个从理论方程到Zemax模型的全过程，我最深的体会是：光学设计没有捷径，每一个参数背后都有其物理意义。扎实的理论推导能让你在优化时心中有数，知道该动哪个变量；而熟练的软件操作则能让你高效地探索设计空间，把理论上的可能性变成现实。这个过程就像解一道复杂的几何题，一开始可能毫无头绪，但当你亲手把一个个参数算出来、输进去、看着系统性能一点点变好，那种成就感，是任何现成设计库都无法给予的。