CN112256055A - 一种采用不动点优化的三体对抗防御预测制导方法 - Google Patents

一种采用不动点优化的三体对抗防御预测制导方法 Download PDF

Info

Publication number
CN112256055A
CN112256055A CN202011086110.6A CN202011086110A CN112256055A CN 112256055 A CN112256055 A CN 112256055A CN 202011086110 A CN202011086110 A CN 202011086110A CN 112256055 A CN112256055 A CN 112256055A
Authority
CN
China
Prior art keywords
defense
missile
guidance
point
prediction
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
CN202011086110.6A
Other languages
English (en)
Other versions
CN112256055B (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 CN202011086110.6A priority Critical patent/CN112256055B/zh
Publication of CN112256055A publication Critical patent/CN112256055A/zh
Application granted granted Critical
Publication of CN112256055B publication Critical patent/CN112256055B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/107—Simultaneous control of position or course in three dimensions specially adapted for missiles
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41—WEAPONS
    • F41G—WEAPON SIGHTS; AIMING
    • F41G3/00—Aiming or laying means
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41—WEAPONS
    • F41H—ARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
    • F41H11/00—Defence installations; Defence devices
    • F41H11/02—Anti-aircraft or anti-guided missile or anti-torpedo defence installations or systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明涉及航空航天技术领域，公开了一种采用不动点优化的三体对抗防御预测制导方法，包括：建立三体对抗防御的运动预测模型，建立预测拦截点位置的积分方程并将预测拦截点转换至以防御导弹为原点的极坐标系，建立防御导弹向预测拦截点制导的时间函数，选取优化区间，使用不动点迭代法求解时间函数的不动点得到最优剩余拦截时间，采用运动预测模型数值积分得到预测拦截点位置，以预测位置为目标生成防御导弹的制导指令。本方法给出了防御导弹向预测拦截点制导时间的解析解表达式，从而能够快速优化在线高效生成制导指令，且能够令防御制导的能量最优，能够有效提升三体对抗防御任务的拦截概率。

Description

一种采用不动点优化的三体对抗防御预测制导方法
技术领域
本发明涉及航空航天技术领域，特别涉及一种采用不动点优化的三体对抗防御预测制导方法，该方法适用于航空航天器的制导与控制。
背景技术
当飞机受到精确制导导弹攻击的时候，传统的对抗方法通常是通过机动规避或投放诱饵等被动的方式保障自身安全，但随着电子技术各方面的进步，这些被动手段已经不足以实现有效防御。在此情境下，一种主动的防御方式应运而生，飞机可以通过主动发射防御性的拦截器的方式实现自我防御。飞机、来袭导弹，以及防御导弹三方飞行器构成了相互耦合的对抗问题，此问题也被称为三体对抗防御问题。
由于拦截高速高机动来袭导弹的难度很大，如何设计防御制导与控制策略在近年来受到了国内外的广泛关注。目前，学者们主要通过最优控制、微分对策控制，以及几何分析等手段推导制导策略。很多方法采用了线性化模型开展分析，虽能求解出最优制导律但与实际应用相差较远；而针对非线性模型所提出的方法通常计算量较大，难以实现在线的制导律运算。本发明基于非线性的三体对抗模型，旨在提出一种能够快速在线高效生成制导指令的预测制导方法。
发明内容
针对传统的三体对抗制导方法存在的上述问题，本发明提出了一种采用不动点优化的三体对抗防御预测制导方法，包括如下步骤：
S1：建立三体对抗防御的运动预测模型：
Figure BDA0002720389540000021
其中下标T和M分别表示载机和来袭导弹，(xT,yT)和(xM,yM)分别表示载机和来袭导弹在惯性坐标系中的位置；γT、γM分别表示二者的速度方向角；aT和aM分别表示二者的侧向加速度，其中载机的侧向加速度假设为已知参数，由载机飞行策略给定；假设来袭导弹采用比例导引律(PN)或扩展比例导引律(APN)来对载机进行制导，来袭导弹的侧向加速度如下：
Figure BDA0002720389540000022
其中，δ表示来袭导弹采取不同制导律时的不同取值，取0时表示使用PN，取1时表示使用APN；VMT和λMT分别表示来袭导弹与载机的接近速度与视线角速度，接近速度可表示为：
Figure BDA0002720389540000023
来袭导弹与载机的视线角速度及其导数为：
Figure BDA0002720389540000024
Figure BDA0002720389540000025
S2：根据S1所建立的运动预测模型，以从当前时刻起的预测时长tg为参数，建立预测拦截点位置的积分方程，并将预测拦截点位置转换至以防御导弹为原点的极坐标系。
以t0时刻来袭导弹的位置(xM,0,yM,0)为初始值，对运动预测模型数值积分求得tg时刻之后来袭导弹在惯性坐标系中的位置
Figure BDA0002720389540000026
即为预测拦截点；
将其转换至防御导弹的极坐标系：
Figure BDA0002720389540000027
Figure BDA0002720389540000031
其中(xD,0,yD,0)和γD,0分别为t0时刻防御导弹的位置和速度方向角。
S3：根据预测拦截点在防御导弹极坐标系中的位置，建立防御导弹向预测拦截点制导的时间函数
Figure BDA0002720389540000032
Figure BDA0002720389540000033
其中k＝1-N，N为比例制导律的比例系数。
S4：选取优化区间，使用不动点迭代法求解函数
Figure BDA0002720389540000034
的不动点，即方程F[r(tg),θ(tg)]＝tg的解，即为防御导弹的最优剩余拦截时间；
Figure BDA0002720389540000035
其中，rDM为防御导弹与来袭导弹的相对距离，rMT来袭导弹与载机的相对距离，VD和VM分别为防御导弹与来袭导弹的速度。
S5：根据S4所得到的最优剩余拦截时间，采用S2中的运动预测模型进行数值积分得到预测拦截点位置。
S6：令防御导弹使用比例导引制导律以S5所得到的预测拦截点位置为目标，生成防御导弹的制导指令：
Figure BDA0002720389540000036
其中，
Figure BDA0002720389540000037
为防御导弹的最大制导过载，Vc,DP为防御导弹与预测拦截点的接近速度：
Figure BDA0002720389540000038
λDP为防御导弹相对预测拦截点的视线角：
Figure BDA0002720389540000041
Figure BDA0002720389540000042
为防御导弹相对预测拦截点的视线角速度：
Figure BDA0002720389540000043
本发明的优点在于：
(1)基于非线性的三体对抗模型采用模型预测控制，将复杂的三体制导问题转化为实函数不动点求解问题，有效降低了制导指令求解的计算量，提高了制导效率；
(2)所提出的预测制导方法相对于传统三体对抗制导方法，消耗的制导能量更少，拦截脱靶量更小，拦截概率更高。
附图说明
图1：采用不动点优化的三体对抗防御预测制导方法的流程图；
图2：三体对抗防御的运动模型示意图；
图3：采用不动点优化的三体对抗防御预测制导方法与传统制导律的拦截弹道对比图；
图4：采用不动点优化的三体对抗防御预测制导方法与传统制导律的制导过载对比图；
图5：初始时刻防御导弹向预测拦截点制导时间函数F(tg)以及其不动点的示意图。
具体实施方式
下文将结合具体附图和实例详细描述本发明具体实施例。
本发明一种采用不动点优化的三体对抗防御预测制导方法，具体步骤如下：
S1：建立三体对抗防御的运动预测模型，如图2所示：
Figure BDA0002720389540000051
其中下标T和M分别表示载机和来袭导弹，(xT,yT)和(xM,yM)分别表示载机和来袭导弹在惯性坐标系中的位置；γT、γM分别表示二者的速度方向角；aT和aM分别表示二者的侧向加速度，其中载机的侧向加速度假设为已知参数，由载机飞行策略给定；假设来袭导弹采用比例导引律(PN)或扩展比例导引律(APN)来对载机进行制导，来袭导弹的侧向加速度如下：
Figure BDA0002720389540000052
其中，δ表示来袭导弹采取不同制导律时的不同取值，取0时表示使用PN，取1时表示使用APN；VMT和λMT分别表示来袭导弹与载机的接近速度与视线角速度，接近速度可表示为：
Figure BDA0002720389540000053
来袭导弹与载机的视线角速度及其导数为：
Figure BDA0002720389540000054
Figure BDA0002720389540000055
S2：根据S1所建立的运动预测模型，以从当前时刻起的预测时长tg为参数，建立预测拦截点位置的积分方程，并将预测拦截点位置转换至以防御导弹为原点的极坐标系。
以t0时刻来袭导弹的位置(xM,0,yM,0)为初始值，对运动预测模型数值积分求得tg时刻之后来袭导弹在惯性坐标系中的位置
Figure BDA0002720389540000056
即为预测拦截点；
将其转换至防御导弹的极坐标系：
Figure BDA0002720389540000057
Figure BDA0002720389540000061
其中(xD,0,yD,0)和γD,0分别为t0时刻防御导弹的位置和速度方向角。
S3：根据预测拦截点在防御导弹极坐标系中的位置，建立防御导弹向预测拦截点制导的时间函数
Figure BDA0002720389540000062
Figure BDA0002720389540000063
其中k＝1-N，N为比例制导律的比例系数。
S4：选取优化区间，使用不动点迭代法求解函数F(tg)的不动点，即方程F[r(tg),θ(tg)]＝tg的解，即为防御导弹的最优剩余拦截时间；
Figure BDA0002720389540000064
其中，rDM为防御导弹与来袭导弹的相对距离，rMT来袭导弹与载机的相对距离，VD和VM分别为防御导弹与来袭导弹的速度。
S5：根据S4所得到的最优剩余拦截时间，采用S2中的运动预测模型进行数值积分得到预测拦截点位置。
S6：令防御导弹使用比例导引制导律以S5所得到的预测拦截点位置为目标，生成防御导弹的制导指令：
Figure BDA0002720389540000065
其中，
Figure BDA0002720389540000066
为防御导弹的最大制导过载，Vc,DP为防御导弹与预测拦截点的接近速度：
Figure BDA0002720389540000067
λDP为防御导弹相对预测拦截点的视线角：
Figure BDA0002720389540000071
Figure BDA0002720389540000073
为防御导弹相对预测拦截点的视线角速度：
Figure BDA0002720389540000072
实施例：
在初始时刻，飞机、来袭导弹的位置分别为：
xT＝0,yT＝0，
xM＝9848.1m,yM＝1736.5m；
飞机、来袭导弹、防御导弹的速度分别为
VT＝300m/s,VM＝800m/s,VD＝500m/s；
制导律中的比例系数取值为N＝3。
令来袭导弹采用比例导引制导律(PN)向飞机实施制导，防御导弹分别采用比例导引制导律(PN)、扩展比例导引制导律(APN)、三点法制导律(CLOS)，以及本发明所提出的的预测制导方法(PG)对来袭导弹实施制导。拦截的弹道对比图如图3所示，防御导弹的制导过载曲线如图4所示。图5为初始时刻防御导弹向预测拦截点制导的时间函数F(tg)，以及其不动点的示意图。采用本发明所提出的的制导方法，其制导拦截的脱靶量、制导能量、制导需求过载均优于其它方法，且本发明提出的方法每帧制导指令生成是时间在100ms以内，能够满足当前弹载计算机需求。

Claims (6)

1.一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，包括如下步骤：
S1：建立三体对抗防御的运动预测模型；
S2：根据S1所建立的运动预测模型，以从当前时刻起的预测时长tg为参数，建立预测拦截点位置的积分方程，并将预测拦截点位置转换至以防御导弹为原点的极坐标系；
S3：根据预测拦截点在防御导弹极坐标系中的位置，建立防御导弹向预测拦截点制导的时间函数
Figure FDA0002720389530000013
S4：选取优化区间，使用不动点迭代法求解函数
Figure FDA0002720389530000012
的不动点，即方程F[r(tg),θ(tg)]＝tg的解，即为防御导弹的最优剩余拦截时间；
S5：根据S4所得到的最优剩余拦截时间，采用S2中的运动预测模型进行数值积分得到预测拦截点位置；
S6：令防御导弹使用比例导引制导律以S5所得到的预测拦截点位置为目标，生成防御导弹的制导指令。
2.如权利要求1所述的一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，S1中所述三体对抗防御的运动预测模型的建立方法为：
Figure FDA0002720389530000011
其中下标T和M分别表示载机和来袭导弹，(xT,yT)和(xM,yM)分别表示载机和来袭导弹在惯性坐标系中的位置；γT、γM分别表示二者的速度方向角；aT和aM分别表示二者的侧向加速度，其中载机的侧向加速度假设为已知参数，由载机飞行策略给定；假设来袭导弹采用比例导引律(PN)或扩展比例导引律(APN)来对载机进行制导，来袭导弹的侧向加速度如下：
Figure FDA0002720389530000021
其中，δ表示来袭导弹采取不同制导律时的不同取值，取0时表示使用PN，取1时表示使用APN；VMT和λMT分别表示来袭导弹与载机的接近速度与视线角速度，接近速度可表示为：
Figure FDA0002720389530000022
来袭导弹与载机的视线角速度及其导数为：
Figure FDA0002720389530000023
Figure FDA0002720389530000024
3.如权利要求1所述的一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，S2中所述的将预测拦截点位置转换至以防御导弹为原点的极坐标系，具体方法为：
以t0时刻来袭导弹的位置(xM,0,yM,0)为初始值，对运动预测模型数值积分求得tg时刻之后来袭导弹在惯性坐标系中的位置
Figure FDA0002720389530000025
即为预测拦截点；
将其转换至防御导弹的极坐标系：
Figure FDA0002720389530000026
Figure FDA0002720389530000027
其中(xD,0,yD,0)和γD,0分别为t0时刻防御导弹的位置和速度方向角。
4.如权利要求1所述的一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，S3中所述的防御导弹向预测拦截点制导的时间函数
Figure FDA0002720389530000028
具体方法为：
Figure FDA0002720389530000029
其中k＝1-N，N为比例制导律的比例系数。
5.如权利要求1所述的一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，S4中所述的选取优化区间，具体选取方法为：
Figure FDA0002720389530000031
其中，rDM为防御导弹与来袭导弹的相对距离，rMT来袭导弹与载机的相对距离，VD和VM分别为防御导弹与来袭导弹的速度。
6.如权利要求1所述的一种采用不动点优化的三体对抗防御预测制导方法，其特征在于，S6中所述的生成防御导弹的制导指令，具体方法为：
Figure FDA0002720389530000032
其中，
Figure FDA0002720389530000033
为防御导弹的最大制导过载，Vc,DP为防御导弹与预测拦截点的接近速度：
Figure FDA0002720389530000034
λDP为防御导弹相对预测拦截点的视线角：
Figure FDA0002720389530000035
Figure FDA0002720389530000036
为防御导弹相对预测拦截点的视线角速度：
Figure FDA0002720389530000037
CN202011086110.6A 2020-10-12 2020-10-12 一种采用不动点优化的三体对抗防御预测制导方法 Active CN112256055B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011086110.6A CN112256055B (zh) 2020-10-12 2020-10-12 一种采用不动点优化的三体对抗防御预测制导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011086110.6A CN112256055B (zh) 2020-10-12 2020-10-12 一种采用不动点优化的三体对抗防御预测制导方法

Publications (2)

Publication Number Publication Date
CN112256055A true CN112256055A (zh) 2021-01-22
CN112256055B CN112256055B (zh) 2021-11-05

Family

ID=74242764

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011086110.6A Active CN112256055B (zh) 2020-10-12 2020-10-12 一种采用不动点优化的三体对抗防御预测制导方法

Country Status (1)

Country Link
CN (1) CN112256055B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102020006465A1 (de) 2020-10-21 2022-04-21 Diehl Defence Gmbh & Co. Kg Abfangflugkörper und Verfahren zu dessen Lenkung
CN116222310A (zh) * 2023-04-13 2023-06-06 哈尔滨工业大学 三维空间中基于rbf_g的二对一同步区域覆盖拦截方法
DE102022001285A1 (de) 2022-04-13 2023-10-19 Diehl Defence Gmbh & Co. Kg Verfahren zum Lenken eines Flugkörpers
DE102022001287A1 (de) 2022-04-13 2023-10-19 Diehl Defence Gmbh & Co. Kg Verfahren zum Lenken eines Abfangflugkörpers

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107917646A (zh) * 2017-01-10 2018-04-17 北京航空航天大学 一种基于目标终端可达区预测的红外空空导弹抗干扰制导方法
CN108534614A (zh) * 2018-03-23 2018-09-14 清华大学 一种三维全向实时预测制导方法
CN109857145A (zh) * 2018-11-27 2019-06-07 北京航空航天大学 一种基于迭代预测命中点的增程型拦截弹预测制导方法
US10317852B1 (en) * 2015-10-29 2019-06-11 National Technology & Engineering Solutions Of Sandia, Llc Predictive guidance flight
CN110017729A (zh) * 2019-04-18 2019-07-16 西安交通大学 一种带碰撞角约束的多导弹时间协同制导方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10317852B1 (en) * 2015-10-29 2019-06-11 National Technology & Engineering Solutions Of Sandia, Llc Predictive guidance flight
CN107917646A (zh) * 2017-01-10 2018-04-17 北京航空航天大学 一种基于目标终端可达区预测的红外空空导弹抗干扰制导方法
CN108534614A (zh) * 2018-03-23 2018-09-14 清华大学 一种三维全向实时预测制导方法
CN109857145A (zh) * 2018-11-27 2019-06-07 北京航空航天大学 一种基于迭代预测命中点的增程型拦截弹预测制导方法
CN110017729A (zh) * 2019-04-18 2019-07-16 西安交通大学 一种带碰撞角约束的多导弹时间协同制导方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
史恒 等: ""主动防御的最优预测协同制导律研究"", 《空间控制技术与应用》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102020006465A1 (de) 2020-10-21 2022-04-21 Diehl Defence Gmbh & Co. Kg Abfangflugkörper und Verfahren zu dessen Lenkung
DE102020006465B4 (de) 2020-10-21 2022-06-30 Diehl Defence Gmbh & Co. Kg Abfangflugkörper und Verfahren zu dessen Lenkung
DE102022001285A1 (de) 2022-04-13 2023-10-19 Diehl Defence Gmbh & Co. Kg Verfahren zum Lenken eines Flugkörpers
DE102022001287A1 (de) 2022-04-13 2023-10-19 Diehl Defence Gmbh & Co. Kg Verfahren zum Lenken eines Abfangflugkörpers
CN116222310A (zh) * 2023-04-13 2023-06-06 哈尔滨工业大学 三维空间中基于rbf_g的二对一同步区域覆盖拦截方法
CN116222310B (zh) * 2023-04-13 2024-04-26 哈尔滨工业大学 三维空间中基于rbf_g的二对一同步区域覆盖拦截方法

Also Published As

Publication number Publication date
CN112256055B (zh) 2021-11-05

Similar Documents

Publication Publication Date Title
CN112256055B (zh) 一种采用不动点优化的三体对抗防御预测制导方法
CN108534614B (zh) 一种三维全向实时预测制导方法
Li et al. A UAV path planning method based on deep reinforcement learning
Yu et al. Guidance law with circular no-fly zone constraint
CN112648886B (zh) 一种组合制导目标拦截方法及系统
Zhao et al. Multiple missiles cooperative guidance based on leader-follower strategy
CN114020021A (zh) 一种多导弹分布式协同制导律的设计方法及系统
CN115951695A (zh) 空战模拟环境中基于三方博弈的动态战术控制域解算方法
CN113911398B (zh) 飞行器单脉冲躲避策略确定方法及系统
CN112033235B (zh) 一种采用黄金分割搜索的空战三体对抗最优预测制导方法
CN114035616A (zh) 一种飞行器对移动目标打击控制方法及系统
CN116974208B (zh) 基于捷联导引头的旋翼无人机目标打击控制方法及系统
CN110657713B (zh) 采用预测制导的主动防御导弹最优发射时机搜索方法
CN115129085B (zh) 一种多群组机器人协同执行任务的方法
Lin et al. Control strategies for target-attacker-defender games of USVs
CN114819055A (zh) 一种基于lstm网络的导弹拦截点预测方法
Dong et al. An integrated scheme of a smart net capturer for MUAVs
Lee et al. Impact time and angle control guidance with rendezvous concept
Chauhan et al. A Computing Based Simulation Model for Missile Guidance in Planar Domain
Wang et al. The application of proportional navigation in the process of UAV air combat guidance and optimization of proportional parameter
CN113051706B (zh) 一种采用虚拟点预测的三体对抗防御制导方法
Wang et al. An online proportional guidance midcourse guidance method for near space targets based on air-breathing hypersonic platform
Singh et al. Impact angle constrained guidance with collision avoidance from multiple moving obstacles
Li et al. Multi-target Strike Planning in Unknown Dynamic Environment
Hanghang et al. A method to solve the unreachable zone of mid-range air-to-air missile

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