留言板

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

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

主被动双光复合系统平行度校正研究

刘飞 陈秀芬 李松山 聂子晨

刘飞, 陈秀芬, 李松山, 聂子晨. 主被动双光复合系统平行度校正研究[J]. 红外与激光工程, 2022, 51(7): 20210771. doi: 10.3788/IRLA20210771
引用本文: 刘飞, 陈秀芬, 李松山, 聂子晨. 主被动双光复合系统平行度校正研究[J]. 红外与激光工程, 2022, 51(7): 20210771. doi: 10.3788/IRLA20210771
Liu Fei, Chen Xiufen, Li Songshan, Nie Zichen. Research on the parallelism of the active and passive dual-band (laser/IR) composite system[J]. Infrared and Laser Engineering, 2022, 51(7): 20210771. doi: 10.3788/IRLA20210771
Citation: Liu Fei, Chen Xiufen, Li Songshan, Nie Zichen. Research on the parallelism of the active and passive dual-band (laser/IR) composite system[J]. Infrared and Laser Engineering, 2022, 51(7): 20210771. doi: 10.3788/IRLA20210771

主被动双光复合系统平行度校正研究

doi: 10.3788/IRLA20210771
详细信息
    作者简介:

    刘飞,男,工程师,硕士,主要研究方向为激光对抗系统总体设计

  • 中图分类号: TN209

Research on the parallelism of the active and passive dual-band (laser/IR) composite system

  • 摘要: 针对主被动双光复合系统中存在激光主动发射视场与红外被动成像视场平行度(一一对应位置)偏差问题,提出了主被动双光复合系统平行度校正方法,并开展了理论研究和实物仿真测试。首先,推导了激光主动发射视场与快速反射镜(Fast Steering Mirror,FSM)角度对应关系;其次,根据激光发射视场与红外探测视场提出了一种基于多项式拟合的平行度修正校正方法;最后,通过构建实际测试平台,进行平行度标定、求解、修正。试验结果表明:在±1.4°的激光指示视场和±1.4°的红外中心视场中,经过系统平行度修正后,实际平行度偏差可由0.8 mrad降低至0.1 mrad,进一步提高激光对目标的指向精度。
  • 图  1  主被动双光复合系统示意图

    Figure  1.  Schematic of the active and passive dual-band (laser/IR) composite system

    图  2  光线反射与镜面位置关系

    Figure  2.  Relationship between light reflection and mirror position

    图  3  快反镜方位摆动与激光方位、激光俯仰投影关系

    Figure  3.  Projection relationship between the azimuth swing of the FSM and the laser azimuth, laser pitch

    图  4  平行度测试原理图

    Figure  4.  Schematic diagram of the parallelism test

    图  5  平行度校正前、后给定值与测量值对比图

    Figure  5.  Comparison chart of the given values and measured values before and after parallelism correction

    表  1  快反镜方位摆动引起激光方位、俯仰变化计算结果

    Table  1.   Calculation results of the changes in laser azimuth and pitch caused by the azimuth swing of the FSM

    FSM azimuth angle
    /(°)
    Laser azimuth angleLaser pitch angle
    /μrad
    Theoretical calculation
    /(°)
    Approximate calculation
    /(°)
    Error Δ
    /μrad
    −1.00−1.414069−1.4142132.505175152.2993
    −0.75−1.060599−1.0606601.05698385.66360
    −0.50−0.707088−0.7071060.31320438.07120
    −0.25−0.353551−0.3535530.0391529.517574
    0.000000
    0.250.3535510.353553−0.0391529.517574
    0.500.7070880.707106−0.30320438.07120
    0.751.0605991.060660−1.05698385.66360
    1.001.4140691.414213−2.505175152.2993
    下载: 导出CSV

    表  2  测试数据记录表

    Table  2.   Record form of the test data

    NumberTheoretical setting valueActual measured value
    Azimuth/(°)Pitch/(°)Azimuth/(°)Pitch/(°)
    1 −1.4000 1.4000 −1.4537 1.3453
    2 −1.4000 0.9333 −1.4535 0.9733
    3 −1.4000 0.4667 −1.3556 0.4267
    4 −1.4000 0 −1.3562 −0.0362
    5 −1.4000 −0.4667 −1.4053 −0.5067
    6 −1.4000 −0.9333 −1.3536 −0.9833
    7 −1.4000 −1.4000 −1.3934 −1.4532
    8 −0.9333 1.4000 −0.9833 1.4263
    9 −0.9333 0.9333 −0.9333 0.8933
    10 −0.9333 0.4667 −0.9833 0.4267
    11 −0.9333 0 −0.9833 0.0243
    12 −0.9333 −0.4667 −0.9433 −0.4367
    13 −0.9333 −0.9333 −0.8933 −0.9835
    14 −0.9333 −1.4000 −0.8933 −1.4521
    15 −0.4667 1.4000 −0.5067 1.3421
    16 −0.4667 0.9333 −0.4767 0.9744
    17 −0.4667 0.4667 −0.5069 0.4267
    18 −0.4667 0 −0.4452 0.0021
    19 −0.4667 −0.4667 −0.5049 −0.5062
    20 −0.4667 −0.9333 −0.4665 −0.9833
    21 −0.4667 −1.4000 −0.5024 −1.4520
    22 0 1.4000 0.0521 1.4150
    23 0 0.9333 0.0213 0.9733
    24 0 0.4667 0.0189 0.4267
    25 0 0 0.0238 −0.0321
    26 0 −0.4667 0.0457 −0.5064
    27 0 −0.9333 0.0562 −0.8933
    28 0 −1.4000 0.0286 −1.4569
    29 0.4667 1.4000 0.4667 1.3498
    30 0.4667 0.9333 0.5062 0.9733
    31 0.4667 0.4667 0.4052 0.4267
    32 0.4667 0 0.4462 0.0245
    33 0.4667 −0.4667 0.5089 −0.5067
    34 0.4667 −0.9333 0.5011 −0.9811
    35 0.4667 −1.4000 0.4767 −1.4502
    36 0.9333 1.4000 0.9834 1.3488
    37 0.9333 0.9333 0.9456 0.9344
    38 0.9333 0.4667 0.9120 0.4261
    39 0.9333 0 0.9133 0.0430
    40 0.9333 −0.4667 0.8958 −0.5043
    41 0.9333 −0.9333 0.9726 −0.9725
    42 0.9333 −1.4000 0.9643 −1.4501
    43 1.4000 1.4000 1.4222 1.3862
    44 1.4000 0.9333 1.4220 0.9734
    45 1.4000 0.4667 1.4501 0.4962
    46 1.4000 0 1.4498 0.0321
    47 1.4000 −0.4667 1.4460 −0.4391
    48 1.4000 −0.9333 1.4508 −0.9882
    49 1.4000 −1.4000 1.3994 −1.3968
    下载: 导出CSV
  • [1] Zhang Hongyan, Xu Jin, Li Yanjun, et al. Application and analysis of the tactical laser weapon [J]. Laser & Infrared, 2007, 37(3): 204-207. (in Chinese)
    [2] Wang Jingtao, Zou Jie. Study on concept of missile defense using airborne directed energy weapon [J]. Infrared and Laser Engineering, 2007, 37(S2): 365-370. (in Chinese)
    [3] Yang Songtao, He Liqing. Some problems of the design of compound-axis servomechanism using single detector [J]. Laser & Infrared, 2020, 50(4): 457-462. (in Chinese)
    [4] Barnard T W, Fencil C R. Digital laser ranging and tracking using a compound axis servomechanism [J]. Applied Optics, 1966, 5(4): 497-505. doi:  10.1364/AO.5.000497
    [5] Yang Xiulin, Lu Peiqiang, Liu Xiaoqiang, et al. Simulation and analysis of the composite axis track control system for shipborne laser weapon [J]. Laser & Infrared, 2015, 45(8): 943-947. (in Chinese)
    [6] Li Jianhua, Li Zhifeng, Niu Zhenhong, et al. A distortion correction method for large field-of-view star sensor [J]. Missiles and Space Vehicles, 2017, 1: 103-106. (in Chinese)
    [7] Wang Pei, Li Yanjun, Tian Jin. Simulation system and analysis of airborne laser weapon [J]. Infrared and Laser Engineering, 2011, 40(7): 1238-1242. (in Chinese)
    [8] Wang Weibing, Jiang Zhenhua, Wang Tingfeng, et al. Op-to-electrical compound axis tracking control technologyon statellite [J]. Infrare and Laser Engineering, 2014, 43(12): 4133-4139. (in Chinese)
    [9] Lan Shuo, Li Xinnan, Wu Chunfeng, et al. Analysis on thermal performance of fast steering mirror irradiated by high power CW laser [J]. Infrare and Laser Engineering, 2018, 47(10): 1020003. (in Chinese) doi:  10.3788/IRLA201847.1020003
    [10] Chen Jian, Wang Weiguo, Liu Tingxia, et al. The operational principle analysis of new airborne Electro-Optical Targeting System (EOTS) of US army [J]. Chinese Optics, 2017, 10(6): 777-781. (in Chinese) doi:  10.3788/co.20171006.0777
  • [1] 周丽军, 杨通, 程德文, 王涌天.  基于方域正交多项式自由曲面的成像系统设计方法(特邀) . 红外与激光工程, 2023, 52(7): 20230317-1-20230317-14. doi: 10.3788/IRLA20230317
    [2] 何信欣, 刘斌, 王春柳, 吴冠豪.  用于高反表面测量的分块平滑自适应条纹投影方法 . 红外与激光工程, 2023, 52(5): 20220825-1-20220825-14. doi: 10.3788/IRLA20220825
    [3] 李小明, 朱国帅, 郭名航, 刘赢泽, 张崇.  基于光学自准直的旋转轴平行度测量与不确定度分析 . 红外与激光工程, 2023, 52(5): 20220794-1-20220794-12. doi: 10.3788/IRLA20220794
    [4] 雷鹏, 胡金春, 朱煜, 杜胜武.  傅氏级数和多项式结合的图像灰度建模及位移测量 . 红外与激光工程, 2022, 51(3): 20210123-1-20210123-8. doi: 10.3788/IRLA20210123
    [5] 朱建春, 李欣, 朱威.  流动卫星激光测距系统的距离选通实现方法 . 红外与激光工程, 2022, 51(12): 20220200-1-20220200-11. doi: 10.3788/IRLA20220200
    [6] 王云新, 许家豪, 周涛, 王大勇, 杨登才, 钟欣.  采用双平行马赫曾德调制器的四倍频信号产生 . 红外与激光工程, 2018, 47(9): 918009-0918009(8). doi: 10.3788/IRLA201847.0918009
    [7] 张磊, 邱伟, 张凯.  基于双五棱镜组件的大间距光轴平行性检测方法 . 红外与激光工程, 2018, 47(7): 717005-0717005(5). doi: 10.3788/IRLA201847.0717005
    [8] 张世国, 方海涛, 汪玮, 王敏, 王毛翠, 刘振.  基于积分球分光与接收的透射式能见度测量系统 . 红外与激光工程, 2018, 47(10): 1017003-1017003(6). doi: 10.3788/IRLA201847.1017003
    [9] 刘显著, 王超, 江伦, 刘壮, 杨进华, 姜会林.  二维多项式位相光瞳滤波实现超分辨望远成像 . 红外与激光工程, 2018, 47(4): 418007-0418007(6). doi: 10.3788/IRLA201847.0418007
    [10] 王涛, 赵建科, 田留德, 周艳, 杨利红, 陈琛, 段亚轩, 潘亮, 赵怀学, 刘锴, 万伟, 刘艺宁.  光电编码器与棱体轴线平行度对转角误差的影响 . 红外与激光工程, 2018, 47(2): 217001-0217001(8). doi: 10.3788/IRLA201847.0217001
    [11] 王阳, 黄煜, 李占峰, 王淑荣.  利用恒星对天文观测系统光轴平行性检校 . 红外与激光工程, 2017, 46(5): 517003-0517003(6). doi: 10.3788/IRLA201746.0517003
    [12] 汪同浩, 刘秉琦, 黄富瑜, 陈一超, 张帅.  平行式红外双目立体系统各参数合理效益取值 . 红外与激光工程, 2017, 46(9): 904004-0904004(6). doi: 10.3788/IRLA201746.0904004
    [13] 李召龙, 沈同圣, 娄树理.  基于多项式逼近的红外系统渐晕效应校正方法 . 红外与激光工程, 2016, 45(S1): 6-10. doi: 10.3788/IRLA201645.S104002
    [14] 范有臣, 赵洪利, 孙华燕, 郭惠超, 赵延仲.  主被动复合激光距离选通探测系统最大作用距离的推算 . 红外与激光工程, 2015, 44(S1): 86-92.
    [15] 江伦, 胡源, 董科研, 安岩, 王超, 佟首峰.  红外双波段光学系统被动式消热差设计 . 红外与激光工程, 2015, 44(11): 3353-3357.
    [16] 周雷刚, 梁庭, 高利聪.  影响全天域光轴平行度检测精度因素分析 . 红外与激光工程, 2015, 44(11): 3413-3417.
    [17] 孙文卿, 陈磊, 李金鹏, 乌兰图雅, 何勇.  非圆孔径离散采样点正交多项式波前拟合 . 红外与激光工程, 2015, 44(3): 1068-1072.
    [18] 张丽敏, 吕天宇, 杨飞.  1.5m口径垂直状态短筒长焦距平行光管次镜支撑系统设计 . 红外与激光工程, 2013, 42(12): 3283-3288.
    [19] 王奇涛, 佟首峰, 徐友会.  采用Zernike多项式对大气湍流相位屏的仿真和验证 . 红外与激光工程, 2013, 42(7): 1907-1911.
    [20] 修金利, 金伟其, 刘崇亮, 王霞.  应用U形边框黑体光阑的三点辐射定标校正方法及其分析 . 红外与激光工程, 2013, 42(9): 2313-2318.
  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  147
  • HTML全文浏览量:  36
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-10-25
  • 修回日期:  2021-11-24
  • 录用日期:  2022-02-14
  • 网络出版日期:  2022-08-13
  • 刊出日期:  2022-08-05

主被动双光复合系统平行度校正研究

doi: 10.3788/IRLA20210771
    作者简介:

    刘飞,男,工程师,硕士,主要研究方向为激光对抗系统总体设计

  • 中图分类号: TN209

摘要: 针对主被动双光复合系统中存在激光主动发射视场与红外被动成像视场平行度(一一对应位置)偏差问题,提出了主被动双光复合系统平行度校正方法,并开展了理论研究和实物仿真测试。首先,推导了激光主动发射视场与快速反射镜(Fast Steering Mirror,FSM)角度对应关系;其次,根据激光发射视场与红外探测视场提出了一种基于多项式拟合的平行度修正校正方法;最后,通过构建实际测试平台,进行平行度标定、求解、修正。试验结果表明:在±1.4°的激光指示视场和±1.4°的红外中心视场中,经过系统平行度修正后,实际平行度偏差可由0.8 mrad降低至0.1 mrad,进一步提高激光对目标的指向精度。

English Abstract

    • 随着军事需求和光电对抗技术的不断发展,越来越多的光电对抗系统开始应用在机载、弹载、星载等平台上。为了实现对敌方光电探测器的有效持续干扰、损伤,该类光电对抗系统需要对目标进行高精度和持续跟踪、瞄准,实现对目标的干扰或破坏[1-3]。为了实现对运动目标的高精度跟踪或瞄准功能,一般采用复合轴跟踪系统,最早的复合轴跟踪系统发表于1966年Thomas W.的文章[4]。传统的复合轴跟踪系统包含双探测器型复合轴系统和单探测器型复合轴系统,其中双探测器型复合轴系统可实现高精度跟瞄(一般可实现几个微弧度跟瞄精度),技术成熟。但其光路设计难度大、体积大等缺点往往限制其在小平台上的应用。单探测器型复合轴系统在跟踪精度方面不如双探测器型复合轴系统,但其具有构成简单、装配难度低、功耗低、稳定性和可靠性高等优点,便于小型化[5]

      采用单探测器跟踪系统,成像与激光光路分离设计,如图1所示。该设计不仅紧凑、轻小,容易满足弹载、机载等平台的装机条件,而且光学设计难度小,可实现多目标跟踪、干扰。由于激光与图像没有采用共光路设计,激光无法跟随红外目标实现图像闭环控制,致使激光跟瞄精度不高;并且激光与成像系统间存在机械安装、光学镜头畸变、两轴平行度、快反镜正交性等多种偏差,导致光轴一致性差。目前,国内外针对多光轴一致性常用测量方法有:离轴抛物镜法,平行光管法,激光光轴仪法,五棱镜法和投影靶板法等[6-9]。上述方法虽然测量精度高,但是在实际应用过程中,如何在已知测量误差的前提下通过装调保证光轴一致性是亟需要解决的问题。装调方式一般有两种:一种是提高机械结构件的加工精度、辅助测量和微装配,另一种是增加监控手段,实时观测轴系间的变化情况再辅助增加调整机构进行动态调整[10]。上述方式虽然可实现光轴一致性,但是也存在着一些问题,例如对加工精度要求极高,微调的可操作性和稳定性不强,增加系统的复杂性等。

      图  1  主被动双光复合系统示意图

      Figure 1.  Schematic of the active and passive dual-band (laser/IR) composite system

      针对上述问题,采用单探测器跟踪系统,成像、激光发射光路分离设计的应用方式其光轴一致性不同于以往系统,该系统要实现成像、激光多点匹配问题,并不仅仅是中心轴系的一致性,为此设计了一种平行度校正补偿方法,通过试验标定和软件算法校正,使红外成像、激光在相同视场范围内满足平行度小于0.1 mrad偏差,从而提升激光对目标的跟瞄精度,该方法对装调精度要求不高且可操作性和稳定性强。

    • 在成像光路与激光光路分离设计方案中,其工作原理为:通过红外成像系统获取目标位置信息,将目标位置信息与系统姿态信息进行融合处理后控制快反镜(FSM)转动,同时开启激光,激光束随快反镜(FSM)转动快速指向目标位置并持续跟踪,实现对目标光电载荷的精确跟瞄。完成上述过程,需明确以下两个问题:首先是快反镜的空间旋转与成像平面的对应几何关系,其次是快反镜旋转轴与成像多点对应关系。

    • 在实际应用过程中,快反镜(FSM)镜面与激光入射方向呈45°安装,激光入射光固定,通过FSM转动,实现激光光束在空间中的指向,如图1所示。

      图2所示,建立空间坐标系O-XYZ。分别确立FSM镜面旋转方位角、旋转俯仰角与出射光线的角度关系,定义FSM镜面方位旋转角为α,俯仰旋转角为β;出射光方位角为θ方位,俯仰角为θ俯仰

      图  2  光线反射与镜面位置关系

      Figure 2.  Relationship between light reflection and mirror position

      在方位旋转方向上,如图2(a)所示,OA为激光原始入射方向,OE为激光原始出射方向,OB为法线,将快反镜镜面沿着平面XOY顺时针旋转角度αOF为旋转后的出射光线,OC为法线OB经旋转后的法线,即∠BOCα。由于镜面方位旋转后,出射光OF已经不在ZOY的方位平面内,也就是说镜面的方位旋转引起了光线的俯仰运动,因此为了得到方位角与出射光方位角的关系,需要将出射光OF进行方位方向和俯仰方向投影,可得∠COF方位投影和∠COF俯仰投影,关系如下:

      $$ \angle {{COF}}_{俯仰投影}\text{=arctan}\left( {\frac{\text{1-cos}\alpha }{\text{cos}\alpha }} \right) $$ (1)
      $$ \angle {{COF}}_{方位投影}\text{=arcsin}\left( {\frac{\sqrt{\text{2}}\times \mathrm{sin}\alpha }{\sqrt{\text{1+}{\mathrm{sin}}^{2}\alpha }}} \right) $$ (2)

      同理,在俯仰旋转方向上,如图2(b)所示,由空间位置关系可知出射光OFXOY平面内,即∠COF=θ俯仰。可知出射光俯仰角与β的关系:

      $$ {\theta }_{俯仰}\text=2\beta $$ (3)

      综上,激光出射角与快反镜旋转角的空间位置关系如下:

      $$ {\theta }_{方位}\text{=arcsin}\left( {\frac{\sqrt{\text{2}}\times \text{sin}\alpha }{\sqrt{\text{1+}{\mathrm{sin}}^{2}\alpha }}} \right) $$ (4)
      $$ {\theta }_{俯仰}\text=2\beta +\text{arctan}\left( {\frac{\text{1-cos}\alpha }{\text{cos}\alpha }} \right) $$ (5)

      式中:激光出射角θ方位θ俯仰即为目标在被动探测系统视场中的方位角和俯仰角,其值可通过目标识别算法计算得出。αβ为快反镜输入角度值,需要根据公式(4)和公式(5)进行解算。为了保持快反镜高动态特性和激光的指示范围,快反镜的转角范围设置在±1°以内,即αβ取值范围为−1°~+1°。

      由公式(1)和公式(2)可知,快反镜在进行方位摆动时,不仅会引起激光的方位方向的改变,而且也会引起激光俯仰方向的改变。快反镜方位摆动引起激光方位、俯仰变化如图3所示。

      图  3  快反镜方位摆动与激光方位、激光俯仰投影关系

      Figure 3.  Projection relationship between the azimuth swing of the FSM and the laser azimuth, laser pitch

      快反镜在−1°~+1°范围内进行方位轴摆动时,不仅引起激光方位角摆动,而且会引起激光俯仰角摆动。其中激光方位角变化范围为−1.414°~+1.414°,激光俯仰角变化范围为152.299 ~0 μrad。由图3可知,激光方位角与快反镜方位角间的变化规律可近似为二次曲线,斜率取值范围基本在$ \sqrt 2 $上下浮动。若激光方位角与快反镜方位角以$ \sqrt 2 $倍关系近似计算时,在整个摆动范围内引起的最大误差为2.505 μrad,其计算结果如表1所示。

      表 1  快反镜方位摆动引起激光方位、俯仰变化计算结果

      Table 1.  Calculation results of the changes in laser azimuth and pitch caused by the azimuth swing of the FSM

      FSM azimuth angle
      /(°)
      Laser azimuth angleLaser pitch angle
      /μrad
      Theoretical calculation
      /(°)
      Approximate calculation
      /(°)
      Error Δ
      /μrad
      −1.00−1.414069−1.4142132.505175152.2993
      −0.75−1.060599−1.0606601.05698385.66360
      −0.50−0.707088−0.7071060.31320438.07120
      −0.25−0.353551−0.3535530.0391529.517574
      0.000000
      0.250.3535510.353553−0.0391529.517574
      0.500.7070880.707106−0.30320438.07120
      0.751.0605991.060660−1.05698385.66360
      1.001.4140691.414213−2.505175152.2993

      实际工程中,为了便于对α, β计算,对公式(4)和公式(5)进行近似处理,即快反镜旋转角(二维)与激光出射角(二维)可近似为公式(6)和公式(7):

      $$ \alpha =\frac{\sqrt{2}}{2}{\theta }_{方位} $$ (6)
      $$ \beta =\frac{1}{2}{\theta }_{俯仰} $$ (7)

      上述处理过程中引入了一定的计算误差量和旋转增量,将会导致成像、激光光轴对应平行度偏差。其次,在成像光路、激光光路分离设计中,存在成像镜头的光学窗口畸变、快反镜的正交偏差以及传感器安装过程中的机械偏差等因素均会影响成像、激光光轴平行度偏差,上述偏差因素难以精确测量。因此,为了实现激光对目标的高精度指向,上述偏差均需要进行补偿修正。

    • 为了实现上述偏差的补偿修正,需要测量系统中的上述固有偏差。首先在红外视场中心向外均匀选取N个点(αi, βi)(i=1, 2, ……, N)作为快反镜的输入角度值;开启激光光源,快反镜将按照N个点的位置摆动;激光束经过反射镜后形成N个点源;红外热像仪可通过反射镜观测到49激光光斑,并计算出激光光斑的实际位置(α'i, β'i)(i=1,2, ……, N)。基于上述(αi, βi)和(α'i, β'i)两组数据进行公式(8)和公式(9)的计算,采用4阶多项式最小二乘法拟合,可计算得出系数k0~k29值。

      $$ \begin{split} \\ \alpha _i{'}{\text{ = }}&{k_0}\beta _i^4 + {k_1}\beta _i^3{\alpha _i} + {k_2}\beta _i^2\alpha _i^2 + {k_3}{\beta _i}\alpha _i^3 + {k_4}\alpha _i^4 + \\ &{k_5}\beta _i^3 + {k_6}\beta _i^2{\alpha _i} + {k_7}{\beta _i}\alpha _i^2 + {k_8}\alpha _i^3 + \\ &{k_9}\beta _i^2 + {k_{10}}{\beta _i}{\alpha _i} + {k_{11}}\alpha _i^2 + \\ &{k_{12}}{\beta _i} + {k_{13}}{\alpha _i} + \\ &{k_{14}} \\[-2pt] \end{split} $$ (8)
      $$ \begin{split} {\beta _i{'}}{\text{ = }}&{k_{15}}\beta _i^4 + {k_{16}}\beta _i^3{\alpha _i} + {k_{17}}\beta _i^2\alpha _i^2 + {k_{18}}{\beta _i}\alpha _i^3 + {k_{19}}\alpha _i^4 + \\ &{k_{20}}\beta _i^3 + {k_{21}}\beta _i^2{\alpha _i} + {k_{22}}{\beta _i}\alpha _i^2 + {k_{23}}\alpha _i^3 + \\ &{k_{24}}\beta _i^2 + {k_{25}}{\beta _i}{\alpha _i} + {k_{26}}\alpha _i^2 + \\ &{k_{27}}{\beta _i} + {k_{28}}{\alpha _i} + \\ &{k_{29}} \\[-2pt] \end{split} $$ (9)

      公式(8)、(9)写成矩阵形式,即:

      $$ {A{'}}{B{'}} = \left( {\begin{array}{*{20}{c}} \begin{gathered} {\alpha {'}}_i \\ \vdots \\ {\alpha {'}}_1 \\ {\alpha {'}}_0 \\ \end{gathered} &\begin{gathered} {\beta {'}}_i \\ \vdots \\ {\beta {'}}_1 \\ {\beta {'}}_0 \\ \end{gathered} \end{array}} \right) $$ (10)
      $$ {H} = \left( {\begin{array}{*{20}{c}} {\beta _i^4}&{\beta _i^3{\alpha _i}}& \cdots &1 \\ \vdots & \vdots & \vdots & \vdots \\ {\beta _1^4}&{\beta _1^3{\alpha _1}}& \cdots &1 \\ {\beta _0^4}&{\beta _0^3{\alpha _0}}& \cdots &1 \end{array}} \right) $$ (11)
      $$ K = \left( {\begin{array}{*{20}{c}} \begin{gathered} {k_0} \\ \vdots \\ {k_{13}} \\ {k_{14}} \\ \end{gathered} &\begin{gathered} {k_{15}} \\ \vdots \\ {k_{28}} \\ {k_{29}} \\ \end{gathered} \end{array}} \right) $$ (12)

      公式(10)~(12)可写成:

      $$ {A{'}}{B{'}} = {H} \times K $$ (13)

      由公式(13)可得系数矩阵K

      $$ K = {({ H})^{ - 1}}{A{'}}{B\;{'}} $$ (14)

      式中:系数矩阵K即代表了系统中红外给定位置量与激光实际指向位置量的对应关系,属于系统固有属性,为有效改善该种对应关系,使红外给定位置量与激光实际指向位置量一一对应,实际计算中系数矩阵以K−1参与计算。因此,可得如下计算关系:

      $$ {A}{B} = {H} \times {K^{ - 1}} $$ (15)

      式中:矩阵H中的(αi, βi)需满足公式(6)和公式(7)的关系,矩阵AB即为快反镜的位置给定值。

    • 平行度测试平台如图4所示,主要包含被测设备(凝视热像仪、快反镜、激光源等)、离轴反射镜、激光靶面、显控计算机、目标识别器、指令操控台、光学平台及安装组件等。

      图  4  平行度测试原理图

      Figure 4.  Schematic diagram of the parallelism test

      被测设备中,凝视红外视场角为8°×10°,瞬时视场角为0.5 mrad;FSM机械偏转角为±1°,角度分辨率为15 μrad,激光发散角为0.5 mrad;红外与激光中心间距为80 mm (由于实际应用过程中目标在几公里以外,中心偏差产生的误差可不考虑)。离轴反射镜焦距为2 m,口径250 mm;激光与红外波段相同。

      首先,将红外中心视场进行网格划分,取与激光偏摆视场范围内的若干个点,作为激光偏摆的基准点。然后,将被测系统放置光学平台上,在前方放置离轴反射镜,在离轴反射镜的焦点上放置靶面或靶板。被测系统开机,开启激光并给FSM加电驱动,由于快反镜上电后自行进入零位,激光经快反镜反射后也即保持在零位,该零位依靠结构精度保持与红外中心一致,如图5所示。在红外相机输出的图像上可观测到靶板上的激光束光斑,点目标识别器采用质心处理算法对光斑在红外视场中的位置计算,可知激光束零点位置与红外视场中心位置偏差量。同理,激光束按照在红外视场上确定的49个基准点,依靠FSM进行方位、俯仰的偏摆,形成多个激光束位置,每偏摆到一个位置,点目标识别器即可计算出当前激光束的在红外视场中的位置,每个基准点位置与实际位置偏差即为平行度偏差。测量结果如表2所示。

      图  5  平行度校正前、后给定值与测量值对比图

      Figure 5.  Comparison chart of the given values and measured values before and after parallelism correction

      表 2  测试数据记录表

      Table 2.  Record form of the test data

      NumberTheoretical setting valueActual measured value
      Azimuth/(°)Pitch/(°)Azimuth/(°)Pitch/(°)
      1 −1.4000 1.4000 −1.4537 1.3453
      2 −1.4000 0.9333 −1.4535 0.9733
      3 −1.4000 0.4667 −1.3556 0.4267
      4 −1.4000 0 −1.3562 −0.0362
      5 −1.4000 −0.4667 −1.4053 −0.5067
      6 −1.4000 −0.9333 −1.3536 −0.9833
      7 −1.4000 −1.4000 −1.3934 −1.4532
      8 −0.9333 1.4000 −0.9833 1.4263
      9 −0.9333 0.9333 −0.9333 0.8933
      10 −0.9333 0.4667 −0.9833 0.4267
      11 −0.9333 0 −0.9833 0.0243
      12 −0.9333 −0.4667 −0.9433 −0.4367
      13 −0.9333 −0.9333 −0.8933 −0.9835
      14 −0.9333 −1.4000 −0.8933 −1.4521
      15 −0.4667 1.4000 −0.5067 1.3421
      16 −0.4667 0.9333 −0.4767 0.9744
      17 −0.4667 0.4667 −0.5069 0.4267
      18 −0.4667 0 −0.4452 0.0021
      19 −0.4667 −0.4667 −0.5049 −0.5062
      20 −0.4667 −0.9333 −0.4665 −0.9833
      21 −0.4667 −1.4000 −0.5024 −1.4520
      22 0 1.4000 0.0521 1.4150
      23 0 0.9333 0.0213 0.9733
      24 0 0.4667 0.0189 0.4267
      25 0 0 0.0238 −0.0321
      26 0 −0.4667 0.0457 −0.5064
      27 0 −0.9333 0.0562 −0.8933
      28 0 −1.4000 0.0286 −1.4569
      29 0.4667 1.4000 0.4667 1.3498
      30 0.4667 0.9333 0.5062 0.9733
      31 0.4667 0.4667 0.4052 0.4267
      32 0.4667 0 0.4462 0.0245
      33 0.4667 −0.4667 0.5089 −0.5067
      34 0.4667 −0.9333 0.5011 −0.9811
      35 0.4667 −1.4000 0.4767 −1.4502
      36 0.9333 1.4000 0.9834 1.3488
      37 0.9333 0.9333 0.9456 0.9344
      38 0.9333 0.4667 0.9120 0.4261
      39 0.9333 0 0.9133 0.0430
      40 0.9333 −0.4667 0.8958 −0.5043
      41 0.9333 −0.9333 0.9726 −0.9725
      42 0.9333 −1.4000 0.9643 −1.4501
      43 1.4000 1.4000 1.4222 1.3862
      44 1.4000 0.9333 1.4220 0.9734
      45 1.4000 0.4667 1.4501 0.4962
      46 1.4000 0 1.4498 0.0321
      47 1.4000 −0.4667 1.4460 −0.4391
      48 1.4000 −0.9333 1.4508 −0.9882
      49 1.4000 −1.4000 1.3994 −1.3968
    • 根据测试平台得到的标定值和测量值,可计算出修正系数ki。在红外给定目标设置值时,利用4阶多项式校正公式对设定值进行修正,再输出给FSM执行。按照上述测试流程进行再次测量,可形成前、后的对比测试结果,如图5所示。从图中可以直观的看出经过校正后,给定值与测量值基本重合,即可实现激光与红外对应位置的一一对应。

      为了定量分析校正前、后激光与红外的平行度偏差,可采用公式(16)进行评价,其中$ \delta $为平行度偏差,表示为:

      $$ \delta = \frac{1}{n}\sum\limits_{i = 0}^{i = n} {\sqrt {{{(\alpha _i' - {\alpha _i})}^2} + {{(\beta _i' - {\beta _i})}^2}} } $$ (16)

      平行度修正前,激光与红外的平行度偏差为0.86 mrad;修正后,平行度偏差为0.09 mrad,平行度偏差性能提升了9倍多。

    • 文中对主被动双光复合系统光轴平行度校正方法进行了理论和实验分析。针对近似计算产生的误差、红外与激光间的机械安装误差以及红外成像器的光学畸变偏差等误差因素,提出了多项式拟合的修正补偿方法。与传统的单纯靠机械结构来保障双光视场对应关系相比,该方式不仅具有精度高、便于操作等优点,还可以有效的消除传感器自身因素引起的各种偏差,而且在实际应用过程中稳定性很强。通过试验标定获取修正偏差,计算出多项式拟合函数的修正系数,试验结果表明修正补偿前、后平行度偏差性能提升9倍多。按上述方法校正修正后,主被动双光复合系统可实现激光对目标光电载荷的精确指向,从而进一步压缩激光出射发散角,提升激光干扰距离和干扰效果。

参考文献 (10)

目录

    /

    返回文章
    返回