Calculation of Temperature and Radiation Characteristics of Midcourse Ballistic Missiles
-
摘要: 中段飞行弹道导弹目标辐射特性对于探测手段选择、传感器设计等具有重要意义。以采用红外隐身技术的弹道导弹为研究对象,考虑表面隐身涂层特性和进出地影的影响,采用有限体积法计算了弹道导弹在整个中段飞行过程中的温度变化;结合导弹自身辐射与太阳、地球辐射光谱,给出了导弹0~15 μm范围内的辐射特性;研究了导弹表面温度和辐射特性的关系,不同表面涂层在光照与地影状态下的探测差异。结果表明,相比于表面温度,涂层光学参数对导弹辐射特性的影响更大,同一涂层的可见光和红外突防效能存在矛盾,可以采用两种探测方法协作,提高探测能力。Abstract: The radiation characteristics of mid-course ballistic missiles are the basis and premise for their detection and identification. Radiation characteristics have an important guiding meaning in the selection of detection methods, sensor design, etc. Taking a ballistic missile with infrared stealth technology as the research object, and considering the factors of stealth coatings and the influence of earth shadows, the temperature and its variation trend are calculated using the finite-volume method. Combining the radiation of the missile with solar and earth radiation spectra, the radiation characteristics are presented within for missile wavelengths ranging from 0–15 μm. The relationship between the surface temperature and the radiation intensity of the missile is studied. The detection differences of different surface coatings under illumination and shadow conditions are discussed. The results show that the optical parameters of the coating have a greater influence on the radiation characteristics of the missile than the surface temperature. There is disagreement between the penetration effectiveness in the visible and infrared wavebands. Visible and infrared detection methods can be used simultaneously to improve detection capability.
-
Keywords:
- ballistic missiles /
- radiation characteristics /
- stealth coatings
-
0. 引言
近年来,随着红外制冷探测器朝尺寸小、重量轻、功耗低、成本低即低SWaP-C(size, weight and power, cost)方向快速发展,集成小型制冷机的该类中波红外探测器国外已大量应用于武器热瞄镜、便携式手持热像仪、小型无人机、无人车、遥控狙击手和遥控武器站、导弹导引头等空间受限的红外系统[1]。针对该类小型探测器设计一款仅由4片透镜组成的尺寸小、重量轻、成本低的中波红外连续变焦光学系统进而生产结构紧凑、低功耗和低成本的红外变焦热像仪将在手持观瞄具、边防监视系统、小型无人系统等平台得到广泛应用。
目前采用四片式的红外变焦光学多为双视场变焦系统。文献[2]为四片式非制冷双视场,文献[3]为四片制冷型长波双视场系统,文献[4]为4片制冷型中波双视场系统,目前仅采用4片透镜实现连续变焦的制冷型红外光学系统未见报道。
从变焦理论分析,一般常用的机械补偿变焦系统由典型的前固定组、变倍组、补偿组、后固定组4组透镜组成。变倍组一般是负透镜,而补偿组可以是正透镜组也可以是负透镜组,前者为正组补偿系统,后者称为负组补偿系统。从光学系统像差校正难易程度、减少透镜数量、降低光学透镜成本考虑,本文采用无后固定组正组补偿变焦系统,即变焦部分由会聚目标光线的前固定组正透镜、变倍负透镜、补偿正透镜构成。为压缩前固定组物镜口径并满足100%冷屏效率,系统采用二次成像方案,利用单片正光焦度透镜将一次像再次中继成像到探测器焦平面。为压缩轴向尺寸,采用二片平面反射镜将光路U型折转,最终实现仅由4片透镜构成的轻小型中波红外连续变焦光学系统。
1. 连续变焦理论计算模型
机械补偿连续变焦光学系统参数求解就是确定变焦系统在满足像面稳定和焦距在一定范围内连续变化的条件下系统中各组元的焦距、间隔、位移量等参数。通过建立数学模型能方便地计算和分析变焦过程、确定变焦系统高斯光学参数[5]。二组元正组补偿连续变焦系统运动方式如图 1所示。
变焦系统由于只有运动组才产生像面位移,只需抽出变倍组、补偿组加以分析。
因变倍组f2′的移动,引起整个运动组分的像面移动为m32(1-m22)dq,因补偿组f3′的移动,引起整个运动组分的像面移动为(1-m32)dΔ。为达到像面稳定,两个运动组像面移动量的代数和必为零。
$$ m_3^2\left(1-m_2^2\right) \mathrm{d} q+\left(1-m_3^2\right) \mathrm{d} \varDelta=0$$ (1) 而变倍组f2′、补偿组f3′微分移动量dq、dΔ与其倍率变化dm2、dm3之间的关系为:
$$ \mathrm{d} q = \frac{{{f_2}'}}{{m_2^2}}\mathrm{d}{m_2} $$ (2) $$ \mathrm{d} \varDelta=f_3{ }^{\prime} \mathrm{d} m_3 $$ (3) 将(2)、(3)代入(1),经整理得到二组元连续变焦微分方程如下:
$$ \frac{{1 - m_2^2}}{{m_2^2}}{f_2}'\mathrm{d}{m_2} + \frac{{1 - m_3^2}}{{m_3^2}}{f_3}'\mathrm{d}{m_3} = 0 $$ (4) 式(4)是多变量全微分型微分方程,设U(m2, m3)为原函数,则有dU(m2, m3)=0。
其通解为:
$$ U({m_2}, {m_3}) = {f_2}'(\frac{1}{{{m_2}}} + {m_2}) + {f_3}'(\frac{1}{{{m_3}}} + {m_3}) = C $$ (5) 式中:C为常量;设变倍组f2′、补偿组f3′初始状态都处于系统长焦位置,则:
$$ m_{2}=m_{2l}\text{;}m_{3}=m_{3l} $$ (6) 同样有:
$$ {f_2}'(\frac{1}{{{m_{2l}}}} + {m_{2l}}) + {f_3}'(\frac{1}{{{m_{3l}}}} + {m_{3l}}) = C $$ (7) 消去常量C,得到方程的特解:
$$ {f_2}'(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + {f_3}'(\frac{1}{{{m_3}}} - \frac{1}{{{m_{3l}}}} + {m_3} - {m_{3l}}) = 0 $$ (8) 将(8)式整理得到补偿组f3′的倍率m3构成的二次方程:
$$ m_3^2 - b{m_3} + 1 = 0 $$ (9) 其中
$$ b = - \frac{{{f_2}'}}{{{f_3}'}}(\frac{1}{{{m_2}}} - \frac{1}{{{m_{2l}}}} + {m_2} - {m_{2l}}) + (\frac{1}{{{m_{3l}}}} + {m_{3l}}) $$ (10) 解得m3的两根为:
$$ {m_{31}} = \frac{{b + \sqrt {{b^2} - 4} }}{2} $$ (11) $$ {m_{32}} = \frac{{b - \sqrt {{b^2} - 4} }}{2} $$ (12) 系统参数求解过程如下:
1)将(2)式积分并整理得到变倍组f2′的倍率m2;
$$ {m_2} = \frac{1}{{\frac{1}{{{m_{2l}}}} + \frac{{{q_2}}}{{{f_2}'}}}} $$ (13) 2)根据求得的m2按照公式(10)求出系数b,再由(11)、(12)式解得补偿组f3′满足运动方程的两个解m31和m32。
3)根据补偿组的两个解求出满足运动方程补偿像面位移的移动量Δ1和Δ2:
$$ \varDelta_{1}=f_{3}′(m_{31}-m_{3l})$$ (14) $$ \varDelta_{2}=f_{3}′(m_{32}-m_{3l}) $$ (15) 4)求出系统的总变倍比
$$ {\varGamma _1} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{31}}}} $$ (16) $$ {\varGamma _2} = \frac{{{m_{2l}}{m_{3l}}}}{{{m_2}{m_{32}}}} $$ (17) 变倍组f2′每移动q2对应着Γ1和Γ2,变倍组和补偿组一起同步运动直到预定的总变倍比为止,到达变倍比的要求的最终状态为系统短焦位置。根据上述连续变焦微分模型,利用(10)~(17)式,可解得各组分光焦度分配及间隔位置关系。
2. 连续变焦系统设计流程
针对复杂的变焦光学系统,建立模型是至关重要的环节,尽管根据(10)~(17)式可推导求解变焦系统高斯参数,但是当光学系统初始参数选择不合适, 会使得系统光焦度分配不合理导致系统像差校正难度大、间隔不合适导致各组件运动中相互碰撞等情况发生。鉴于连续变焦系统的复杂性,系统建模难度大,因此根据连续变焦理论模型编制连续变焦参数计算程序,辅助建立理想光学模型。建模后,设计工作的重点将放在选型以及评价函数的设置和动态修改上,使得设计系统快速收敛[6]。连续变焦光学系统设计流程如图 2所示。
首先,根据设计指标要求需要建立理想光学模型,确定每个组元的参数;再合理选型选材、设定评价函数,进入优化和全局优化;最后根据设计结果评价成像质量。其中函数优化和像质评价环节反复多次迭代,直至达到设计技术指标要求。
3. 连续变焦光学系统设计
3.1 设计指标
中波制冷型连续变焦光学系统采用昆明物理研究所研制生产的长度方向仅119 mm的小型化制冷型中波红外640×512焦平面探测器组件,该探测器具体参数如表 1所示。连续变焦光学系统主要设计指标见表 2。
表 1 探测器参数Table 1. Parameters of detectorDetector HgCdTe Array size 640×512 Pixel dimension/μm 15 NETD/mK ≤22 Spectral response/μm 3.7−4.8 μm Weight ≤380 g 表 2 光学系统设计指标Table 2. Parameters of optical system designWorking waveband/μm 3.7 to 4.8 Zoom 10:01 Field of view 20°×16° to 2.0°×1.6° F# 4 Focal length/mm 27.0~275 Working temperature/℃ -40 to 60 3.2 设计过程
首先,按照系统指标要求,根据连续变焦理论模型编制连续变焦正组补偿参数计算程序求解变焦光学系统各组元光焦度分配及位置间隔关系。在系统初始参数取值上主要考虑以下几点:
1)系统采用二次成像,既可以压缩前固定组直径又可以满足100%冷屏效率。中继组初始倍率取为-1×,前端变焦核心按照系统实际变焦范围进行取值,无需缩放;
2)系统无后固定组,在理论求解中将后固定组倍率取为1×即将后固定组定为无光焦度的虚拟面;
3)为压缩系统变焦过程中变倍组、补偿组位移量,变倍组长焦初始倍率取较大的倍率值,以符合最速变焦理论;
4)考虑光路U型折转,补偿组和中继组之间需较大的空间安置两片平面反射镜,则补偿组取较大的焦距值;
5)系统处于短焦位置时前固定组与变倍组应留出足够的间隔,使两组透镜不至于相碰,初始值设为0.55,补偿组与无光焦度后固定组虚拟面距离初始设为0.55。
利用计算程序,通过反复调整系统初始参数,观察组元间隔、光焦度分配是否合适,最终确定系统初始值为:
$$ \begin{gathered} f_{2}′=-1、f_{3}′=1.62、m_{2l}=-1.45、\\m_{3l}=-1.34、d_{12d}=0.55、d_{34d}=0.55 \end{gathered} $$ 表 3为连续变焦光学系统参数简易计算程序按照上述初始值计算得到的5个视场位置的变焦间隔参数分配结果。
表 3 变焦系统初始间隔参数Table 3. Initial spacing parameters of optical systemFocal length/mm 275 215 150 78.6 25.9 f1/f2 spacing/mm 83.49 80.06 74.56 59.44 18.90 f2/f3 spacing/mm 13.02 24.01 38.57 68.69 125.94 f3/f4 spacing/mm 67.23 59.67 50.61 35.61 18.90 其次,将上述程序计算的变焦光学系统各组元光焦度及间隔位置参数输入ZEMAX光学设计程序,得到系统近轴光学结构,如图 3所示。通过近轴光学结构分析,系统在各位置的焦距值与计算值吻合、总长一致、间隔布局合理、变焦曲线连续,验证程序计算结果正确,可进入下一步选型工作。
再次,考虑各组元材料与调焦镜的选取,由于系统需满足-40℃~+60℃工作环境下成像清晰要求,按照光学系统无热化设计理论,系统的光焦度分配、材料选取、元件间隔都要满足光焦度、消色差、消热差3个方程[7]。
$$ \sum\limits_{i = 1}^j {{h_i}{\varphi _i}} = \varphi $$ (18) $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\theta _i}} = 0 $$ (19) $$ {\left( {\frac{1}{{{h_1}\varphi }}} \right)^2}\sum\limits_{i = 1}^j {h_i^2{\varphi _i}{\chi _i}} = \sum\limits_{i = 1}^j {{\alpha _i}{L_i}} $$ (20) 式中:hi、ϕi、θi、χi分别为系统各透镜组近轴光线高度、光焦度、色差系数及热差系数;h1为第一透镜近轴光线高;ϕ为系统总光焦度;αi为各透镜间隔镜筒材料的线膨胀系数;Li为各间隔镜筒长度。
该系统采用分步设计技术实现系统连续变焦无热化。首先选取满足上述(18)(19)公式的光焦度和材料分配,实现系统常温状态连续变焦成像清晰及高低温情况下离焦量的线性变化,系统光学材料主要选用硅单晶、锗单晶及硒化锌。其次利用主动补偿技术使系统满足(20)式要求,由于系统透镜数量少,变倍组、补偿组采用凸轮轨道变焦,中继组的轴向移动会产生变倍效果即系统视场焦距发生变化,因此采用前固定组轴向移动来进行主动调焦消热。通过上述无热化设计方法,光学系统在-40℃~+60℃温度范围内保持其性能基本不变。
最后,设置多重结构,合理设计像差评价函数,不断调整优化。为提升连续变焦系统各视场成像清晰度要求,通过设置多个高次非球面和二元衍射面,以提供更多的自由度,有利于球差、色差、像散等各类像差的校正。
3.3 设计结果
轻小型中波红外连续变焦光学系统最终设计结果如图 4所示,从上到下依次为长焦275 mm、中焦100 mm、短焦27 mm的系统图。
系统前固定组采用硅单晶材料用于会聚目标景物光线、压缩变倍组透镜尺寸;变倍组采用锗单晶材料,利用锗单晶高折射率、高色散的特性实现大倍率的变焦;补偿组采用硒化锌材料主要利用其较低的温度折射率系数使其在高低温工作环境只产生较小的离焦量以便实时补偿;中继组选择具有低的温度折射率系数、低色散的硅单晶材料平衡前端变焦核的残留色差。
整个光学系统由4片透镜、两个平面反射镜组成,其中最大透镜为第一透镜其加工直径为71 mm,平面反射镜U型折转后光学系统轴向尺寸长度172 mm,横向尺寸宽度108 mm,光学零件总重量为64 g。该系统光学透镜数量少、重量轻,光路紧凑体积小,冷屏效率100%,适配小型化制冷探测器符合连续变焦光学系统轻小型设计理念。
3.4 二元衍射面分析
由于二元衍射面在消热差、消色差方面的优异特性,系统采用了两个二元衍射面用于减少透镜数量、简化系统设计、提高系统成像质量。
在补偿组透镜硒化锌材料上引入的二元衍射面参数为Norm Radius=15 mm,A1=-33.416,A2=4.580。经计算得到二元衍射面环带深度随透镜径向的变化如图 5所示。硒化锌二元面环带数为4,最大环带深度2.918 μm,最小环带间隔宽度为2.06 mm。该面型环带间隔宽、加工环带数量少,易于单点金刚石车削加工。
在中继组硅单晶材料引入的二元衍射面参数为Norm Radius=16 mm,A1=-60.402,A2=-1.254。其环带深度随透镜径向的变化如图 6所示。硅透镜二元面环带数为9,最大环带深度1.72 μm,最小环带间隔宽度为0.86 mm。该透镜由于材料硬、环带多相对加工难度大,目前昆明物理所光学中心采用单点金刚石车削加工工艺,能制造出满足指标要求的硅基底二元光学元件。
3.5 凸轮曲线计算
本文应用动态光学理论[8],根据像移补偿公式计算补偿组运动曲线。由于变焦组和补偿组均为沿光轴的一维移动,稳像方程为:
$$ {\beta _{2m}}{\beta _2}(1 - {\beta _{1m}}{\beta _1}){q_1} + (1 - {\beta _{2m}}{\beta _2}){q_2} = 0 $$ (21) 式中:β1为变倍组初始位置的垂轴放大率;β1m为变倍组运动后的垂轴放大率;β2为补偿组初始位置的垂轴放大率;β2m补偿组运动后的垂轴放大率;q1为变倍组沿光轴位移量;q2为补偿组沿光轴位移量。
式中:
$$ {\beta _{1m}} = \frac{{{\beta _1}{f_1}'}}{{{f_1}' - {\beta _1}{q_1}}} $$ (22) $$ {\beta _{2m}} = \frac{{{\beta _2}{f_2}'}}{{{f_2}' + (1 - {\beta _1}{\beta _{1m}}){\beta _2}{q_2} - {\beta _2}{q_2}}} $$ (23) 由(21)式~(23)式可得出q1和q2的运动关系,即:
$$ Aq_{2}^{2}+Bq_{2}+C=0$$ (24) 式中:
$$ A=(f_{1}′-β_{1}q_{1})β_{2}\text{;} $$ $$ \begin{align} B=&β_{1}β_{2}q_{1}^{2}+[f_{2}′(1-β_{2}^{2})β_{1}-f_{1}′(1-β_{1}^{2})β_{2}]q_{1}- \\&f_{1}′f_{2}′(1-β_{1}^{2}) \end{align}$$ $$ C=β_{2}^{2}f_{2}′[β_{1}q_{1}-f_{1}′(1-β_{1}^{2})]q_{1} $$ 得到:
$$ {q_2} = \frac{{ - B \pm \sqrt {{B^2} - 4AC} }}{{2A}} $$ (25) 根据上述求解公式计算该变焦系统凸轮曲线如图 7所示。变倍组最大行程为28.6 mm、补偿组最大行程35 mm;补偿组曲线变化平滑,有利于凸轮轨道加工。
4. 系统像质评价
4.1 光学传递函数
系统3个焦距状态下的光学传递函数(MTF)如图 8所示。在3个焦距状态下的系统MTF满足使用要求,光学系统成像质量清晰。
4.2 点列图
系统3个焦距状态下的点列图如图 9所示。各视场RMS(root mean square)均小于一个像素,最大弥散斑RMS半径为12.6 μm,小于像元尺寸;最大弥散斑几何半径为24.5 μm,与系统艾里斑半径20.9 μm相当。系统成像质量良好,满足使用要求。
4.3 畸变
系统畸变情况如图 10所示,在长焦端小视场位置时,最大畸变量为1.2%,在短焦端大视场位置时的最大畸变量为2.4%,该变焦系统畸变对连续成像无明显影响。
4.4 系统高低温成像质量
在高低温工作环境中,系统采用轴向移动前固定组进行主动调焦消热。系统长焦端受环境温度变化,成像质量影响较大,本文主要分析长焦275 mm在高低温下经补偿后的系统成像质量。图 11为系统在高低温下长焦经补偿后的系统调制传递函数。图 12为系统在高低温下长焦经补偿后的系统点列图。从高低温传函图及点列图中看出系统在-40℃~60℃范围内成像质量良好,满足使用要求。
5. 结论
基于小型化制冷中波640×512、像元间距15 μm的焦平面探测器,设计一款具有SWaP-C特征的正组补偿连续变焦光学系统。系统由4片透镜两片平面反射镜组成,F#为4、视场变化范围为20°×16°~2.0°×1.6°,变倍比为10×、U型折叠后系统包络尺寸为172 mm×108 mm、最大物镜口径71 mm、光学零件总重量64 g、零件加工工艺成熟,变焦凸轮曲线平滑,在-40℃~60℃范围内保持较好的成像质量。该轻小型中波红外连续变焦光学系统在导航、搜索、跟踪、警戒、侦察等领域具有广阔的市场前景。
-
表 1 几种涂层材料可见光吸收率及红发射率[17]
Table 1 Solar absorptance and IR emissivity of several coating materials
Numbers Paints αS εIR αS/εIR No.1 White paint 0.19 0.94 0.2 No.2 Polyester resin 0.17 0.5 0.34 No.3 Aluminum silicon paint 0.25 0.28 0.89 No.4 Grey paint 0.87 0.87 1 No.5 Aluminum paint 0.54 0.45 1.2 No.6 Graphite paint 0.782 0.49 1.6 No.7 Aluminum foil 0.192 0.036 5.33 -
[1] 李士刚, 张力争. 弹道导弹突防措施分析[J]. 指挥控制与仿真, 2014, 36(6): 73-76. https://www.cnki.com.cn/Article/CJFDTOTAL-QBZH201406017.htm LI Shigang, ZHANG Lizheng. Analysis on penetration strategies for ballistic missile[J]. Command Control & Simulation, 2014, 36(6): 73-76. https://www.cnki.com.cn/Article/CJFDTOTAL-QBZH201406017.htm
[2] Ender T, Leurck R F, Weaver B, et al. Systems-of-systems analysis of ballistic missile defense architecture effectiveness through surrogate modeling and simulation[J]. IEEE Systems Journal, 2010, 4(2): 156-166. DOI: 10.1109/JSYST.2010.2045541
[3] 丰松江, 聂万胜. 中段弹道导弹和背景红外辐射研究概况[J]. 红外, 2006, 27(11): 19-22. https://www.cnki.com.cn/Article/CJFDTOTAL-HWAI200611007.htm FENG Songjiang, NIE Wansheng. Study on infrared radiation of ballistic missile targets and their background in midcourse[J]. Infrared, 2006, 27(11): 19-22. https://www.cnki.com.cn/Article/CJFDTOTAL-HWAI200611007.htm
[4] 申文涛, 朱定强, 蔡国飙. 中段弹道目标的温度场与红外辐射特性计算[J]. 宇航学报, 2010, 31(9): 2210-2217. DOI: 10.3873/j.issn.1000-1328.2010.09.023 SHEN Wentao, ZHU Dingqiang, CAI Guobiao. Calculation of temperature field and infrared radiation characteristics of midcourse ballistic target[J]. Journal of Astronautics, 2010, 31(9): 2210-2217. DOI: 10.3873/j.issn.1000-1328.2010.09.023
[5] 罗晨星, 盛文. 弹头壁面温度场建模及红外辐射特性研究[J]. 空军预警学院学报, 2014, 28(4): 254-257. DOI: 10.3969/j.issn.2095-5839.2014.04.006 LUO Chenxing, SHENG Wen. Modeling of missile wall temperature field and its infrared radiation characteristics[J]. Journal of Air Force Early Warning Academy, 2014, 28(4): 254-257. DOI: 10.3969/j.issn.2095-5839.2014.04.006
[6] 杨星, 吴晓迪, 陈杰. 包络球与空间气球间温度分布及红外辐射特性的差异[J]. 红外技术, 2018(4): 395-399. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201804015.htm YANG Xing, WU Xiaodi, CHEN Jie. Differences of temperature distribution and infrared radiation feature between envelope balloon and spatial balloon[J]. Infrared Technology, 2018(4): 395-399. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201804015.htm
[7] 孙成明, 袁艳, 赵飞. 空间目标天基成像探测信噪比分析[J]. 红外与激光工程, 2015, 44(5): 1654-1659. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201505047.htm SUN Chengming, YUAN Yan, ZHAO Fei. Analysis of SNR for space-based imaging detection of space object[J]. Infrared and Laser Engineering, 2015, 44(5): 1654-1659. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201505047.htm
[8] 王盈, 黄建明, 魏祥泉. 空间目标在轨红外成像仿真[J]. 红外与激光工程, 2015, 44(9): 2593-2597. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201509010.htm WANG Ying, HUANG Jianming, WEI Xiangquan. Infrared imaging simulation of space target in orbit[J]. Infrared and Laser Engineering, 2015, 44(9): 2593-2597. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201509010.htm
[9] 黄浩, 陶华敏, 陈尚锋. 基于混合融合策略的双波段红外小目标检测方法[J]. 红外与激光工程, 2014, 43(9): 2827-2831. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201409008.htm HUANG Hao, TAO Huamin, CHEN Shangfeng. Dual-band infrared dim target detection based on hybrid fusion algorithm[J]. Infrared and Laser Engineering, 2014, 43(9): 2827-2831. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201409008.htm
[10] Adomeit U, Ebert R. Improved target detection by IR dual-band image fusion[C]//Proc. Of SPIE, 2009, 7481: 74810B.
[11] 杨桄, 童涛, 陆松岩, 等. 基于多特征的红外与可见光图像融合[J]. 光学精密工程, 2014, 22(2): 489-496. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201402034.htm YANG Guang, TONG Tao, LU Songyan, et al. Fusion of infrared and visible images based on multi-features[J]. Optics and Precision Engineering, 2014, 22(2): 489-496. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201402034.htm
[12] 王宇庆, 王索建. 红外与可见光融合图像的质量评价[J]. 中国光学, 2014, 7(3): 396-401. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201403007.htm WANG Yuqing, WANG Suojian. Quality assessment method of IR and visible fusion image[J]. Chinese Journal of Optics, 2014, 7(3): 396-401. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201403007.htm
[13] 王霄, 高思莉, 李范鸣. 深空双波段红外动态场景仿真与目标分析[J]. 激光与红外, 2017, 47(9): 1123-1127. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201709013.htm WANG Xiao, GAO Sili, LI Fanming. Simulation and target analysis of dual-band infrared dynamic scene in deep space[J]. Infrared and Laser Engineering, 2017, 47(9): 1123-1127. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201709013.htm
[14] 李波. 红外隐身技术的应用及发展趋势[J]. 中国光学, 2013, 6(6): 818-823. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201306004.htm LI Bo. Application and development trend of infrared stealth technology[J]. Chinese Journal of Optics, 2013, 6(6): 818-823. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201306004.htm
[15] 周啟航, 张刘, 霍明英, 等. 弹道导弹中段突防弹道设计与验证[J]. 光学精密工程, 2015, 23(9): 2645-2655. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201509028.htm ZHOU Qihang, ZHANG Liu, HUO Mingying, et al. Design and validation of ballistic missile midcourse penetration trajectory[J]. Optics and Precision Engineering, 2015, 23(9): 2645-2655. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201509028.htm
[16] 张骏, 杨华, 凌永顺, 等. 弹道导弹中段弹头表面温度场分布理论分析[J]. 红外与激光工程, 2005, 34(5): 582-586. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200505017.htm ZHANG Jun, YANG Hua, LING Yongshun, et al. Theoretical analysis of temperature field on the surface of ballistic missile warhead in midcourse[J]. Infrared and Laser Engineering, 2005, 34(5): 582-586. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ200505017.htm
[17] 刘石泉. 弹道导弹突防技术导论[M]. 北京: 中国宇航出版社, 2003: 179-184.