-
随着军事需求和光电对抗技术的不断发展,越来越多的光电对抗系统开始应用在机载、弹载、星载等平台上。为了实现对敌方光电探测器的有效持续干扰、损伤,该类光电对抗系统需要对目标进行高精度和持续跟踪、瞄准,实现对目标的干扰或破坏[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(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)所示,由空间位置关系可知出射光OF在XOY平面内,即∠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 angle Laser pitch angle
/μradTheoretical calculation
/(°)Approximate calculation
/(°)Error Δ
/μrad−1.00 −1.414069 −1.414213 2.505175 152.2993 −0.75 −1.060599 −1.060660 1.056983 85.66360 −0.50 −0.707088 −0.707106 0.313204 38.07120 −0.25 −0.353551 −0.353553 0.039152 9.517574 0.00 0 0 0 0 0.25 0.353551 0.353553 −0.039152 9.517574 0.50 0.707088 0.707106 −0.303204 38.07120 0.75 1.060599 1.060660 −1.056983 85.66360 1.00 1.414069 1.414213 −2.505175 152.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所示,主要包含被测设备(凝视热像仪、快反镜、激光源等)、离轴反射镜、激光靶面、显控计算机、目标识别器、指令操控台、光学平台及安装组件等。
被测设备中,凝视红外视场角为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
Number Theoretical setting value Actual 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倍多。按上述方法校正修正后,主被动双光复合系统可实现激光对目标光电载荷的精确指向,从而进一步压缩激光出射发散角,提升激光干扰距离和干扰效果。
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,进一步提高激光对目标的指向精度。Abstract: Aiming at the deviation of the parallelism(on to one correspondence position) between the active laser field of view and the infrared passive field of view in the parallelism correction of the active and passive dual-band (laser/IR) composite system, the parallelism correction method of the active and passive dual-band (laser/IR) composite system was proposed and the theoretical research and physical simulation test were carried out. First, derived the angle correspondence between the active laser emission field of view and the fast steering mirror. Secondly, a parallelism correction method based on polynomial fitting is proposed according to the laser emission field of view and infrared detection field of view. Finally, the actual test platform is constructed to calibrate, solve and correct. The test results showed that the actual maximum deviation can be reduced from 0.8 mrad to less than 0.1 mrad in the ±1.4° indicated field of view of laser and ±1.4° infrared central field of view, which further improve the aiming accuracy of the laser.
-
表 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 angle Laser pitch angle
/μradTheoretical calculation
/(°)Approximate calculation
/(°)Error Δ
/μrad−1.00 −1.414069 −1.414213 2.505175 152.2993 −0.75 −1.060599 −1.060660 1.056983 85.66360 −0.50 −0.707088 −0.707106 0.313204 38.07120 −0.25 −0.353551 −0.353553 0.039152 9.517574 0.00 0 0 0 0 0.25 0.353551 0.353553 −0.039152 9.517574 0.50 0.707088 0.707106 −0.303204 38.07120 0.75 1.060599 1.060660 −1.056983 85.66360 1.00 1.414069 1.414213 −2.505175 152.2993 表 2 测试数据记录表
Table 2. Record form of the test data
Number Theoretical setting value Actual 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 -
[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 -