搜索

文章查询

x

留言板

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

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

基于矢量像差理论的离轴反射光学系统初始结构设计

操超 廖志远 白瑜 范真节 廖胜

基于矢量像差理论的离轴反射光学系统初始结构设计

操超, 廖志远, 白瑜, 范真节, 廖胜
PDF
HTML
导出引用
导出核心图
  • 传统的离轴反射光学系统初始结构设计方法是先求取轴对称反射光学系统结构, 然后通过光瞳离轴、视场离轴或二者结合的方法实现无遮拦设计. 由于同轴光学系统像差分布规律不适用于离轴光学系统, 因此离轴后的反射光学系统结构像差较大, 而且系统无遮拦设计过程复杂. 本文提出了一种基于矢量像差理论的离轴反射光学系统初始结构设计方法, 可以直接获取光瞳离轴、视场离轴或二者结合的无遮拦离轴反射光学系统初始结构. 该方法可以获得较好的离轴反射光学系统初始结构供光学设计软件进一步优化. 针对面阵探测器, 设计了一个长波红外离轴三反光学系统, 通过光瞳离轴和视场离轴实现无遮拦设计, 光学系统成像质量好, 反射镜不存在倾斜和偏心, 光学系统易于装调.
      通信作者: 廖志远, liaozhiyuan1@163.com
    • 基金项目: 国家自然科学基金(批准号: 61501429)资助的课题.
    [1]

    沈本兰, 常军, 王希, 牛亚军, 冯树龙 2014 物理学报 63 144201

    Shen B L, Chang J, Wang X, Niu Y J, Feng S L 2014 Acta Phys. Sin. 63 144201

    [2]

    Jiang Z Y, Li L, Huang Y F 2009 Chin. Phys. B 18 2774

    [3]

    Fuerschbach K, Davis G E, Thompson K P, Rolland J P 2014 Opt. Lett. 39 2896

    [4]

    Meng Q Y, Wang H Y, Wang K J, Wang Y, Ji Z H, Wang D 2016 Appl. Opt. 55 8962

    [5]

    Nie Y F, Gross H, Zhong Y, Duerr F 2017 Appl. Opt. 56 5630

    [6]

    刘军, 刘伟奇, 康玉思, 吕博, 冯睿, 柳华, 魏忠伦 2013 光学学报 33 1022002

    Liu J, Liu W Q, Kang Y S, Lü B, Feng R, Liu H, Wei Z L 2013 Acta Opt. Sin. 33 1022002

    [7]

    徐奉刚, 黄玮, 徐明飞 2016 光学学报 36 1222002

    Xu F G, Huang W, Xu M F 2016 Acta Opt. Sin. 36 1222002

    [8]

    李东熙, 卢振武, 孙强, 刘华, 张云翠 2007 物理学报 56 5766

    Li D X, Lu Z W, Sun Q, Liu H, Zhang Y C 2007 Acta Phys. Sin. 56 5766

    [9]

    夏春秋, 钟兴, 金光 2015 光学学报 35 0922002

    Xia C Q, Zhong X, Jin G 2015 Acta Opt. Sin. 35 0922002

    [10]

    Thompson K P 2005 J. Opt. Soc. Am. A 22 1389

    [11]

    孙金霞, 潘国庆, 刘英 2013 物理学报 62 094203

    Sun J X, Pan G Q, Liu Y 2013 Acta Phys. Sin. 62 094203

    [12]

    Wang Y, Zhang X, Wang L J, Wang C 2014 Chin. Phys. B 23 014202

    [13]

    Schmid T, Rolland J P, Rakich A, Thompson K P 2010 Opt. Express 18 17433

    [14]

    Zhong Y, Gross H 2017 Opt. Express 25 10016

    [15]

    Shi H D, Jiang H L, Zhang X, Wang C, Liu T 2016 Appl. Opt. 55 6782

    [16]

    庞志海, 樊学武, 任国瑞, 丁蛟腾, 徐亮, 凤良杰 2016 红外与激光工程 45 0618002

    Pang Z H, Fan X W, Ren G R, Ding J T, Xu L, Feng L J 2016 Infrared Laser Eng. 45 0618002

    [17]

    Thompson K P 1980 Ph. D. Dissertation (Tucson: University of Arizona)

    [18]

    邱立军, 付霖宇, 董琪, 顾钧元 2018 兵器装备工程学报 39 88

    Qiu L J, Fu L Y, Dong Q, Gu J Y 2018 J. Ordan. Equip. Eng. 39 88

    [19]

    Yan F, Zhang X J 2009 Opt. Express 17 16809

    [20]

    Zhou G Y, Chen Y X, Wang Z G, Song H W 1999 Appl. Opt. 38 4281

  • 图 1  光瞳离轴光学系统示意图

    Fig. 1.  Schematic ofoptical system with an offset aperture stop.

    图 2  光线追迹模型

    Fig. 2.  Rays trace model.

    图 3  设计流程

    Fig. 3.  Flow chart of design process.

    图 4  误差函数收敛曲线

    Fig. 4.  Convergence curve of the error function.

    图 5  离轴三反光学系统初始结构性能分析 (a)系统布局; (b)点列图均方根直径; (c)畸变网格

    Fig. 5.  Initial configurationperformance analysis of the designed system: (a) System layout; (b) RMS spot diameter; (c) distortion grid.

    图 6  离轴三反光学系统性能分析 (a)系统布局; (b)点列图均方根直径; (c)畸变网格; (d)调制传递函数

    Fig. 6.  Performance analysis of the designed system: (a) System layout; (b) RMS spot diameter; (c) distortion grid; (d) MTF.

    表 1  系统参数

    Table 1.  System specifications.

    ParameterSpecification
    Wavelength range/μm8 to 12
    Focal length/mm200
    F-number2.5
    Field of view4° × 6°
    Pixel size/μm30
    下载: 导出CSV

    表 2  参数范围和求解结果

    Table 2.  Parameters ranges and solution of designed system.

    ParameterRangeSolution
    d1/mmd1 = –d2126.4621
    d2/mm[–200, –50]–126.4621
    d3/mm[50, 200]144.1817
    r1/mm[–500, –50]–297.3433
    r2/mm[–500, –50]–451.5485
    r3/mm[–500, –50]–131.0513
    k1[–5, 5]–0.5789
    k2[–5, 5]4.7595
    k3[–5, 5]–0.1263
    下载: 导出CSV
  • [1]

    沈本兰, 常军, 王希, 牛亚军, 冯树龙 2014 物理学报 63 144201

    Shen B L, Chang J, Wang X, Niu Y J, Feng S L 2014 Acta Phys. Sin. 63 144201

    [2]

    Jiang Z Y, Li L, Huang Y F 2009 Chin. Phys. B 18 2774

    [3]

    Fuerschbach K, Davis G E, Thompson K P, Rolland J P 2014 Opt. Lett. 39 2896

    [4]

    Meng Q Y, Wang H Y, Wang K J, Wang Y, Ji Z H, Wang D 2016 Appl. Opt. 55 8962

    [5]

    Nie Y F, Gross H, Zhong Y, Duerr F 2017 Appl. Opt. 56 5630

    [6]

    刘军, 刘伟奇, 康玉思, 吕博, 冯睿, 柳华, 魏忠伦 2013 光学学报 33 1022002

    Liu J, Liu W Q, Kang Y S, Lü B, Feng R, Liu H, Wei Z L 2013 Acta Opt. Sin. 33 1022002

    [7]

    徐奉刚, 黄玮, 徐明飞 2016 光学学报 36 1222002

    Xu F G, Huang W, Xu M F 2016 Acta Opt. Sin. 36 1222002

    [8]

    李东熙, 卢振武, 孙强, 刘华, 张云翠 2007 物理学报 56 5766

    Li D X, Lu Z W, Sun Q, Liu H, Zhang Y C 2007 Acta Phys. Sin. 56 5766

    [9]

    夏春秋, 钟兴, 金光 2015 光学学报 35 0922002

    Xia C Q, Zhong X, Jin G 2015 Acta Opt. Sin. 35 0922002

    [10]

    Thompson K P 2005 J. Opt. Soc. Am. A 22 1389

    [11]

    孙金霞, 潘国庆, 刘英 2013 物理学报 62 094203

    Sun J X, Pan G Q, Liu Y 2013 Acta Phys. Sin. 62 094203

    [12]

    Wang Y, Zhang X, Wang L J, Wang C 2014 Chin. Phys. B 23 014202

    [13]

    Schmid T, Rolland J P, Rakich A, Thompson K P 2010 Opt. Express 18 17433

    [14]

    Zhong Y, Gross H 2017 Opt. Express 25 10016

    [15]

    Shi H D, Jiang H L, Zhang X, Wang C, Liu T 2016 Appl. Opt. 55 6782

    [16]

    庞志海, 樊学武, 任国瑞, 丁蛟腾, 徐亮, 凤良杰 2016 红外与激光工程 45 0618002

    Pang Z H, Fan X W, Ren G R, Ding J T, Xu L, Feng L J 2016 Infrared Laser Eng. 45 0618002

    [17]

    Thompson K P 1980 Ph. D. Dissertation (Tucson: University of Arizona)

    [18]

    邱立军, 付霖宇, 董琪, 顾钧元 2018 兵器装备工程学报 39 88

    Qiu L J, Fu L Y, Dong Q, Gu J Y 2018 J. Ordan. Equip. Eng. 39 88

    [19]

    Yan F, Zhang X J 2009 Opt. Express 17 16809

    [20]

    Zhou G Y, Chen Y X, Wang Z G, Song H W 1999 Appl. Opt. 38 4281

  • [1] 孙金霞, 潘国庆, 刘英. 面对称光学系统的初级波像差理论研究. 物理学报, 2013, 62(9): 094203. doi: 10.7498/aps.62.094203
    [2] 吕向博, 朱菁, 杨宝喜, 黄惠杰. 基于ybar-y图的光学结构计算方法研究. 物理学报, 2015, 64(11): 114201. doi: 10.7498/aps.64.114201
    [3] 吴逢铁, 江新光, 刘彬, 邱振兴. 轴棱锥产生无衍射光束自再现特性的几何光学分析. 物理学报, 2009, 58(5): 3125-3129. doi: 10.7498/aps.58.3125
    [4] 王 方, 朱启华, 蒋东镔, 张清泉, 邓 武, 景 峰. 多程放大系统主放大级光学优化设计. 物理学报, 2006, 55(10): 5277-5282. doi: 10.7498/aps.55.5277
    [5] 李永大, 葛振杰, 孙 强, 张云翠, 王 健, 孙金霞, 刘建卓, 董科研. 折射/衍射混合红外双焦光学系统设计. 物理学报, 2006, 55(9): 4602-4607. doi: 10.7498/aps.55.4602
    [6] 裴琳琳, 吕群波, 王建威, 刘扬阳. 编码孔径成像光谱仪光学系统设计. 物理学报, 2014, 63(21): 210702. doi: 10.7498/aps.63.210702
    [7] 冯帅, 常军, 牛亚军, 穆郁, 刘鑫. 一种非对称双面离轴非球面反射镜检测补偿变焦光路设计方法. 物理学报, 2019, 68(11): 114201. doi: 10.7498/aps.68.20182253
    [8] 沈本兰, 常军, 王希, 牛亚军, 冯树龙. 三反射主动变焦系统设计. 物理学报, 2014, 63(14): 144201. doi: 10.7498/aps.63.144201
    [9] 刘飞, 魏雅喆, 韩平丽, 刘佳维, 邵晓鹏. 基于共心球透镜的多尺度广域高分辨率计算成像系统设计. 物理学报, 2019, 68(8): 084201. doi: 10.7498/aps.68.20182229
    [10] 严雄伟, 王振国, 蒋新颖, 郑建刚, 李敏, 荆玉峰. 基于微透镜阵列匀束的激光二极管面阵抽运耦合系统分析. 物理学报, 2018, 67(18): 184201. doi: 10.7498/aps.67.20172473
    [11] 张晓晖, 张爽, 孙春生. 粗糙海面对高斯分布激光光束的反射模型推导. 物理学报, 2016, 65(14): 144204. doi: 10.7498/aps.65.144204
    [12] 王肇圻, 母国光, 于 斌, 卢振武, 孙 强. 谐衍射双波段红外超光谱探测系统研究. 物理学报, 2004, 53(3): 756-761. doi: 10.7498/aps.53.756
    [13] 任洪亮. 有限远共轭显微镜光镊设计和误差分析. 物理学报, 2013, 62(10): 100701. doi: 10.7498/aps.62.100701
    [14] 徐平, 杨伟, 张旭琳, 罗统政, 黄燕燕. 集成化导光板下表面微棱镜二维分布设计. 物理学报, 2019, 68(3): 038502. doi: 10.7498/aps.68.20181684
    [15] 张书赫, 邵梦, 张盛昭, 周金华. 傅里叶域中的光线. 物理学报, 2019, 68(21): 214202. doi: 10.7498/aps.68.20190839
    [16] 张书赫, 邵梦, 周金华. 光线庞加莱球法构建的结构光场及其传输特性研究. 物理学报, 2018, 67(22): 224204. doi: 10.7498/aps.67.20180918
    [17] 温昌礼, 季家镕, 窦文华, 冯向华, 徐蓉, 门涛, 刘长海. 制备多模聚硅氧烷光波导关键技术的改进. 物理学报, 2012, 61(9): 094202. doi: 10.7498/aps.61.094202
    [18] 刘岩, 李健军, 高冬阳, 翟文超, 胡友勃, 郭园园, 夏茂鹏, 郑小兵. I类自发参量下转换相关光子圆环的时间相关特性研究. 物理学报, 2016, 65(19): 194211. doi: 10.7498/aps.65.194211
    [19] 张旭琳, 杨伟, 罗统政, 黄燕燕, 雷蕾, 李贵君, 徐平. 集成化导光板下表面微棱镜二维分布公式探究. 物理学报, 2019, 68(21): 218501. doi: 10.7498/aps.68.20190854
    [20] 张娟, 袁广才. 微腔效应对顶发射串联蓝光有机电致发光器件性能的影响. 物理学报, 2020, (): . doi: 10.7498/aps.69.20191576
  • 引用本文:
    Citation:
计量
  • 文章访问数:  276
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-05
  • 修回日期:  2019-04-11
  • 上网日期:  2019-08-16
  • 刊出日期:  2019-07-01

基于矢量像差理论的离轴反射光学系统初始结构设计

  • 1. 中国科学院光电技术研究所, 成都 610209
  • 2. 中国科学院大学, 北京 100049
  • 通信作者: 廖志远, liaozhiyuan1@163.com
    基金项目: 国家自然科学基金(批准号: 61501429)资助的课题.

摘要: 传统的离轴反射光学系统初始结构设计方法是先求取轴对称反射光学系统结构, 然后通过光瞳离轴、视场离轴或二者结合的方法实现无遮拦设计. 由于同轴光学系统像差分布规律不适用于离轴光学系统, 因此离轴后的反射光学系统结构像差较大, 而且系统无遮拦设计过程复杂. 本文提出了一种基于矢量像差理论的离轴反射光学系统初始结构设计方法, 可以直接获取光瞳离轴、视场离轴或二者结合的无遮拦离轴反射光学系统初始结构. 该方法可以获得较好的离轴反射光学系统初始结构供光学设计软件进一步优化. 针对面阵探测器, 设计了一个长波红外离轴三反光学系统, 通过光瞳离轴和视场离轴实现无遮拦设计, 光学系统成像质量好, 反射镜不存在倾斜和偏心, 光学系统易于装调.

English Abstract

    • 反射光学系统具有无色差、重量轻、热稳定性好和大口径等优点, 而且相对于同轴反射系统, 离轴反射系统不存在中心遮拦, 能量利用率高, 因此离轴反射光学系统被广泛应用于空间光学系统中[13]. 离轴反射光学系统由于非旋转对称, 像差分布形式比较复杂, 另外还需要避免光线发生遮挡, 因此离轴反射光学系统设计是光学设计领域的难点. 通过轴对称光学系统的赛德尔像差理论求解同轴反射光学系统结构再进行离轴设计是比较常用的离轴反射光学系统设计方法[46]; 另外, 文献[7,8]提出根据等光程原理和正弦条件求取同轴反射光学系统结构, 该方法可以直接获得非球面同轴反射光学系统结构, 之后再进行离轴设计; 夏春秋等[9]提出了一种基于微分方程的离轴光学系统设计方法, 该方法也是先求取同轴反射光学系统结构再进行离轴设计. 传统的离轴反射光学系统初始结构设计主要是先求取轴对称反射光学系统结构, 然后通过光瞳离轴、视场离轴或二者结合的方法实现无遮拦设计. 由于同轴光学系统像差分布规律不适用于离轴光学系统, 因此离轴后得到的反射光学系统结构像差较大, 与最终优化结果偏差较大, 而且光学系统无遮拦设计过程复杂.

      矢量像差理论由Thompson等提出并不断完善, 可以用来描述非旋转对称光学系统的像差特性[10,11]. 近年来, 矢量像差理论得到很大发展, 被用于分析非旋转对称光学系统像差特性和面型误差对光学系统性能的影响[12,13], 另外矢量像差理论还被用于指导光学系统装调、设计和优化[14,15]等.

      本文提出了一种基于矢量像差理论的离轴反射光学系统初始结构设计方法, 可以直接获得光瞳离轴、视场离轴或二者结合的无遮拦离轴反射光学系统初始结构. 基于矢量像差理论和高斯括号法, 推导了光阑位于任意位置处的离轴反射光学系统的初级矢量波像差表达式. 针对目前主流的光学设计软件均采用阻尼最小二乘法作为优化算法, 该算法是一种局部优化算法, 容易陷于局部最优解, 因此选择一个好的初始结构是光学设计非常重要的一步. 本文通过遗传算法求取离轴反射光学系统的初始结构. 遗传算法是一种全局优化算法, 可以获得全局最优解. 针对面阵探测器, 设计了一个长波红外离轴三反光学系统, 通过光瞳离轴和视场离轴实现无遮拦设计, 系统成像质量好, 反射镜不存在倾斜和偏心, 光学系统易于装调.

    • 光瞳离轴、视场离轴或二者结合的方法是实现无遮拦设计的主要方式. 对于视场离轴的光学系统, 只是在分析时取轴对称系统的轴外视场部分, 因此视场离轴光学系统像差表现形式与轴对称光学系统一致; 对于光瞳离轴的光学系统, 由于光瞳存在偏移, 像差表现形式变得十分复杂, 光瞳离轴光学系统可以看作是轴对称系统的光瞳截取轴外一部分形成的子光学系统, 如图1所示.

      图  1  光瞳离轴光学系统示意图

      Figure 1.  Schematic ofoptical system with an offset aperture stop.

      光瞳离轴反射系统孔径及其对应轴对称系统孔径的关系如(1)式所示[16], ρ'为轴对称反射系统孔径归一化矢量, ρ为离轴反射系统孔径归一化矢量, P1为孔径偏移归一化矢量, B为离轴反射系统与轴对称反射系统的孔径缩放比.

      $\left\{{ \begin{aligned} & B = {R_2}/{R_1} \\ & {{{\rho}}'} = {{{r}}_1}/{R_1} \\ & {{\rho}} = {{{r}}_2}/{R_2} \\ &{{{P}}_1} = {{P}}/{R_1} \\ &{{{r}}_1} = {{{r}}_2} + {{P}} \\ &{{{\rho}}'} = B{{\rho}} + {{{P}}_1} \\ \end{aligned} }\right.,$

      其中R1为轴对称反射系统的孔径半径, R2为离轴反射系统的孔径半径, P为光瞳偏移矢量.

      轴对称反射光学系统的初级矢量波像差为[17]

      $\begin{split} 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) + \frac{1}{2}\sum\limits_j {{W_{222j}}} \left( {{{{H}}^2} \cdot {{{{\rho}} '}^2}} \right) \\ & + \sum\limits_j {{W_{220Mj}}} \left( {{{H}} \cdot {{H}}} \right)\left( {{{\rho}} ' \cdot {{\rho}} '} \right) \\ &+ \sum\limits_j {{W_{311j}}} \left( {{{H}} \cdot {{H}}} \right)\left( {{{H}} \cdot {{\rho}} '} \right), \end{split}$

      其中H为轴对称反射系统视场归一化矢量, Wj为不同表面的初级波像差系数.

      根据(1)和(2)式可得离轴反射光学系统的初级矢量波像差为

      $\left\{ \begin{split} &W = \sum\limits_j {{W_{040j}}} {\left( {\left( {B{{\rho}} + {{{P}}_1}} \right) \cdot \left( {B{{\rho}} + {{{P}}_1}} \right)} \right)^2} + \sum\limits_j {{W_{311j}}} \left( {{{H}} \cdot {{H}}} \right)\left( {{{H}} \cdot \left( {B{{\rho}} + {{{P}}_1}} \right)} \right)\\ &\;\;\;\;\;\;\;\; + \frac{1}{2}\sum\limits_j {{W_{222j}}} \left( {{{{H}}^2} \cdot {{\left( {B{{\rho}} + {{{P}}_1}} \right)}^2}} \right){\kern 1pt} {\kern 1pt} + \sum\limits_j {{W_{220Mj}}} \left( {{{H}} \cdot {{H}}} \right)\left( {\left( {B{{\rho}} + {{{P}}_1}} \right) \cdot \left( {B{{\rho}} + {{{P}}_1}} \right)} \right)\\ &\;\;\;\;\;\;\;\; + \sum\limits_j {{W_{131j}}} \left( {{{H}} \cdot \left( {B{{\rho}} + {{{P}}_1}} \right)} \right)\left( {\left( {B{{\rho}} + {{{P}}_1}} \right) \cdot \left( {B{{\rho}} + {{{P}}_1}} \right)} \right),\\ &{W_{220Mj}} = {W_{220j}} + \frac{1}{2}{W_{222j}}. \end{split} \right.$

      将(3)式进行展开并忽略其中的常数项, 可得离轴光学系统的初级矢量波像差为

      $\begin{aligned} &W = {B^4}\sum\limits_j {{W_{040j}}} {\left( {{{\rho}} \cdot {{\rho}} } \right)^2} + {B^3}\left( {4\sum\limits_j {{W_{040j}}} {{{P}}_1} + \sum\limits_j {{W_{131j}}} {{H}}} \right) \cdot {{\rho}} \left( { {{\rho}} \cdot {{\rho}} } \right)\\ &\;\;\;\;\;\;\; + {B^2}\left( {2\sum\limits_j {{W_{040j}}} {{{P}}_1}^2 + \sum\limits_j {{W_{131j}}} {{{P}}_1}{{H}} + \frac{1}{2}\sum\limits_j {{W_{222j}}} {{{H}}^2}} \right) \cdot {{{\rho}} ^2}\\ &\;\;\;\;\;\;\; + {B^2}\left( {4\sum\limits_j {{W_{040j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right) + 2\sum\limits_j {{W_{131j}}} \left( {{{{P}}_1} \cdot {{H}}} \right) + \sum\limits_j {{W_{220Mj}}\left( {{{H}} \cdot {{H}}} \right)} } \right)\left( {{{\rho}} \cdot {{\rho}} } \right)\\ &\;\;\;\;\;\;\; + B\left( {4\sum\limits_j {{W_{040j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right){{{P}}_1} + 2\sum\limits_j {{W_{131j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right){{H}} + \sum\limits_j {{W_{131j}}} {{{P}}_1}^2{{{H}}^*}} \right) \cdot {{\rho}} \\ &\;\;\;\;\;\;\; + B\left( {\sum\limits_j {{W_{222j}}} {{{H}}^2}{{{P}}_1}^* + 2\sum\limits_j {{W_{220Mj}}} \left( {{{H}} \cdot {{H}}} \right){{{P}}_1} + \sum\limits_j {{W_{311j}}} \left( {{{H}} \cdot {{H}}} \right){{H}}} \right) \cdot {{\rho}} , \end{aligned}$

      其中P12, H2, ρ2表示向量的乘法[17]; H*, P1*表示向量的共轭.

      初级矢量波像差系数不受光学元件倾斜和偏心的影响, 因此可以通过对轴对称光学系统主光线和边缘光线追迹获得初级矢量波像差系数, 本文通过高斯括号法推导光瞳位于任意位置处的离轴反射系统初级矢量波像差系数, 如图2所示.

      图  2  光线追迹模型

      Figure 2.  Rays trace model.

      根据高斯括号法原理, 近轴光线从第i个光学表面传播到第j个光学表面可以用四个高斯常数来描述, 高斯常数如(5)式所示[14]:

      $\left\{ \begin{split} & ^i{A_j} = \left[ {{\varPhi _i}, - {e_i},{\varPhi _{i + 1}}, \cdots,{\varPhi _{j - 1}}, - {e_{j - 1}}} \right], \\ & ^i{B_j} = \left[ { - {e_i},{\varPhi _{i + 1}}, - {e_{i + 1}}, \cdots,{\varPhi _{j - 1}}, - {e_{j - 1}}} \right], \\ & ^i{C_j} = \left[ {{\varPhi _i}, - {e_i},{\varPhi _{i + 1}}, - {e_{i + 1}}, \cdots,{\varPhi _{j - 1}}, - {e_{j - 1}},{\varPhi _j}} \right], \\ & ^i{D_j} = \left[ { - {e_i},{\varPhi _{i + 1}}, - {e_{i + 1}}, \cdots,{\varPhi _{j - 1}}, - {e_{j - 1}},{\varPhi _j}} \right], \end{split} \right.$

      其中ei = di/ni, di为不同光学表面间的距离, ni为折射率; Фi表示第i个表面的光焦度.

      近轴光线从第i个光学表面传播到第j个光学表面, 光线的入射高度hj, hi和孔径角ui, ui', uj, uj' 变化关系可以用四个高斯常数来表示, 如(6)式所示:

      $\left\{ \begin{split} & \left( {\begin{array}{*{20}{c}} {{h_j}} \\ {{n_j}{u_j}'} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{}^i{A_j}}&{{}^i{B_j}} \\ {{}^i{C_j}}&{{}^i{D_j}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{h_i}} \\ {{n_{i - 1}}{u_i}} \end{array}} \right), \\ & \left( {\begin{array}{*{20}{c}} {{h_j}} \\ {{n_{j - 1}}{u_j}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{}^i{A_j}}&{{}^i{B_j}} \\ {{}^i{C_{j - 1}}}&{{}^i{D_{j - 1}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{h_i}} \\ {{n_{i - 1}}{u_i}} \end{array}} \right), \\ & \left( {\begin{array}{*{20}{c}} {{h_j}} \\ {{n_j}{u_j}'} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{}^{i + 1}{A_j}}&{{}^i{B_j}} \\ {{}^{i + 1}{C_j}}&{{}^i{D_j}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{h_i}} \\ {{n_i}{u_i}'} \end{array}} \right), \\ & \left( {\begin{array}{*{20}{c}} {{h_j}} \\ {{n_{j - 1}}{u_j}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{}^{i + 1}{A_j}}&{{}^i{B_j}} \\ {{}^{i + 1}{C_{j - 1}}}&{{}^i{D_{j - 1}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{h_i}} \\ {{n_i}{u_i}'} \end{array}} \right). \\ \end{split} \right.$

      根据(6)式可以得到主光线和边缘光线的近轴光线追迹数据, 由初级波像差系数和几何像差系数的关系, 轴对称系统的初级波像差系数如(7)式所示:

      $ \left\{ \begin{split} & {W_{040j}} = \frac{1}{8}{S_{Ij}} = - \frac{1}{8}{A_j}^2{h_j}\left( {\frac{{{u_j}'}}{{{n_j}}} - \frac{{{u_j}}}{{{n_{j - 1}}}}} \right) + \tau, \\ & {W_{131j}} = \frac{1}{2}{S_{IIj}} = - \frac{1}{2}{A_j}\overline {{A_j}} {h_j}\left( {\frac{{{u_j}'}}{{{n_j}}} - \frac{{{u_j}}}{{{n_{j - 1}}}}} \right) + \frac{{\overline {{h_j}} }}{{{h_j}}}\tau, \\ & {W_{222j}} = \frac{1}{2}{S_{IIIj}} = - \frac{1}{2}{\overline {{A_j}} ^2}{h_j}\left( {\frac{{{u_j}'}}{{{n_j}}} - \frac{{{u_j}}}{{{n_{j - 1}}}}} \right) + {\left( {\frac{{\overline {{h_j}} }}{{{h_j}}}} \right)^2}\tau, \\ & {W_{220j}} = \frac{1}{4}{S_{IVj}} = - \frac{1}{4}{H_j}^2{c_j}\left( {\frac{1}{{{n_j}}} - \frac{1}{{{n_{j - 1}}}}} \right), \\ &{W_{311j}} = \frac{1}{2}{S_{Vj}} = - \frac{1}{2}\left[ \frac{{\overline {{A_j}} }}{{{A_j}}}{H_j}^2{c_j}\left( {\frac{1}{{{n_j}}} - \frac{1}{{{n_{j - 1}}}}} \right) \right.+ \left. \frac{{{{\overline {{A_j}} }^3}}}{{{A_j}}}{h_j}\left( {\frac{{{u_j}'}}{{{n_j}}} - \frac{{{u_j}}}{{{n_{j - 1}}}}} \right) \right] + {\left( {\frac{{\overline {{h_j}} }}{{{h_j}}}} \right)^3}\tau, \\ & \tau = {c_j}^3\left( {{n_j} - {n_{j - 1}}} \right){k_j}{h_j}^4, \end{split} \right. $

      其中: ${A_j} = \left( {u_j^\prime - {u_j}} \right)/\left( {1/{n_j} - 1/{n_{j - 1}}} \right)$, $\overline {A_ j} =$$ ( - {H_j} + \overline {{{\bar h}_j}}{A_j} )/{h_j} $, Hj为光学系统拉赫不变量, kj为光学表面二次常数, cj为光学表面曲率半径.

      根据(4)式和(6), (7)式可以获得光瞳位于任意位置处的离轴光学系统初级矢量波像差. 初级矢量波像差系数如(8)式所示, 其中球差系数C1和场曲系数C4为标量, 彗差系数C2、像散系数C3和畸变系数C5为矢量.

      $\left\{ \begin{aligned} &{C_1} = {B^4}\sum\limits_j {{W_{040j}},} \\ &{{{C}}_2} = {B^3}\left( {4\sum\limits_j {{W_{040j}}} {{{P}}_1} + \displaystyle\sum\limits_j {{W_{131j}}} {{H}}} \right),\\ &{{{C}}_3} = {B^2}\left( {2\sum\limits_j {{W_{040j}}} {{{P}}_1}^2 + \displaystyle\sum\limits_j {{W_{131j}}} {{{P}}_1}{{H}} + \frac{1}{2}\sum\limits_j {{W_{222j}}} {{{H}}^2}} \right),\\ &{C_4} = {B^2}\left( {4\displaystyle\sum\limits_j {{W_{040j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right) + 2\sum\limits_j {{W_{131j}}} \left( {{{{P}}_1} \cdot {{H}}} \right) + \sum\limits_j {{W_{220Mj}}\left( {{{H}} \cdot {{H}}} \right)} } \right),\\ &{{{C}}_5} = B\left( {4\displaystyle\sum\limits_j {{W_{040j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right){{{P}}_1} + 2\sum\limits_j {{W_{131j}}} \left( {{{{P}}_1} \cdot {{{P}}_1}} \right){{H}} + \sum\limits_j {{W_{131j}}} {{{P}}_1}^2{{{H}}^*}} \right.\\ &\;\;\;\;\;\;\;\; \sum\limits_j {{W_{222j}}} {{{H}}^2}{{{P}}_1}^* + 2\displaystyle\sum\limits_j {{W_{220Mj}}} \left( {{{H}} \cdot {{H}}} \right){{{P}}_1} + \left. {\sum\limits_j {{W_{311j}}} \left( {{{H}} \cdot {{H}}} \right){{H}}} \right). \end{aligned} \right.$

    • 建立了评价离轴反射系统性能的误差函数, 误差函数由中心视场的初级像差系数和焦距约束条件组成, 也可以根据系统的需求添加其他约束条件, 误差函数如(9)式所示:

      $\begin{split} G =\;& {\omega _1}{\left\| {{C_1}} \right\|_1} \!+\! {\omega _2}{\left\| {{{{C}}_2}} \right\|_1} \!+ \!{\omega _3}{\left\| {{{{C}}_3}} \right\|_1} \!+\! {\omega _4}{\left\| {{C_4}} \right\|_1} \\ &+ {\omega _5}{\left\| {{{{C}}_5}} \right\|_1} + {\omega _6}{\left\| {{f'} - {h_1}/{u_k}'} \right\|_1}, \end{split}$

      其中ωi(i =1, 2,···, 6)为不同项的权重, f’为光学系统焦距, ${\left\| {} \right\|_1}$表示1-范数, uk' 表示边缘光线在第k个表面的出射角度.

      根据(4)式和(6)—(9)式, 误差函数可以表示为光瞳偏移矢量、中心视场偏移矢量和光学系统结构参数的函数, 如(10)式所示:

      $\begin{split} G &= G\left( {B,{{{P}}_1},{{H}},{r_j},{d_j},{k_j}} \right) \\ & = {\omega _1}{\left\| {{C_1}} \right\|_1} + {\omega _2}{\left\| {{{{C}}_2}} \right\|_1} + {\omega _3}{\left\| {{{{C}}_3}} \right\|_1} \\ &\quad + {\omega _4}{\left\| {{C_4}} \right\|_1} + {\omega _5}{\left\| {{{{C}}_5}} \right\|_1}\\ &\quad + {\omega _6}{\left\| {{f'} - {h_1}/{u_k}'} \right\|_1}. \end{split} $

      误差函数G是光学系统性能的综合反映, 其值越小表示光学系统的初始结构越好. 由于光学设计软件采用的是一种局部优化算法, 容易陷入局部最优解, 因此好的初始结构可以获得更优的设计结果. 本文通过遗传算法[18]求解误差函数G的全局最优解, 遗传算法是一种高度并行、随机和自适应的全局优化算法, 已在光学系统设计中得到广泛应用[19,20]. 遗传算法优化误差函数的设计流程如图3所示: 首先进行编码和随机产生初始种群; 然后计算误差函数值, 函数值可以用来评价个体的优劣, 个体的函数值越小被选择的概率越大; 再通过选择、交叉和变异等操作产生新一代种群; 重复以上操作, 直到满足迭代终止条件, 解码并输出优化结果. 通过遗传算法可以直接获得无遮拦的离轴反射光学系统初始结构, 之后通过光学设计软件对初始结构进行进一步优化.

      图  3  设计流程

      Figure 3.  Flow chart of design process.

    • 根据第二部分提出的离轴反射光学系统初始结构设计方法, 针对长波红外面阵探测器, 设计了一个无遮拦离轴三反光学系统, 光学系统在x方向的视场角为6°, 在y方向的视场角为4°, 光学系统参数如表1所列.

      ParameterSpecification
      Wavelength range/μm8 to 12
      Focal length/mm200
      F-number2.5
      Field of view4° × 6°
      Pixel size/μm30

      表 1  系统参数

      Table 1.  System specifications.

    • 首先建立了光学系统的误差函数G, 通过光瞳离轴和视场离轴(光瞳离轴或视场离轴可视为该方式的特例)实现无遮拦设计, 光瞳偏移量沿 + Y方向80 mm, 孔径缩小比B为1/3, 视场偏移量为2°, 因此光学系统的误差函数G如(11)式所示:

      $\begin{split} G &= G\left( {B,{{{P}}_1},{{H}},{r_j},{d_j},{k_j}} \right) \\ & = {\left\| {{C_1}} \right\|_1} + {\left\| {{{{C}}_2}} \right\|_1} + {\left\| {{{{C}}_3}} \right\|_1} + {\left\| {{C_4}} \right\|_1}\\ &\quad + {\left\| {{{{C}}_5}} \right\|_1} + 5{\left\| {{f'} - {h_1}/{u_4}'} \right\|_1}, \end{split} $

      其中: P1=(0, 2/3), H=(0, –0.5).

      通过遗传算法在结构参数边界范围内求取误差函数的最优解, 设置合适的结构参数范围避免光线发生遮拦, 参数范围和求解结果如表2所列, 误差函数的收敛曲线如图4所示.

      ParameterRangeSolution
      d1/mmd1 = –d2126.4621
      d2/mm[–200, –50]–126.4621
      d3/mm[50, 200]144.1817
      r1/mm[–500, –50]–297.3433
      r2/mm[–500, –50]–451.5485
      r3/mm[–500, –50]–131.0513
      k1[–5, 5]–0.5789
      k2[–5, 5]4.7595
      k3[–5, 5]–0.1263

      表 2  参数范围和求解结果

      Table 2.  Parameters ranges and solution of designed system.

      图  4  误差函数收敛曲线

      Figure 4.  Convergence curve of the error function.

      离轴三反光学系统的初始结构如图5(a)所示, 光学系统光线无遮拦, 系统是一个二次成像光学系统, 在中间像面处可以放置视场光阑, 有效地抑制杂散光; 初始结构点列图如图5(b)所示, 不同视场点列图均方根(root mean square, RMS)直径关于中心视场对称分布, 表明初始结构具有较好的像差分布特性; 初始结构畸变网格如图5(c)所示, 初始结构畸变较小, 该方法可以直接获得成像质量较好的无遮拦离轴反射光学系统初始结构.

      图  5  离轴三反光学系统初始结构性能分析 (a)系统布局; (b)点列图均方根直径; (c)畸变网格

      Figure 5.  Initial configurationperformance analysis of the designed system: (a) System layout; (b) RMS spot diameter; (c) distortion grid.

      为了提高初始结构的成像质量, 通过光学设计软件对初始结构进行进一步优化, 光学系统最终结构如图6(a)所示, 在优化过程中反射镜的倾斜和偏心不作为优化变量, 通过使用非球面提高成像质量, 因此反射镜不存在倾斜和偏心, 离轴反射光学系统易于装调; 光学系统点列图如图6(b)所示, 点列图RMS直径小于探测器像元尺寸30 µm, 满足使用要求; 光学系统畸变网格如图6(c)所示, 光学系统最大畸变小于0.6%; 光学系统调制传递函数(modulation transfer function, MTF)如图6(d)所示, MTF接近衍射极限, 表明光学系统具有好的成像质量.

      图  6  离轴三反光学系统性能分析 (a)系统布局; (b)点列图均方根直径; (c)畸变网格; (d)调制传递函数

      Figure 6.  Performance analysis of the designed system: (a) System layout; (b) RMS spot diameter; (c) distortion grid; (d) MTF.

    • 本文提出了一种基于矢量像差理论的离轴反射光学系统初始结构设计方法, 通过设置合适的参数范围可以直接获取光瞳离轴、视场离轴或二者结合的离轴反射光学系统初始结构. 推导了光瞳位于任意位置处的离轴反射光学系统初级矢量波像差表达式, 建立了评价光学系统综合性能的误差函数, 针对光学设计软件通常采用的是一种局部优化算法, 初始结构性能会影响光学设计优化结果, 通过遗传算法对误差函数进行求解, 可以获得误差函数的全局最优解, 该方法可以获得较好的离轴反射光学系统初始结构. 针对长波红外面阵探测器, 设计了一个离轴三反光学系统, 视场角为4° × 6°, 光学系统成像质量好, 反射镜不存在倾斜和偏心, 光学系统易于装调.

参考文献 (20)

目录

    /

    返回文章
    返回