CN112504016A - 一种适应协同任务规划的目标不可逃逸区可靠预测方法 - Google Patents

一种适应协同任务规划的目标不可逃逸区可靠预测方法 Download PDF

Info

Publication number
CN112504016A
CN112504016A CN202010993510.9A CN202010993510A CN112504016A CN 112504016 A CN112504016 A CN 112504016A CN 202010993510 A CN202010993510 A CN 202010993510A CN 112504016 A CN112504016 A CN 112504016A
Authority
CN
China
Prior art keywords
target
missile
escape
escape area
mti
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
CN202010993510.9A
Other languages
English (en)
Other versions
CN112504016B (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN202010993510.9A priority Critical patent/CN112504016B/zh
Publication of CN112504016A publication Critical patent/CN112504016A/zh
Application granted granted Critical
Publication of CN112504016B publication Critical patent/CN112504016B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41HARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
    • F41H11/00Defence installations; Defence devices
    • F41H11/02Anti-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)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明提供了一种适应协同任务规划的目标不可逃逸区可靠预测方法,其特征在于,包含以下步骤:步骤1:选择导弹飞行包络中的至少2个特征点;步骤2:根据目标的位置和速度,设置目标在空域中的空间分布;步骤3:以导弹飞行包络中的特征点为起始点,目标为终点,建立目标不可逃逸区,并建立导弹当前状态即特征点的初始状态与目标不可逃逸区的映射关系;步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力。本发明根据导弹和目标的飞行状态信息,建立导弹攻击目标的不可逃逸区模型,有效预测目标逃逸代价,为导弹协同任务规划提供支撑。

Description

一种适应协同任务规划的目标不可逃逸区可靠预测方法
技术领域
本发明涉及飞行器协同任务规划技术领域,具体涉及一种适应协同任务规划的目标不可逃逸区可靠预测方法。
本发明的应用领域为战术武器的智能化,其特征是解决战术武器进行协同任务规划时,难以准确判断当前飞行状态能否有效命中目标的问题,能够为战术武器提供自主决策依据。
背景技术
随着未来作战模式体系化、网络化程度加深,导弹的作战空域不断扩大,有人机、无人机、导弹等高低价值目标混合,目标机动能力强,作战环境更加复杂。为了能够在高动态作战环境下实现多对多对抗需求,导弹需要进一步提高自主决策和任务规划能力,在满足任务需求的同时能够保证准确命中目标。一般任务规划方法主要通过考虑航程、任务时序、燃料消耗、目标价值等因素建立描述飞行器执行某一任务的代价模型,飞行器对任务目标的可达能力通过人为规定。实际多弹协同攻击目标时弹目相对运动快、飞行时间短,导弹攻击目标的能力无法通过人为规定,因此上述方法无法解决多弹协同攻击时目标的逃逸概率问题,需要采用更加精确的方法预测目标不可逃逸区,保证任务规划的可靠性。
发明内容
本发明的目的在于提供一种适应协同任务规划的目标不可逃逸区可靠预测方法,根据导弹和目标的飞行状态信息,建立导弹攻击目标的不可逃逸区模型,有效预测目标逃逸代价,为导弹协同任务规划提供支撑。
为了达到上述目的,本发明通过以下技术方案实现:一种适应协同任务规划的目标不可逃逸区可靠预测方法,其特点是,包含以下步骤:
步骤1:选择导弹飞行包络中的至少2个特征点;
步骤2:根据目标的位置和速度,设置目标在空域中的空间分布;
步骤3:以导弹飞行包络中的特征点为起始点,目标为终点,建立目标不可逃逸区,并建立导弹当前状态即特征点的初始状态与目标不可逃逸区的映射关系;
步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力。
优选地,所述的步骤1中,飞行包络内的特征点初始状态信息按照如下特征向量形式确定:
Ui=[Hi,Viivi,Pi]T
Figure RE-GDA0002933818450000021
Figure RE-GDA0002933818450000022
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure RE-GDA0002933818450000023
为高度上界,Vi为导弹速度,V为速度下界,
Figure RE-GDA0002933818450000024
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure RE-GDA0002933818450000025
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure RE-GDA0002933818450000026
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure RE-GDA0002933818450000027
为推力上界,i为特征点序号。
优选地,所述的步骤2中目标空间分布的表示方式为:
Umti=[Xmti,Ymti,Zmti,Vmtxi,Vmtyi,Vmtzi]T
式中,Umti表示目标空间分布,Xmti,Ymti,Zmti为目标相对导弹位置在惯性坐标系X、Y、Z方向的表示,Vmtxi,Vmtyi,Vmtzi为目标相对导弹速度在惯性坐标系X、Y、Z方向的表示。
优选地,所述的步骤3中导弹当前状态与目标不可逃逸区的映射关系建立方法为:
步骤31:选择导弹的每一个特征点,通过蒙特卡洛仿真在六自由度仿真模型中循环攻击设置好的空间分布目标;
步骤32:记录导弹每一个特征点的初始状态,当该特征点所攻击的目标脱靶量在设置的最大脱靶量范围内时,记录该目标相对导弹位置、相对导弹速度;
步骤33:对于每一个特征点,统计所记录的目标相对导弹位置的上下界,并以特征点初始状态、目标相对导弹速度为输入,通过如下插值函数形式建立导弹当前状态与目标不可逃逸区的映射关系;
Figure RE-GDA0002933818450000028
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure RE-GDA0002933818450000031
为目标相对导弹X方向的上下界,
Figure RE-GDA0002933818450000032
为目标相对导弹Y方向的上下界,
Figure RE-GDA0002933818450000033
为目标相对导弹Z方向的上下界。
优选地,所述的步骤4中目标逃逸代价模型为:
Figure RE-GDA0002933818450000034
Figure RE-GDA0002933818450000035
Figure RE-GDA0002933818450000036
J=ζx×ζy×ζz
式中,ζxyz分别表示目标在惯性坐标系下X、Y、Z三个方向的逃逸代价,J为目标逃逸总代价。
本发明一种适应协同任务规划的目标不可逃逸区可靠预测方法具有以下优点:采用以当前飞行状态进行插值的方式求解目标不可逃逸区能够避免复杂的算法计算,满足导弹高速飞行下快速在线任务规划需求;通过对目标逃逸代价的计算,能够扩展一般任务规划模型中对于导弹命中目标能力的描述,兼顾最优任务分配和有效命中目标。本方法通过建立导弹攻击目标不可逃逸区模型,解决了一般任务规划模型只考虑任务执行代价,难以适用于多弹协同任务规划的问题,能够为多弹协同任务规划的应用提供技术支撑。
附图说明
图1为本发明一种适应协同任务规划的目标不可逃逸区可靠预测方法的流程图;
图2为目标不可逃逸区建模示意图;
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,一种适应协同任务规划的目标不可逃逸区可靠预测方法,包含以下步骤:
步骤1:选择导弹飞行包络中的至少2个特征点;
根据已知导弹飞行包络,按照高度、速度、弹道倾角、弹道偏角、推力变化区间选取导弹飞行空域中的N个特征点,并以导弹位置矢量、速度矢量的形式进行记录。
特征点初始状态信息按照如下特征向量确定:
Ui=[Hi,Viivi,Pi]T
Figure RE-GDA0002933818450000041
Figure RE-GDA0002933818450000042
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure RE-GDA0002933818450000043
为高度上界,Vi为导弹速度,V为速度下界,
Figure RE-GDA0002933818450000044
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure RE-GDA0002933818450000045
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure RE-GDA0002933818450000046
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure RE-GDA0002933818450000047
为推力上界,i为特征点序号。
步骤2:根据目标的位置和速度,设置目标在空域中的空间分布;
根据实际目标的速度范围和飞行空域范围,以目标相对导弹位置矢量、目标相对导弹速度矢量的形式设置目标在空间的分布。目标空间分布的表示方式如下:
Umti=[Xmti,Ymti,Zmti,Vmtxi,Vmtyi,Vmtzi]T
式中,Umti表示目标空间分布,Xmti,Ymti,Zmti为目标相对导弹位置在惯性坐标系X、Y、Z方向的表示,Vmtxi,Vmtyi,Vmtzi为目标相对导弹速度在惯性坐标系X、Y、Z方向的表示。
步骤3:以导弹飞行包络中的特征点为起始点,目标为终点,建立目标不可逃逸区,并建立导弹当前状态即特征点的初始状态与目标不可逃逸区的映射关系;
以各特征点对应的高度、速度、弹道倾角、弹道偏角、推力和目标相对导弹速度为输入,以各特征点对应的可攻击目标相对位置的上下界为输出,以插值函数的方式建立导弹当前状态与目标不可逃逸区的映射关系。
导弹当前状态与目标不可逃逸区的映射关系建立方法为:
步骤31:选择导弹的每一个特征点,通过蒙特卡洛仿真在六自由度仿真模型中循环攻击设置好的空间分布目标;
步骤32:记录导弹每一个特征点的初始状态,当该特征点所攻击的目标脱靶量在设置的最大脱靶量范围内时,记录该目标相对导弹位置、相对导弹速度;
步骤33:对于每一个特征点,统计所记录的目标相对导弹位置的上下界,并以特征点初始状态、目标相对导弹速度为输入,通过如下插值函数形式建立导弹当前状态与目标不可逃逸区的映射关系;
Figure RE-GDA0002933818450000051
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure RE-GDA0002933818450000052
为目标相对导弹X方向的上下界,
Figure RE-GDA0002933818450000053
为目标相对导弹Y方向的上下界,
Figure RE-GDA0002933818450000054
为目标相对导弹Z方向的上下界。
步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力。
根据导弹和目标当前状态插值计算目标不可逃逸区,结合目标当前位置计算目标逃逸代价。目标逃逸代价模型为:
Figure RE-GDA0002933818450000055
Figure RE-GDA0002933818450000056
Figure RE-GDA0002933818450000057
J=ζx×ζy×ζz
式中,ζxyz分别表示目标在惯性坐标系下X、Y、Z三个方向的逃逸代价,J为目标逃逸总代价。
基于上述详细步骤,即可完成目标不可逃逸区预测。
具体应用:假设导弹和目标在纵向平面内运动,导弹飞行包络上特征点为:Ui=[Hi,Vi,0,0,0]T,且Hi的取值范围为[4000,11500],Vi的取值范围为[675,825]。目标相对导弹的空间分布为:Umti=[30000,Ymti,Zmti,-400,50,0]T,且 Ymti=Hi,Zmti的取值范围为[-8000,8000]。
按照本发明中步骤1-3建立导弹攻击目标不可逃逸区映射关系,并将导弹攻击目标不可逃逸区用三维图表示,如图2所示,图中,三维曲面为特征点上的导弹对其左右两侧目标的最大攻击包络,Z向最大攻击范围为 -5280m~5280m。通过设置表1所示的导弹初始状态和目标初始状态,以表中 Hi,Vi,Vmti为输入插值得到不可逃逸区Z mti
Figure RE-GDA0002933818450000061
并按照步骤4计算出目标逃逸代价J,如表1所示。
表1目标不可逃逸区预测结果
Figure RE-GDA0002933818450000062
可以看出,本方法能够准确根据导弹和目标的飞行状态,预测目标不可逃逸区,并计算目标逃逸代价,为协同任务规划提供了支撑。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
本发明未详细说明部分属于本领域技术人员公知常识。

Claims (5)

1.一种适应协同任务规划的目标不可逃逸区可靠预测方法,其特征在于,包含以下步骤:
步骤1:选择导弹飞行包络中的至少2个特征点;
步骤2:根据目标的位置和速度,设置目标在空域中的空间分布;
步骤3:以导弹飞行包络中的特征点为起始点,目标为终点,建立目标不可逃逸区,并建立导弹当前状态即特征点的初始状态与目标不可逃逸区的映射关系;
步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力。
2.如权利要求1所述的目标不可逃逸区可靠预测方法,其特征在于,所述的步骤1中,飞行包络内的特征点的初始状态信息按照如下特征向量形式确定:
Ui=[Hi,Viivi,Pi]T
Figure RE-FDA0002933818440000011
Figure RE-FDA0002933818440000012
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure RE-FDA0002933818440000013
为高度上界,Vi为导弹速度,V为速度下界,
Figure RE-FDA0002933818440000014
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure RE-FDA0002933818440000015
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure RE-FDA0002933818440000016
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure RE-FDA0002933818440000017
为推力上界,i为特征点序号。
3.如权利要求1所述的的目标不可逃逸区可靠预测方法,其特征在于,所述的步骤2中目标在空域中的空间分布的表示方式为:
Umti=[Xmti,Ymti,Zmti,Vmtxi,Vmtyi,Vmtzi]T
式中,Umti表示目标空间分布,Xmti,Ymti,Zmti为目标相对导弹位置在惯性坐标系X、Y、Z方向的表示,Vmtxi,Vmtyi,Vmtzi为目标相对导弹速度在惯性坐标系X、Y、Z方向的表示。
4.如权利要求1所述的目标不可逃逸区可靠预测方法,其特征在于,所述的步骤3中导弹当前状态与目标不可逃逸区的映射关系建立方法为:
步骤31:选择导弹的每一个特征点,通过蒙特卡洛仿真在六自由度仿真模型中循环攻击设置好的空间分布目标;
步骤32:记录导弹每一个特征点的初始状态,当该特征点所攻击的目标脱靶量在设置的最大脱靶量范围内时,记录该目标相对导弹位置、相对导弹速度;
步骤33:对于每一个特征点,统计所记录的目标相对导弹位置的上下界,并以特征点初始状态、目标相对导弹速度为输入,通过如下插值函数形式建立导弹当前状态与目标不可逃逸区的映射关系;
Figure RE-FDA0002933818440000021
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure RE-FDA0002933818440000022
为目标相对导弹X方向的上下界,
Figure RE-FDA0002933818440000023
为目标相对导弹Y方向的上下界,
Figure RE-FDA0002933818440000024
为目标相对导弹Z方向的上下界。
5.如权利要求1所述的目标不可逃逸区可靠预测方法,其特征在于,所述的步骤4中目标逃逸代价模型为:
Figure RE-FDA0002933818440000025
Figure RE-FDA0002933818440000026
Figure RE-FDA0002933818440000027
J=ζx×ζy×ζz
式中,ζxyz分别表示目标在惯性坐标系下X、Y、Z三个方向的逃逸代价,J为目标逃逸总代价。
CN202010993510.9A 2020-09-21 2020-09-21 一种适应协同任务规划的目标不可逃逸区可靠预测方法 Active CN112504016B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010993510.9A CN112504016B (zh) 2020-09-21 2020-09-21 一种适应协同任务规划的目标不可逃逸区可靠预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010993510.9A CN112504016B (zh) 2020-09-21 2020-09-21 一种适应协同任务规划的目标不可逃逸区可靠预测方法

Publications (2)

Publication Number Publication Date
CN112504016A true CN112504016A (zh) 2021-03-16
CN112504016B CN112504016B (zh) 2022-11-18

Family

ID=74953527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010993510.9A Active CN112504016B (zh) 2020-09-21 2020-09-21 一种适应协同任务规划的目标不可逃逸区可靠预测方法

Country Status (1)

Country Link
CN (1) CN112504016B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0119627D0 (en) * 2001-08-13 2001-10-03 Zabrama Michael Aquarius
WO2003054647A1 (fr) * 2001-12-18 2003-07-03 Thales Procede d'optimisation de l'interception de cibles mobiles, systeme d'interception et missile mettant en oeuvre un tel procede
US20040050240A1 (en) * 2000-10-17 2004-03-18 Greene Ben A. Autonomous weapon system
CN105976061A (zh) * 2016-05-11 2016-09-28 黄安祥 一种提高虚拟智能对手飞机生存概率的方法
US20160320165A1 (en) * 2014-01-01 2016-11-03 Israel Aerospace Industries Ltd. Interception missle and warhead therefor
GB201706440D0 (en) * 2016-04-25 2017-06-07 Bae Systems Plc System integration
CN111091273A (zh) * 2019-11-29 2020-05-01 上海航天控制技术研究所 一种基于能力预测的多弹协同任务规划方法
CN111191359A (zh) * 2019-12-24 2020-05-22 中国人民解放军海军工程大学 一种基于少量实际数据防空导弹空间杀伤区模型构建和射击诸元解算方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040050240A1 (en) * 2000-10-17 2004-03-18 Greene Ben A. Autonomous weapon system
GB0119627D0 (en) * 2001-08-13 2001-10-03 Zabrama Michael Aquarius
WO2003054647A1 (fr) * 2001-12-18 2003-07-03 Thales Procede d'optimisation de l'interception de cibles mobiles, systeme d'interception et missile mettant en oeuvre un tel procede
US20160320165A1 (en) * 2014-01-01 2016-11-03 Israel Aerospace Industries Ltd. Interception missle and warhead therefor
GB201706440D0 (en) * 2016-04-25 2017-06-07 Bae Systems Plc System integration
CN105976061A (zh) * 2016-05-11 2016-09-28 黄安祥 一种提高虚拟智能对手飞机生存概率的方法
CN111091273A (zh) * 2019-11-29 2020-05-01 上海航天控制技术研究所 一种基于能力预测的多弹协同任务规划方法
CN111191359A (zh) * 2019-12-24 2020-05-22 中国人民解放军海军工程大学 一种基于少量实际数据防空导弹空间杀伤区模型构建和射击诸元解算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
隋起胜 等: "《太空安全防御技术》", 30 June 2018 *

Also Published As

Publication number Publication date
CN112504016B (zh) 2022-11-18

Similar Documents

Publication Publication Date Title
CN111580556B (zh) 一种时空约束下多无人机协同路径规划与制导方法
Kumar et al. Impact time guidance for large heading errors using sliding mode control
CN108153328B (zh) 一种基于分段贝塞尔曲线的多导弹协同航迹规划方法
CN111027206B (zh) 具有规定性能的拦截机动目标自适应滑模控制方法
CN109871031A (zh) 一种固定翼无人机的轨迹规划方法
CN113093733B (zh) 一种无人艇集群对海打击方法
CN110187640A (zh) 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法
CN112824998A (zh) 马尔可夫决策过程的多无人机协同航路规划方法和装置
CN115951695A (zh) 空战模拟环境中基于三方博弈的动态战术控制域解算方法
CN109827478B (zh) 一种带落角约束与过载约束的制导方法
CN115033024A (zh) 基于攻击时间和角度共同约束的多导弹三维协同制导方法
CN114740883A (zh) 一种协同点侦察任务规划跨层联合优化方法
CN113093789B (zh) 一种基于路径点优选的飞行器禁飞区规避轨迹规划方法
CN112504016B (zh) 一种适应协同任务规划的目标不可逃逸区可靠预测方法
CN114035616A (zh) 一种飞行器对移动目标打击控制方法及系统
CN109543135A (zh) 一种基于信息最大化的多auv协同同时导引方法
Shi et al. Cooperative prediction guidance law in target-attacker-defender scenario
CN112379599B (zh) 一种针对空中机动目标的三维协同中制导方法
CN113110576B (zh) 一种可连续切换的自适应固定时间收敛协同制导方法
CN115145295A (zh) 一种动态环境下无人机在线自主航迹优化控制方法
Zhao et al. Real-time path planning for multi-UAVs with share of threats information
CN114756052B (zh) 一种基于无人机群的多目标协同追踪方法
Shin et al. Membership set-based mid course guidance: application to manoeuvreing target interception
Liu et al. Research on cooperative search technology of heterogeneous UAVs in complex environments
Ni et al. Optimized Application Of Multi-UAV Electromagnetic Deception Jamming

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