ZEMAX优化函数结构浅探.doc

ZEMAX优化函数结构浅探 公安部第一研究所 许正光 各种光学自动软件最终都归结到优化函数结构和优化过程算法的问题。最近，本人在使用ZEMAX过程中，仔细分析了一下构成ZEMAX软件的优化函数构成以及优化过程算法，有些心得，留给入门的朋友们共享。 优化函数结构 凡使用过SOD88软件(北京理工大学光电工程系开发)或者ZEMAX、CODEV的设计人员都知道，优化的参数包括以下几个种类：光学特性参数，例如焦距、入瞳距离、成像尺寸或者物高、物距，镜片间空气间距、镜片厚度等等；像质参数，例如畸变、场曲、彗差等等。ZEMAX将所有这些要求达到的目标都作为一个优化元附加一定的权重系数组成一个优化函数，并且通过改变结构参数使得这个优化函数趋向最小。 数学表达式为： ，其中为各个优化元的权重系数，为系统结构参数光学追迹出来的各个优化元，例如焦距、畸变、彗差等等，为该优化元的目标值。 优化过程有局部优化和全局优化两种。局部优化是指，通过改变系统结构参数的数值(半径、厚度、光学玻璃材料)计算出各个优化元的数值，然后构成整个优化函数的值的计算过程。该过程的思路是解决当前状态已经处于“U”型中的某个位置，迫使其落到“U”中间的最小位置。 全局优化和局部优化不同的是，优化过程类似于一个搜索过程，这个搜索过程在结构参数限定的某个区域内进行优化，优化函数可能经历若干过波峰和波谷(多个极值之间)进行。由于采用的方法不同，构成了多种全局优化算法。全局优化能够避开某个局部极值寻找到更加优良的结构形式，使得光学设计距离完全自动化更进了一步。当然，目前的各种算法都还有一定局限性，例如搜索能力强度、计算复杂程度，由此影响计算速度、计算资源需求量以及误差累计造成的准确度等等问题。 但是不管怎么样，现有的几种光学设计软件基于现有的高度发达计算机水平、光学设计发展水平和数学优化算法等，已经能够很好的满足具有一定光学设计经验知识的设计者们。 ZEMAX的缺省优化函数结构 入门的光学设计者通常知道在进行结构优化时选用default缺省的优化函数，然后加入少量的优化目标例如焦距来进行优化分析，但是对于这个缺省结构怎么构成的常常缺乏深入分析，这在一定程度上限制了我们进一步充分利用软件优化能力的水平发挥。 实际上，缺省函数的构成结构并不复杂，它和ZEMAX提供给设计人员的“Default Merit Function(缺省优化函数)”紧密相关。如下图所示。 第一行中的“Optimization Function and Reference (优化函数和参数方式)”。他的主要思想是：从某个视场代表物点发出若干条光线，在像面上有一个分布形式，按照各个象差的定义进行象差数值计算。第一个框中有“RMS(方均根)”和“PTV(峰谷差)”，第二个框中有差值计算的各种依据：WaveFront(波前)、Spot Radius(像点尺寸)、Spot X(X方向度量尺寸)、Spot Y(Y方向度量尺寸)、Spot X+Y(X和Y方向平均度量尺寸)。第三个框中有“Cetriod(重心点)”、“Chief Ray(主光线)”、“Mean（平均值）”。 举例说明各自配合的含义，缺省状况一般是：峰谷值＋波前＋质心点，优化的目标是：通过重心的光线到达像面时的相位和其他光线该相位状态时所在的位置有一个位置差别，将这些位置差别减小到最少，即认为各种像差都可能趋近于零。从定义而言，这个组合适合于象差不是很大的场合，对于小象差系统其优化能力非常明显。对于较大象差系统，采用“RMS + Spot Radius + Chief Ray”效果更为明显，后者以到达像面上的各条光线与主光线的差别来进行度量。 当然，这两种方法对于畸变的校正能力是有限的，因此对于畸变要求较高的系统需要加入畸变这个优化目标元进行优化校正。 那么，优化所使用的光线如何获得呢？其实这些光线的分布形式也是从缺省函数中定义的，如左图： 软件提供了两种分布形式，一种是“高斯积分”形式：将轴对称的入瞳面分为数个环和扇面，每一个扇面中选择中心光线作为代表进行计算，每一个视场的所选用的光线数目都是两者的乘积；后一种模式是“矩阵模式”，它对于轴对称的入瞳按照正方形进行各种密度的抽样，一般而言，矩阵模式因为具有去渐晕的能力，在实际的设计中可能更为贴近实际效果。不过由于设计人员的习惯性，一般采用传统的按照孔径直径进行选择的前种方法。通常，光线的数量越多，抽样越密，计算得精确度越高，这主要取决于设计人员的实际系统和计算机资源。 以上的设置就构成了每一个视场对应点的优化结构，这些光线的具体分布与视场大小或者物面尺寸、入瞳直径或者相对孔径有关。显然，每一条光线的优化目标都是和所选择的参考光线的差别为零，因此在优化函数列表中可以看到“target 值”都为零。这样构成的结构在优