引　言

对光学系统而言，像差是影响其成像质量的主要因素^[1]，然而光学系统的像差通常会随着光束打在光学面上的入射角的增大变得愈发严重。随着非球面在各类光学系统中的广泛应用，应用非球面可以有效控制像差、改善光学系统成像性能也成为了人们的普遍认识^[2-3]。Miks等^[4]提出了一种包含一个或两个非球面光学表面的镜头的简单设计方法，在指定物、像位置上可以有效校正球差并减小彗差的贡献。非球面不仅在轴对称光学系统中被广泛应用，在众多的鱼眼镜头优化设计^[5-6]过程中也被广泛应用。张宝龙等采用非球面技术，用仅四片透镜就实现视场角为210°的超广角鱼眼镜头设计要求^[7]。然而，应用像差分析的方法在光学系统中选择最适非球面位置及确定非球面系数初值的方法鲜有讨论。

对于鱼眼镜头这类超大视场光学系统，其在很多场合例如无人驾驶、安全监控等^[8-12]方面都有着广泛的应用。但由于在其工作环境中，物点发出的光束会以很大的入射角进入光学系统，这种掠入射光束经光学系统成像后，在子午和弧矢平面内的聚焦位置和波阵面参数可能完全不一致，波前的形状将严重偏离球面，具有平面对称光学系统的成像特性。赛德尔理论只能用于常规的轴对称光学系统的近轴像差分析，不再适合此类系统。近年来，LU等采用超环面作为参考波阵面，发展了适用于平面对称光学系统的波像差理论^[13-14]，可以被用来分析包括非球面系数在内的光学系统的结构参数对成像性能的影响，从而为寻找最适非球面奠定理论基础。

本文首先阐述了LU的平面对称光学系统像差理论的理论基础，接着阐述鱼眼镜头光学系统中各光学面上各类波像差分布的计算方法，然后针对一个全球面鱼眼镜头实例，根据所建立的系统成像性能评价函数随非球面系数改变的变化规律，确定非球面位置以及非球面系数的初值，随后重新计算光学系统的波像差分布，并与原先系统的波像差分布图进行对比，最后应用基于调制传递函数（MTF）作为评价函数的自适应归一化实数编码遗传算法，优化该鱼眼镜头实例。

1 理论基础 1.1 LU的波像差理论

LU应用波像差方法是针对平面对称光学系统在子午和弧矢方向上焦点可能完全分离的特点，以超环面作为参考波阵面，发展出的适用于平面对称多元件光栅系统的像差理论，并根据光线几何，通过多项式拟合的方法求得光线在光学元件表面和入射及像差波阵面之间的映射关系，导出精确的平面对称光学系统的波像差表达式，再由波像差求得光学系统的像差表达式。该理论可以精确处理一般的平面对称多元件光学系统的像差，并且可以应用于研究各种面形的反射镜、光栅和像平面在任意位置上的光学系统的成像。

由文献[13]可以知道，平面对称光栅系统的波像差 $ W $ 与物方空间波像差 $ {W}_{\mathrm{o}\mathrm{b}\mathrm{j}} $ 、像方空间波像差 $ {W}_{\mathrm{i}\mathrm{m}\mathrm{a}} $ 以及受光栅刻槽影响的光程差 $ S $ 相关，具体表达式如下：

$ W = {W_{{\text{obj}}}} + {W_{{\text{ima}}}} - S $

(1)

式中： $ S=nm\lambda $ ， $ \lambda $ 为光的波长， $ n $ 为光栅刻槽的槽数； $ m $ 为光栅的衍射级， $ m=0,\pm 1,\pm 2,\cdots $ 。可将式（1）右边写成四阶波像差系数 $ {w}_{ijk} $ 的形式：

$ W = \sum\limits_{ijk}^4 {{w_{ijk}}} {x^i}{y^j}{u^k},\left( {i + j + k \leqslant 4} \right) $

(2)

式中： $ (x,y) $ 为光线打在光学面的坐标； $ u $ 为物点到光学面连线与光轴的夹角。关于波像差系数 $ {w}_{ijk} $ 以及具体应用可参见2.1节内容。

LU波像差理论只涉及到孔径像差如彗差波像差、球差波像差的计算，而对场曲波像差及色差波像差的计算方法未有深入研究。该理论在近年被FAN^[15]等补充，使得该理论可以利用解析的方法精确计算出光学系统中除畸变外的所有像差，为寻找最适非球面奠定了理论基础。

1.2 光学表面的面型系数

在超大视场光学系统（如鱼眼镜头）中将光学元件表面定义为二次圆锥表面，应用平面对称像差理论计算波像差时需要以主光线打在光学元件表面的位置为坐标原点O，这样就需要转换坐标。图1所示为一条主光线经过一个二次圆锥光学面反射后的光路图：将坐标方程从以 ${O}{\rm{'}}$ 为坐标原点的 ${x}{{\rm{'}}}{y}{\rm{'}}{z}{\rm{' }}$ 坐标系转换至以 $ O $ 点为坐标原点 $ xyz $ 坐标轴^[16]； $ y $ 轴和 ${y}{\rm{'}}$ 轴均垂直于纸面向外； $ \alpha $ 和 $ \beta $ 分别表示入射角和反射角； $ {\omega }_{i-1} $ 和 $ {\omega }_{i} $ 分别为第i光学面的入射和出射光线与光轴的夹角，即视场角，其中， $ i=\mathrm{1,2},\cdots $ 。

以 $ {O}' $ 为坐标原点的坐标系 $ {x}'{y}'{z}' $ 绕光轴 $ {z}' $ 旋转的二次圆锥曲面方程为：

$ {{x}^{\prime }}^{2}+{{y}^{\prime }}^{2}={a}_{1}{z}^{\prime }+{a}_{2}{{z}^{\prime }}^{2},\; {a}_{1}=2{R}_{0} $

(3)

式中： $ {R}_{0} $ 是二次圆锥曲面在坐标原点 $ {O}' $ 处的曲率半径，即 $ {R}_{0}=\overline{{C}_{1}{O}'} $ ； $ {a}_{2} $ 是确定二次圆锥曲面类型的系数，即非球面系数，当 $ {a}_{2} ＞ 0 $ 时为双曲面，当 $ {a}_{2}=0 $ 时为抛物面，当 $ -1{ ＜ a}_{2} ＜ 0 $ 时为长椭球面， $ {a}_{2}=-1 $ 时为球面， $ {a}_{2} ＜ -1 $ 为扁椭球面。 $ {a}_{2} $ 与Zemax软件中表示非球面特性的圆锥系数 $ k $ 的转换关系为：

$ {a_2} = - 1 - k $

(4)

设 $ O $ 点在平面坐标系 $ {x}'{O}'{z}'$ 下的坐标为 $ \left({x}_{0}^{\mathrm{*}},{z}_{0}^{\mathrm{*}}\right) $ ， $ x $ 轴为主光线相交于光学元件表面 $ O $ 处的切线， $ z $ 轴为 $ O $ 处的法线。根据图1所示，将坐标系 $ xyz $ 转换为 $ {x}'{y}'{z}'$ 的表达式为：

$ \left\{ \begin{gathered} x' = x\cos \theta + z\sin \theta + x_0^* \hfill \\ y' = y \hfill \\ z' = - x\sin \theta + z\cos \theta + z_0^* \hfill \\ \end{gathered} \right. $

(5)

同时，从图1的几何关系可以得到：

$ \left\{ \begin{gathered} \cos \theta = - \frac{{{a_1} + 2{a_2}z_0^ * }}{{\sqrt {{{\left( {{a_1} + 2{a_2}z_0^ * } \right)}^2} + 4x_0^{ * {\text{2}}}} }} \hfill \\ \sin \theta = \frac{{2x_0^ * }}{{\sqrt {{{\left( {{a_1} + 2{a_2}z_0^ * } \right)}^2} + 4x_0^{ * {\text{2}}}} }} \hfill \\ \end{gathered} \right. $

(6)

将式（5）和（6）带入式（3），并作四阶泰勒级数展开，就可以得到适用于四阶波像差理论的二次圆锥曲面的方程：

$ \begin{split} z & = {c_{2,0}}{x^2} + {c_{0,2}}{y^2} + {c_{3,0}}{x^3} + {c_{1,2}}x{y^2} + {c_{4,0}}{x^4}+ \\ & {c_{0,4}}{y^4} + {c_{2,2}}{x^2}{y^2} \end{split}$

(7)

式中， $ {c}_{i,j} $ 为面型系数，具体可参见文献[16]。

由图1的几何关系可知 $ O $ 点处光学元件表面的子午方向的曲率半径 ${\varGamma}$ 和弧矢方向的曲率半径 $ \rho $ 分别为：

$ {\varGamma} {\text{ = }} - \overline {{C_1}O'} = z_0^ * - \frac{{{x^ * }}}{{\tan \theta }} = \frac{{{a_1}}}{2} + \left( {1 + {a_2}} \right)z_0^ * $

(8)

$ \rho {\text{ = }} - \overline {{C_1}O} = - \frac{{\sqrt {4x{{_0^*}^2}(1 + {a_2}) + {a_1}^2} }}{2} = \frac{1}{{2{c_{0,2}}}} $

(9)

对于有由多个光学面组成的光学系统，从第 $ i $ 光学面到 $ i+1 $ 光学面的坐标传递方程 $ {s}_{0} $ 的表达式为：

$ {s_{0\left( {i + 1} \right)}} = z_{0\left( i \right)}^ * - {d_i} - \frac{{x_{0\left( i \right)}^*}}{{\tan {\omega _{0\left( {i + 1} \right)}}}} $

(10)

式中 $ {d}_{i} $ 是从第 $ i $ 光学面到第 $ i+1 $ 光学面在沿光轴的光学间隔。

通过以上推导内容可将一般平面对称二次圆锥曲面系数转换为LU的平面对称光学系统像差理论中的面型系数形式。

1.3 主光线传输方程

应用LU的平面对称光学系统波像差理论计算光学系统波像差时还需要追迹出各光学面的主光线参数，包括入射角 $ {\alpha }_{i} $ ，折射或反射角 $ {\beta }_{i} $ ，视场角 $ {\omega }_{i} $ 等参数。图2所示为一条主光线 $ \overline{{A}_{i}{O}_{i}{O}_{i+1}} $ 在子午平面内的传输示意图：依次穿过光学介质的折射率分别为 $ {n}_{i-1},{n}_{i} $ 和 $ {n}_{i+1} $ ； $ {\omega }_{i} $ 的符号与主光线相对光轴的斜率一致； $ {\alpha }_{i} $ , $ {\beta }_{i} $ 的符号与主光线相对于入射坐标系的法线斜率一致；第 $ i $ 到第 $ i+1 $ 光学面的光学间隔为 $ \overline{{D}_{i}{D}_{i+1}}={d}_{i} $ ；光学面右侧为正，反之为负；光学面曲率半径的符号规定为圆心在光学面右侧为正，反之为负。

由图2可得，主光线经过任一光学面的传输方程为^[16-17]：

$ \sin {\alpha _{i + 1}} = \frac{{{\varGamma _{i + 1}} + {d_i} - {\varGamma _i}}}{{{\rho _{i + 1}}}}\sin {\omega _i} + \frac{{{\rho _i}}}{{{\rho _{i + 1}}}}\sin {\beta _i} $

(11)

$ {\omega _i} = {\omega _{i - 1}} + {\beta _i} - {\alpha _i} = {\omega _0} + \sum\limits_{i = 1}^i {\left( {{\beta _i} - {\alpha _i}} \right)} $

(12)

$ {\beta _{i + 1}} = {\sin ^{ - 1}}\left( {\frac{{{n_i}}}{{{n_{i + 1}}}}\sin {\alpha _{i + 1}}} \right) $

(13)

式中： ${\rho _i} = \overline {{o_i}{c_i}} $ ； ${\varGamma _i} = - \overline {{C_i}{D_i}} $ ； ${\varGamma _{i + 1}} = \overline {{D_{i + 1}}{C_{i + 1}}} $ 。

上述理论及推导为计算鱼眼镜头中的波像差提供了理论基础。

2 波像差计算 2.1 孔径像差波像差

对于由g个光学面组成的超大视场平面对称光学系统，应用式（2）即可满足其计算精度，由于没有光栅，其表达式为：

$ \begin{array}{*{20}{c}} {W = \displaystyle\sum\limits_1^g {\displaystyle\sum\limits_{ij0}^4 {{w_{ij0}}} {x^i}{y^j}} ,}&{\left( {i + j \leqslant 4} \right)} \end{array} $

(14)

其中：

$ {w_{ij0}} = n{M_{ij0}}\left( {\alpha ,{r_{\text{M}}},{r_{\text{S}}},0} \right) + n'{M_{ij0}}\left( {\beta ,{{r'}_{\text{M}}},{{r'}_{\text{S}}},0} \right) $

(15)

其中： $ {M}_{ij0} $ 为波像差系数，参考文献[18]中有详细讨论； $ {r}_{\mathrm{M}\left(i\right)}'$ , $ {r}_{\mathrm{S}\left(i\right)}'$ 表示第 $ i $ 光学面在像空间的子午和弧矢焦距，根据准基理论 $ {w}_{200}=0 $ 和 $ {w}_{020}=0 $ 计算而来^[13,19]，表达式为：

$ \begin{gathered} {{r'}_{{\text{M}}\left( i \right)}} = \frac{{{{n'}_i}{{\cos }^2}{\beta _i}}}{{2{c_{2,0}}\left( {{n_i}\cos {\alpha _i} + {{n'}_i}\cos {\beta _i}} \right) - \dfrac{{{n_i}{{\cos }^2}{\alpha _i}}}{{{r_{{\text{M}}\left( i \right)}}}}}} \hfill \\ {{r'}_{{\text{S}}\left( i \right)}} = \frac{{{{n'}_i}}}{{2{c_{0,2}}\left( {{n_i}\cos {\alpha _i} + {{n'}_i}\cos {\beta _i}} \right) - \dfrac{{{n_i}}}{{{r_{{\text{s}}\left( i \right)}}}}}} \hfill \\ \end{gathered} $

(16)

第 $ i+1 $ 光学面沿主光线在子午方向的物方焦距 $ {r}_{\mathrm{M}(i+1)} $ 和弧矢方向的物方焦距 $ {r}_{\mathrm{S}(i+1)} $ 分别为：

$ \begin{array}{cc}{r}_{\text{M}\left(i+1\right)}={\overline{d}}_{i}-{{r}^{\prime }}_{\text{M}\left(i\right)}\text{，}& {r}_{\text{S}\left(i+1\right)}={\overline{d}}_{i}-{{r}^{\prime }}_{\text{S}\left(i\right)}\end{array} $

(17)

其中： $ \overline{{d}_{i}} $ 为第 $ i $ 到第 $ i+1 $ 光学面沿主光线的光学间隔，表达式为：

$ {\overline {{d_i}}} = \frac{{{\rho _i}\sin \left( {{\omega _{i - 1}} - {\alpha _i}} \right) - {\rho _{i + 1}}\sin \left( {{\omega _i} - {\alpha _{i + 1}}} \right)}}{{\sin {\omega _i}}} $

(18)

在四阶波像差理论中，任一光学面球差波像差 $ {W}_{\mathrm{s}\mathrm{p}\mathrm{h}} $ 和彗差波像差 $ {W}_{\mathrm{c}\mathrm{o}\mathrm{m}\mathrm{a}} $ 的计算表达式为：

$ \left\{ \begin{gathered} {W_{{\text{sph}}}} = {w_{400}}{x^4} + {w_{220}}{x^2}{y^2} + {w_{040}}{y^4} \hfill \\ {W_{{\text{coma}}}} = {w_{300}}{x^3} + {w_{120}}x{y^2} \hfill \\ \end{gathered} \right. $

(19)

2.2 主光线像差波像差

由g个光学面组成的光学系统在子午方向和弧矢方向的场曲波像差 ${W}_{\mathrm{M}\left({g}\right)}^{\mathrm{*}},{W}_{\mathrm{S}\left({g}\right)}^{\mathrm{*}}$ 分别由以下公式计算而来^[15]：

$ \left\{ \begin{gathered} W_{{\text{M}}\left( g \right)}^ * = \frac{{{{n'}_g}{{\cos }^2}{\beta _g}}}{2}\left( {\frac{1}{{{{r'}_{{\text{M}}\left( g \right)}}}} - \frac{1}{{{{r'}_{{\text{C}}\left( g \right)}}}}} \right)x_g^2 \hfill \\ W_{{\text{S}}\left( g \right)}^ * = \frac{{{{n'}_g}}}{2}\left( {\frac{1}{{{{r'}_{{\text{S}}\left( g \right)}}}} - \frac{1}{{{{r'}_{{\text{C}}\left( g \right)}}}}} \right)y_g^2 \hfill \\ \end{gathered} \right. $

(20)

式中： $\left({x}_{{g}},{y}_{{g}}\right)$ 为光线打在第g光学面上的坐标； ${r}_{\mathrm{M}\left({g}\right)}'$ 和 ${r}_{\mathrm{S}\left({g}\right)}'$ 分别为子午和弧矢方向上的像方焦距； ${r}_{\mathrm{C}\left({g}\right)}'$ 为第g光学面到高斯像面之间沿主光线的距离； ${n}_{{g}}'$ 为像方折射率。

由g个光学面组成光学系统的轴向和垂轴色差波像差 ${W}_{\mathrm{C}\mathrm{L}\left({g}\right)}^{\mathrm{*}},{W}_{\mathrm{C}\mathrm{T}\left({g}\right)}^{\mathrm{*}}$ ，计算分别如下^[15]：

$ W_{{\text{CL}}\left( g \right)}^ * = \frac{{{{n'}_g}}}{2}\left( {\frac{1}{{{{r'}_{{\text{F}}\left( g \right)}}}} - \frac{1}{{{{r'}_{{\text{C}}\left( g \right)}}}}} \right)x_{g\_0}^2 $

(21)

$ W_{{\text{CT}}\left( g \right)}^ * = \frac{{{{n'}_g}{x_{g\_0}}}}{{{{r'}_{0\left( g \right)}}}}\left( {{{y'}_{{\text{F}}\left( g \right)}}\cos {\omega _{{\text{F}}\left( g \right)}} - {{y'}_{{\text{C}}\left( g \right)}}\cos {\omega _{{\text{C}}\left( g \right)}}} \right) $

(22)

其中： ${r}_{\mathrm{F}\left({g}\right)}'$ 和 ${r}_{\mathrm{C}\left({g}\right)}'$ 分别为正入射时第g光学面F光和C光的像方焦距，即 $ \alpha =180 $ °和 $ \beta =0 $ °时应用式(18)求得； ${x}_{{g}\_0}$ 为光线在主光线与光学面交点的切面上的 $ x $ 坐标，即光线在光学面上坐标的线性近似； ${\omega }_{\mathrm{F}\left({g}\right)}$ 和 ${\omega }_{\mathrm{C}\left({g}\right)}$ 分别为F光和C光在第g光学面上的入射视场角； ${y}_{\mathrm{F}\left({g}\right)}'$ 和 ${y}_{\mathrm{C}\left({g}\right)}'$ 为F光和C光在像面上的投影高度可由下式计算^[6,20]：

$ {y'_g} = \tan {\omega _g}\left( {\frac{{{\rho _g}\sin {\beta _g}}}{{\sin {\omega _g}}} - {\Gamma _g} + {{r}'_{0\left( g \right)}}} \right) $

(23)

其中， ${r}_{0\left({g}\right)}'$ 表示第g光学面的高斯像距。

在式（6）～（23）中，( $ {x}_{i},{y}_{i} $ )为光束投影到第 $ i $ 光学面上的孔径坐标，由孔径光阑半径逐一线性传递所得：

$ \begin{array}{*{20}{c}} {{x_i} = {A_i}{x_{i + 1}},}&{{y_i} = {B_i}{y_{i + 1}}} \end{array} $

(24)

其中：

$ \begin{array}{*{20}{c}} {{A_i} = \dfrac{{{{r'}_{{\text{M}}\left( i \right)}}\cos {\alpha _{i + 1}}}}{{{r_{{\text{M}}\left( {i + 1} \right)}}\cos {\beta _i}}},}&{{B_i} = - \dfrac{{{{r'}_{{\text{S}}\left( i \right)}}}}{{{r_{{\text{S}}\left( {i + 1} \right)}}}}} \end{array} $

(25)

由波像差的叠加原理可知第g光学面的主光线波像差为：

$ \left\{\begin{array}{l}{W}_{\text{M}\left(g\right)}={W}_{\text{M}\left(g\right)}^{\ast }-{W}_{\text{M}\left(g-1\right)}^{\ast }\\ {W}_{\text{S}\left(g\right)}={W}_{\text{S}\left(g\right)}^{\ast }-{W}_{\text{S}\left(g-1\right)}^{\ast }\\ {W}_{\text{CL}\left(g\right)}={W}_{\text{CL}\left(g\right)}^{\ast }-{W}_{\text{CL}\left(g-1\right)}^{\ast }\\ {W}_{\text{CT}\left(g\right)}\text={W}_{\text{CT}\left(g\right)}^{\ast }-{W}_{\text{CT}\left(g-1\right)}^{\ast }\end{array} \right.$

(26)

应用式（26）可以计算光学系统中每一光学面的各类波像差的分布，并且可以得到系统波像差随非球面系数的变化规律，从而选出对非球面系数变化较为敏感的光学面。

3 利用非球面优化鱼眼镜头 3.1 定义评价函数

在视场范围内取 $ \mu $ 个物点对鱼眼镜头光学系统成像性能进行数值评价，该评价函数是将第2节所计算的波像差转换到像面上的像差 $ Q $ ，表达式为^[15,20]：

$ Q = \sum\limits_{k = 1}^\mu {{\varepsilon _k}\left( {Q_{x\left( k \right)}^2 + Q_{y\left( k \right)}^2 + {\eta _k}Q_{\eta \left( k \right)}^2 + {\mu _k}Q_{{\text{C}}\left( k \right)}^2} \right)} $

(27)

式中： $ {\varepsilon }_{k} $ 、 $ {\eta }_{k} $ 和 $ {\mu }_{k} $ 分别表示各自的权重系数，是根据具体光学系统像差分布决定的； $ {Q}_{x\left(k\right)} $ 和 $ {Q}_{y\left(k\right)} $ 表示第k视场角的孔径像差分量 $ k=\mathrm{0,} 1,2,\cdots $ 。分别在光束覆盖的区域范围内进行积分求和：

$ \left\{ \begin{gathered} Q_{x\left( k \right)}^2 = \frac{8}{{{\text{π }}{W_{\rm{q}}}L}}\int_{ - {W_q}/2}^{{W_{\rm{q}}}/2} {\int_0^{\frac{{\sqrt {1 - 4{x^2}/{W_{\rm{q}}}^2} }}{{L/2}}} {{{(x' - \overline {x'} )}^2}} {\text{d}}x{\text{d}}y} \hfill \\ Q_{y\left( k \right)}^2 = \frac{8}{{{\text{π }}{W_{\rm{q}}}L}}\int_{ - {W_{\rm{q}}}/2}^{{W_{\rm{q}}}/2} {\int_0^{\frac{{\sqrt {1 - 4{x^2}/{W_{\rm{q}}}^2} }}{{L/2}}} {y{'^2}{\text{d}}x{\text{d}}y} } \hfill \\ \end{gathered} \right. $

(28)

$ \overline {x'} = \frac{8}{{{\text{π }}{W_{\rm{q}}}L}}\int_{ - {W_{\rm{q}}}/2}^{{W_{\rm{q}}}/2} {\int_0^{\frac{{\sqrt {1 - 4{x^2}/{W_{\rm{q}}}^2} }}{{L/2}}} {x'{\text{d}}x{\text{d}}y} } $

(29)

式中： ${W}_{{\rm{q}}}$ 和 $ L $ 分别表示孔径光束在最后光学面上沿子午 $ x $ 和弧矢 $ y $ 方向的投影长度，由于光束在子午和弧矢平面内的焦点位置不一致，圆形的孔径光阑传递至各光学面的投影均为椭圆形，故上式中的积分区域为椭圆域； $ {x}'$ 和 $ {y}'$ 表示光线在像面上的像差表达式，分别计算为：

$ \left\{ \begin{gathered} x' = \frac{1}{{\cos {\omega _g}}}\left( {{d_{100}}{x_g} + {d_{200}}x_g^2 + {d_{020}}y_g^2 + {d_{300}}x_g^3 + {d_{120}}{x_g}y_g^2} \right) \hfill \\ y' = {h_{010}}{y_g} + {h_{110}}{x_g}{y_g} + {h_{210}}x_g^2{y_g} + {h_{030}}y_g^3 \hfill \\ \end{gathered} \right. $

(30)

其中： $ {d}_{ij0} $ 和 $ {h}_{ij0} $ 表示像差系数，详细在文献[13]中讨论，其中已包含离焦像差，即子午和弧矢方向场曲对成像性能的贡献。

$ {Q}_{\eta \left(k\right)} $ 和 $ {Q}_{\mathrm{c}\left(k\right)} $ 分别表示第k视场角光学系统的轴向色差和垂轴色差的评价函数，即轴向色差和垂轴色差在成像面上的贡献分量，计算表达式为：

$ {Q_{\eta \left( k \right)}} = \frac{{2{x_{g\_0}}\left| {{{r'}_{{\text{F}}\left( g \right)}} - {{r'}_{{\text{C}}\left( g \right)}}} \right|}}{{{{r'}_{0\left( g \right)}}}} $

(31)

$ {Q_{{\text{C}}\left( k \right)}} = \left| {{{y'}_{{\text{F}}\left( g \right)}} - {{y'}_{{\text{C}}\left( g \right)}}} \right| $

(32)

3.2 单个非球面优化

在光学系统优化中加入非球面，就能在传统全球面优化的基础上多出一个或多个优化自由度，从而能得到更好的结果。由于非球面制造及成本原因，如何选择正确的光学面作为非球面是一个非常值得研究的问题。

为研究如何寻找适合系统优化的最适非球面，需要通过前述理论计算研究光学系统中各光学面的波像差分布以及波像差随非球面系数变化的规律，下面以一个实际光学系统为例研究变化规律。该系统的视场角为 $ 2\omega ={220}$ °， $ F=2.8 $ ，该光学系统来自于文献[21]，光路示意图如图3所示，相关系统参量见表1 Original栏。

由于在优化过程中，为了用较短的时间得到优化结果，参数优化范围并不是特别大，因此优化后的光学系统结构及光路较原系统变化较为微小，反映在图像中是无法进行区别的，并且非球面系数的变化也无法在图中显示出来，因此本文中其余优化结果的结构及光路均可参照图3。

对于一个光学系统而言，非球面应当选择对整个系统影响较大的那个面，若其作为非球面对系统的影响可以忽略不记（即MTF曲线在其作为非球面和球面时变化不大），即可认定其不适合作为非球面。据此，以光学面对非球面系数的敏感度（系统成像性能随其非球面系数变化的剧烈程度）作为选择非球面的依据，具体过程如下。

首先可由第1、2节关于波像差的计算公式得到系统中各个光学面的波像差以及系统总体的波像差分布规律，根据式（19）～（22）以及式（26）可得到每个光学面的波像差分布特点，结果如图4，图中所示为视场角 $ {\omega }_{0}={110} $ ° (最大视场物点) 时各光学面的波像差分布。

其中： $ {W}_{\mathrm{M}} $ 和 $ {W}_{\mathrm{S}} $ 分别为场曲波像差在子午和弧失方向上的分量； $ {W}_{\mathrm{C}\mathrm{L}} $ 和 $ {W}_{\mathrm{C}\mathrm{T}} $ 分别为轴向和垂轴色差波像差； $ {W}_{\mathrm{s}\mathrm{p}\mathrm{h}} $ 和 $ {W}_{\mathrm{c}\mathrm{o}\mathrm{m}\mathrm{a}} $ 分别为球差波像差和彗差波像差。

从图4中可以看到虽然系统中有几个面的场曲波像差（子午）、球差波像差和彗差波像差较为显著，但系统本身已经对其做出了矫正，最终像面上的像差还在合理范围内。根据波像差分布情况可以决定Q值权重系数，由于该光学系统球差、彗差等孔径像差较为突出，而垂轴色差及轴向色差较小，因此可令 $ {\varepsilon }_{k}=2 $ 、 $ {\eta }_{k}=0.5 $ 和 $ {\mu }_{k}=0.5 $ 。

接下来应用式（27）～（32）依次分析光学系统成像性能Q随各光学面的非球面系数 $ {a}_{2} $ 变化的规律。 $ {a}_{2} $ 的遍历范围为 $ \left(-\mathrm{3,1}\right) $ ，变化步距为0.001；取 $ {0} $ °、 $ {10} $ °、 $ {20} $ °、 $ {30} $ °、 $ {40} $ °、 $ {50} $ °、 $ {60} $ °、 $ {70} $ °、 $ {80} $ °、 $ {90} $ °、 $ {100} $ °、 $ {110} $ °这12个视场物点；Q值随非球面系数变化的规律如图5所示，在 $ {a}_{2} $ 的遍历范围内，第6、8和9面为非球面时系统的Q值变化幅度较大。达到极小值时对应的 $ {a}_{2} $ 值见表2。

图6是第9面作为非球面，非球面系数 $ {a}_{2}\left(9\right)= -1.109 $ 时对应的波像差分布。

从表2中可以清楚的看出，在将第9面作为非球面后，评价函数Q的值下降较为明显。对比图6和图4，可以验证系统总的像差是降低的。因此，通过改变光学元件的面型，找到合适的非球面系数，能够优化光路传播、调整系统波像差分布情况，使最终成像性能得到提升。在上述步骤中，利用本文提出的方法可以找到合适的非球面位置及其初始非球面系数，之后需要结合其余光学参量进行优化，确定结果是否符合期望。

原始光学系统视场角虽然为 $ {110} $ °，但其原始MTF曲线未能显示 $ {90} $ °以后的曲线，因此先进行了一轮全球面优化。所得到的MTF曲线图见图7，优化后参数见表1 Optimize Ⅰ栏。在图7中，MTF曲线图是在Zemax中应用快速傅里叶变换(FFT)方式得到的；在子午（Mer.）和弧矢（Sag.）方向空间频率分别取10 lp/mm 和30 lp/mm；经过全球面优化后各个视场下的MTF曲线均有所提高且 $ {90} $ °之后的曲线图也能显示出来，说明在全球面优化后的系统已能在 $ {110} $ °视场内有良好的性能。

我们应用以MTF作为评价函数的自适应归一化实数编码遗传算法^[22-23]优化该实例镜头。此次优化遵循先进行大范围多代数粗略搜索参数如下：种群规模设为100；交叉概率为0.8；初始变异概率为0.1；求解精度为0.001；优化800代；范围设置为初始参数的±50%；优化20次。从结果中选择优化较好的5个系统的参数值作为小范围精确搜索的初始值，此时，其余参数不变，只将优化代数缩小为200代；搜索范围为±(10%～20%)；每组优化10次；选取MTF最高的参数作为最终系统^[16]。以第9面为例进行单个非球面的优化，优化后的MTF曲线如图8所示，各光学面参数见表3 Optimized Ⅱ。表3相较于表2增加了非球面系数（Aspheric coefficient）一项，并且由于材料与折射率均未改变，因此未在表3中列出。

图8相较于图7，第9面作为非球面优化后，MTF曲线较全球面优化后提升较多，视场角在90°之前曲线平滑，下降缓慢，成像质量优良，在90°之后出现骤降，说明在大视场下成像质量不是很好，但相较于全球面优化仍然有所提升，说明本文提出的方法对于非球面的选取是有效的，对于大视场下优化结果较差的情况，可以考虑通过增加非球面个数尝试解决。

3.3 同时进行两个非球面优化

由3.2节的结论可以看出，在对光学系统进行优化时选择对非球面系数变化较为敏感的光学面作为待优化的非球面能够有效降低系统像差，并得到较好的MTF曲线。这对单个非球面是有效的，因此可以考虑将之运用于多个非球面优化系统中，以下将采用两个非球面优化该系统。

采用多个光学面作为非球面时，由于前一个光学面的非球面系数改变将会导致在它之后光路发生变化，光线到达下一个非球面时的位置与理论不符，因此，多个非球面参与优化时，最终优化结果并不是单个非球面影响的叠加，需要得到多个非球面共同作用下对系统整体的影响。图9为几组两个非球面组合时系统Q值的变化规律图。

依据图5，在本文中将待选非球面分为两类：一类是对非球面系数变化敏感度较低的光学面（下简称为非敏感面），如第11、12面；另一类为对非球面系数变化敏感度较高的光学面（下简称敏感面），如第6、8面。在对两个非球面组合时得到的Q值变化图中，发现当两个非敏感度面组合时，如图9（a），其系统Q值变化微小；两个敏感面组合时，如图9（c），其系统Q值变化较为剧烈；而敏感面与非敏感面组合时，如图9（b），系统Q值变化规律趋向于与敏感面的变化相同。

由此，在多个光学面作为非球面参与优化系统时，仍然是敏感面对改变系统像差的贡献较大，因此可以确定在用多个非球面优化系统时仍应采用对系统Q值影响较大的光学面作为待选非球面。在本文中将第6、8和9面两两组合得到三种优化组合，表4列出了三种方案的最小Q值及对应的非球面系数，图10为第8、9面为非球面时对应的波像差分布图。

对比表2和表4，相较于加入单个非球面，加入两个非球面时，系统Q值进一步降低，比单个非球面时下降了大约23%。从图10中可以看出，第8面和第9面作为非球面后其波像差下降较为明显，而图6中第9面单独作为非球面时对自身波像差影响较小，说明多个非球面组合可以有效的抑制自身波像差，优化像差分布，进一步降低系统像差，提升成像质量。从表2可以得出，在加入单个非球面时系统Q值降低大约为50%，而在此基础上再加入非球面系统Q值下降幅度变小，说明随着非球面个数的增加，其对系统优化的能力会越来越低，最终可能导致再增加非球面个数也无法提高成像质量。

之后，仍以 MTF作为评价函数的自适应归一化实数编码遗传算法优化该镜头，步骤及参数与单个非球面优化时相同。经优化得到最终结果为第8面非球面系数为 $ {a}_{2}\left(8\right)=-1.646 $ ，第9面非球面系数为 $ {a}_{2}\left(9\right)=-1.092 $ 时系统的MTF曲线最好。优化后的系统参数见表3 Optimized Ⅲ栏，优化后的MTF曲线如图11所示。

从图11中可以看出：第8和第9面作为非球面优化后，光学系统的MTF曲线较为平滑，未出现陡降，平均MTF值在0.5以上，最低处也大于0.3；与图7、8相比MTF曲线整体上移，90°视场角后的曲线也较好，下降幅度不大，说明两个面作为非球面优化更能兼顾大视场成像。因此相较于单个非球面优化，多个非球面优化结果更好。

4 结果讨论

所采用的LU的四阶波像差理论，其对于大视场小孔径光学系统波像差计算较为精确，但当孔径较大时计算精度有所下降，结果可能偏离实际。目前六阶波像差理论尚未完善，无法用于本文中的理论计算，可以期待六阶理论完善时，结合本文提出的方法对大孔径光学系统进行更精确的计算研究。对于第6、8和9面三个光学面组合优化，本文尚未进行探讨是否会有更好的结果。对于在光学系统中加入非球面进行优化虽然能有效改善系统成像质量，但随着非球面的增多，再增加非球面个数最终可能会使成像质量的提升微乎其微甚至下降，因此优化时选择非球面的个数也是值得研究的。

5 结　论

本文基于LU的平面对称光学系统的像差理论，发展了一种有效确定非球面位置的方法，从而可以为超大视场光学系统的设计优化提供指导方法。该方法通过寻找非球面系数对系统成像质量影响较大的光学面，采用单个对非球面系数敏感的光学面或多个敏感面之间的组合作为非球面参与光学系统的优化可以有效的降低系统像差，得到较好的MTF曲线图。

本文通过实例验证了该方法的有效性，并且在优化光学系统时通过增加非球面的个数可以有效提升优化效果，并且多个非球面组合时，其对超大视场的优化明显优于单个非球面优化的结果，通过两个非球面组合优化得到了比单个非球面更好的结果，确定了该方法可以应用于复数非球面系统的优化，进一步增加了可以参加系统优化的变量，为光学系统的优化提供了新的思路。