CN112033235A - 一种采用黄金分割搜索的空战三体对抗最优预测制导方法 - Google Patents

一种采用黄金分割搜索的空战三体对抗最优预测制导方法 Download PDF

Info

Publication number
CN112033235A
CN112033235A CN202010923657.0A CN202010923657A CN112033235A CN 112033235 A CN112033235 A CN 112033235A CN 202010923657 A CN202010923657 A CN 202010923657A CN 112033235 A CN112033235 A CN 112033235A
Authority
CN
China
Prior art keywords
missile
defense
incoming
coordinate system
guidance
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
CN202010923657.0A
Other languages
English (en)
Other versions
CN112033235B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN202010923657.0A priority Critical patent/CN112033235B/zh
Publication of CN112033235A publication Critical patent/CN112033235A/zh
Application granted granted Critical
Publication of CN112033235B publication Critical patent/CN112033235B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B15/00Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
    • F42B15/01Arrangements thereon for guidance or control
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Chemical & Material Sciences (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明涉及航空航天技术领域,公开了一种采用黄金分割搜索的空战三体对抗最优预测制导方法,包括:建立空战三体对抗交战三方的运动学模型,根据交战三方实时运动状态选取优化指标函数,采用黄金分割法优化求解得到最优剩余拦截时间,将来袭导弹运动状态积分至剩余拦截时刻并选取为虚拟目标点,将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。本方法基于预测思想,将飞机在空战三体对抗场景中的耦合状态考虑在制导律设计中,选取未来相对距离为优化指标,针对非线性耦合运动模型进行优化,所设计制导律有效提升了防御导弹拦截概率,降低了脱靶量,且需求制导过载小于传统制导律,能够使低性能防御导弹完成对高速高机动来袭导弹的拦截。

Description

一种采用黄金分割搜索的空战三体对抗最优预测制导方法
技术领域
本发明涉及反导技术领域,特别涉及一种采用黄金分割搜索的空战三体对抗最优预测制导方法,该方法适用于航空航天器的制导与控制。
背景技术
空战三体对抗问题由飞机、来袭导弹、防御导弹三方飞行器对象共同构成,来袭导弹瞄准飞机进行制导攻击,飞机需要发射防御导弹对来袭导弹实施制导拦截,从而完成自我防御。由于该问题既包含了防御导弹对来袭导弹的拦截制导,又包含了飞机对来袭导弹机动规避的耦合关系,因而防御制导律的设计难度很大,近年来受到了国内外的广泛关注。
在传统的制导律研究中,主要是针对一对一情况设计的制导律,而在空战三体对抗场景下,防御导弹与来袭导弹的相对速度较大,由于未能考虑飞机在制导任务中的耦合关系,制导产生的脱靶量较大、精度较低,制导需求过载较高,制导能量消耗也较大。针对此问题,预测制导方法的提出将是一种高效和最优的解决方案。
发明内容
针对传统制导方法存在的上述问题,本发明提出了一种采用黄金分割搜索的空战三体对抗最优预测制导方法,包括如下步骤:
S1:建立空战三体对抗交战三方的三维运动学模型。
首先选取为每个飞行器的状态矢量为:
Figure BDA0002667565260000011
其中,xi、yi、zi分别表示飞行器i在地面坐标系中三轴的位置分量,θi
Figure BDA0002667565260000012
分别表示弹道倾角与弹道偏角,Vi表示飞行器的速度标量,T、A、D分别代表飞机、来袭导弹与防御导弹;
交战三方飞行器的三维运动学模型为:
Figure BDA0002667565260000021
其中,ayi、azi分别为飞行器i的机动加速度在本体坐标系中y轴与z轴方向上的分量,K1、K2分别为速度衰减系数。
S2:根据交战三方实时的运动状态,选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离。
选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离:
J(tg)=||rD(t+tg)-rA(t+tg)||
其中rD(t+tg)、rA(t+tg)分别为t+tg时刻防御导弹与来袭导弹的位置矢量,其求法如下:
t+tg时刻来袭导弹的位置矢量可假设来袭导弹对飞机使用比例制导律,通过运动学模型由t时刻来袭导弹和飞机的状态积分至t+tg时刻得到,运动学模型中
ayA=-N·Vc,TA·ΩzA
azA=-N·Vc,TA·ΩyA
N表示制导比例系数,Vc,TA为防御导弹与飞机的接近速度:
Figure BDA0002667565260000022
ΩyA和ΩzA为来袭导弹针对飞机的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
Figure BDA0002667565260000023
其中
Figure BDA0002667565260000024
为来袭导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure BDA0002667565260000031
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
A]E=rTA×vTA/|rTA|2
其中rTA、vTA分别为来袭导弹与飞机的相对位置与相对速度矢量:
Figure BDA0002667565260000032
Figure BDA0002667565260000033
t+tg时刻防御导弹的位置矢量可假设防御导弹向固定目标点rA(t+tg)使用比例制导律,通过运动学模型由t时刻防御导弹的状态积分至t+tg时刻得到,运动学模型中
ayD=-N·Vc,DA·ΩzD
azD=-N·Vc,DA·ΩyD
N表示制导比例系数,Vc,DA为防御导弹与固定目标点的接近速度:
Figure BDA0002667565260000034
ΩyD和ΩzD为防御导弹针对固定目标点的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
Figure BDA0002667565260000036
其中
Figure BDA0002667565260000037
为防御导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure BDA0002667565260000035
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
Figure BDA0002667565260000041
S3:采用黄金分割法优化求解得到防御导弹相对来袭导弹的最优剩余拦截时间。
首先使用线性化方法估计得到当前t时刻防御导弹相对来袭导弹的线性剩余拦截时间
Figure BDA0002667565260000042
Figure BDA0002667565260000043
其中rDA为防御导弹与来袭导弹的相对距离,Vc,DA为防御导弹与来袭导弹的接近速度;
选定待优化的最优剩余拦截时间初始搜索区间
Figure BDA0002667565260000046
根据黄金分割率确定两个测试边界值:
Figure BDA00026675652600000415
Figure BDA0002667565260000045
Figure BDA0002667565260000047
Figure BDA0002667565260000048
分别代入优化指标函数J(tg),比较大小后确定新的搜索区间边界;若
Figure BDA0002667565260000049
则取
Figure BDA00026675652600000410
同时保留右边界tR,若
Figure BDA00026675652600000411
则取
Figure BDA00026675652600000412
同时保留左边界tL,继续进行下轮迭代搜索;直至||tR-tL||<0.01时结束搜索,所得tg *=tL即为所求最优剩余拦截时间。
S4:将来袭导弹的运动学模型以当前的运动状态为初始状态,积分至最优剩余拦截时间时刻,选取积分后的位置为虚拟目标点。
在t时刻,将来袭导弹的运动学模型由
Figure BDA00026675652600000413
积分至t+tg *时刻,得到末状态为
Figure BDA00026675652600000414
从而得到虚拟目标点位置为:
Figure BDA0002667565260000044
S5:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
在t时刻,防御导弹与虚拟目标点的相对位置为:
Figure BDA0002667565260000051
防御导弹与虚拟目标点的接近速度为:
Figure BDA0002667565260000052
防御导弹相对于虚拟目标点在地面坐标系下的视线角速度为:
Figure BDA0002667565260000053
防御导弹的预测制导需求过载为:
Figure BDA0002667565260000054
Figure BDA0002667565260000055
其中,ΩEy和ΩEz为[ΩD]E在y轴和z轴上的分量,τ为一阶惯性环节时间常数,g为重力加速度。
本发明的优点在于:
(1)使用预测的思想设计制导律,使得拦截的有效性大大提升,有效降低了制导拦截的脱靶量;
(2)预测制导律相对于传统制导律的需求过载更小,消耗能量更少,能够使得低过载低性能防御导弹完成对高速高机动来袭导弹的拦截。
附图说明
图1:采用黄金分割搜索的空战三体对抗最优预测制导方法的流程图;
图2:空战三体对抗模型飞行器运动关系示意图;
图3:采用黄金分割搜索的空战三体对抗最优预测制导方法与比例制导律的拦截弹道对比图;
图4:采用黄金分割搜索的空战三体对抗最优预测制导方法与比例制导律的制导过载对比图。
具体实施方式
下文将结合具体附图和实例详细描述本发明具体实施例。
本发明一种采用黄金分割搜索的空战三体对抗最优预测制导方法,具体步骤如下:
S1:建立空战三体对抗交战三方的三维运动学模型。
空战三体对抗模型飞行器运动关系示意图如图2所示;
首先选取为每个飞行器的状态矢量为:
Figure BDA0002667565260000062
其中,xi、yi、zi分别表示飞行器i在地面坐标系中三轴的位置分量,θi
Figure BDA0002667565260000063
分别表示弹道倾角与弹道偏角,Vi表示飞行器的速度标量,T、A、D分别代表飞机、来袭导弹与防御导弹;
交战三方飞行器的三维运动学模型为:
Figure BDA0002667565260000061
其中,ayi、azi分别为飞行器i的机动加速度在本体坐标系中y轴与z轴方向上的分量,K1、K2分别为速度衰减系数。
S2:根据交战三方实时的运动状态,选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离。
选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离:
J(tg)=||rD(t+tg)-rA(t+tg)||
其中rD(t+tg)、rA(t+tg)分别为t+tg时刻防御导弹与来袭导弹的位置矢量,其求法如下:
t+tg时刻来袭导弹的位置矢量可假设来袭导弹对飞机使用比例制导律,通过运动学模型由t时刻来袭导弹和飞机的状态积分至t+tg时刻得到,运动学模型中
ayA=-N·Vc,TA·ΩzA
azA=-N·Vc,TA·ΩyA
N表示制导比例系数,Vc,TA为防御导弹与飞机的接近速度:
Figure BDA0002667565260000071
ΩyA和ΩzA为来袭导弹针对飞机的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
Figure BDA0002667565260000076
其中
Figure BDA0002667565260000077
为来袭导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure BDA0002667565260000072
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
A]E=rTA×vTA/|rTA|2
其中rTA、vTA分别为来袭导弹与飞机的相对位置与相对速度矢量:
Figure BDA0002667565260000073
Figure BDA0002667565260000074
t+tg时刻防御导弹的位置矢量可假设防御导弹向固定目标点rA(t+tg)使用比例制导律,通过运动学模型由t时刻防御导弹的状态积分至t+tg时刻得到,运动学模型中
ayD=-N·Vc,DA·ΩzD
azD=-N·Vc,DA·ΩyD
N表示制导比例系数,Vc,DA为防御导弹与固定目标点的接近速度:
Figure BDA0002667565260000075
ΩyD和ΩzD为防御导弹针对固定目标点的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
Figure BDA0002667565260000085
其中
Figure BDA0002667565260000086
为防御导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure BDA0002667565260000081
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
Figure BDA0002667565260000082
S3:采用黄金分割法优化求解得到防御导弹相对来袭导弹的最优剩余拦截时间。
首先使用线性化方法估计得到当前t时刻防御导弹相对来袭导弹的线性剩余拦截时间
Figure BDA0002667565260000083
Figure BDA0002667565260000084
其中rDA为防御导弹与来袭导弹的相对距离,Vc,DA为防御导弹与来袭导弹的接近速度;
选定待优化的最优剩余拦截时间初始搜索区间
Figure BDA00026675652600000813
根据黄金分割率确定两个测试边界值:
Figure BDA00026675652600000814
Figure BDA00026675652600000815
Figure BDA0002667565260000087
Figure BDA0002667565260000088
分别代入优化指标函数J(tg),比较大小后确定新的搜索区间边界;若
Figure BDA0002667565260000089
则取
Figure BDA00026675652600000810
同时保留右边界tR,若J
Figure BDA00026675652600000811
则取
Figure BDA00026675652600000812
同时保留左边界tL,继续进行下轮迭代搜索;直至||tR-tL||<0.01时结束搜索,所得tg *=tL即为所求最优剩余拦截时间。
S4:将来袭导弹的运动学模型以当前的运动状态为初始状态,积分至最优剩余拦截时间时刻,选取积分后的位置为虚拟目标点。
在t时刻,将来袭导弹的运动学模型由
Figure BDA0002667565260000097
积分至t+tg *时刻,得到末状态为
Figure BDA0002667565260000098
从而得到虚拟目标点位置为:
Figure BDA0002667565260000091
S5:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
在t时刻,防御导弹与虚拟目标点的相对位置为:
Figure BDA0002667565260000092
防御导弹与虚拟目标点的接近速度为:
Figure BDA0002667565260000093
防御导弹相对于虚拟目标点在地面坐标系下的视线角速度为:
Figure BDA0002667565260000094
防御导弹的预测制导需求过载为:
Figure BDA0002667565260000095
Figure BDA0002667565260000096
其中,ΩEy和ΩEz为[ΩD]E在y轴和z轴上的分量,τ为一阶惯性环节时间常数,g为重力加速度。
实施例:
初始时刻来袭导弹从飞机迎头方向来袭,飞机、来袭导弹、防御导弹的初始状态分别为:
XT0=[0 5000 0 0 0 270]T
XA0=[25000 5000 2000 0.087 3.054 800]T
XD0=[0 5000 0 0 0 500]T
将防御导弹的制导方法使用本文所设计的空战三体对抗最优预测制导方法的仿真结果与经典的比例导引法(PN)进行了对比。使用本文制导方法与比例导引法的拦截弹道对比图如图3所示,两种方法下的防御导弹的制导过载曲线如图4所示。可以看到,采用本文最优预测制导方法在空战三体对抗场景中实施的拦截,其脱靶量远小于传统的比例制导方法,且对防御导弹的需求过载也小于比例制导。

Claims (4)

1.一种采用黄金分割搜索的空战三体对抗最优预测制导方法,其特征在于,包括如下步骤:
S1:建立空战三体对抗交战三方的三维运动学模型;
S2:根据交战三方实时的运动状态,选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离;
S3:采用黄金分割法优化求解得到防御导弹相对来袭导弹的最优剩余拦截时间;
S4:将来袭导弹的运动学模型以当前的运动状态为初始状态,积分至最优剩余拦截时间时刻,选取积分后的位置为虚拟目标点;
S5:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
2.如权利要求1所述的一种采用黄金分割搜索的空战三体对抗最优预测制导方法,其特征在于,S1中所述空战三体对抗交战三方的三维运动学模型的建立方法为:
首先选取为每个飞行器的状态矢量为:
Figure FDA0002667565250000011
其中,xi、yi、zi分别表示飞行器i在地面坐标系中三轴的位置分量,θi
Figure FDA0002667565250000012
分别表示弹道倾角与弹道偏角,Vi表示飞行器的速度标量,T、A、D分别代表飞机、来袭导弹与防御导弹;
交战三方飞行器的三维运动学模型为:
Figure FDA0002667565250000013
其中,ayi、azi分别为飞行器i的机动加速度在本体坐标系中y轴与z轴方向上的分量,K1、K2分别为速度衰减系数。
3.如权利要求1所述的一种采用黄金分割搜索的空战三体对抗最优预测制导方法,其特征在于,S2中所述优化指标函数的选取即tg时刻后防御导弹与来袭导弹之间的预测距离的计算方法为:
选取优化指标函数为tg时刻后防御导弹与来袭导弹之间的预测距离:
J(tg)=||rD(t+tg)-A(t+tg)||
其中rD(t+tg)、rA(t+tg)分别为t+tg时刻防御导弹与来袭导弹的位置矢量,其求法如下:
t+tg时刻来袭导弹的位置矢量可假设来袭导弹对飞机使用比例制导律,通过运动学模型由t时刻来袭导弹和飞机的状态积分至t+tg时刻得到,运动学模型中
ayA=-N·Vc,TA·ΩzA
azA=-N·Vc,TA·ΩyA
N表示制导比例系数,Vc,TA为防御导弹与飞机的接近速度;ΩyA和ΩzA为来袭导弹针对飞机的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
Figure FDA0002667565250000021
其中
Figure FDA0002667565250000022
为来袭导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure FDA0002667565250000023
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
A]E=rTA×vTA/|rTA|2
其中rTA、vTA分别为来袭导弹与飞机的相对位置与相对速度矢量;
t+tg时刻防御导弹的位置矢量可假设防御导弹向固定目标点rA(t+tg)使用比例制导律,通过运动学模型由t时刻防御导弹的状态积分至t+tg时刻得到,运动学模型中
ayD=-N·Vc,DA·ΩzD
azD=-N·Vc,DA·ΩyD
N表示制导比例系数,Vc,DA为防御导弹与固定目标点的接近速度;ΩyD和ΩzD为防御导弹针对固定目标点的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
Figure FDA0002667565250000031
其中
Figure FDA0002667565250000032
为防御导弹弹道坐标系至地面坐标系的方向余弦矩阵
Figure FDA0002667565250000033
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
Figure FDA0002667565250000034
4.如权利要求1所述的一种采用黄金分割搜索的空战三体对抗最优预测制导方法,其特征在于,所述S3中防御导弹相对来袭导弹的最优剩余拦截时间的优化求解方法为:
首先使用线性化方法估计得到当前t时刻防御导弹相对来袭导弹的线性剩余拦截时间
Figure FDA0002667565250000035
Figure FDA0002667565250000036
其中rDA为防御导弹与来袭导弹的相对距离,vc,DA为防御导弹与来袭导弹的接近速度;
选定待优化的最优剩余拦截时间初始搜索区间
Figure FDA0002667565250000037
根据黄金分割率确定两个测试边界值:
Figure FDA0002667565250000038
Figure FDA0002667565250000039
Figure FDA00026675652500000310
Figure FDA00026675652500000311
分别代入优化指标函数J(tg),比较大小后确定新的搜索区间边界;若
Figure FDA00026675652500000312
则取
Figure FDA00026675652500000313
同时保留右边界tR,若
Figure FDA00026675652500000314
则取
Figure FDA00026675652500000315
同时保留左边界tL,继续进行下轮迭代搜索;直至‖tR-tL‖<0.01时结束搜索,所得tL即为所求最优剩余拦截时间。
CN202010923657.0A 2020-09-04 2020-09-04 一种采用黄金分割搜索的空战三体对抗最优预测制导方法 Active CN112033235B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010923657.0A CN112033235B (zh) 2020-09-04 2020-09-04 一种采用黄金分割搜索的空战三体对抗最优预测制导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010923657.0A CN112033235B (zh) 2020-09-04 2020-09-04 一种采用黄金分割搜索的空战三体对抗最优预测制导方法

Publications (2)

Publication Number Publication Date
CN112033235A true CN112033235A (zh) 2020-12-04
CN112033235B CN112033235B (zh) 2021-07-27

Family

ID=73591516

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010923657.0A Active CN112033235B (zh) 2020-09-04 2020-09-04 一种采用黄金分割搜索的空战三体对抗最优预测制导方法

Country Status (1)

Country Link
CN (1) CN112033235B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859918A (zh) * 2021-01-18 2021-05-28 中国人民解放军火箭军工程大学 一种多智能体多模自适应协同制导律确定方法及系统
CN115468454A (zh) * 2022-08-30 2022-12-13 南京理工大学 虚拟力法与引战配合相结合的多弹协同对抗策略

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105759612A (zh) * 2016-03-22 2016-07-13 北京航空航天大学 带落角约束的微分对策反拦截机动突防/精确打击导引方法
CN106843272A (zh) * 2017-02-28 2017-06-13 北京航空航天大学 一种具有终端速度、弹道倾角和过载约束的显式制导律
CN108363299A (zh) * 2018-01-26 2018-08-03 北京航空航天大学 一种外大气层拦截最优末制导方法
CN108534614A (zh) * 2018-03-23 2018-09-14 清华大学 一种三维全向实时预测制导方法
CN109115035A (zh) * 2018-08-06 2019-01-01 哈尔滨工业大学 一种基于弹道成形的运载火箭子级返回段精确制导方法
CN110425943A (zh) * 2019-08-06 2019-11-08 西北工业大学 面向变质心飞行器的工程化再入制导方法
CN110657713A (zh) * 2019-08-12 2020-01-07 清华大学 采用预测制导的主动防御导弹最优发射时机搜索方法
CN111551080A (zh) * 2020-04-03 2020-08-18 北京航空航天大学 一种攻击时间控制的制导方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105759612A (zh) * 2016-03-22 2016-07-13 北京航空航天大学 带落角约束的微分对策反拦截机动突防/精确打击导引方法
CN106843272A (zh) * 2017-02-28 2017-06-13 北京航空航天大学 一种具有终端速度、弹道倾角和过载约束的显式制导律
CN108363299A (zh) * 2018-01-26 2018-08-03 北京航空航天大学 一种外大气层拦截最优末制导方法
CN108534614A (zh) * 2018-03-23 2018-09-14 清华大学 一种三维全向实时预测制导方法
CN109115035A (zh) * 2018-08-06 2019-01-01 哈尔滨工业大学 一种基于弹道成形的运载火箭子级返回段精确制导方法
CN110425943A (zh) * 2019-08-06 2019-11-08 西北工业大学 面向变质心飞行器的工程化再入制导方法
CN110657713A (zh) * 2019-08-12 2020-01-07 清华大学 采用预测制导的主动防御导弹最优发射时机搜索方法
CN111551080A (zh) * 2020-04-03 2020-08-18 北京航空航天大学 一种攻击时间控制的制导方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859918A (zh) * 2021-01-18 2021-05-28 中国人民解放军火箭军工程大学 一种多智能体多模自适应协同制导律确定方法及系统
CN112859918B (zh) * 2021-01-18 2023-10-27 中国人民解放军火箭军工程大学 一种多智能体多模自适应协同制导律确定方法及系统
CN115468454A (zh) * 2022-08-30 2022-12-13 南京理工大学 虚拟力法与引战配合相结合的多弹协同对抗策略

Also Published As

Publication number Publication date
CN112033235B (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
CN108534614B (zh) 一种三维全向实时预测制导方法
CN112033235B (zh) 一种采用黄金分割搜索的空战三体对抗最优预测制导方法
CN112256055B (zh) 一种采用不动点优化的三体对抗防御预测制导方法
CN112558631A (zh) 基于测量信息的带大落角约束变参数制导方法
CN114020021B (zh) 一种多导弹分布式协同制导律的设计方法及系统
CN112859921B (zh) 一种三维同时攻击鲁棒协同制导律设计方法
CN112648886B (zh) 一种组合制导目标拦截方法及系统
CN117171877A (zh) 基于时机博弈的高超声速飞行器机动突防策略设计方法
CN114035616A (zh) 一种飞行器对移动目标打击控制方法及系统
CN113911398A (zh) 飞行器单脉冲躲避策略确定方法及系统
CN116360489A (zh) 一种基于成型制导律的协同中制导律设计方法
CN114153143B (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN112729015B (zh) 一种基于拦截弹导引律辨识的机动突防方法
Awad et al. Evaluation and enhancing missile performance via real time flight simulation model
CN113051706B (zh) 一种采用虚拟点预测的三体对抗防御制导方法
Wang et al. Finite-time Simultaneous-Arrival-To-Origin Sliding Model Control for Multi-missile Systems
Zhou et al. Finite-time sliding mode guidance law design with impact angle constraints
Wang et al. Design of missile combined terminal guidance law
CN114216368B (zh) 一种变速条件下弹群协同制导方法及系统
Hanghang et al. A method to solve the unreachable zone of mid-range air-to-air missile
Wang et al. The application of proportional navigation in the process of UAV air combat guidance and optimization of proportional parameter
CN117870460A (zh) 一种基于滑模控制的进攻弹俯冲段制导方法及系统
Li et al. SDRE Optimal Sliding Mode Guidance Law Design with Attack Angle Constraint
Marković et al. Engagement areas of missiles in the proportional navigated flight powered by air breathing engines
Tahk et al. Suboptimal Guidance Based on Pursuit and Impact Angle Control for Long-Range Air-to-Air Missiles

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