CN108227184B - 自由曲面成像光学***的设计方法 - Google Patents

Abstract

本发明涉及自由曲面成像光学***的设计方法，包括：建立初始***；选取特征光线，求解自由曲面a上的特征数据点，拟合得到自由曲面a；以此类推，得到第一自由曲面成像光学***；将第一自由曲面成像光学***作为初始***，多次迭代，得到第二自由曲面成像光学***；将第二自由曲面成像光学***作为参考***，选一个曲面位置自由度，在该曲面位置自由度的值附近选取多个值，改变对应曲面的位置，得到多个第三自由曲面成像光学***；将每个第三自由曲面成像光学***作为初始***，重复迭代步骤，得到多个第四自由曲面成像光学***；挑选出满足约束条件且在评价指标下最优的第四自由曲面成像光学***；以此类推，直到所有曲面位置自由度均被使用过。

Description

自由曲面成像光学***的设计方法

技术领域

本发明涉及光学设计领域，尤其涉及一种自由曲面成像光学***的自动化设计方法。

背景技术

现有的成像光学***的设计方法中，设计流程一般是先找到一个好的初始结构，然后进行计算机辅助的优化。

然而，这种设计方法并不是十分智能，功能也不够强大。其一，此种设计方法高度依赖初始结构，但可用的初始结构一般是有限的，如果初始结构的结构与参数与设计要求相差很远，可能需要进行长时间的优化。而且对于一些指标先进，结构特殊的成像光学***，很可能无法找到合适的初始结构。如此，设计者很可能花费大量的时间进行优化，甚至失败。其二，优化过程是一个“技术”与“思维”的有机结合，软件优化算法是很重要的，而人对优化过程的指导与干预以及人的经验与技能等也很重要，因此整个设计过程高度依赖人力参与。

发明内容

综上所述，确有必要提供一种自由曲面成像光学***的自动化设计方法。

一种自由曲面成像光学***的设计方法，其包括以下步骤：

步骤S1，建立一初始***，该初始***包括多个初始曲面，且该初始***中的一个初始曲面对应待设计自由曲面成像光学***中一个自由曲面；

步骤S2，根据待设计自由曲面成像光学***的设计要求，建立自由曲面成像光学***的约束条件；

步骤S3，将待设计自由曲面成像光学***中的一个自由曲面定义为自由曲面a，保持所述多个初始曲面不变，选取K条特征光线R_i(i＝1，2…K)，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面a上的多个交点，进而得到自由曲面a上的多个特征数据点，将该自由曲面a上多个特征数据点进行曲面拟合，得到所述自由曲面a；

步骤S4，将待设计自由曲面成像光学***中的另一个自由曲面定义为自由曲面b，保持所述自由曲面a以及自由曲面a对应的初始曲面之外的其它初始曲面不变，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面b的多个交点，进而得到所述自由曲面b上的多个特征数据点，将该自由曲面b上的多个特征数据点进行曲面拟合，得到所述自由曲面b；

步骤S5，以此类推，直到获得待设计自由曲面成像光学***中所有自由曲面，进而得到一第一自由曲面成像光学***；

步骤S6，将所述第一自由曲面成像光学***再次作为初始***，经过多次迭代，得到一第二自由曲面成像光学***；

步骤S7，将所述第二自由曲面成像光学***作为一第一参考***，选定该第一参考***中的一个曲面位置自由度，并将其定义为第一曲面位置自由度，在第一曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第三自由曲面成像光学***；将每个“有扰动”的第三自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第四自由曲面成像光学***；挑选出满足所述约束条件且在评价指标下最优的第四自由曲面成像光学***；

S8：选取一第二参考***，选定第二参考***中的一个曲面位置自由度，并将其定义为第二曲面位置自由度，该第二曲面位置自由度与所述第一曲面位置自由度不同；在第二曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第五自由曲面成像光学***；将每个“有扰动”的第五自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第六自由曲面成像光学***；挑选出满足所述约束条件，且在评价指标下最优的第六自由曲面成像光学***；以及

S9：以此类推，直到所有曲面位置自由度均被使用过，得到所述待设计自由曲面成像光学***。

相较于现有技术，本发明提供的自由曲面成像光学***的设计方法，可以仅使用由简单平面组成的初始***作为设计的输入，在给定预期的物像关系或光线映射关系的情况下，给定约束条件以后，实现高成像质量自由曲面***的自动设计，克服了现有的光学***设计方法中初始结构匮乏的问题，而且方法不需使用光学设计软件优化，大大减少了人力。

附图说明

图1为本发明实施例提供的自由曲面成像光学***设计方法的流程图。

图2为本发明实施例提供的自由曲面成像光学***中每个视场中特征光线选择方法示意图。

图3为本发明实施例提供的求解特征数据点时特征光线起点与终点示意图。

图4为本发明实施例提供的自由曲面成像光学***的迭代方式。

具体实施方式

下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。

请参阅图1，本发明提供一种自由曲面成像光学***的设计方法，其包括以下步骤：

步骤S1，建立一初始***，该初始***包括多个初始曲面，且该初始***中的一个初始曲面对应待设计自由曲面成像光学***中一个自由曲面；

步骤S2，根据待设计自由曲面成像光学***的设计要求，建立自由曲面成像光学***的约束条件；

步骤S3，将待设计自由曲面成像光学***中的一个自由曲面定义为自由曲面a，保持所述多个初始曲面不变，选取K条特征光线R_i(i＝1，2…K)，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面a上的多个交点，进而得到自由曲面a上的多个特征数据点，将该自由曲面a上多个特征数据点进行曲面拟合，得到所述自由曲面a；

步骤S4，将待设计自由曲面成像光学***中的另一个自由曲面定义为自由曲面b，保持所述自由曲面a以及自由曲面a对应的初始曲面之外的其它初始曲面不变，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面b的多个交点，进而得到所述自由曲面b上的多个特征数据点，将该自由曲面b上的多个特征数据点进行曲面拟合，得到所述自由曲面b；

步骤S5，以此类推，直到获得待设计自由曲面成像光学***中所有自由曲面，进而得到一第一自由曲面成像光学***；

步骤S6，将所述第一自由曲面成像光学***再次作为初始***，经过多次迭代，得到一第二自由曲面成像光学***；

步骤S7，将所述第二自由曲面成像光学***作为一第一参考***，选定该第一参考***中的一个曲面位置自由度，并将其定义为第一曲面位置自由度，在第一曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第三自由曲面成像光学***；将每个“有扰动”的第三自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第四自由曲面成像光学***；挑选出满足所述约束条件且在评价指标下最优的第四自由曲面成像光学***；

S8：选取一第二参考***，选定第二参考***中的一个曲面位置自由度，并将其定义为第二曲面位置自由度，该第二曲面位置自由度与所述第一曲面位置自由度不同；在第二曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第五自由曲面成像光学***；将每个“有扰动”的第五自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第六自由曲面成像光学***；挑选出满足所述约束条件且在评价指标下最优的第六自由曲面成像光学***；以及

S9：以此类推，直到所有曲面位置自由度均被使用过，得到所述待设计自由曲面成像光学***。

步骤S1中，所述多个初始曲面可以为平面、球面等。所述多个初始曲面的具***置根据待设计自由曲面成像光学***的实际需要进行选择。

步骤S2中，所述约束条件可以为控制光线遮拦，消除曲面间位置干涉，控制自由曲面之间的位置关系等。

步骤S3中，所述K条特征光线的选取可通过以下方法进行：

根据需求选取M个视场，并将每个视场的孔径分成N等份，并从每一等份中选取不同孔径位置的P条特征光线，这样一共选取了K＝M×N×P条对应不同视场不同孔径位置的特征光线。所述孔径可以为圆形、长方形、正方形、椭圆形或其他规则或不规则的形状。请参阅图2，优选的，所述视场孔径为圆形孔径，将每个视场的圆形孔径等分成N个角度，间隔为φ，因此有

沿着每个角度的半径方向取P个不同的孔径位置，那么一共取K＝M×N×P条对应不同视场不同孔径位置的特征光线。

请一并参阅图3，为了得到一个待求的自由曲面Ω上的所有特征数据点P_i(i＝1，2…K)，将借助特征光线R_i(i＝1，2…K)与待求的自由曲面Ω的前一个曲面Ω′及后一个曲面Ω″的交点。在求解每条特征光线R_i(i＝1，2…K)对应的待求的自由曲面Ω上的特征数据点P_i(i＝1，2…K)时，将该特征光线R_i与前一个曲面Ω′的交点定义为该特征光线的起点S_i，特征光线R_i与后一个曲面Ω″的交点定义为该特征光线的终点E_i。当待设计的***与特征光线确定后，该特征光线R_i的起点S_i是确定的，且易于通过光线追迹即物像关系得到，特征光线的终点E_i可通过物像关系求解。在理想状态下，特征光线R_i从Ω′上的S_i射出后，经过Ω上的P_i，交于Ω″上的E_i，并最终交目标面于其理想目标点T_i，ideal。如果Ω″是目标面，特征光线的终点就是其目标点T_i；如果在待求自由曲面Ω和目标面之间还有其他面，特征光线的终点是Ω″上从特征数据点到其对应的目标点的光程的一阶变分为零的点，即

所述自由曲面a上特征数据点P_i(i＝1，2…K)可以通过以下两种计算方法获得。

第一种计算方法包括以下步骤：

步骤a，取定一条特征光线R₁与所述自由曲面a对应的初始曲面的第一交点为特征数据点P₁；

步骤b，在得到i(1≤i≤K-1)个特征数据点P_i后，根据斯涅尔定律的矢量形式求解该特征数据点P_i处的单位法向量

步骤c，过所述i(1≤i≤K-1)个特征数据点P_i分别做一第一切平面，得到i个第一切平面，该i个第一切平面与其余K-i条特征光线相交得到i×(K-i)个第二交点，从该i×(K-i)个第二交点中选取出与所述i(1≤i≤K-1)个特征数据点P_i距离最近的一个第二交点，作为所述自由曲面a的下一个特征数据点P_i+1；

步骤d，重复步骤b和c，直到计算得到自由曲面a上的所有特征数据点P_i(i＝1，2…K)，通过曲面拟合可以得到所述自由曲面a的方程式。

步骤b中，每个特征数据点P_i处的单位法向量

可以根据斯涅尔(Snell)定律的矢量形式求解。当待求的自由曲面Ω为折射面时，则每个特征数据点P_i(i＝1，2…K)处的单位法向量满足：

其中，

分别是沿着光线入射与出射方向的单位矢量，n，n′分别为待求的自由曲面Ω前后两种介质的折射率。

类似的，当待求的自由曲面Ω为反射面时，则每个特征数据点P_i(i＝1，2…K)处的单位法向量

满足：

由于，所述特征数据点P_i(i＝1，2…K)处的单位法向量

与所述特征数据点P_i(i＝1，2…K)处的切平面垂直。故，可以得到特征数据点P_i(i＝1，2…K)处的切平面。

所述第一种算法的计算复杂度为

当设计中采用的特征光线的数量较多时，该方法需要较长的计算时间。

第二种计算方法包括以下步骤：

步骤a′，取定一第一条特征光线R₁与所述自由曲面a对应的初始曲面的第一交点为特征数据点P₁；

步骤b′，在得到第i(1≤i≤K-1)个特征数据点P_i后，根据斯涅尔定律的矢量形式求解第i个特征数据点P_i处的单位法向量

进而求得P_i处的单位切向量

步骤c′，仅过所述第i(1≤i≤K-1)个特征数据点P_i做一第一切平面并与其余K-i条特征光线相交，得到K-i个第二交点，从该K-i个第二交点中选取出与所述第i个特征数据点P_i距离最短的第二交点Q_i+1，并将其对应的特征光线及与所述第i个特征数据点P_i的最短距离分别定义为R_i+1和D；

步骤d′，过特征数据点P_i(1≤i≤K-1)之前已求得的i-1个第一特征数据点分别做一第二切平面，得到i-1个第二切平面，该i-1个第二切平面与所述特征光线R_i+1相交得到i-1个第三交点，在每一第二切平面上每一第三交点与其所对应的特征数据点P_i形成一交点对，在所述交点对中，选出交点对中距离最短的一对，并将距离最短的交点对的第三交点和最短距离分别定义为Q_(i+1)′和D_i′；

步骤e′，比较D_i与D_i′，如果D_i≤D_i′，则把Q_i+1取为下一个特征数据点P_i+1，反之，则把Q_(i+1)′取为下一个特征数据点P_i+1；以及

步骤f′，重复步骤b′到e′，直到计算得到自由曲面a上的所有特征数据点P_i(i＝1，2…K)，通过曲面拟合可以得到所述自由曲面a对应的的方程式。

步骤b′中，每个特征数据点P_i处的单位法向量的计算方法与所述第一种算法步骤b中相同。

所述第二种算法的计算复杂度为

当设计中个采用的特征光线的数量较多时，所述第二种算法明显比第一种算法的计算复杂度小。优选的，采用所述第二种算法逐点求解特征数据点P_i(i＝1，2…K)。

所述自由曲面a的方程式包括一二次曲面项和一自由曲面项，该自由曲面a的方程式可表达为：

其中，

是二次曲面项，c是二次曲面顶点处的曲率，k是二次曲面系数；

是自由曲面项，Aj为每项对应的系数，该自由曲面项可以为xy多项式项，Chebyshev多项式，或Zernike多项式项等。

将所述初始***所在的空间定义一第一三维直角坐标系。优选的，将光束传播方向定义为z轴，垂直于z轴方向的平面为xoy平面。

将所述自由曲面a上的多个特征数据点P_i(i＝1，2…K)进行曲面拟合包括以下步骤：

S31：将所述多个特征数据点P_i(i＝1，2…K)在所述第一直角坐标系中拟合成一球面，得到所述球面的曲率c及其对应的曲率中心(x_c，y_c，z_c)，

S32：将中心采样视场主光线对应的特征数据点(x_o，y_o，z_o)定义为球面的顶点，并以该球面的顶点为原点，过曲率中心与球面顶点的直线为z轴，建立一第二三维直角坐标系；

S33：将所述多个特征数据点P_i(i＝1，2…K)在第一三维直角坐标系中的坐标(x_i，yi，z_i)及其对应的法向量(α_i，β_i，γ_i)变换为第二三维直角坐标系中的坐标(x′_i，y′_i，z′_i)及其法向量(α′_i，β′_i，γ′_i)；

S34：根据所述多个特征数据点P_i(i＝1，2…K)在第二三维直角坐标系中的坐标(x′_i，y′_i，z′_i)，以及步骤S42中求得的球面的曲率c，将特征数据点P_i(i＝1，2…K)在第二三维直角坐标系中拟合成一个二次曲面，得到二次曲面系数k；以及

S35：将所述多个特征数据点P_i(i＝1，2…K)在第二三维直角坐标系中的二次曲面上的坐标与法向量分别从坐标(x′_i，y′_i，z′_i)和法向量(α′_i，β′_i，γ′_i)中排除掉，得到残余坐标与残余法向量，将该残余坐标与残余法向量拟合成一个自由曲面，该自由曲面的方程式与所述二次曲面的方程式相加即可得到所述自由曲面a的方程式。

通常的，光学***关于三维直角坐标系的yoz面对称，因此，步骤S41中，所述球面在第二三维直角坐标系yoz面内相对于在第一三维直角坐标系yoz面内的倾斜角θ为：

所述多个特征数据点P_i(i＝1，2…K)在第二三维直角坐标系中的坐标(x′_i，y′_i，z′_i)与法向量(α′_i，β′_i，γ′_i)与在第一三维直角坐标系中的坐标(x_i，y_i，z_i)和法向量(α_i，β_i，γ_i)的关系式分别为：

和

在第二三维直角坐标系中，将在二次曲面上的特征数据点的坐标与法向量分别定义为(x′_i，y′_i，z′_is)和(α′_is，β′_is，γ′_is)。将法向的z′分量归一化为-1，将在二次曲面上的特征数据点坐标(x′_i，y′_i，z′_is)与法向量(α′_is，β′_is，γ′_is)分别从坐标(x′_i，y′_i，z′_i)和法向量(α′_i，β′_i，γ′_i)排除掉之后，得到残余坐标(x″_i，y″_i，z″_i)与残余法向量(α″_i，β″_i，-1)分别为：

(x_i″，y_i″，Z_i″)＝(x_i′，y_i′，Z_i′-Z_is′)和

步骤S35中，所述将残余坐标与残余法向量拟合得到自由曲面的步骤包括：

S351：在所述第二三维直角坐标系中，用排除掉二次曲面项的多项式曲面作为待构建自由曲面的面形描述，即

其中，g_j(x，y)为多项式的某一项，P＝(p₁，p₂，L，p_J)^T为系数集合；

S352：获得残余坐标拟合误差d₁(P)，即所述残余坐标值(x″_i，y″_i，z″_i)(i＝1，2，…，n)与所述自由曲面在z′轴方向残余坐标差值的平方和；以及残余法向矢量拟合误差d₂(P)，即所述残余法向量N_i＝(α″_i，β″_i，-1)(i＝1，2，…，n)与所述自由曲面法向量的矢量差的模值的平方和，

其中，

Z＝(z₁，z₂，L，z_I)^T，U＝(u₁，u₂，L，u_I)^T，V＝(v₁，v₂，L，v_I)^T，

S353：获得评价函数其中，w为权重且大于0；

S354：选择不同的权重w，并令所述评价函数f(P)的梯度

从而获得多组不同的P及其对应的多个自由曲面面形z＝f(x，y；P)；以及

S355：获得具有最佳的成像质量的最终自由曲面面形。

步骤S4中，自由曲面b上多个特征数据点的求解方法与步骤S3中自由曲面a上多个特征数据点的求解方法相同，将自由曲面b上多个特征数据点进行曲面拟合的方法与步骤S3中将自由曲面a上多个特征数据点进行曲面拟合的方法也相同。自由曲面b的方程式与自由曲面a的方程式也相同。

步骤S5中，其余的待求的反射镜的反射面的求解方法也与步骤S3中自由曲面a的求解方法相同。将其余的待求的反射镜的反射面上多个特征数据点进行曲面拟合的方法与步骤S3中将自由曲面a上多个特征数据点进行曲面拟合的方法也相同。其余的待求的反射镜的反射面的方程式与自由曲面a的方程式也相同。

所述待求的反射镜的反射面的数量可以与步骤S1中初始曲面的数量相同，也可以小于步骤S1中初始曲面的数量。本实施例中，所述待求的反射镜的反射面的数量为3个。

步骤S6中，所述迭代过程可以通过以下两种方法进行：

第一种方法包括以下步骤：

S61：将所述K条特征光线与步骤S3中得到的自由曲面a的多个交点作为特征数据点，分别求解每个特征数据点处的法向量，按照步骤S3中的方法将该多个特征数据点进行曲面拟合，得到一个自由曲面a′；

S62：将所述K条特征光线与步骤S4中得到的自由曲面b的多个交点作为特征数据点，分别求解每个特征数据点处的法向量，按照步骤S4中的方法将该多个特征数据点进行曲面拟合，得到一个自由曲面b′；

S63：以此类推，直至步骤S5中所有的待求自由曲面均重新求解完成，完成第一次迭代；

S64：以第一次迭代之后得到的自由曲面成像光学***再次作为初始***，重复步骤S61到S63，经过多次迭代后，所述K条特征光线与目标面的实际交点与理想目标点比较接近，迭代过程完成。

第二种方法包括以下步骤：将步骤S5中得到的自由曲面成像光学***作为初始***，重复步骤S3到S5，完成第一次迭代；将第一次迭代中得到的自由曲面成像光学***再次作为初始***，重复步骤S3到S5，完成第二次迭代；以此类推，经过多次迭代后，所述K条特征光线与目标面的实际交点与理想目标点比较接近，迭代过程完成。

优选的，所述迭代过程采用所述第一种方法进行，该方法比较简单，且节省时间。

请参阅图4，所述迭代可以采用正常迭代，负反馈，或逐步逼近等迭代方式。当采用正常迭代时，目标点T_i为理想目标点T_i，ideal，即T_i＝T_i，ideal。

当采用负反馈时，负反馈函数为：

其中，ε是负反馈系数，ε＞0；T_i ^*是本次迭代之前(上次迭代之后)特征光线与目标面的实际交点；Δ为负反馈阈值，Δ＞0，采用负反馈阀值有利于避免因T_i ^*和T_i，ideal差别太大，使T_i变化较大而引起的迭代不稳定；负反馈迭代方式速度比较快。

当采用逐步逼近迭代方式时T_i＝T_i ^*+ρ(T_i，ideal-T_i ^*)，其中，ρ是逼近系数，ρ＞0。逐步逼近迭代方式比较稳定。

所述自由曲面成像光学***的设计方法中待求的自由曲面的求解顺序不限，可以根据实际需要进行调换。

步骤S7中，将上述得到的第二自由曲面成像光学***作为一第一参考***，选定该第一参考***中的一个曲面位置自由度，所述第一曲面位置自由度可以为第一参考***中某个自由曲面或像面的x方向位移，y方向位移，z方向位移，绕x轴旋转角度，绕y轴旋转角度，绕z轴旋转角度中的任意一个。

在第一曲面位置自由度的值的附近选取多个值，改变对应曲面(包括自由曲面或像面)的位置，得到多个相对于第一参考***有一定变化的***。将这些***定义为一系列“有扰动”的第三自由曲面成像光学***。

如果将每个“有扰动”的第三自由曲面成像光学***作为初始***，仅重复进行步骤S6时，第一轮迭代不能使用负反馈模式。

所述评价指标可以是特征光线与目标面的实际交点与理想目标点的偏差的均方根(RMS)值σ_RMS，该σ_RMS的表达式为：

其中，σ_m是第m条特征光线与目标面的实际交点与理想目标点的距离。

步骤S8中，如果满足所述结构约束条件，且在评价指标下最优的第四自由曲面成像光学***的σ_RMS小于所述第一参考***的σ_RMS，该最优的第四自由曲面成像光学***作为所述第二参考***；反之，第一参考***作为第二参考***。

在第二曲面位置自由度的值的附近选取多个值，改变对应曲面(包括自由曲面或像面)的位置，得到多个相对于第二参考***有一定变化的***。将这些***定义为一系列“有扰动”的第五自由曲面成像光学***。

如果将每个“有扰动”的第五自由曲面成像光学***作为初始***，仅重复进行步骤S6时，第一轮迭代不能使用负反馈模式。

可以理解，在S7-S9的整个过程中，所述约束条件也可以变化。例如：在早期过程中，约束条件可以限制的松一些；在后期过程中，约束条件逐渐变得严格。在某一实例中，约束条件在步骤S7和S8中比较松，在步骤S9中比较严格。

在S7-S9的后期过程中，可以使用比步骤S3中更多的特征光线，即特征光线的数量大于K。自由曲面拟合使用的方程式的最高次数也可以提升。在某一实例中，在步骤S9中，特征光线的数量大于K，且自由曲面拟合使用的方程式的最高次数也进行了提升。

在整个设计过程前期，自由曲面随设计过程面形变化比较大。如果自由曲面成像光学***的孔径光阑不在该光学***的最前端，自由曲面面形的较大变形可能会引起入瞳位置发生较大的变化，进而使参与成像的光线发生变化。因此，可以在整个设计过程前期，将自由曲面成像光学***的入瞳直径选取的适当大一些，来减小入瞳位置变化的影响。在整个设计过程后期，由于自由曲面的面形变化不大，自由曲面成像光学***入瞳直径可恢复正常。在某一实施例中，自由曲面成像光学***的孔径光阑不在该光学***的最前端，在步骤S1-S8中，将自由曲面成像光学***的入瞳直径选取的适当大一些，在步骤S9中，自由曲面成像光学***入瞳直径恢复正常。

对于一般的自由曲面成像光学***，对***像质的要求更高，对***畸变的要求不高。因此，在S7-S9的后期过程中，可以选用“像质优先”的模式。在“像质优先”模式下，每条特征光线在构建与迭代的过程中的理想目标点不再是预设的理想像点，而是自由曲面成像光学***每轮迭代之前对应的视场的主光线与目标面的实际交点。此时，可以用评价指标σ^* _RMS来评价构建与迭代过程的效果，

其中，

是第m条特征光线与目标面的实际交点与自由曲面成像光学***每轮迭代之前对应的视场的主光线与目标面的实际交点的距离。

所述自由曲面成像光学***的设计方法可进一步包括一对步骤S9中得到的自由曲面成像光学***进行进一步设计的步骤。具体地，可将S9中得到的待设计自由曲面成像光学***作为新的参考***，至少一次重复进行S7-S9的步骤，得到最终的自由曲面成像光学***。可以理解，该对步骤S9中得到的进行进一步设计的步骤并不是必需的，可以根据实际需要设计。

本发明提供的自由曲面成像光学***设计方法具有以下优点：其一，在将特征数据点进行曲面拟合时，同时考虑了特征数据点的坐标和法向量，使得到的自由曲面的面形更加准确，且提高了自由曲面成像光学***的光学性能；其二，通过建立两个三维直角坐标系，在第一三维直角坐标系中拟合得到二次曲面项后，将特征数据点的坐标和法向量变换到第二直角坐标系中的坐标和法向量，进而拟合得到自由曲面项，该拟合方法比较精确；其三，将逐点构建得到的自由曲面成像光学***再次作为初始***，经过多次迭代，直到特征光线与目标面的实际交点与理想目标点接近，进而提高自由曲面成像光学***的成像质量，降低自由曲面成像光学***的畸变和平均RMS弥散斑直径，方法简单。其四，可以仅使用由简单平面组成的***作为设计的输入，在给定预期的物像关系或光线映射关系，以及一定的结构约束条件以后，可以实现高成像质量自由曲面成像***的自动设计，克服了传统光学设计初始结构匮乏的问题。而且该方法不需使用光学设计软件优化，大大减少了人力。

另外，本领域技术人员还可在本发明精神内作其它变化，当然这些依据本发明精神所作的变化，都应包含在本发明所要求保护的范围内。

Claims (10)

1.一种自由曲面成像光学***的设计方法，其包括以下步骤：

步骤S1，建立一初始***，该初始***包括多个初始曲面，且该初始***中的一个初始曲面对应待设计自由曲面成像光学***中一个自由曲面；

步骤S2，根据待设计自由曲面成像光学***的设计要求，建立自由曲面成像光学***的约束条件；

步骤S3，将待设计自由曲面成像光学***中的一个自由曲面定义为自由曲面a，保持所述多个初始曲面不变，选取K条特征光线R_i(i＝1，2…K)，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面a上的多个交点，进而得到自由曲面a上的多个特征数据点，将该自由曲面a上多个特征数据点进行曲面拟合，得到所述自由曲面a；

步骤S4，将待设计自由曲面成像光学***中的另一个自由曲面定义为自由曲面b，保持所述自由曲面a以及自由曲面a对应的初始曲面之外的其它初始曲面不变，根据物像关系及斯涅尔定律逐点求解所述K条特征光线与自由曲面b的多个交点，进而得到所述自由曲面b上的多个特征数据点，将该自由曲面b上的多个特征数据点进行曲面拟合，得到所述自由曲面b；

步骤S5，以此类推，直到获得待设计自由曲面成像光学***中所有自由曲面，进而得到一第一自由曲面成像光学***；

步骤S6，将所述第一自由曲面成像光学***再次作为初始***，经过多次迭代，得到一第二自由曲面成像光学***；

步骤S7，将所述第二自由曲面成像光学***作为一第一参考***，选定该第一参考***中的一个曲面位置自由度，并将其定义为第一曲面位置自由度，在第一曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第三自由曲面成像光学***；将每个“有扰动”的第三自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第四自由曲面成像光学***；挑选出满足所述约束条件且在评价指标下最优的第四自由曲面成像光学***；

S8：选取一第二参考***，选定第二参考***中的一个曲面位置自由度，并将其定义为第二曲面位置自由度，该第二曲面位置自由度与所述第一曲面位置自由度不同；在第二曲面位置自由度的值的附近选取多个值，改变对应曲面的位置，得到多个“有扰动”的第五自由曲面成像光学***；将每个“有扰动”的第五自由曲面成像光学***作为初始***，重复进行步骤S3-S6或步骤S6，得到多个第六自由曲面成像光学***；挑选出满足所述约束条件且在评价指标下最优的第六自由曲面成像光学***；以及

S9：以此类推，直到所有曲面位置自由度均被使用过，得到所述待设计自由曲面成像光学***。

2.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，所述约束条件包括控制光线遮拦、消除曲面间位置干涉以及控制自由曲面之间的位置关系。

3.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，所述评价指标为特征光线R_i(i＝1，2…K)与目标面的实际交点与理想目标点的偏差的均方根值σ_RMS，该σ_RMS的表达式为：

其中，σ_m是第m条特征光线与目标面的实际交点与理想目标点的距离。

4.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，所述第一曲面位置自由度为第一参考***中自由曲面或像面的x方向位移，y方向位移，z方向位移，绕x轴旋转角度，绕y轴旋转角度，绕z轴旋转角度中的任意一个。

5.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，步骤S7中，如果将每个“有扰动”的第三自由曲面成像光学***作为初始***，仅重复进行步骤S6时，第一轮迭代不能使用负反馈模式。

6.如权利要求3所述的自由曲面成像光学***的设计方法，其特征在于，步骤S8中，满足所述约束条件，且在评价指标下最优的第四自由曲面成像光学***的σ_RMS小于所述第一参考***的σ_RMS，该最优的第四自由曲面成像光学***作为所述第二参考***。

7.如权利要求3所述的自由曲面成像光学***的设计方法，其特征在于，满足所述约束条件，且在评价指标下最优的第四自由曲面成像光学***的σ_RMS大于等于所述第一参考***的σ_RMS，第一参考***作为第二参考***。

8.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，进一步包括一对步骤S9中得到的自由曲面成像光学***进行进一步设计的步骤，具体地，将S9中得到的待设计自由曲面成像光学***作为新的参考***，至少一次重复S7-S9的步骤。

9.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，所述自由曲面a上的多个特征数据点的求解包括以下步骤：

步骤a′，取定一第一条特征光线R₁与所述自由曲面a对应的初始曲面的第一交点为特征数据点P₁；

步骤b′，在得到第i(1≤i≤K-1)个特征数据点P_i后，根据斯涅尔定律的矢量形式求解第i个特征数据点P_i处的单位法向量

进而求得P_i处的单位切向量

步骤c′，仅过所述第i(1≤i≤K-1)个特征数据点P_i做一第一切平面并与其余K-i条特征光线相交，得到K-i个第二交点，从该K-i个第二交点中选取出与所述第i个特征数据点P_i距离最短的第二交点Q_i+1，并将其对应的特征光线及与所述第i个特征数据点P_i的最短距离分别定义为R_i+1和D_i；

步骤d′，过特征数据点P_i(1≤i≤K-1)之前已求得的i-1个第一特征数据点分别做一第二切平面，得到i-1个第二切平面，该i-1个第二切平面与所述特征光线R_i+1相交得到i-1个第三交点，在每一第二切平面上每一第三交点与其所对应的特征数据点P_i形成一交点对，在所述交点对中，选出交点对中距离最短的一对，并将距离最短的交点对的第三交点和最短距离分别定义为Q_(i+1)′和D_i′；

步骤e′，比较D_i与D_i′，如果D_i≤D_i′，则把Q_i+1取为下一个特征数据点P_i+1，反之，则把Q_(i+1)′取为下一个特征数据点P_i+1；以及

步骤f′，重复步骤b′到e′，直到计算得到自由曲面a上的所有特征数据点P_i(i＝1，2…K)。

10.如权利要求1所述的自由曲面成像光学***的设计方法，其特征在于，所述进行多次迭代过程包括以下步骤：

S61：将所述K条特征光线与步骤S3中得到的自由曲面a的多个交点作为特征数据点，分别求解每个特征数据点处的法向量，按照步骤S3中的方法将该多个特征数据点进行曲面拟合，得到一个自由曲面a′；

S62：将所述K条特征光线与步骤S4中得到的自由曲面b的多个交点作为特征数据点，分别求解每个特征数据点处的法向量，按照步骤S4中的方法将该多个特征数据点进行曲面拟合，得到一个自由曲面b′；

S63：以此类推，直至步骤S5中所有的自由曲面均重新求解完成，完成第一次迭代；

S64：以第一次迭代之后得到的自由曲面成像光学***再次作为初始***，重复步骤S61到S63，经过多次迭代后，得到所述第二自由曲面成像光学***，迭代过程完成。

Images