留言板

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

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

基于量子遗传算法的自由曲面离轴反射光学系统设计

缪麟 田博宇 孙年春 张彬

缪麟, 田博宇, 孙年春, 张彬. 基于量子遗传算法的自由曲面离轴反射光学系统设计[J]. 红外与激光工程, 2021, 50(12): 20210365. doi: 10.3788/IRLA20210365
引用本文: 缪麟, 田博宇, 孙年春, 张彬. 基于量子遗传算法的自由曲面离轴反射光学系统设计[J]. 红外与激光工程, 2021, 50(12): 20210365. doi: 10.3788/IRLA20210365
Miao Lin, Tian Boyu, Sun Nianchun, Zhang Bin. Design of freeform surface off-axis reflective optical systems based on quantum genetic algorithm[J]. Infrared and Laser Engineering, 2021, 50(12): 20210365. doi: 10.3788/IRLA20210365
Citation: Miao Lin, Tian Boyu, Sun Nianchun, Zhang Bin. Design of freeform surface off-axis reflective optical systems based on quantum genetic algorithm[J]. Infrared and Laser Engineering, 2021, 50(12): 20210365. doi: 10.3788/IRLA20210365

基于量子遗传算法的自由曲面离轴反射光学系统设计

doi: 10.3788/IRLA20210365
基金项目: 四川省科技计划(2018JY0553);中国科学院自适应光学重点实验室基金(LAOF1801)
详细信息
    作者简介:

    缪麟,男,硕士生,主要从事光学设计方面的的研究工作

  • 中图分类号: O436

Design of freeform surface off-axis reflective optical systems based on quantum genetic algorithm

  • 摘要: 针对自由曲面离轴反射光学系统的初始结构较少、优化过程复杂的问题,提出了一种基于量子遗传算法的离轴反射光学系统设计方法,即采用量子遗传算法结合像差理论得到高像质的初始光学系统结构,再对其离轴化处理并引入自由曲面,进而设计得到自由曲面离轴反射光学系统。该方法可以方便、快捷地为自由曲面离轴反射光学系统提供像质良好、特定布局的初始结构。以离轴三反系统为例,采用该方法设计了一套工作在可见光波段、视场大小为10°×4°、F数为2的主三镜一体化的自由曲面离轴三反系统。仿真结果表明,该离轴三反系统视场较大、结构紧凑、成像质量良好,其最大RMS波像差为0.0485 λ (λ=0.6328 μm),全视场的调制传递函数(MTF)在50 lp/mm处大于0.9,且一体化集成结构有利于简化系统的对准及装调。
  • 图  1  同轴三反初始结构示意图

    Figure  1.  Schematic diagram of the initial structure of the on-axis three-mirror system

    图  2  设计流程图

    Figure  2.  Flowchart of the design process

    图  3  误差函数的收敛曲线

    Figure  3.  Convergence curve of error function

    图  4  初始结构的(a) 布局及(b) 调制传递函数

    Figure  4.  (a) Layout and (b) MTF of the initial structure of the on-axis three-mirror system

    图  5  (a) 间隔约束和(b) 视场拓展示意图

    Figure  5.  Schematic diagram of (a) layout constraint and (b) expansion the field of view

    图  6  光学系统:(a) 结构布局;(b) M1~M3元件集成化仿真

    Figure  6.  Optical system: (a) Structure layout;(b) Integrated simulation of M1-M3 components

    图  7  自由曲面离轴三反系统性能评价参数:(a)视场内能量集中度;(b) RMS波前误差;(c) 光学调制传递函数

    Figure  7.  Performance evaluation of off-axis three-mirror system with freeform surface. (a) Encircle energy in the field of view; (b) RMS wavefront error; (c) MTF

    表  1  光学系统参数

    Table  1.   Optical system parameters

    ParameterSpecification
    Field of view/(°) 10×4
    F-number 2
    Wavelength/nm 400-760
    Effective focal length/mm 120
    MTF/lp·mm−1 >0.9@50
    Pixel size/μm 10
    下载: 导出CSV

    表  2  参数取值范围

    Table  2.   Ranges of parameters

    Parameterα1α2$ - e_1^2$$ - e_2^2$$ - e_3^2$
    Range[0.1,0.5][1.5,2.1][−5,5][−5,5][−5,5]
    下载: 导出CSV

    表  3  量子门旋转更新策略

    Table  3.   Quantum gate rotation update strategy

    xibestif(x)>f(besti)Δθis(Ai,Bi)
    AiBi>0AiBi<0Ai=0Bi=0
    0 0 False 0 0 0 0 0
    0 0 True 0 0 0 0 0
    0 1 False 0.01π 1 −1 0 ±1
    0 1 True 0.01π −1 1 ±1 0
    1 0 False 0.005π −1 1 ±1 0
    1 1 True 0.005π 1 −1 0 ±1
    1 1 False 0 0 0 0 0
    1 1 True 0 0 0 0 0
    下载: 导出CSV

    表  4  同轴三反系统的初始结构参数

    Table  4.   Initial configuration parameters of the on-axis three-mirror system

    MirrorRadius/mmThickness/mmConic
    M1 −467.19 −141.68 0.117
    M2 −109.60 141.68 2.052
    M3 −143.19 −96.51 −0.044
    下载: 导出CSV

    表  5  三反系统自由曲面多项式系数值

    Table  5.   Polynomial coefficients of freeform surface for three-mirror system

    ItemM1M2M3
    y −5.825E-003 −1.148E-003 −5.519E-004
    x2 −1.352E-003 2.242E-004 −1.851E-004
    y2 −1.336E-003 2.196E-004 −1.895E-004
    x2y 1.374E-008 −4.828E-008 −5.616E-009
    y3 −2.060E-008 −1.371E-007 −6.372E-009
    x4 7.533E-010 1.052E-010 −3.402E-010
    x2y2 1.136E-009 −2.277E-010 −7.071E-010
    y4 6.287E-010 −1.3E-010 −3.559E-010
    x4y −2.124E-012 3.946E-014
    x2y 3 1.941E-013 6.297E-014
    y5 −7.222E-015 8.015E-014
    x6 −3.345E-015 −1.147E-015
    x4y2 8.395E-015 −4.474E-015
    x2y4 −2.407E-015 −5.253E-015
    y6 −7.531E-016 −1.681E-015
    下载: 导出CSV

    表  6  离轴三反系统的最终结构参数

    Table  6.   Final configuration parameters of the off-axis three-mirror system

    MirrorRadius/mmThickness/mmConic
    M1 −436.92 −218.86 −1.611
    M2 −243.10 218.86 3.461
    M3 −329.40 −240.47 0.548
    下载: 导出CSV
  • [1] Zhu J, Zhang B, Hou W, et al. Design of an oblique camera based on a field-dependent parameter [J]. Applied Optics, 2019, 58(21): 5650. doi:  10.1364/AO.58.005650
    [2] Seo So Hyun, Ryu Jae Myung, Choi Hojong. Focus-adjustable head mounted display with off-axis system [J]. Applied Sciences, 2020, 10(21): 7931-7931. doi:  10.3390/app10217931
    [3] Chen Bingxu, Liao Zhiyuan, Cao Chao, et al. Design of the freeform imaging system with large field of view and large relative aperture [J]. Infrared and Laser Engineering, 2020, 49(8): 20200005. (in Chinese)
    [4] Sang B. The EnMAP hyperspectral imaging spectrometer: instrument concept, calibration and technologies[C]//SPIE, 2008, 7089: 350-365.
    [5] Meng Qingyu, Zhang Wei, Long Funian. Analysis on detection ability of space-based space target visible camera [J]. Infrared and Laser Engineering, 2012, 41(8): 2079-2084. (in Chinese) doi:  10.3969/j.issn.1007-2276.2012.08.022
    [6] Mill, J D, O’Neil R R, Price S. Midcourse space experiment: introduction to the spacecraft, instruments, and scientific objectives [J]. Journal of Spacecraft and Rockets, 1994, 31(5): 900-907. doi:  10.2514/3.55673
    [7] Yang Tong, Duan Yingzhe, Cheng Dewen, et al. Freeform imaging optical system, design: Theories, development and applications [J]. Acta Optica Sinica, 2021, 41(1): 0108001. (in Chinese)
    [8] Xiao Xisheng, Yu Qinghua, Chen Guilin, et al. Locating optimal freeform surfaces for off-axis optical systems [J]. Optics Communications, 2020, 467: 125757. doi:  10.1016/j.optcom.2020.125757
    [9] Liu Jun, Wei Huang, Fan Hongjie. A novel method for finding the initial structure parameters of optical systems via a genetic algorithm [J]. Optics Communications, 2016, 361: 28-35. doi:  10.1016/j.optcom.2015.10.036
    [10] Li Na, Huang Ying. Design of large-field and low-distortion freeform space optical system with 3D construction method [J]. Acta Optica Sinica, 2016, 36(3): 0322001. (in Chinese)
    [11] 王超. 自由曲面表征函数及其应用研究[D]. 中国科学院研究生院(长春光学精密机械与物理研究所), 2014.

    Wang Chao. Research on characterization function and application of freeform[D]. Changchun: University of Chinese Academy of Sciences (Changchun Institute of Optics, Fine Mechanics and Physics), 2014. (in Chinese)
    [12] Li Hua, Fan Zhenjie, Sun Wei, et al. Analysis and test of transmit-receive isolation of off-axis three-mirror anastigmat optical antenna applied to laser communication [J]. Journal of Modern Optics, 2020, 67(13): 1217-1225. doi:  10.1080/09500340.2020.1823505
    [13] Li Junyang, Li Wenqiang, Long You, et al. Design of off-axis three-mirror optical system based on freeform surface [J]. Laser and Optoelectronics Progress, 2020, 57(9): 092202. (in Chinese)
    [14] Xu Fenggang, Huang Wei, Xu Mingfei. Design of off-axis three-mirror optical system based on Wassermann-Wolf equations [J]. Acta Optica Sinica, 2016, 36(12): 1222002. (in Chinese)
    [15] Wang Yunqi, Liu Weiqi, Zhang Daliang, et al. Design of off-axis three-mirror optical system with wide field of view based on transfer matrix [J]. Infrared and Laser Engineering, 2016, 45(4): 0418003. (in Chinese)
    [16] Jun Zhu, Wei Hou, Xiao Dongzhang, 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
    [17] 潘君骅. 光学非球面的设计、加工与检验[M]. 苏州大学: 苏州大学出版社, 2004.

    Pan Junhua. The Design, Manufacture and Test of the Aspherical Optical Surfaces[M]. Suzhou: Soochow University Press, 2001. (in Chinese)
    [18] Xu Fenggang, Huang Wei. Application of genetic algorithm in the design of off-axis four-mirror optical system [J]. Optics and Precision Enginerring, 2017, 25(8): 2076-2082. (in Chinese) doi:  10.3788/OPE.20172508.2076
    [19] Babenko B, Yang M H, Belongie S. Visual tracking with online multiple instance learning [C]//Computer Vision and Pattern Recognition, IEEE Conference on, IEEE, 2009, 33(8): 983-990.
    [20] 魏雨. 基于自由曲面的大视场推扫式光学系统设计[D]. 西安工业大学, 2020.

    Wei Yu. Design of large field of view pushover optical system based on freeform surface[D]. Xi'an: Xi'an Technological University, 2020. (in Chinese)
    [21] Han K H, Kim J H. Genetic quantum algorithm and its application to combinatorial optimization problem [J]. Congress on Evolutionary Computation, 2003, 2(2): 1354-1360.
  • [1] 廖清明, 冯泽心.  面向光束整形的自由曲面衍射光学设计方法(特邀) . 红外与激光工程, 2023, 52(7): 20230430-1-20230430-11. doi: 10.3788/IRLA20230430
    [2] 周丽军, 杨通, 程德文, 王涌天.  基于方域正交多项式自由曲面的成像系统设计方法(特邀) . 红外与激光工程, 2023, 52(7): 20230317-1-20230317-14. doi: 10.3788/IRLA20230317
    [3] 许宁晏, 高志山, 陈露, 黄静, 邹宇通, 袁群.  应用自由曲面的紧凑型长焦手机镜头设计(特邀) . 红外与激光工程, 2023, 52(7): 20230322-1-20230322-10. doi: 10.3788/IRLA20230322
    [4] 任成明, 孟庆宇, 秦子长.  大型自由曲面离轴三反光学系统降敏设计(特邀) . 红外与激光工程, 2023, 52(7): 20230287-1-20230287-11. doi: 10.3788/IRLA20230287
    [5] 王合龙, 陈建发, 黄浩阳, 崔泽曜.  基于自由曲面的离轴三反光学系统研制 . 红外与激光工程, 2023, 52(3): 20220523-1-20220523-8. doi: 10.3788/IRLA20220523
    [6] 钱壮, 莫言, 樊润东, 谈昊, 冀慧茹, 马冬林.  制冷型大面阵自由曲面离轴三反光学系统设计(特邀) . 红外与激光工程, 2023, 52(7): 20230339-1-20230339-9. doi: 10.3788/IRLA20230339
    [7] 芮丛珊, 曾春梅, 冯志强, 夏成樑, 洪洋.  离轴反射式头戴显示光学系统的自由曲面设计方法 . 红外与激光工程, 2022, 51(10): 20211119-1-20211119-11. doi: 10.3788/IRLA20211119
    [8] 许宁晏, 陈露, 黄静, 邹宇通, 袁群, 高志山.  自由曲面成像光学系统的初始结构设计方法 . 红外与激光工程, 2022, 51(2): 20210852-1-20210852-12. doi: 10.3788/IRLA20210852
    [9] 蒋婷婷, 冯华君, 李奇.  自由曲面变焦的内调焦式光学系统设计 . 红外与激光工程, 2021, 50(4): 20200290-1-20200290-8. doi: 10.3788/IRLA20200290
    [10] 毛姗姗, 李艳秋, 刘克, 刘丽辉, 郑猛, 闫旭.  高数值孔径自由曲面极紫外光刻物镜光学设计 . 红外与激光工程, 2019, 48(8): 814002-0814002(7). doi: 10.3788/IRLA201948.0814002
    [11] 姚艳霞, 袁群, 陈露, 窦沂蒙, 殷慧敏, 高志山.  结合面型和视场优化策略的自由曲面设计方法 . 红外与激光工程, 2018, 47(10): 1018001-1018001(8). doi: 10.3788/IRLA201847.1018001
    [12] 赵宇宸, 何欣, 张凯, 刘强, 崔永鹏, 孟庆宇.  轻小型大视场自由曲面离轴光学系统设计 . 红外与激光工程, 2018, 47(12): 1218001-1218001(7). doi: 10.3788/IRLA201847.1218001
    [13] 姜晰文, 贾学志, 丛杉珊.  自由曲面在制冷型离轴三反光学系统的应用 . 红外与激光工程, 2018, 47(9): 918004-0918004(7). doi: 10.3788/IRLA201847.0918004
    [14] 关姝, 王超, 佟首峰, 姜会林, 常帅, 范雪冰.  空间激光通信离轴两镜反射望远镜自由曲面光学天线设计 . 红外与激光工程, 2017, 46(12): 1222003-1222003(8). doi: 10.3788/IRLA201746.1222003
    [15] 刘军, 黄玮.  反射式自由曲面头盔显示器光学系统设计 . 红外与激光工程, 2016, 45(10): 1018001-1018001(6). doi: 10.3788/IRLA201645.1018001
    [16] 孟祥翔, 刘伟奇, 张大亮, 姜国华, 朱秀庆, 杨建明.  双自由曲面大视场头盔显示光学系统设计 . 红外与激光工程, 2016, 45(4): 418004-0418004(6). doi: 10.3788/IRLA201645.0418004
    [17] 孟庆宇, 汪洪源, 王严, 纪振华, 王栋.  大线视场自由曲面离轴三反光学系统设计 . 红外与激光工程, 2016, 45(10): 1018002-1018002(8). doi: 10.3788/IRLA201645.1018002
    [18] 王洪, 陈赞吉, 吴衡, 葛鹏.  自由曲面LED汽车前照灯光学透镜设计方法 . 红外与激光工程, 2014, 43(5): 1529-1534.
    [19] 陈颖聪, 文尚胜, 罗婉霞, 陈津桥, 谢嘉宁.  自由曲面底板的LED 光学设计 . 红外与激光工程, 2014, 43(9): 2947-2953.
    [20] 孙秀辉, 杜惊雷, 尹韶云, 邓启凌, 杜春雷.  非球面内表面调制的自由曲面配光透镜设计方法 . 红外与激光工程, 2013, 42(1): 163-166.
  • 加载中
图(7) / 表(6)
计量
  • 文章访问数:  411
  • HTML全文浏览量:  161
  • PDF下载量:  66
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-06-02
  • 修回日期:  2021-08-05
  • 网络出版日期:  2022-01-06
  • 刊出日期:  2021-12-31

基于量子遗传算法的自由曲面离轴反射光学系统设计

doi: 10.3788/IRLA20210365
    作者简介:

    缪麟,男,硕士生,主要从事光学设计方面的的研究工作

基金项目:  四川省科技计划(2018JY0553);中国科学院自适应光学重点实验室基金(LAOF1801)
  • 中图分类号: O436

摘要: 针对自由曲面离轴反射光学系统的初始结构较少、优化过程复杂的问题,提出了一种基于量子遗传算法的离轴反射光学系统设计方法,即采用量子遗传算法结合像差理论得到高像质的初始光学系统结构,再对其离轴化处理并引入自由曲面,进而设计得到自由曲面离轴反射光学系统。该方法可以方便、快捷地为自由曲面离轴反射光学系统提供像质良好、特定布局的初始结构。以离轴三反系统为例,采用该方法设计了一套工作在可见光波段、视场大小为10°×4°、F数为2的主三镜一体化的自由曲面离轴三反系统。仿真结果表明,该离轴三反系统视场较大、结构紧凑、成像质量良好,其最大RMS波像差为0.0485 λ (λ=0.6328 μm),全视场的调制传递函数(MTF)在50 lp/mm处大于0.9,且一体化集成结构有利于简化系统的对准及装调。

English Abstract

    • 离轴反射系统具有可消除中心遮拦、无色差、光路可折叠、能量利用率高等特点[1-3],而离轴三反系统除了具有反射式光学系统的特点外,还具有可实现大视场、镜面数较少等优点,被广泛应用于空间光学系统中[4-6],但也因其失去了旋转对称性而不可避免地会引入离轴像差。自由曲面具有很强的像差校正能力及较多的设计自由度,可有效应用于离轴反射系统中离轴像差的校正[7-8]。然而,自由曲面离轴反射光学系统初始结构较少,其优化过程复杂,致使其设计难度较大。同时,在离轴三反系统装调过程中,若以主镜为基准,依次调节次镜和三镜的偏心与倾斜,装调自由度共12个,由此带来的装调误差和制造成本较大,因而降低装调难度很有必要。若将主镜和三镜集成一体,会使得装调自由度减少六个,使得系统结构紧凑和易于装调。

      现有的空间光学系统设计软件大多依赖于初始结构的选择[9],直接获取一个与目标要求接近且初始像差较小的初始结构,不仅可节约后续优化的时间,而且可降低对高水平设计经验的要求。因此,选择合理的初始结构是离轴光学系统设计极为重要的一步。现有的自由曲面离轴三反系统的初始结构获取方法主要有两种:一种方法是在专利库里寻找初始结构,进而不断优化得到满足设计要求的自由曲面离轴光学系统[10-11]。然而,在大多数情况下,尤其是存在特定布局要求时,很难在专利库中找到相匹配的初始系统结构。因此,另一种方法是通过理论求解得到同轴三反系统初始结构,再经过离轴和优化得到最终的光学系统设计[12-13]。徐奉刚等[14]、王蕴琦等[15]和朱钧等[16]分别基于Wassermann-Wolf微分方程、光学传递矩阵和光焦度分配求解出同轴三反初始结构,进而设计得到离轴反射系统。然而,理论推导需要给定足够的初始条件才能求解,而若要满足特定布局结构要求,其求解过程就变得尤为复杂和耗时。因此,需要寻求快速方法以获取初始结构,通过像差理论[17]建立由同轴遮拦比αi(i=1,2)、放大率βi(i=1,2)和圆锥常数$ - e_j^2$(j=1,2,3)七个变量组成的用于评价光学系统初始结构性能的误差约束函数,进而将初始结构参数问题转化为利用优化算法寻找误差约束函数极小值问题,从而获取高性能的同轴初始结构[18]。然而,误差约束函数变量通常较多,搜索空间极大,冗余信息过多,从而造成计算耗时长。同时,遗传算法等经典优化算法在解决高维问题时容易陷于局部最优。

      文中提出了一种自由曲面离轴三反系统初始结构的快速设计方法。该方法通过减少误差约束函数中的变量数目以及引入量子遗传算法(Quantum genetic algorithm,QGA),无需给定初始条件即可快速得到初始结构。量子遗传算法是一种基于量子计算原理、结合经典遗传算法的全局优化算法[19]。相比于传统的遗传算法,量子遗传算法具有更好的种群多样性、更快的收敛速度和更强的全局搜索能力,更容易跳出局部最优解而获得全局最优解。文中首先基于像差理论得到主三镜一体化同轴三反系统的像差表达式,进而建立以同轴遮拦比和圆锥常数为变量的评价光学系统性能的误差约束函数;其次,采用量子遗传算法求解误差约束函数最小值的最优解,并将最优解作为初始结构,进而通过合理的离轴设置实现系统无遮拦设计,最终得到满足设计要求的主三镜一体化自由曲面离轴三反系统。

    • 离轴三反系统是在同轴三反系统基础上进行离轴化得到的,其常见布局如图1所示。该系统主要由主镜M1、次镜M2、三镜M3组成。每个反射镜面的圆锥系数分别为$ - e_1^2$$ - e_2^2$$ - e_3^2$M1M2的遮拦比分别为α1α2M2M3的放大率分别为β1β2

      图  1  同轴三反初始结构示意图

      Figure 1.  Schematic diagram of the initial structure of the on-axis three-mirror system

      由近轴光学理论,可根据遮拦比和放大率得到反射镜的曲率半径和相应距离的表达式分别为:

      $$\left\{ {\begin{array}{l} {{r_1} = \dfrac{2}{{{\beta _{\rm{1}}}{\beta _{\rm{2}}}}}{f^\prime }\;\;{r_2} = \dfrac{{2{\alpha _{\rm{1}}}}}{{\left( {1 + {\beta _{\rm{1}}}} \right){\beta _{\rm{2}}}}}{f^\prime }\;\;{r_3} = \dfrac{{2{\alpha _{\rm{1}}}{\alpha _2}}}{{1 + {\beta _{\rm{2}}}}}{f^\prime }} \\ {{d_1} = \dfrac{{1 - {\alpha _{\rm{1}}}}}{{{\beta _{\rm{1}}}{\beta _{\rm{2}}}}}{f^\prime }\;\;{d_2} = \dfrac{{{\alpha _{\rm{1}}}\left( {1 - {\alpha _2}} \right)}}{{{\beta _{\rm{2}}}}}{f^\prime }\;\;{d_3} = {\alpha _1}{\alpha _2}{f^\prime }} \end{array}} \right.$$ (1)

      式中:f '为三反光学系统总焦距。

      通过追踪边缘光线和主光线,结合像差理论,得到轮廓参数α1α2β1β2$ - e_1^2$$ - e_2^2$$ - e_3^2$表示的单色像差,其像差系数SSSSS是关于α1α2β1β2$ - e_1^2$$ - e_2^2$$ - e_3^2$的函数[20]

      为了降低装调难度,设计主镜和三镜一体化结构。于是,M1M3M2之间距离应相等,即d1=−d2,可得到:

      $${\beta _{\rm{1}}}{\rm{ = }}\dfrac{{1 - {\alpha _{\rm{1}}}}}{{{\alpha _{\rm{1}}}\left( {{\alpha _{\rm{2}}} - 1} \right)}}$$ (2)

      在平像场条件下,根据像差理论可得到SIV=0,即:

      $${\beta _{\rm{2}}}{\rm{ = }}\dfrac{1}{{{\alpha _2}{\beta _{\rm{1}}}\left( {1 - {\alpha _1}} \right) + {\alpha _2} - 1}}$$ (3)

      主三镜一体化结构的像差系数SSSS是关于αi(i=1,2)、$ - e_j^2$(j=1,2,3)的复杂的函数,其隐式表达为:

      $$\left\{ {\begin{array}{l} {{S_{\rm{I}}} = {S_{\rm{I}}}\left( {{\alpha _i}, - e_j^2} \right)} \\ {{S_{{\rm{II}}}} = {S_{{\rm{II}}}}\left( {{\alpha _i}, - e_j^2} \right)} \\ {{S_{{\rm{III}}}} = {S_{{\rm{III}}}}\left( {{\alpha _i}, - e_j^2} \right)} \\ {{S_{\rm{V}}} = {S_{\rm{V}}}\left( {{\alpha _i}, - e_j^2} \right)} \end{array}} \right.$$ (4)

      光学系统初始结构设计是光学设计过程中的关键。由公式(4)可知,选择合适的轮廓参数α1α2$ - e_1^2$$ - e_2^2$$ - e_3^2$,使得SSSS趋于0,即可获得小像差的系统初始结构布局。

      基于以上像差分析,文中建立了评价光学系统初始结构性能的误差约束函数Q。该误差函数由主三镜一体化结构的各个像差系数组成,也可以根据系统其他布局要求添加约束条件,其表达式为:

      $$ \begin{array}{c}Q=Q({\alpha }_{1}{\text{、}}{\alpha }_{2}{\text{、}}{e}_{1}^{2}{\text{、}}{e}_{2}^{2}{\text{、}}{e}_{3}^{2}) =\left|{S}_{\rm{I}}\right|+\left|{S}_{\rm{II}}\right|+\left|{S}_{\rm{III}}\right|+\left|{S}_{\rm{V}}\right|\end{array}$$ (5)

      误差约束函数Q的值越小,代表光学系统的成像质量越好。采用量子遗传算法对误差约束函数Q进行多次迭代求出其最小值的最优解,从而确定系统初始结构布局,再通过光学设计软件对初始结构进一步优化。量子遗传算法是一种全局优化算法,可以获得全局最优解[21],其原理是建立在量子位和量子叠加态的基础上。一个量子位可以表示处于1态(表示|1〉)、0态(表示|0〉)或者两者相叠加状态的信息。量子位可以表示为:

      $$\left| \varphi \right\rangle {\rm{ = }}A\left| 0 \right\rangle + B\left| 1 \right\rangle $$ (6)

      式中:AB分别表示量子状态被发现为0态和1态的概率幅度值。将量子比特的概率幅值表示应用于染色体的编码,使得一条染色体可以表示多个态的叠加,并利用量子旋转门向着最优个体方向进化,达到了种群更新和完成进化的目的。基于量子遗传算法搜索误差约束函数最小值的最优解的流程如图2右侧所示。

      图  2  设计流程图

      Figure 2.  Flowchart of the design process

      图2左侧为主要设计流程图,其具体设计过程:第一步,根据像差分析建立评价光学系统性能的误差约束函数Q;第二步,采用量子遗传算法求解误差约束函数最优解,获取初始结构遮拦比αi(i=1,2)和圆锥常数$ - e_j^2$(j=1,2,3);第三步,基于第二步求解得到的参数值,根据公式(1)~(3)计算结构参数r1r2r3d1d2d3;第四步,将获得的初始结构输入到光学设计软件中进一步优化,继而得到离轴系统设计结构。在第二步中,初始结构遮拦比αi(i=1,2)和圆锥常数$ - e_j^2$(j=1,2,3)的合理求解至关重要。首先,使用量子位将遮拦比αi(i=1,2)、圆锥常数$ - e_j^2$(j=1,2,3)编码成一条染色体,将量子比特概率幅值表示应用于染色体的编码,初始化种群。然后,对种群进行测量,计算种群中个体适应度,进而保留种群中最优个体和适应度。最后,判断是否满足条件,如果不满足则进入下一代迭代,用当前最优个体诱导量子旋转产生新的种群,使得种群向着最优个体进化,再次返回测量步骤,如此循环下去,直至满足条件,迭代出误差约束函数最小值的最优解,进而通过解码即可得到光学系统初始遮拦比αi(i=1,2)和圆锥常数$ - e_j^2$(j=1,2,3)。

    • 在合理限定参数范围之后,利用量子遗传算法得到误差约束函数最小值的最优解,从而得到主三镜一体化的同轴三反系统的初始结构。在此基础上,对初始系统进行离轴化处理,并引入自由曲面替代二次曲面校正离轴像差,最终实现成像质量良好、结构紧凑的自由曲面主三镜一体化离轴三反系统设计。表1给出了所设计的光学系统视场角范围为X方向[−5°,5°]、Y方向[−14°,−10°]的具体参数。

      表 1  光学系统参数

      Table 1.  Optical system parameters

      ParameterSpecification
      Field of view/(°) 10×4
      F-number 2
      Wavelength/nm 400-760
      Effective focal length/mm 120
      MTF/lp·mm−1 >0.9@50
      Pixel size/μm 10

      一个合理的同轴三反初始结构布局须满足一定的边界条件:d1<0,d2>0和d3<0。因此,在一次成像系统中满足焦距为负和光焦度分配的原则,结合公式(1)可知,参数α1α2应满足:0<α1<1,α2>1。在上述条件下,结合对系统体积要求及加工难易程度等因素的综合考虑,设置了如表2所示的参数范围。

      表 2  参数取值范围

      Table 2.  Ranges of parameters

      Parameterα1α2$ - e_1^2$$ - e_2^2$$ - e_3^2$
      Range[0.1,0.5][1.5,2.1][−5,5][−5,5][−5,5]

      在量子遗传算法中,种群更新是通过量子旋转门完成的。文中选择的量子旋转门更新策略如表3所示,其中,旋转角度θi=s(Ai,Bi)×Δθis(Ai,Bi)表示旋转方向,Δθi为旋转角的角度增量。

      表 3  量子门旋转更新策略

      Table 3.  Quantum gate rotation update strategy

      xibestif(x)>f(besti)Δθis(Ai,Bi)
      AiBi>0AiBi<0Ai=0Bi=0
      0 0 False 0 0 0 0 0
      0 0 True 0 0 0 0 0
      0 1 False 0.01π 1 −1 0 ±1
      0 1 True 0.01π −1 1 ±1 0
      1 0 False 0.005π −1 1 ±1 0
      1 1 True 0.005π 1 −1 0 ±1
      1 1 False 0 0 0 0 0
      1 1 True 0 0 0 0 0

      根据公式(5)所定义的误差约束函数和表2列出的参数范围,可通过量子遗传算法搜索误差约束函数最小值的最优解,从而获取具有合理布局的初始结构。误差约束函数的收敛曲线如图3所示,从中可得到误差约束函数Q的收敛结果以及相应的轮廓参数取值,其圆锥常数$ - e_j^2$(j=1,2,3)在表4中给出。

      图  3  误差函数的收敛曲线

      Figure 3.  Convergence curve of error function

      表 4  同轴三反系统的初始结构参数

      Table 4.  Initial configuration parameters of the on-axis three-mirror system

      MirrorRadius/mmThickness/mmConic
      M1 −467.19 −141.68 0.117
      M2 −109.60 141.68 2.052
      M3 −143.19 −96.51 −0.044

      将误差约束函数收敛结果(即轮廓参数)代入公式(1)~(3)中,计算得到如表4所示的同轴三反系统参数,并可将其作为离轴三反系统的初始结构。图4(a)(b)分别为初始结构的布局和调制传递函数(MTF)曲线。从图4中可以看出,该初始结构紧凑,且MTF接近衍射极限,成像质量良好。上述结果表明,基于上述量子遗传算法可以搜索得到满足条件的初始结构。

      图  4  初始结构的(a) 布局及(b) 调制传递函数

      Figure 4.  (a) Layout and (b) MTF of the initial structure of the on-axis three-mirror system

      以上述方法将得到的初始结构作为优化起点,为了增大光学系统视场,将光阑置于次镜,并对同轴三反初始结构采用视场离轴和孔径离轴相结合的综合离轴方式,以避免光线遮拦。同时,根据光线追迹原理,为了给反射镜留有一定的支撑结构,需要对结构布局间隔进行约束。

      图5(a)图5(b)分别为间隔约束和视场扩展示意图。其中,图5(a)给出了间隔的约束范围,在优化过程中,为保证主三镜在Z轴位置一致,需满足|d1|=|d2|。如图5(b)所示,在增大视场过程中,为了在整个视场范围内获得比较好的成像质量,选取一半视场范围内的九个视场点:(0°,−10°),(0°,−12°),(0°,−14°),(2.5°,−10°),(2.5°,−12°),(2.5°,−14°),(5°,−10°),(5°,−12°),(5°,−14°)。由于边缘视场的像差要比中心视场大,在优化设计过程中,需要根据优化的多次中间结果不断调节各个视场的优化权重,以平衡各视场在像面上的成像质量。由于自由曲面可以提供更多的设计自由度,因而引入自由曲面替代二次曲面校正离轴带来的像差,以大幅提高成像性能。采用渐进式逐步逼近最优策略,从低阶项开始逐步引入高阶项作为变量进行优化。

      图  5  (a) 间隔约束和(b) 视场拓展示意图

      Figure 5.  Schematic diagram of (a) layout constraint and (b) expansion the field of view

      在表征自由曲面时选择XY多项式,其表达式与数控加工使用的表达式相同,适合加工。此外,为了保证光学系统关于YOZ平面对称,只使用XY多项式中x的偶数项,即先对只含有y项的参数进行优化,再选择含有x项的偶次项进行优化。XY多项式的面型矢高定义为在二次基底基础上加上泰勒展开式,即:

      $${\textit{z}}\left( {x,y} \right) = \dfrac{{c\left( {{x^2} + {y^2}} \right)}}{{1 + \sqrt {1 - (k + 1){c^2}\left( {{x^2} + {y^2}} \right)} }} + \sum\limits_{i = 1}^N {{A_i}{E_i}\left( {x,y} \right)} $$ (7)
      $$i = \dfrac{{\left( {{{\left( {m + n} \right)}^2} + m + 3n} \right)}}{2} + 1$$ (8)

      式中:c为顶点处曲率;k为圆锥系数;AiXY多项式第i项的系数;Ei(x,y)为XY多项式的第i项。

      经过对初始结构优化后,最终得到X方向上视场为10°(−5°~5°)、Y方向上视场为4°(−14°~−10°)的离轴三反系统,其自由曲面的多项式系数如表5所示,系统的最终布局如图6(a)所示,其结构参数列于表6。比较表4表6可知,虽然最终结构参数相比于初始结构参数有一定变化,但仍接近于初始结构参数,从而说明了合理的初始结构设计有利于快速、有效地实现最终设计目标,从而大幅降低后续优化难度。在图6(b)中,主镜M1和三镜M3在垂直方向上近似相切,且M1~M3元件的中心部分不接收任何光线,使得在一块基板上制造两个反射镜变得可能,进而可以由一个完整的加工过程中制造得到。

      表 5  三反系统自由曲面多项式系数值

      Table 5.  Polynomial coefficients of freeform surface for three-mirror system

      ItemM1M2M3
      y −5.825E-003 −1.148E-003 −5.519E-004
      x2 −1.352E-003 2.242E-004 −1.851E-004
      y2 −1.336E-003 2.196E-004 −1.895E-004
      x2y 1.374E-008 −4.828E-008 −5.616E-009
      y3 −2.060E-008 −1.371E-007 −6.372E-009
      x4 7.533E-010 1.052E-010 −3.402E-010
      x2y2 1.136E-009 −2.277E-010 −7.071E-010
      y4 6.287E-010 −1.3E-010 −3.559E-010
      x4y −2.124E-012 3.946E-014
      x2y 3 1.941E-013 6.297E-014
      y5 −7.222E-015 8.015E-014
      x6 −3.345E-015 −1.147E-015
      x4y2 8.395E-015 −4.474E-015
      x2y4 −2.407E-015 −5.253E-015
      y6 −7.531E-016 −1.681E-015

      图  6  光学系统:(a) 结构布局;(b) M1~M3元件集成化仿真

      Figure 6.  Optical system: (a) Structure layout;(b) Integrated simulation of M1-M3 components

      表 6  离轴三反系统的最终结构参数

      Table 6.  Final configuration parameters of the off-axis three-mirror system

      MirrorRadius/mmThickness/mmConic
      M1 −436.92 −218.86 −1.611
      M2 −243.10 218.86 3.461
      M3 −329.40 −240.47 0.548

      自由曲面离轴三反系统的视场内能量集中度、RMS波前误差和调制传递函数MTF分别如图7(a)~(c)所示。从图7(a)中可以看出,在视场内80%的能量均集中在艾里斑(RAiry=1.558 μm)内,其能量集中度较高。分析图7(b)可知,系统的最大RMS波前误差为0.0485λ(参考波长λ=0.6328 μm),平均RMS波前像差为0.0377λ,可见其系统的波前误差较小。图7(c)表明,各视场的MTF均接近衍射极限,在50 lp/mm处大于0.9,成像质量良好。

      图  7  自由曲面离轴三反系统性能评价参数:(a)视场内能量集中度;(b) RMS波前误差;(c) 光学调制传递函数

      Figure 7.  Performance evaluation of off-axis three-mirror system with freeform surface. (a) Encircle energy in the field of view; (b) RMS wavefront error; (c) MTF

    • 文中提出了一种基于量子遗传算法的自由曲面离轴三反光学系统设计方法。该方法采用量子遗传算法,通过使得构造的评价光学系统性能的误差约束函数最小化,进而确定合理的初始结构参数。结果表明,基于量子遗传算法获得的初始结构可以为三反系统提供很好的设计起点,大幅降低了后续优化难度。首先,从同轴三反系统的像差理论出发,建立了评价主三镜一体结构光学性能的误差约束函数。其次,利用量子遗传算法全局搜索能力强的特点,通过迭代方式求解误差约束函数最小值的最优解,从而获得像差较小的二次曲面初始结构。在此基础上,将该初始结构作为设计起点,可降低后续优化设计的难度。最后,通过对初始结构离轴化处理,并引入自由曲面校正离轴像差,最终完成了自由曲面主三镜一体离轴三反系统设计。该系统视场较大、F数较小,且全视场的MTF在50 lp/mm处均大于0.9,接近衍射极限,具有结构紧凑、视场较大、成像质量良好等突出特点。所提出的基于量子遗传算法的自由曲面离轴三反光学系统设计方法不仅可降低自由曲面主三镜一体化离轴三反系统的设计难度,而且可推广适用于其他类型的离轴反射式系统设计。

参考文献 (21)

目录

    /

    返回文章
    返回