CN112033235A - 一种采用黄金分割搜索的空战三体对抗最优预测制导方法 - Google Patents
一种采用黄金分割搜索的空战三体对抗最优预测制导方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting 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:建立空战三体对抗交战三方的三维运动学模型。
首先选取为每个飞行器的状态矢量为:
交战三方飞行器的三维运动学模型为:
其中,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为防御导弹与飞机的接近速度:
ΩyA和ΩzA为来袭导弹针对飞机的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
[Ω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轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
S3:采用黄金分割法优化求解得到防御导弹相对来袭导弹的最优剩余拦截时间。
其中rDA为防御导弹与来袭导弹的相对距离,Vc,DA为防御导弹与来袭导弹的接近速度;
根据黄金分割率确定两个测试边界值:
将和分别代入优化指标函数J(tg),比较大小后确定新的搜索区间边界;若则取同时保留右边界tR,若则取同时保留左边界tL,继续进行下轮迭代搜索;直至||tR-tL||<0.01时结束搜索,所得tg *=tL即为所求最优剩余拦截时间。
S4:将来袭导弹的运动学模型以当前的运动状态为初始状态,积分至最优剩余拦截时间时刻,选取积分后的位置为虚拟目标点。
S5:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
在t时刻,防御导弹与虚拟目标点的相对位置为:
防御导弹与虚拟目标点的接近速度为:
防御导弹相对于虚拟目标点在地面坐标系下的视线角速度为:
防御导弹的预测制导需求过载为:
其中,ΩEy和ΩEz为[ΩD]E在y轴和z轴上的分量,τ为一阶惯性环节时间常数,g为重力加速度。
本发明的优点在于:
(1)使用预测的思想设计制导律,使得拦截的有效性大大提升,有效降低了制导拦截的脱靶量;
(2)预测制导律相对于传统制导律的需求过载更小,消耗能量更少,能够使得低过载低性能防御导弹完成对高速高机动来袭导弹的拦截。
附图说明
图1:采用黄金分割搜索的空战三体对抗最优预测制导方法的流程图;
图2:空战三体对抗模型飞行器运动关系示意图;
图3:采用黄金分割搜索的空战三体对抗最优预测制导方法与比例制导律的拦截弹道对比图;
图4:采用黄金分割搜索的空战三体对抗最优预测制导方法与比例制导律的制导过载对比图。
具体实施方式
下文将结合具体附图和实例详细描述本发明具体实施例。
本发明一种采用黄金分割搜索的空战三体对抗最优预测制导方法,具体步骤如下:
S1:建立空战三体对抗交战三方的三维运动学模型。
空战三体对抗模型飞行器运动关系示意图如图2所示;
首先选取为每个飞行器的状态矢量为:
交战三方飞行器的三维运动学模型为:
其中,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为防御导弹与飞机的接近速度:
ΩyA和ΩzA为来袭导弹针对飞机的视线角速度在地面坐标系中y轴与z轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
[Ω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轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
S3:采用黄金分割法优化求解得到防御导弹相对来袭导弹的最优剩余拦截时间。
其中rDA为防御导弹与来袭导弹的相对距离,Vc,DA为防御导弹与来袭导弹的接近速度;
根据黄金分割率确定两个测试边界值:
将和分别代入优化指标函数J(tg),比较大小后确定新的搜索区间边界;若则取同时保留右边界tR,若J则取同时保留左边界tL,继续进行下轮迭代搜索;直至||tR-tL||<0.01时结束搜索,所得tg *=tL即为所求最优剩余拦截时间。
S4:将来袭导弹的运动学模型以当前的运动状态为初始状态,积分至最优剩余拦截时间时刻,选取积分后的位置为虚拟目标点。
S5:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
在t时刻,防御导弹与虚拟目标点的相对位置为:
防御导弹与虚拟目标点的接近速度为:
防御导弹相对于虚拟目标点在地面坐标系下的视线角速度为:
防御导弹的预测制导需求过载为:
其中,Ω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:将防御导弹以虚拟目标点为目标,采用比例导引律进行制导。
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轴方向上的分量,由弹道坐标系中来袭导弹针对飞机的视线角速度矢量转换得到:
来袭导弹针对飞机在地面坐标系中的视线角速度矢量为:
[Ω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轴方向上的分量,由弹道坐标系中防御导弹针对固定目标点的视线角速度矢量转换得到:
防御导弹针对固定目标点在地面坐标系中的视线角速度矢量为:
4.如权利要求1所述的一种采用黄金分割搜索的空战三体对抗最优预测制导方法,其特征在于,所述S3中防御导弹相对来袭导弹的最优剩余拦截时间的优化求解方法为:
其中rDA为防御导弹与来袭导弹的相对距离,vc,DA为防御导弹与来袭导弹的接近速度;
根据黄金分割率确定两个测试边界值:
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)
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)
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 | 北京航空航天大学 | 一种攻击时间控制的制导方法 |
-
2020
- 2020-09-04 CN CN202010923657.0A patent/CN112033235B/zh active Active
Patent Citations (8)
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)
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 | |
Shin | UAV Evasive Maneuver Algorithm based on Moth Natural Behavior |
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 |