留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

大视场大相对孔径自由曲面成像系统设计

陈炳旭 廖志远 操超 白瑜 牟达

陈炳旭, 廖志远, 操超, 白瑜, 牟达. 大视场大相对孔径自由曲面成像系统设计[J]. 红外与激光工程. doi: 10.3788/IRLA202049.20200005
引用本文: 陈炳旭, 廖志远, 操超, 白瑜, 牟达. 大视场大相对孔径自由曲面成像系统设计[J]. 红外与激光工程. doi: 10.3788/IRLA202049.20200005
Chen Bingxu, Liao Zhiyuan, Cao Chao, Bai Yu, Mu Da. Design of the freeform imaging system with large field of view and large relative aperture[J]. Infrared and Laser Engineering. doi: 10.3788/IRLA202049.20200005
Citation: Chen Bingxu, Liao Zhiyuan, Cao Chao, Bai Yu, Mu Da. Design of the freeform imaging system with large field of view and large relative aperture[J]. Infrared and Laser Engineering. doi: 10.3788/IRLA202049.20200005

大视场大相对孔径自由曲面成像系统设计

doi: 10.3788/IRLA202049.20200005
基金项目: 国家自然科学基金(61501429)
详细信息
    作者简介:

    陈炳旭(1994-),女,吉林长春人,硕士研究生,主要从事光学设计方面的研究工作。Email:18143096012@163.com

    牟达(1979-),女,吉林长春人,副教授,硕士生导师,博士,主要从事现代光学设计及测量方面的研究。Email:md-79@126.com

  • 中图分类号: 0436

Design of the freeform imaging system with large field of view and large relative aperture

  • 摘要: 近年来,自由曲面在高性能离轴反射成像系统中应用越来越广泛。文中基于矢量像差理论和遗传算法求解初始结构,采用视场离轴和孔径离轴相结合的方式,利用XY多项式描述的自由曲面增大自由度校正离轴像差,设计了一个大视场大相对孔径的制冷型离轴反射光学系统。系统工作波段为长波红外8~12 μm,焦距400 mm,F数2,视场角$ 8^{\circ} \times 5^{\circ} $,平均RMS波像差为0.037054λ(λ=9 μm)。探测器冷光阑和光学系统出瞳相匹配,实现100%冷光阑效率,系统能量集中度较高,成像质量良好。
  • 图  1  设计过程

    Figure  1.  Flow chart of design process.

    图  2  点列图均方根直径变化趋势

    Figure  2.  The change trend of RMS diameter

    图  3  间隔约束范围

    Figure  3.  Interval constraint range

    图  4  光学系统初始结构的设计过程和最终结果。(a)初始系统。(b) φ(x1,y1)结果。(c)最终结果。

    Figure  4.  In-progress and final results for design of the initial configuration. (a) Initial system. (b) Results for system construction using Φ(X1, Y1)(c) Final resul

    图  6  最终系统的均方根波前误差

    Figure  6.  Field map of RMS wavefront error for final system

    图  7  最终系统的调制传递函数

    Figure  7.  MTF for the final system

    图  8  全视场下能量包围分布图

    Figure  8.  Encircled energy distribution of all field.

    图  5  (a)PM三维自由曲面面型 (b)SM三维自由曲面面型 (c)TM三维自由曲面面型

    Figure  5.  (a) PM 3D Free-form surface (b) SM 3D Free-form surface (c) TM 3D Free-form surface

    图  9  公差概率曲线

    Figure  9.  Tolerance probability curve

    表  1  Zernike多项式与XY多项式对应像差关系

    Table  1.   Zernike polynomials and XY polynomials correspond to aberration relations

    TermZernike
    Polynomial
    XY PolynomialCorresponding Aberration type
    1$1$$1$Piston(constant)
    ${Z_1}$${T_1}$0th order in wavefront
    2$\rho \cos \theta $$X$Distortion-Tilt(x-axis)
    ${Z_2}$${T_2}$2nd order in wavefront
    3$\rho \sin \theta $$Y$Distortion-Tilt(y-axis)
    ${Z_3}$${T_3}$2nd order in wavefront
    4$2{\rho ^2} - 1$$2\left( {{x^2} + {y^2}} \right) - 1$Defocus
    ${Z_4}$${T_1} + {T_4} + {T_6}$2nd order in wavefront
    5${\rho ^2}\cos \left( {2\theta } \right)$${x^2} - {y^2}$Astigmatism,Primary(axis at 0°or 90°)
    ${Z_5}$${T_4} + {T_6}$4th order in wavefront
    6${\rho ^2}\sin \left( {2\theta } \right)$$2xy$Astigmatism,Primary(axis at ±45°)
    ${Z_6}$${T_5}$4th order in wavefront
    7$\left( {3{\rho ^2} - 2\rho } \right)\cos \theta $$3{x^3} + 3x{y^3} - 2x$Coma,Primary(x-axis)
    ${Z_7}$${T_2} + {T_7} + {T_9}$4th order in wavefront
    8$\left( {3{\rho ^2} - 2\rho } \right)\sin \theta $$3{x^2}y + 2{y^3} - 2y$Coma,Primary(y-axis)
    ${Z_8}$${T_3} + {T_8} + {T_{10}}$4th order in wavefront
    9$6{\rho ^4} - 6{\rho ^2} + 1$$6{x^4} + 12{x^2}{y^2} - 6{x^2} + 6{y^4} - 6{y^2} + 1$Spherical aberration,Primary
    ${Z_9}$${T_1} + {T_4} + {T_6} + {T_{11}} + {T_{13}} + {T_{15}}$4th order in wavefront
    下载: 导出CSV

    表  2  光学系统参数

    Table  2.   Optical system specification

    ParameterSpecification
    Wavelength range/μm8-12
    Focal length/mm400
    F-number2
    Field of view$8^{\circ }\times 5^{\circ} $
    Pixel size/μm30
    下载: 导出CSV

    表  3  解空间结构参数范围和求解结果

    Table  3.   Ranges of configuration parameters and optimization result

    ParameterRangeResult
    d1/mm[−250,−100]−201.4271
    d2/mm[100,250]210.6225
    r1/mm[−800,−50]−780.0523
    r2/mm[−800,−50]−115.2459
    r3/mm[−800,−50]−279.4206
    k1[−10,10]9.0266
    k2[−10,10]2.5360
    k3[−10,10]0.1883
    α1[10,30]15°
    α2[10,30]25.0181°
    α3[10,30]10.0006°
    下载: 导出CSV

    表  4  系统公差分配表

    Table  4.   Tolerances distribution of the system

    Surface△R(mm)△KRMS surface error(λ=9 μm)DisplacementTilt
    (X/μm)(Y/μm)(Z/μm) (X/’’)(Y/’’)(Z/’’)
    PM0.020.10%1/50λ202025303020
    SM0.050.05%1/40λ404020202020
    TM0.050.10%1/50λ202020202030
    下载: 导出CSV

    表  5  全视场公差变动表现

    Table  5.   Tolerances performance of all field

    Cumulative probabilityChange in MTF
    F1F2F3F4F5F6F7F8F9
    50.0%0.002 40.004 30.004 60.010 50.006 40.012 80.015 60.007 60.005 4
    84.1%0.012 70.012 50.017 40.036 80.017 70.024 50.038 90.026 30.024 7
    97.7%0.018 50.025 60.029 80.073 70.028 80.039 70.062 20.045 00.031 6
    99.9%0.022 40.037 80.052 50.103 00.048 40.068 50.098 40.063 70.046 7
    下载: 导出CSV
  • [1] 赵宇宸, 何欣, 张凯, 等. 轻小型大视场自由曲面离轴光学系统设计[J]. 红外与激光工程, 2018, 47(12): 225−231.

    Zhao Y C, He X, Zhang K, et al. Optical design of mini-aturized and large field of view off-axis optical system based on freeform surface [J]. Infrared and Laser Engineering, 2018, 47(12): 225−231. (in Chinese)
    [2] 王伟. 自由曲面离轴反射式光学系统设计[D]. 南京理工大学, 2016.

    Wang W. The optical design of off-axis reflective system with freeform surface[D]. Nanjing University of Science and Technology, 2016.
    [3] Fuerschbach K, Davis G E, Thompson K P, et al. Asse-mbly of a freeform off-axis optical system employing three φ-polynomial Zernike mirrors [J]. Optics Letters, 2014, 39(10): 2896−2899. doi:  10.1364/OL.39.002896
    [4] Zhu J, Hou W, et al. Design of a low F-number freeform off-axis three mirror system with rectangular field-of-view [J]. Journal of Optics, 2015, 17(1): 015605. doi:  10.1088/2040-8978/17/1/015605
    [5] 孟庆宇, 汪洪源, 王严, 等. 大线视场自由曲面离轴三反光学系统设计[J]. 红外与激光工程, 2016, 45(10): 1018002. doi:  10.3788/IRLA201645.1018002

    Meng Q Y, Wang H Y, Wang Y, et al. Off-axis threem-irror freeform optical system with large linear field of view [J]. Infrared and Laser Engineering, 2016, 45(10): 1018002. (in Chinese) doi:  10.3788/IRLA201645.1018002
    [6] 姜晰文. 致冷型离轴反射光学系统研究[D]. 长春理工大学, 2016.

    Jiang X W. Studies on cooled off-axis reflective optical systems[D]. Changchun University of Science and Tech-nology, 2016.
    [7] 张庭成, 廖志波. 6倍制冷型中红外连续变焦光学系统设计[J]. 光学学报, 2012, 32(11): 1122004. doi:  10.3788/AOS201232.1122004

    Zhang T C, Liao Z B. Design of 6× cooled thermal im-aging middle infrared continuous zoom optical system [J]. Acta Optica Sinica, 2012, 32(11): 1122004. (in Chinese) doi:  10.3788/AOS201232.1122004
    [8] Thompson K. Description of the third-order optical ab-errations of near-circular pupil optical systems without symmetry [J]. Journal of the Optical Society of America A, 2005, 22(7): 1389−1401. doi:  10.1364/JOSAA.22.001389
    [9] 庞志海, 樊学武, 任国瑞, 等. 离轴反射光学系统像差特性研究[J]. 红外与激光工程, 2016, 45(9): 0618002.

    Pang Z H, Fan X W, Ren G R, et al. Study of aberration characteristic of off-axis reflective system [J]. Infrared and Laser Engineering, 2016, 45(9): 0618002. (in Chinese)
    [10] 操超, 廖胜, 廖志远, 等. 制冷型自由曲面离轴反射光学系统设计[J]. 光学学报, 2019, 39(11): 1122001. doi:  10.3788/AOS201939.1122001

    Cao C, Liao S, Liao Z Y, et al. Design of Cooled Free-form Off-axis Reflective Optical System [J]. Acta Optica Sinica, 2019, 39(11): 1122001. (in Chinese) doi:  10.3788/AOS201939.1122001
    [11] Thompson K P, Schmid T, Cakmakci O, et al. Real-ray-based method for locating individual surface aberration field centers in imaging optical systems without rotational sy-mmetry [J]. Journal of the Optical Society of America A, 2009, 26(6): 1503−1517. doi:  10.1364/JOSAA.26.001503
    [12] Cao C, Liao S, Liao Z, et al. Initial configuration design method for off-axis reflective optical systems using nodal aberration theory and genetic algorithm [J]. Optical En-gineering, 2019, 58(10): 105101.
    [13] 徐奉刚, 黄玮. 遗传算法在离轴四反光学系统设计中的应用[J]. 光学精密工程, 2017, 25(8): 2076−2082. doi:  10.3788/OPE.20172508.2076

    Xu F G, Huang W. Application of genetic algorithm in the design of off-axis four-mirror optical system [J]. Optics and Precision Engineering, 2017, 25(8): 2076−2082. (in Chinese) doi:  10.3788/OPE.20172508.2076
    [14] 程颖. 光学自由曲面设计方法及应用研究[D]. 天津大学, 2013.

    Cheng Y. Study on Design and Application of Freeform Optics[D]. Tianjin University, 2013.
    [15] 李旭阳, 倪栋伟, 杨明洋, 等. 基于自由曲面的大视场空间相机光学系统设计[J]. 光子学报, 2018, 47(9): 922003. doi:  10.3788/gzxb20184709.0922003

    Li X Y, Ni D W, Yang M Y, et al. Design of Large Field of View Space Camera Optical System Based on Freeform Surface [J]. Acta Photonica Sinica, 2018, 47(9): 922003. (in Chinese) doi:  10.3788/gzxb20184709.0922003
    [16] 张磊, 刘东, 师途, 等. 光学自由曲面面形检测技术[J]. 中国光学, 2017, 10(3): 283−199. doi:  10.3788/co.20171003.0283

    Zhang L, Liu D, Shi T, et al. Optical freeform surfaces testing technologies [J]. Chinese Optics, 2017, 10(3): 283−199. (in Chinese) doi:  10.3788/co.20171003.0283
  • [1] 尚超, 陈宝国.  弹载大视场红外线扫图像几何畸变仿真研究 . 红外与激光工程, doi: 10.3788/IRLA202049.0204002
    [2] 胡雪蕾, 高明, 陈阳.  大视场曲面仿生复眼光学系统设计 . 红外与激光工程, doi: 10.3788/IRLA202049.0114002
    [3] 刘智颖, 高柳絮, 黄蕴涵.  offner型连续变焦中波红外光谱成像系统设计 . 红外与激光工程, doi: 10.3788/IRLA201948.0718003
    [4] 孟祥月, 王洋, 张磊, 付跃刚, 顾志远.  大相对孔径宽光谱星敏感器光学镜头设计 . 红外与激光工程, doi: 10.3788/IRLA201948.0718005
    [5] 王鑫, 王向军.  大视场双目视觉定位系统中多目标稀疏匹配 . 红外与激光工程, doi: 10.3788/IRLA201847.0726001
    [6] 赵宇宸, 何欣, 张凯, 刘强, 崔永鹏, 孟庆宇.  轻小型大视场自由曲面离轴光学系统设计 . 红外与激光工程, doi: 10.3788/IRLA201847.1218001
    [7] 王向军, 于潼, 张佳丽, 刘峰, 王越.  DEM和单相机的弹落点坐标测量方法 . 红外与激光工程, doi: 10.3788/IRLA201847.0917005
    [8] 李重阳, 董欣, 岳丽清, 马丽娜, 张继友.  航天大视场遥感相机畸变测试方法 . 红外与激光工程, doi: 10.3788/IRLA201847.1117003
    [9] 袁影, 王晓蕊, 吴雄雄, 穆江浩, 张艳.  多孔径压缩编码超分辨率大视场成像方法 . 红外与激光工程, doi: 10.3788/IRLA201746.0824001
    [10] 张伟, 张合, 张祥金.  小型大瞬时视场光学探测系统优化设计 . 红外与激光工程, doi: 10.3788/IRLA201645.0518002
    [11] 刘巍, 李肖, 马鑫, 贾振元, 陈玲, 刘惟肖.  采用复合式靶标的近景大视场相机标定方法 . 红外与激光工程, doi: 10.3788/IRLA201645.0717005
    [12] 王臣臣, 邹刚毅, 庞志海, 李瑞昌, 樊学武.  大视场可见红外一体化光学系统设计 . 红外与激光工程, doi: 10.3788/IRLA201645.1018003
    [13] 孟祥翔, 刘伟奇, 张大亮, 姜国华, 朱秀庆, 杨建明.  双自由曲面大视场头盔显示光学系统设计 . 红外与激光工程, doi: 10.3788/IRLA201645.0418004
    [14] 程少园, 张丽, 高卫军, 王劲强.  大视场空间相机侧摆成像时几何参数分析 . 红外与激光工程,
    [15] 王虎, 罗建军.  新型紧凑型大相对孔径可见光光学系统 . 红外与激光工程,
    [16] 时光, 张海洋, 耿仲博, 杨苏辉, 赵长明.  大视场激光主动探测系统实验 . 红外与激光工程,
    [17] 霍炬, 杨宁, 杨明, 董文博.  基于小靶标拼接的大视场摄像机标定方法 . 红外与激光工程,
    [18] 庞志海, 樊学武, 邹刚毅, 赵惠.  新型大视场无遮拦三反光学系统设计 . 红外与激光工程,
    [19] 马洪涛, 张晓辉, 韩冰.  宽光谱、大视场小畸变望远系统设计 . 红外与激光工程,
    [20] 明名, 杨飞, 赵金宇, 张丽敏, 吴小霞.  折反式大口径、大视场、宽光谱光学系统 . 红外与激光工程,
  • 加载中
计量
  • 文章访问数:  913
  • HTML全文浏览量:  811
  • 被引次数: 0
出版历程
  • 网络出版日期:  2020-04-30

大视场大相对孔径自由曲面成像系统设计

doi: 10.3788/IRLA202049.20200005
    作者简介:

    陈炳旭(1994-),女,吉林长春人,硕士研究生,主要从事光学设计方面的研究工作。Email:18143096012@163.com

    牟达(1979-),女,吉林长春人,副教授,硕士生导师,博士,主要从事现代光学设计及测量方面的研究。Email:md-79@126.com

基金项目:  国家自然科学基金(61501429)
  • 中图分类号: 0436

摘要: 近年来,自由曲面在高性能离轴反射成像系统中应用越来越广泛。文中基于矢量像差理论和遗传算法求解初始结构,采用视场离轴和孔径离轴相结合的方式,利用XY多项式描述的自由曲面增大自由度校正离轴像差,设计了一个大视场大相对孔径的制冷型离轴反射光学系统。系统工作波段为长波红外8~12 μm,焦距400 mm,F数2,视场角$ 8^{\circ} \times 5^{\circ} $,平均RMS波像差为0.037054λ(λ=9 μm)。探测器冷光阑和光学系统出瞳相匹配,实现100%冷光阑效率,系统能量集中度较高,成像质量良好。

English Abstract

    • 离轴反射系统因其显著的宽光谱、高透过率、无色差的优点,被广泛应用于诸多成像领域[1]。为避免光线遮拦,离轴系统通常缺失旋转对称性,进而引入离轴像差。常规的旋转对称表面,如球面、非球面,校正这类像差的能力较弱,而光学自由曲面具有非旋转对称和多自由度性能,适用于校正此类像差。近几年来,自由曲面广泛应用于各类离轴反射成像系统[2],美国罗切斯特大学设计了一个非制冷型短焦红外的紧凑型自由曲面离轴三反光学系统[3],清华大学提出了具有小F数自由曲面离轴三反系统的设计方法[4],中国科学院长春光学精密机械与物理研究所孟庆宇采用XY多项式自由曲面设计了一款$30^{\circ }\times 1^{\circ} $的大线视场角离轴三反系统[5],长春理工大学姜晰文设计了一款$3.663^{\circ }\times 2.931^{\circ} $的制冷型离轴三反系统[6],使用了Zernike多项式自由曲面面型。与非制冷型系统相比,制冷型系统的灵敏度更高,响应速度更快,探测距离更远[7],但针对制冷型探测器,冷光阑限制了系统视场角的扩大。在制冷型离轴反射光学系统中,同时实现大视场和大相对孔径比较困难。

      光学系统设计中,直接获取一个接近目标的初始结构是很有必要的,不但可以降低对高水平设计经验的要求,而且能节约大量的软件优化设计时间。文中基于矢量像差理论和高斯矩阵光线追迹方法,利用遗传算法迭代出最小误差函数,采用光瞳离轴和视场离轴,获得接近目标的初始结构,实现了大视场大相对孔径的成像光学系统设计。最终设计的自由曲面制冷型离轴反射光学系统具有良好的成像质量,满足目前光学系统的加工、检测、装调能力,对制冷型离轴反射光学系统的设计有一定的参考价值。

    • 矢量像差理论由Shack和Thompon提出并由Thompson[8]完善,可用来描述非旋转对称系统的像差场分布规律,适用于离轴反射系统的像差关系。根据矢量像差理论,光瞳离轴、视场离轴与其对应的轴对称光学系统的关系[9]如公式(1)、(2)所示:

      $${H_{Aj}} = H - {\sigma _j}$$ (1)
      $${\rho ^{'}} = B\rho + {P_1}$$ (2)

      ${H_{Aj}}$为倾斜偏心后表面上的有效视场,$H$为轴对称光学系统的视场矢量,${\sigma _j}$为第$j$面上的偏移矢量。${\rho ^{'}}$为轴对称光学系统的孔径矢量,$\rho $为离轴光学系统的孔径矢量,${P_1}$为光瞳的偏心矢量,$B$为光瞳在离轴光学系统和同轴光学系统中的缩放比,$B = {{{R_2}} /{{R_1}}},{P_1} = {P /{{R_1}}}$${R_1}$为轴对称光学系统的光瞳半径,${R_2}$为离轴光学系统的光瞳半径。

      同轴反射光学系统的矢量波像差系数可以表示为:

      $$\left\{ \begin{aligned} & W = {\sum\limits_{{j}} {{W_{040j}}\left( {{\rho ^{'}} \cdot {\rho ^{'}}} \right)} ^2} + \sum\limits_j {{W_{131j}}\left( {H \cdot {\rho ^{'}}} \right) \times \left( {{\rho ^{'}} \cdot {\rho ^{'}}} \right)} \\ & \quad\quad + \frac{1}{2}\sum\limits_j {{W_{222j}}\left( {{H^2} \cdot {\rho ^{'}}^2} \right)} {\rm{ + }}\sum\limits_{{j}} {{W_{220Mj}}\left( {H \cdot H} \right)\left( {{\rho ^{'}} \cdot {\rho ^{'}}} \right)} \\ & \quad\quad + \sum\limits_j {{W_{131j}}\left( {H \cdot H} \right)\left( {{\rho ^{'}} \cdot {\rho ^{'}}} \right)} \\ & {W_{220Mj}} = {W_{220j}} + \frac{1}{2}{W_{222j}} \end{aligned} \right.$$ (3)

      ${W_{040j}},{W_{131j}},{W_{222j}},{W_{220j}},{W_{311j}}$依次是同轴光学系统的初级球差、慧差、像散、场曲、畸变。

      根据(1)(2)(3)式,离轴反射系统的初级波像差系数由(4)式表示,球差系数${C_1}$和场区系数${C_4}$为标量,慧差系数${C_2}$,像散系数${C_3}$和畸变系数${C_5}$为矢量[10]

      $$\begin{split} & {{C}_{1}}={{B}^{4}}\sum\limits_{j}{{{W}_{040j}}} \\ & {{C}_{2}}={{B}^{3}}\left( 4\sum\limits_{j}{{{W}_{040j}}{{P}_{1}}+\sum\limits_{j}{{{W}_{131j}}{{H}_{Aj}}}} \right) \\ & {{C}_{3}}={{B}^{2}}\left( 2\sum\limits_{j}{{{W}_{040j}}{{P}_{1}}^{2}+\sum\limits_{j}{{{W}_{131j}}{{P}_{1}}{{H}_{Aj}}}} \right.\left. +\frac{1}{2}\sum\limits_{j}{{{W}_{222j}}{{H}^{2}}_{Aj}} \right) \\ & {{C}_{4}}={{B}^{2}}\left( 4\sum\limits_{j}{{{W}_{040j}}\left( {{P}_{1}}\centerdot {{P}_{1}} \right)+2\sum\limits_{j}{{{W}_{131j}}\left( {{P}_{1}}\centerdot {{H}_{Aj}} \right)}} \right.\left. +\sum\limits_{j}{{{W}_{220Mj}}\left( {{H}_{Aj}}\centerdot {{H}_{Aj}} \right)} \right) \\ & {{C}_{5}}=B\left( 4\sum\limits_{j}{{{W}_{040j}}\left( {{P}_{1}}\centerdot {{P}_{1}} \right){{P}_{1}}+2\sum\limits_{j}{{{W}_{131j}}\left( {{P}_{1}}\centerdot {{P}_{1}} \right){{H}_{Aj}}}} \right.+\sum\limits_{j}{{{W}_{131j}}{{P}_{1}}^{2}{{H}^{*}}_{Aj}} \\ & \quad\;\;\left. +\sum\limits_{j}{{{W}_{222j}}{{H}^{2}}_{Aj}{{P}_{1}}^{*}}+2\sum\limits_{j}{{{W}_{220Mj}}\left( {{H}_{Aj}}\centerdot {{H}_{Aj}} \right){{P}_{1}}+}\sum\limits_{j}{{{W}_{311j}}\left( {{H}_{Aj}}\centerdot {{H}_{Aj}} \right){{H}_{Aj}}} \right) \end{split}$$ (4)

      离轴反射光学系统的初级波像差是同轴反射光学系统的初级波像差和位移矢量的函数,根据反射矩阵和传递矩阵获取主光线和边缘光线的近轴光线追迹数据,进而获得同轴光学系统的初级波像差系数[11-12]

      使(4)式中各像差系数的绝对值之和最小,计算出离轴光学系统的初级波像差系数,由于离轴光学系统是非旋转对称系统,单独的中心视场点的成像质量不能代表整个系统的成像性能,所以构造评价离轴反射系统性能的误差函数,由全局视场的初级像差系数组成,误差函数如式所示:

      $$ \left\{ \begin{aligned} & Q = \int\limits_{\rm{c}}^d {\int\limits_a^b {Q\left( H \right)d{H_x}{H_y}\left( {a \le {H_x} \le b,c \le {H_y} \le d} \right)} } \\ & Q\left( H \right) = {\omega _1}\left| {{C_1}} \right| + {\omega _2}{\left\| {{C_2}} \right\|_1} + {\omega _3}{\left\| {{C_3}} \right\|_1} + {\omega _4}\left| {{C_4}} \right| + {\omega _5}{\left\| {{C_5}} \right\|_1} \end{aligned} \right. $$ (5)

      其中,${\omega _i}$表示各个波像差系数的权重因子${\left\| {} \right\|_1}$为1-范数。

      公式(6)将公式(5)表达成离散点拟合的形式,更方便计算误差函数,其中N为选取的视场数量。

      $$Q{\rm{ = }}\frac{1}{N}\sum\limits_{i = 1}^N {Q\left( {{H_i}} \right)} $$ (6)

      根据近轴光线追迹数据,可以计算系统的焦距${f^{'}}$如式(7),后截距${d_k}$如式(8):

      $${f^{'}} = {{ - {h_0}} /{{u^{'}}_k}}$$ (7)
      $${d_k} = {{ - {h_k}} /{{u^{'}}_k}}$$ (8)

      光线离开光学系统的最后一个表面与系统孔径光阑的距离L为:

      $$ L = {{{ - \overline {{h_k}} } /{\overline {{u_k}} }}^{'}} $$ (9)

      为满足制冷型离轴反射系统冷光阑效率达到100%,需要使探测器冷屏和光学系统出瞳大小和位置完全相同,则可以在误差评价函数中构造以下条件:

      $$s = {d_k} - L$$ (10)

      最终构造的误差评价函数包括所有选取的视场像差系数、焦距大小和冷光阑匹配条件,如式(11):

      $$ Q{\rm{ = }}\frac{1}{N}\sum\limits_{i = 1}^N {Q\left( {{H_i}} \right)} {\rm{ + }}{\omega _6}\left| {{f^{'}} + {{{h_0}} /{{u^{'}}_k}}} \right| + {\omega _7}\left| {{d_k} - L - s} \right| $$ (11)

      误差函数Q越小,代表光学系统总体成像质量越高,文中选用遗传算法多次迭代计算得出最小值。遗传算法是一种高效率的全局搜索计算方法,通过选择、交叉、变异,逐渐搜索到接近问题最优解的区域。与传统优化方法相比,其最大的区别在于遗传算法并不是从单一解或局部区域开始迭代,而是在问题可能的解空间内进行全局搜索与计算,覆盖范围广,有利于择取最优解[13]。通过遗传算法对误差函数的优化过程如图1所示,首先根据系统各项像差系数、焦距以及冷光阑匹配条件建立误差函数,编码解空间范围并产生初始种群,计算多组误差函数数值,经过选择交叉变异循环优化,从多组解空间多次迭代计算出误差函数最小的最优解,解码即可得到光学系统的初始结构参数。将其输入光学设计软件CodeV中,做为离轴反射光学系统的初始结构的起点进行下一步优化。

      图  1  设计过程

      Figure 1.  Flow chart of design process.

    • 自由曲面的描述方法分为参数描述法和多项式描述法,参数法不能用多项式描述,如贝塞尔曲面、B样条曲面[14],因为其描述的面型没有固定的表达形式,在加工和检测过程中会有诸多问题。多项式法描述的面型比较平滑,而且精度较高,目前发展较为成熟,可以加工和检测并投入使用。在成像光学系统中,通常使用Zernike多项式和XY多项式表达自由曲面,其中,Zernike多项式在单位圆内正交,各项像差互不影响,XY多项式与Zernike多项式之间可以互相转换,即可表达各项像差,又能直接对应目前的数控车床加工形式,使用XY多项式,增大光学系统设计自由度,可以有效扩大视场角,校正离轴像差,易于加工和生产。

      XY多项式的表达形式:

      $$z = \frac{{c{r^2}}}{{1 + \sqrt {1 - \left( {1 + k} \right){c^2}{r^2}} }} + \sum\limits_{j = 2}^N {{C_j}{x^m}{y^n}} $$ (12)

      其中,

      $$j = \frac{{{{\left( {m + n} \right)}^2} + m + 3n}}{2} + 1$$ (13)

      式中$c$为圆锥曲率,$k$为圆锥常数,${C_j}{x^m}{y^n}$为的XY多项式的第$j$项,$m$$n$为非负整数,分别代表$x$$y$的阶数,在光学设计软件中,XY多项式阶数通常取十阶以内,用来表达自由曲面面型,即$m + n < 10$[15]

      在使用自由曲面优化离轴反射光学系统时,要根据光学系统所包含像差情况,优化自由曲面中能平衡像差的表达式系数。Zernike多项式可以直观描述光学系统像差,将XY多项式与Zernike多项式对应起来,可以在使用XY多项式作为优化变量时准确校正系统像差。现整理出Zernike多项式和XY多项式分别在极坐标系和笛卡尔坐标系下表现形式,并在表1中描述其与Seidel像差的对应关系。表1${Z_i}$为Zernike多项式系数,${T_i}$为XY多项式系数。

      表 1  Zernike多项式与XY多项式对应像差关系

      Table 1.  Zernike polynomials and XY polynomials correspond to aberration relations

      TermZernike
      Polynomial
      XY PolynomialCorresponding Aberration type
      1$1$$1$Piston(constant)
      ${Z_1}$${T_1}$0th order in wavefront
      2$\rho \cos \theta $$X$Distortion-Tilt(x-axis)
      ${Z_2}$${T_2}$2nd order in wavefront
      3$\rho \sin \theta $$Y$Distortion-Tilt(y-axis)
      ${Z_3}$${T_3}$2nd order in wavefront
      4$2{\rho ^2} - 1$$2\left( {{x^2} + {y^2}} \right) - 1$Defocus
      ${Z_4}$${T_1} + {T_4} + {T_6}$2nd order in wavefront
      5${\rho ^2}\cos \left( {2\theta } \right)$${x^2} - {y^2}$Astigmatism,Primary(axis at 0°or 90°)
      ${Z_5}$${T_4} + {T_6}$4th order in wavefront
      6${\rho ^2}\sin \left( {2\theta } \right)$$2xy$Astigmatism,Primary(axis at ±45°)
      ${Z_6}$${T_5}$4th order in wavefront
      7$\left( {3{\rho ^2} - 2\rho } \right)\cos \theta $$3{x^3} + 3x{y^3} - 2x$Coma,Primary(x-axis)
      ${Z_7}$${T_2} + {T_7} + {T_9}$4th order in wavefront
      8$\left( {3{\rho ^2} - 2\rho } \right)\sin \theta $$3{x^2}y + 2{y^3} - 2y$Coma,Primary(y-axis)
      ${Z_8}$${T_3} + {T_8} + {T_{10}}$4th order in wavefront
      9$6{\rho ^4} - 6{\rho ^2} + 1$$6{x^4} + 12{x^2}{y^2} - 6{x^2} + 6{y^4} - 6{y^2} + 1$Spherical aberration,Primary
      ${Z_9}$${T_1} + {T_4} + {T_6} + {T_{11}} + {T_{13}} + {T_{15}}$4th order in wavefront
    • 为实现100%冷光阑匹配,将制冷型探测器孔径光阑设置于光学系统的出瞳上,最终设计的光学系统视场角为X方向[−4°,+4°],Y方向[2.5°,7.5°]的,其具体参数如表2

      避免离轴反射系统结构发生光线遮拦,在构造误差评价函数时需合理配置视场和孔径的离轴量。求解的初始结构视场角为X方向[−2.5°,+2.5°],Y方向为[2.5°,7.5°]的二次曲面离轴反射光学系统,为了在整个视场范围内获得比较好成像质量,选取一半视场范围内的9个视场点:(0°, 2.5°), (0°, 5°), (0°, 7.5°), (1.5°, 2.5°), (1.5°, 5°), (1.5°, 7.5°), (2.5°, 2.5°), (2.5°, 5°), (2.5°, 7.5°)。利用遗传算法对构造的误差函数进行迭代计算,得到最小值。解空间的参数范围和求解结果如表3,表中d为镜片间隔,r为曲率半径,k为二次项系数,$\alpha $为倾斜角度。求解系统的点列图均方根直径变化趋势如图2所示。

      表 2  光学系统参数

      Table 2.  Optical system specification

      ParameterSpecification
      Wavelength range/μm8-12
      Focal length/mm400
      F-number2
      Field of view$8^{\circ }\times 5^{\circ} $
      Pixel size/μm30

      表 3  解空间结构参数范围和求解结果

      Table 3.  Ranges of configuration parameters and optimization result

      ParameterRangeResult
      d1/mm[−250,−100]−201.4271
      d2/mm[100,250]210.6225
      r1/mm[−800,−50]−780.0523
      r2/mm[−800,−50]−115.2459
      r3/mm[−800,−50]−279.4206
      k1[−10,10]9.0266
      k2[−10,10]2.5360
      k3[−10,10]0.1883
      α1[10,30]15°
      α2[10,30]25.0181°
      α3[10,30]10.0006°

      初始的离轴反射光学系统结构,具有比较好的成像质量,利用光学设计软件对初始光学系统进行优化,为避免光线遮拦和杂散光,须严格控制系统中可能出现的光线遮挡或产生干涉的位置距离,图3详细说明了间隔的约束范围。使用XY多项式描述的自由曲面对系统进行优化时,保证系统关于YOZ面对称,只使用含有x偶次项的多项式,先选取只含有y的参数项,再加上与x偶次幂相关的参数项作为变量,优化过程中同时扩大X方向视场角,最终实现$8^{\circ }\times 5^{\circ} $的大视场制冷型离轴反射系统,扩大视场过程中结构和像质变化情况如图4。从图4(c)中可以看出,最终系统结构紧凑,点列图中各视场弥散斑直径最大的是24.779 μm,小于一个像元尺寸。

      图  2  点列图均方根直径变化趋势

      Figure 2.  The change trend of RMS diameter

      图  3  间隔约束范围

      Figure 3.  Interval constraint range

      图  4  光学系统初始结构的设计过程和最终结果。(a)初始系统。(b) φ(x1,y1)结果。(c)最终结果。

      Figure 4.  In-progress and final results for design of the initial configuration. (a) Initial system. (b) Results for system construction using Φ(X1, Y1)(c) Final resul

      离轴反射系统的最大RMS波像差为0.056λ(λ=9 μm),平均RMS波像差为0.037λ(λ=9 μm),如图6所示。光学调制传递函数如图7所示,调制传递函数接近衍射极限,光学调制传递函数在奈奎斯特频率(16 lp/mm)大于0.33,有良好的成像质量。图8中全视场范围内80%的能量都集中在直径为0.06 mm的包围圆内,表明光学系统80%以上的能量都集中在探测器2个像元内,系统的能量集中度较高。

      图  6  最终系统的均方根波前误差

      Figure 6.  Field map of RMS wavefront error for final system

      图  7  最终系统的调制传递函数

      Figure 7.  MTF for the final system

      图  8  全视场下能量包围分布图

      Figure 8.  Encircled energy distribution of all field.

      根据以上的分析结果可以看出,自由曲面的应用有效校正了大视场离轴三反系统的像差,在$8^{\circ }\times 5^{\circ} $的大视场下,系统能量集中度高,有良好的成像质量。

      为了更直观的描述出自由曲面面型,根据光学系统最终的XY多项式面型系数,结合CodeV软件中光线追迹数据与Matlab软件连接,绘出系统中自由曲面面型,如图5。由图可见,系统中三个自由曲面面型变化比较平缓,没有局部凸起,适于目前的加工水平。对自由曲面的检测通常采用非球面检测的方法,其中计算全息检测法CGH (Computer-Generated Hologram)得到了广泛的应用[16],针对本系统中的三个自由曲面面型,适用于计算全息检测方法,加工难度和误差较小,可以达到较高的检测精度。

      图  5  (a)PM三维自由曲面面型 (b)SM三维自由曲面面型 (c)TM三维自由曲面面型

      Figure 5.  (a) PM 3D Free-form surface (b) SM 3D Free-form surface (c) TM 3D Free-form surface

    • 反射系统的公差主要包括制造公差和装调公差,制造公差包括顶点曲率半径的加工公差(△R)、二次曲面常数项公差(△K),RMS面型公差,装调公差包括在X、Y、Z三个方向的位移公差、绕X轴和Y轴的旋转公差。为保证系统在有效视场内有良好得成像质量,根据现有的加工、装调能力,用蒙特卡罗方法做1 000次公差分析概率统计,以奈奎斯特频(16 lp/mm)处的MTF作为衡量系统性能标准。公差分配结果见表4表5为各概率下每个视场的详细MTF变动表现,系统在奈奎斯特频率(16 lp/mm)处的 MTF变化曲线如图9所示,从图9公差概率曲线可以直观看出,全视场MTF有80%的概率大于0.3。在装调过程中可以通过补偿器和计算机辅助安装,自由曲面面型利用CGH方法检测,易于实现高精度自由曲面光学系统的应用。

      表 4  系统公差分配表

      Table 4.  Tolerances distribution of the system

      Surface△R(mm)△KRMS surface error(λ=9 μm)DisplacementTilt
      (X/μm)(Y/μm)(Z/μm) (X/’’)(Y/’’)(Z/’’)
      PM0.020.10%1/50λ202025303020
      SM0.050.05%1/40λ404020202020
      TM0.050.10%1/50λ202020202030

      表 5  全视场公差变动表现

      Table 5.  Tolerances performance of all field

      Cumulative probabilityChange in MTF
      F1F2F3F4F5F6F7F8F9
      50.0%0.002 40.004 30.004 60.010 50.006 40.012 80.015 60.007 60.005 4
      84.1%0.012 70.012 50.017 40.036 80.017 70.024 50.038 90.026 30.024 7
      97.7%0.018 50.025 60.029 80.073 70.028 80.039 70.062 20.045 00.031 6
      99.9%0.022 40.037 80.052 50.103 00.048 40.068 50.098 40.063 70.046 7

      图  9  公差概率曲线

      Figure 9.  Tolerance probability curve

    • 文中通过使用视场离轴和孔径离轴,约束系统焦距和冷光阑位置,结合矢量像差理论和遗传算法,采用XY多项式描述的自由曲面面型设计了一款大视场大相对孔径的制冷型离轴反射成像光学系统。迭代计算的初始结构为二次曲面面型,有较好的成像质量,采用XY多项式描述的自由曲面有针对性的校正系统像差,优化过程中同时扩大视场,最终得到应用于长波红外制冷型探测器的自由曲面离轴三反光学系统,系统成像质量良好,焦距400 mm,F数2,视场角$8^{\circ }\times 5^{\circ} $,平均RMS波像差为0.037054λ(λ=9 μm),能量集中度高,实现100%冷光阑效率,光学系统公差合理,对高分辨率、大视场的制冷型空间光学系统设计有一定的参考价值。

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回