CN110686564A - 一种红外半捷联导引头制导方法及系统 - Google Patents

一种红外半捷联导引头制导方法及系统 Download PDF

Info

Publication number
CN110686564A
CN110686564A CN201910977292.7A CN201910977292A CN110686564A CN 110686564 A CN110686564 A CN 110686564A CN 201910977292 A CN201910977292 A CN 201910977292A CN 110686564 A CN110686564 A CN 110686564A
Authority
CN
China
Prior art keywords
representing
bullet
missile
information
target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910977292.7A
Other languages
English (en)
Other versions
CN110686564B (zh
Inventor
任章
梁源
李清东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Beijing University of Aeronautics and Astronautics
Original Assignee
Beijing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Aeronautics and Astronautics filed Critical Beijing University of Aeronautics and Astronautics
Priority to CN201910977292.7A priority Critical patent/CN110686564B/zh
Publication of CN110686564A publication Critical patent/CN110686564A/zh
Application granted granted Critical
Publication of CN110686564B publication Critical patent/CN110686564B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41—WEAPONS
    • F41G—WEAPON SIGHTS; AIMING
    • F41G3/00—Aiming or laying means
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明公开了一种红外半捷联导引头制导方法及系统。所述方法包括：获取攻击目标的失调角信息；将失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值；依据弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令；依据法向过载控制指令，控制导弹对攻击目标进行攻击。本发明能够提高导弹的制导精度。

Description

一种红外半捷联导引头制导方法及系统
技术领域
本发明涉及导弹制导控制技术领域，特别是涉及一种红外半捷联导引头制导方法及系统。
背景技术
随着导弹技术的发展，精确制导武器早已成为现代陆、海、空战中的常规武器，其凭借着命中精度高、可实施远程精确打击等特点，成为现代战争中的主角。早在20世纪70年代中期，美国在越南战争中使用了大量的精确制导导弹，由于它具备精确的制导装置，因此在战争中取得了惊人的作战效果，精确制导武器也逐渐引起人们的注意，据资料统计，在北约对南联盟的空袭中，精确制导武器的使用率高达98％；更有西方的专家认为：精确制导武器是一种能够代替核武器，并且对战争的胜负起到决定性作用的新型武器。
导引头作为导弹的“眼睛”，用于识别、锁定、跟踪敌方目标，是决定导弹杀伤性能的核心部件。传统的导引头无法很好地对多目标进行识别，全向视场能力差，抗干扰能力弱。目前，在导引头领域，国内外的研究热点为红外半捷联导引头，这种导引头取消了传统的万向支架，利用具有瞬时大视场的红外导引头直接结合弹体姿态的调整完成对目标的捕获和追踪。美国的陆基中段导弹防御系统GMD、THAAD导弹都采用了红外半捷联导引头。虽然红外半捷联导引头具有小型化、低成本的优势，但是也带来了滤波方程可观测性弱的问题，利用其对目标进行被动跟踪的滤波系统不仅存在高度的非线性，而且其可观性的强弱与导弹的飞行弹道有关，如何提升滤波方程的可观测性进而提升导弹打击精度成为了目前亟待解决的问题。
发明内容
基于此，有必要提供一种红外半捷联导引头制导方法及系统，以解决红外半捷联导引头存在的滤波方程可观测性弱的问题，进而提高导弹的制导精度。
为实现上述目的，本发明提供了如下方案：
一种红外半捷联导引头制导方法，包括：
获取攻击目标的失调角信息；
将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值；
依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令；
依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
可选的，所述将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息，具体包括：
构建红外半捷联导引头目标跟踪系统模型；
将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
可选的，所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程；
所述状态方程为
Figure BDA0002234051960000021
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure BDA0002234051960000031
表示r的导数，
Figure BDA0002234051960000032
表示vr的导数，
Figure BDA0002234051960000033
表示q的导数，
Figure BDA0002234051960000034
表示ω的导数，
Figure BDA0002234051960000035
表示ε的导数，
Figure BDA0002234051960000036
表示atx的导数，
Figure BDA0002234051960000037
表示aty的导数；
所述测量方程为Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
可选的，所述法向过载控制指令，具体为：
ωr＝f1(t)sin(bt)；
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，
Figure BDA0002234051960000039
表示弹目相对距离估计值，
Figure BDA00022340519600000310
表示弹目接近速度估计值，
Figure BDA00022340519600000311
表示弹目视线角估计值，
Figure BDA00022340519600000312
表示目标在X轴上的加速度投影估计值，
Figure BDA00022340519600000313
表示目标在Y轴上的加速度投影估计值，表示弹目视线角估计值的导数，f1(t)表示基准函数，
Figure BDA00022340519600000315
t表示当前制导系统的时间，a、b、c表示参数，
Figure BDA00022340519600000316
ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
本发明还提供了一种红外半捷联导引头制导系统，包括：
信息获取模块，用于获取攻击目标的失调角信息；
滤波模块，用于将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值；
过载控制指令产生模块，用于依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令；
攻击模块，用于依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
可选的，所述滤波模块，具体包括：
模型构建单元，用于构建红外半捷联导引头目标跟踪系统模型；
估计单元，用于将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
可选的，所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程；
所述状态方程为
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure BDA0002234051960000042
表示r的导数，
Figure BDA0002234051960000043
表示vr的导数，
Figure BDA0002234051960000044
表示q的导数，
Figure BDA0002234051960000045
表示ω的导数，
Figure BDA0002234051960000046
表示ε的导数，
Figure BDA0002234051960000047
表示atx的导数，
Figure BDA0002234051960000048
表示aty的导数；
所述测量方程为Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
可选的，所述法向过载控制指令，具体为：
Figure BDA0002234051960000051
ωr＝f1(t)sin(bt)；
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，
Figure BDA0002234051960000052
表示弹目相对距离估计值，
Figure BDA0002234051960000053
表示弹目接近速度估计值，
Figure BDA0002234051960000054
表示弹目视线角估计值，
Figure BDA0002234051960000055
表示目标在X轴上的加速度投影估计值，表示目标在Y轴上的加速度投影估计值，
Figure BDA0002234051960000057
表示弹目视线角估计值的导数，f1(t)表示基准函数，
Figure BDA0002234051960000058
t表示当前制导系统的时间，a、b、c表示参数，
Figure BDA0002234051960000059
ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
与现有技术相比，本发明的有益效果是：
本发明提出了一种红外半捷联导引头制导方法及系统，其中，采用弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数设计制导律，从而得到导弹的法向过载控制指令，该方法提高了弹目相对运动的机动性，使得导弹在接近目标的过程中能够同时保持弹目视线角处于频繁变化的状态，从而提升红外半捷联导引头目标跟踪系统的目标跟踪精度，进而提升导弹的制导精度，实现导弹对攻击目标的准确攻击。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例一种红外半捷联导引头制导方法的流程图；
图2为本发明实施例参考视线角速率随时间的变化曲线图；
图3为本发明实施例红外半捷联导引头制导方法所得到的导弹弹道与传统比例导引法所得到的弹道的比较结果图；
图4为本发明实施例红外半捷联导引头制导方法所得到的弹目接近速度估计误差与传统比例导引法所得到的弹目接近速度估计误差的比较结果图；
图5为本发明实施例一种红外半捷联导引头制导系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本实施例的红外半捷联导引头制导方法，通过优化导弹飞行弹道，提升红外半捷联导引头目标跟踪系统的可观测性，进而有效提升制导信息精度，最终实现提高导弹制导精度的目的。
本实施例的红外半捷联导引头制导方法的核心思想如下：由于导弹的制导精度很大程度上取决于红外半捷联导引头目标跟踪系统所提供的目标信息的准确程度，所以从某种意义上来说，想要提高制导精度，就需要确保红外半捷联导引头目标跟踪系统所提供的目标信息足够准确；而红外半捷联导引头目标跟踪系统所提供的目标信息的准确程度很大程度上正比于系统本身的可观测性，可观测性强弱又主要由导弹的飞行弹道决定。由此可知，通过优化导弹飞行弹道，进而提高红外半捷联导引头目标跟踪系统的可观测性的方式可以有效的提升制导精度。而通过可观测性相关分析可知，在制导过程中通过增大视线角速率，提高弹目相对运动，提升视线角变化速率，能够有效改善并提高红外半捷联导引头目标跟踪系统的可观测性。
图1为本发明实施例一种红外半捷联导引头制导方法的流程图。参见图1，本实施例的红外半捷联导引头制导方法具体如下：
实施例的红外半捷联导引头制导方法，包括：
步骤S1：获取攻击目标的失调角信息。
具体的，导弹发射后，由红外半捷联导引头目标跟踪系统对攻击目标进行跟踪，获取失调角信息。
步骤S2：将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息。
所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值。
所述步骤S2，具体包括：
21)构建红外半捷联导引头目标跟踪系统模型。所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程。
所述状态方程为
Figure BDA0002234051960000071
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure BDA0002234051960000072
表示r的导数，表示vr的导数，
Figure BDA0002234051960000074
表示q的导数，
Figure BDA0002234051960000075
表示ω的导数，
Figure BDA0002234051960000076
表示ε的导数，
Figure BDA0002234051960000077
表示atx的导数，
Figure BDA0002234051960000078
表示aty的导数。
由于红外半捷联导引头是被动导引头，无法直接测得弹目相对距离，只能获得失调角信息，因此测量方程为Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
在红外半捷联导引头目标跟踪系统模型中，ωa,amx,amy为事先准确已知量，wtx,wty视作过程噪声处理，ε为滤波输入，r,vr,q,atx,aty为滤波输出(待估计变量)。
22)将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波(ExtendedKalmanFilter，EKF)算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
经过上述步骤21)和22)，可得到r,vr,q,atx,aty的估计值，依次为：
Figure BDA0002234051960000081
步骤S3：依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令。
所述法向过载控制指令，具体为：
Figure BDA0002234051960000082
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，表示弹目相对距离估计值，
Figure BDA0002234051960000084
表示弹目接近速度估计值，
Figure BDA0002234051960000085
表示弹目视线角估计值，
Figure BDA0002234051960000086
表示目标在X轴上的加速度投影估计值，
Figure BDA0002234051960000087
表示目标在Y轴上的加速度投影估计值，表示弹目视线角估计值的导数。
步骤S4：依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
作为一种可选的实施方式，预设比例系数F和参考视线角速率ωr的确定过程为：根据上述步骤S3中的法向过载控制指令的形式，如果在制导过程中弹目视线角和参考视线角速率都呈现出振荡的状态，则能增强导弹的机动性能，即增强红外半捷联目标跟踪系统的可观测性。但是在参数的选择上，若F选得过小，则振荡项带来的影响过低，振幅较小，可观测性不会有明显提高；若过大，则振荡项
Figure BDA00022340519600000810
带来的影响过大，振幅较大，可观测性明显提高，但是，这样导致的导弹的可用法向过载小于导弹加速度最大值，视线角速率发散会导致弹目拦截的脱靶量增大，因此会降低打击精度。
因此，结合滤波模型可观测性和脱靶量精度的综合考虑，本实施方式中采用下面的原则来设计参考视线角速率：在弹目拦截进入末制导的时候，初期可以设计不断振荡的参考视线角速率，这样可以增大弹目机动，从而提高滤波模型的可观测性；但是在末期，将参考视线角速率设计为趋近于0。这样能抑制脱靶量，则整个末制导阶段，参考视线角速率曲线是一边振荡一边收敛的。综合实际工程经验，F选择为0.005；参考视线角速率ωr为：
ωr＝f1(t)sin(bt)；
其中，f1(t)表示基准函数，
Figure BDA0002234051960000091
t表示当前制导系统的时间，制导系统的时间的起点为制导开始时刻，也可以根据需要，将其调整为当前制导系统的时间的倍数，a、b、c表示参数，
Figure BDA0002234051960000092
Figure BDA0002234051960000093
ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
各参数的意义如下：a的大小正比于振幅的大小，b的大小正比于震荡的频率，c的大小正比于收敛的速度。在假设弹目之间的初始接近速度为Vr＝1200m/s，末制导阶段初始弹目相对距离为5400m，导航比N＝5，拦截能力ad＝40g的前提下，ωr随时间变化的曲线图如图2所示。从图2中可以看出，ωr的选择满足一边振荡一边收敛的要求。
在设定好ωr与F的基础上，可以得到导弹的法向过载控制指令，进而实现导弹对目标的准确攻击。
下面对上述实施例中的红外半捷联导引头制导方法进行验证，具体的，分别应用上述红外半捷联导引头制导方法与传统的比例导引法得到导弹弹道，二者所得到的弹道的比较如图3所示。其中，PNG代表传统的比例导引法的结果，APNG代表本实施例的红外半捷联导引头制导方法的结果。从图3中可以看出，本实施例的红外半捷联导引头制导方法，可以有效的提升弹目之间相对运动的机动性。采用比例导引法和本实施例的红外半捷联导引头制导方法得到的弹目接近速度vr的估计误差比较图如图4所示，其中，PNG代表传统的比例导引法的结果，APNG代表本实施例的红外半捷联导引头制导方法的结果。从图4中可以看出，本实施例的红外半捷联导引头制导方法，可以有效提升弹目接近速度vr的估计精度。
本实施例的红外半捷联导引头制导方法具有以下优点：
1)采用弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数设计制导律，从而得到导弹的法向过载控制指令，能够提升弹目相对运动的机动性，增强红外半捷联导引头目标跟踪系统可观测性，进而提升弹目相对位姿的估计精度，最终达到提升导弹制导精度的目的。
2)通过采用提升参考视线角速率的方式，提升红外半捷联导引头目标跟踪系统可观测性，此外，根据末制导不同阶段的不同特点，选择不同参考视线角速率：在末制导初期，选用不断振荡的参考视线角速率，从而增强红外半捷联导引头目标跟踪系统的可观测性；在末制导末期，选用收敛趋近于0的参考视线角速率，达到抑制脱靶量，提升制导精度的目的；通过上述参考角速率的选择方式，使得制导方法在增强目标跟踪系统的可观测增强的同时，提高制导精度。
本发明还提供了一种红外半捷联导引头制导系统，图5为本发明实施例一种红外半捷联导引头制导系统的结构示意图。
参见图5，实施例的红外半捷联导引头制导系统，包括：
信息获取模块501，用于获取攻击目标的失调角信息。
滤波模块502，用于将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值。
过载控制指令产生模块503，用于依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令。
攻击模块504，用于依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
作为一种可选的实施方式，所述滤波模块502，具体包括：
模型构建单元，用于构建红外半捷联导引头目标跟踪系统模型。
估计单元，用于将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
作为一种可选的实施方式，所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程。
所述状态方程为
Figure BDA0002234051960000111
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure BDA0002234051960000112
表示r的导数，
Figure BDA0002234051960000113
表示vr的导数，
Figure BDA0002234051960000114
表示q的导数，
Figure BDA0002234051960000115
表示ω的导数，
Figure BDA0002234051960000116
表示ε的导数，
Figure BDA0002234051960000117
表示atx的导数，
Figure BDA0002234051960000118
表示aty的导数。
所述测量方程为Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
作为一种可选的实施方式，所述法向过载控制指令，具体为：
Figure BDA0002234051960000119
ωr＝f1(t)sin(bt)；
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，
Figure BDA0002234051960000121
表示弹目相对距离估计值，
Figure BDA0002234051960000122
表示弹目接近速度估计值，
Figure BDA0002234051960000123
表示弹目视线角估计值，
Figure BDA0002234051960000124
表示目标在X轴上的加速度投影估计值，
Figure BDA0002234051960000125
表示目标在Y轴上的加速度投影估计值，
Figure BDA0002234051960000126
表示弹目视线角估计值的导数，f1(t)表示基准函数，t表示当前制导系统的时间，a、b、c表示参数，ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
本实施例的红外半捷联导引头制导系统，设置有过载控制指令产生模块，提高了弹目相对运动的机动性，使得导弹在接近目标的过程中能够同时保持弹目视线角处于频繁变化的状态，从而提升红外半捷联导引头目标跟踪系统的目标跟踪精度，进而提升导弹的制导精度，实现导弹对攻击目标的准确攻击。
对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种红外半捷联导引头制导方法，其特征在于，包括：
获取攻击目标的失调角信息；
将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值；
依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令；
依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
2.根据权利要求1所述的一种红外半捷联导引头制导方法，其特征在于，所述将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息，具体包括：
构建红外半捷联导引头目标跟踪系统模型；
将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
3.根据权利要求1所述的一种红外半捷联导引头制导方法，其特征在于，所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程；
所述状态方程为
Figure FDA0002234051950000011
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure FDA0002234051950000021
表示r的导数，
Figure FDA0002234051950000022
表示vr的导数，
Figure FDA0002234051950000023
表示q的导数，
Figure FDA0002234051950000024
表示ω的导数，
Figure FDA0002234051950000025
表示ε的导数，表示atx的导数，
Figure FDA0002234051950000027
表示aty的导数；
所述测量方程为
Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
4.根据权利要求1所述的一种红外半捷联导引头制导方法，其特征在于，所述法向过载控制指令，具体为：
Figure FDA0002234051950000028
ωr＝f1(t)sin(bt)；
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，
Figure FDA0002234051950000029
表示弹目相对距离估计值，
Figure FDA00022340519500000210
表示弹目接近速度估计值，
Figure FDA00022340519500000211
表示弹目视线角估计值，
Figure FDA00022340519500000212
表示目标在X轴上的加速度投影估计值，
Figure FDA00022340519500000213
表示目标在Y轴上的加速度投影估计值，
Figure FDA00022340519500000214
表示弹目视线角估计值的导数，f1(t)表示基准函数，
Figure FDA00022340519500000215
t表示当前制导系统的时间，a、b、c表示参数，
Figure FDA00022340519500000216
ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
5.一种红外半捷联导引头制导系统，其特征在于，包括：
信息获取模块，用于获取攻击目标的失调角信息；
滤波模块，用于将所述失调角信息输入至红外半捷联导引头目标跟踪系统模型中，采用滤波算法得到弹目相对位姿估计信息；所述弹目相对位姿估计信息包括弹目相对距离估计值、弹目接近速度估计值、弹目视线角估计值和目标加速度投影估计值；
过载控制指令产生模块，用于依据所述弹目相对位姿估计信息、导航比、参考视线角速率和预设比例系数得到导弹的法向过载控制指令；
攻击模块，用于依据所述法向过载控制指令，控制所述导弹对所述攻击目标进行攻击。
6.根据权利要求5所述的一种红外半捷联导引头制导系统，其特征在于，所述滤波模块，具体包括：
模型构建单元，用于构建红外半捷联导引头目标跟踪系统模型；
估计单元，用于将所述失调角信息输入至所述红外半捷联导引头目标跟踪系统模型中，采用扩展卡尔曼滤波算法对弹目相对位姿信息进行估计，得到弹目相对位姿估计信息；所述弹目相对位姿信息包括弹目相对距离、弹目接近速度、弹目视线角和目标加速度投影。
7.根据权利要求5所述的一种红外半捷联导引头制导系统，其特征在于，所述红外半捷联导引头目标跟踪系统模型包括状态方程和测量方程；
所述状态方程为
Figure FDA0002234051950000031
其中，r表示弹目相对距离，vr表示弹目接近速度，q表示弹目视线角，ω表示弹目视线角速率，ε表示去噪后的失调角信息，atx表示目标在X轴上的加速度投影，aty表示目标在Y轴上的加速度投影，amx表示导弹在X轴上的加速度投影，amy表示导弹在Y轴上的加速度投影，ωa表示红外导引头扫描角速率，α表示目标运动的相关时间，wtx表示目标在X轴上运动的等效急动度，wty表示目标在Y轴上运动的等效急动度，
Figure FDA0002234051950000032
表示r的导数，
Figure FDA0002234051950000033
表示vr的导数，
Figure FDA0002234051950000034
表示q的导数，
Figure FDA0002234051950000035
表示ω的导数，
Figure FDA0002234051950000036
表示ε的导数，
Figure FDA0002234051950000037
表示atx的导数，
Figure FDA0002234051950000041
表示aty的导数；
所述测量方程为
Z＝ε+V；
其中，Z表示失调角信息，V表示均值为零的高斯白噪声。
8.根据权利要求5所述的一种红外半捷联导引头制导系统，其特征在于，所述法向过载控制指令，具体为：
ωr＝f1(t)sin(bt)；
其中，am表示法向过载控制指令，N表示导航比，F表示预设比例系数，ωr表示参考视线角速率，
Figure FDA0002234051950000043
表示弹目相对距离估计值，
Figure FDA0002234051950000044
表示弹目接近速度估计值，
Figure FDA0002234051950000045
表示弹目视线角估计值，
Figure FDA0002234051950000046
表示目标在X轴上的加速度投影估计值，
Figure FDA0002234051950000047
表示目标在Y轴上的加速度投影估计值，表示弹目视线角估计值的导数，f1(t)表示基准函数，
Figure FDA0002234051950000049
t表示当前制导系统的时间，a、b、c表示参数，
Figure FDA00022340519500000410
ri表示开始制导时的弹目相对距离，vri表示开始制导时的弹目接近速度，ad表示导弹的拦截能力，logg表示求取以g为底的对数值，ln表示求取以e为底的对数值，g表示重力加速度，e表示自然常数。
CN201910977292.7A 2019-10-15 2019-10-15 一种红外半捷联导引头制导方法及系统 Active CN110686564B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910977292.7A CN110686564B (zh) 2019-10-15 2019-10-15 一种红外半捷联导引头制导方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910977292.7A CN110686564B (zh) 2019-10-15 2019-10-15 一种红外半捷联导引头制导方法及系统

Publications (2)

Publication Number Publication Date
CN110686564A true CN110686564A (zh) 2020-01-14
CN110686564B CN110686564B (zh) 2020-08-11

Family

ID=69112687

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910977292.7A Active CN110686564B (zh) 2019-10-15 2019-10-15 一种红外半捷联导引头制导方法及系统

Country Status (1)

Country Link
CN (1) CN110686564B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111707148A (zh) * 2020-06-28 2020-09-25 北京航空航天大学 一种时变时延条件下多导弹协同制导方法及系统
CN113587740A (zh) * 2021-07-22 2021-11-02 北京航天飞腾装备技术有限责任公司 一种基于弹目视线角的被动反辐射导引方法及系统
CN114061380A (zh) * 2021-09-29 2022-02-18 湖北航天飞行器研究所 一种概略瞄准制导方法
CN115342815A (zh) * 2022-08-26 2022-11-15 哈尔滨工业大学 反大气层内或临近空间机动目标视线角速率估计方法
CN116625166A (zh) * 2023-07-20 2023-08-22 西安松果电子科技有限公司 激光半主动引导头的结构解旋方法与装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4223531A1 (de) * 1992-07-17 1994-01-20 Bodenseewerk Geraetetech Einrichtung zur Stützung und Re-Initialisierung eines Trägheitsreferenzsystems in einem Flugkörper
US20100030520A1 (en) * 2008-07-31 2010-02-04 Collier Jarrell D System for Real-Time Object Detection and Interception
US8946606B1 (en) * 2008-03-26 2015-02-03 Arete Associates Determining angular rate for line-of-sight to a moving object, with a body-fixed imaging sensor
CN105446352A (zh) * 2015-11-23 2016-03-30 哈尔滨工业大学 一种比例导引制导律辨识滤波方法
CN105865272A (zh) * 2016-05-27 2016-08-17 北京航空航天大学 一种用于半捷联制导导弹的一体化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4223531A1 (de) * 1992-07-17 1994-01-20 Bodenseewerk Geraetetech Einrichtung zur Stützung und Re-Initialisierung eines Trägheitsreferenzsystems in einem Flugkörper
US8946606B1 (en) * 2008-03-26 2015-02-03 Arete Associates Determining angular rate for line-of-sight to a moving object, with a body-fixed imaging sensor
US20100030520A1 (en) * 2008-07-31 2010-02-04 Collier Jarrell D System for Real-Time Object Detection and Interception
CN105446352A (zh) * 2015-11-23 2016-03-30 哈尔滨工业大学 一种比例导引制导律辨识滤波方法
CN105865272A (zh) * 2016-05-27 2016-08-17 北京航空航天大学 一种用于半捷联制导导弹的一体化控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YINGRONG YU ET AL.: "Guidance Information Estimation of the Semi-strapdown Infrared Imaging Seeker", 《第36届中国控制会议论文集》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111707148A (zh) * 2020-06-28 2020-09-25 北京航空航天大学 一种时变时延条件下多导弹协同制导方法及系统
CN113587740A (zh) * 2021-07-22 2021-11-02 北京航天飞腾装备技术有限责任公司 一种基于弹目视线角的被动反辐射导引方法及系统
CN114061380A (zh) * 2021-09-29 2022-02-18 湖北航天飞行器研究所 一种概略瞄准制导方法
CN114061380B (zh) * 2021-09-29 2023-11-03 湖北航天飞行器研究所 一种概略瞄准制导方法
CN115342815A (zh) * 2022-08-26 2022-11-15 哈尔滨工业大学 反大气层内或临近空间机动目标视线角速率估计方法
CN115342815B (zh) * 2022-08-26 2024-04-26 哈尔滨工业大学 反大气层内或临近空间机动目标视线角速率估计方法
CN116625166A (zh) * 2023-07-20 2023-08-22 西安松果电子科技有限公司 激光半主动引导头的结构解旋方法与装置
CN116625166B (zh) * 2023-07-20 2023-10-03 西安松果电子科技有限公司 激光半主动引导头的结构解旋方法与装置

Also Published As

Publication number Publication date
CN110686564B (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
CN110686564B (zh) 一种红外半捷联导引头制导方法及系统
CN108168381B (zh) 一种多枚导弹协同作战的控制方法
CN103090728B (zh) 一种基于滑模控制的带末角约束制导方法
CN105021092B (zh) 一种捷联寻的导引头的制导信息提取方法
Oshman et al. Differential-game-based guidance law using target orientation observations
CN105486307B (zh) 针对机动目标的视线角速率估计方法
KR102140097B1 (ko) 총포 기반 대공 방어용 사격 통제 방법
CN111351401B (zh) 应用于捷联导引头制导飞行器的防侧偏制导方法
CN112069605B (zh) 一种带有攻击时间约束的比例导引律设计方法
WO2022007509A1 (zh) 一种图像制导飞行器延时补偿方法及系统
CN109827478B (zh) 一种带落角约束与过载约束的制导方法
CN106091816B (zh) 一种基于滑模变结构理论的半捷联空空导弹制导方法
CN107726921B (zh) 一种适用于目标机动条件下的主动雷达角制导方法
CN117171877A (zh) 基于时机博弈的高超声速飞行器机动突防策略设计方法
CN114153143B (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN114489101B (zh) 一种无人机的末端制导控制方法及系统
CN113639586B (zh) 一种抗雷达关机的制导方法、系统及介质
CN110471283A (zh) 一种带碰撞角约束的三维鲁棒制导律构建方法
Oshman et al. Enhanced air-to-air missile tracking using target orientation observations
CN109780933A (zh) 一种单兵制导火箭动态目标预测导引方法
CN111026137B (zh) 一种攻击角约束下同时攻击目标的三维分布式协同导引方法
CN111879178A (zh) 一种提高卫星制导弹药制导精度的方法
CN115617063A (zh) 一种带落角约束的飞行器制导控制装置及方法
Alqudsi et al. Modeling of target tracking system for homing missiles and air defense systems
CN105987652A (zh) 姿态角速率估算系统及应用其的弹药

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant