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

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

Info

Publication number
CN112504016B
CN112504016B CN202010993510.9A CN202010993510A CN112504016B CN 112504016 B CN112504016 B CN 112504016B CN 202010993510 A CN202010993510 A CN 202010993510A CN 112504016 B CN112504016 B CN 112504016B
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.)
Active
Application number
CN202010993510.9A
Other languages
English (en)
Other versions
CN112504016A (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 GDA0002933818450000021
Figure GDA0002933818450000022
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure GDA0002933818450000023
为高度上界,Vi为导弹速度,V为速度下界,
Figure GDA0002933818450000024
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure GDA0002933818450000025
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure GDA0002933818450000026
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure 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 GDA0002933818450000028
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure GDA0002933818450000031
为目标相对导弹X方向的上下界,
Figure GDA0002933818450000032
为目标相对导弹Y方向的上下界,
Figure GDA0002933818450000033
为目标相对导弹Z方向的上下界。
优选地,所述的步骤4中目标逃逸代价模型为:
Figure GDA0002933818450000034
Figure GDA0002933818450000035
Figure GDA0002933818450000036
J=ζx×ζy×ζz
式中,ζxyz分别表示目标在惯性坐标系下X、Y、Z三个方向的逃逸代价,J为目标逃逸总代价。
本发明一种适应协同任务规划的目标不可逃逸区可靠预测方法具有以下优点:采用以当前飞行状态进行插值的方式求解目标不可逃逸区能够避免复杂的算法计算,满足导弹高速飞行下快速在线任务规划需求;通过对目标逃逸代价的计算,能够扩展一般任务规划模型中对于导弹命中目标能力的描述,兼顾最优任务分配和有效命中目标。本方法通过建立导弹攻击目标不可逃逸区模型,解决了一般任务规划模型只考虑任务执行代价,难以适用于多弹协同任务规划的问题,能够为多弹协同任务规划的应用提供技术支撑。
附图说明
图1为本发明一种适应协同任务规划的目标不可逃逸区可靠预测方法的流程图;
图2为目标不可逃逸区建模示意图;
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,一种适应协同任务规划的目标不可逃逸区可靠预测方法,包含以下步骤:
步骤1:选择导弹飞行包络中的至少2个特征点;
根据已知导弹飞行包络,按照高度、速度、弹道倾角、弹道偏角、推力变化区间选取导弹飞行空域中的N个特征点,并以导弹位置矢量、速度矢量的形式进行记录。
特征点初始状态信息按照如下特征向量确定:
Ui=[Hi,Viivi,Pi]T
Figure GDA0002933818450000041
Figure GDA0002933818450000042
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure GDA0002933818450000043
为高度上界,Vi为导弹速度,V为速度下界,
Figure GDA0002933818450000044
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure GDA0002933818450000045
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure GDA0002933818450000046
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure 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 GDA0002933818450000051
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure GDA0002933818450000052
为目标相对导弹X方向的上下界,
Figure GDA0002933818450000053
为目标相对导弹Y方向的上下界,
Figure GDA0002933818450000054
为目标相对导弹Z方向的上下界。
步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力。
根据导弹和目标当前状态插值计算目标不可逃逸区,结合目标当前位置计算目标逃逸代价。目标逃逸代价模型为:
Figure GDA0002933818450000055
Figure GDA0002933818450000056
Figure 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 GDA0002933818450000061
并按照步骤4计算出目标逃逸代价J,如表1所示。
表1目标不可逃逸区预测结果
Figure GDA0002933818450000062
可以看出,本方法能够准确根据导弹和目标的飞行状态,预测目标不可逃逸区,并计算目标逃逸代价,为协同任务规划提供了支撑。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
本发明未详细说明部分属于本领域技术人员公知常识。

Claims (3)

1.一种适应协同任务规划的目标不可逃逸区可靠预测方法,其特征在于,包含以下步骤:
步骤1:选择导弹飞行包络中的至少2个特征点;
步骤2:根据目标的位置和速度,设置目标在空域中的空间分布;
步骤3:以导弹飞行包络中的特征点为起始点,目标为终点,建立目标不可逃逸区,并建立导弹当前状态即特征点的初始状态与目标不可逃逸区的映射关系;
步骤4:根据导弹当前状态和目标不可逃逸区映射关系获得目标不可逃逸区,建立目标逃逸代价模型,计算当前目标逃逸代价,预测导弹拦截能力;
所述的步骤3中导弹当前状态与目标不可逃逸区的映射关系建立方法为:
步骤31:选择导弹的每一个特征点,通过蒙特卡洛仿真在六自由度仿真模型中循环攻击设置好的空间分布目标;
步骤32:记录导弹每一个特征点的初始状态,当该特征点所攻击的目标脱靶量在设置的最大脱靶量范围内时,记录该目标相对导弹位置、相对导弹速度;
步骤33:对于每一个特征点,统计所记录的目标相对导弹位置的上下界,并以特征点初始状态、目标相对导弹速度为输入,通过如下插值函数形式建立导弹当前状态与目标不可逃逸区的映射关系;
Figure FDA0003840337810000011
式中,fInterp(Ui,Vmtxi,Vmtyi,Vmtzi)为标准线性插值函数,
Figure FDA0003840337810000012
为目标相对导弹X方向的上下界,
Figure FDA0003840337810000013
为目标相对导弹Y方向的上下界,
Figure FDA0003840337810000014
为目标相对导弹Z方向的上下界;
所述的步骤4中目标逃逸代价模型为:
Figure FDA0003840337810000021
Figure FDA0003840337810000022
Figure FDA0003840337810000023
J=ζx×ζy×ζz
式中,ζxyz分别表示目标在惯性坐标系下X、Y、Z三个方向的逃逸代价,J为目标逃逸总代价。
2.如权利要求1所述的目标不可逃逸区可靠预测方法,其特征在于,所述的步骤1中,飞行包络内的特征点的初始状态信息按照如下特征向量形式确定:
Ui=[Hi,Viivi,Pi]T
Figure FDA0003840337810000024
Figure FDA0003840337810000025
式中,Ui表示飞行包络中第i个特征向量,Hi为导弹高度,H为高度下界,
Figure FDA0003840337810000026
为高度上界,Vi为导弹速度,V为速度下界,
Figure FDA0003840337810000027
为速度上界,θi为导弹弹道倾角,θ为弹道倾角下界,
Figure FDA0003840337810000028
为弹道倾角上界,ψvi为导弹弹道偏角,ψ v为弹道偏角下界,
Figure FDA0003840337810000029
为弹道偏角上界,Pi为导弹推力,P为推力下界,
Figure FDA00038403378100000210
为推力上界,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方向的表示。
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 CN112504016A (zh) 2021-03-16
CN112504016B true 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)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AUPR080400A0 (en) * 2000-10-17 2001-01-11 Electro Optic Systems Pty Limited Autonomous weapon system
GB2380244B (en) * 2001-08-13 2006-02-15 Joseph Zabrana Michael Automated Sound Missile and Associated Defence System
FR2833722B1 (fr) * 2001-12-18 2005-06-17 Thales Sa Procede d'optimisation de l'interception de cibles mobiles, systeme d'interception et missile mettant en oeuvre un tel procede
IL230327B (en) * 2014-01-01 2019-11-28 Israel Aerospace Ind Ltd An interceptor missile and a warhead for it
GB2551874B (en) * 2016-04-25 2018-09-05 Bae Systems Plc System integration
CN105976061A (zh) * 2016-05-11 2016-09-28 黄安祥 一种提高虚拟智能对手飞机生存概率的方法
CN111091273B (zh) * 2019-11-29 2023-08-11 上海航天控制技术研究所 一种基于能力预测的多弹协同任务规划方法
CN111191359B (zh) * 2019-12-24 2021-03-23 中国人民解放军海军工程大学 防空导弹空间杀伤区模型构建和射击诸元解算方法

Also Published As

Publication number Publication date
CN112504016A (zh) 2021-03-16

Similar Documents

Publication Publication Date Title
CN111580556B (zh) 一种时空约束下多无人机协同路径规划与制导方法
CN109871031A (zh) 一种固定翼无人机的轨迹规划方法
CN111027206A (zh) 具有规定性能的拦截机动目标自适应滑模控制方法
CN113093733B (zh) 一种无人艇集群对海打击方法
CN115951695A (zh) 空战模拟环境中基于三方博弈的动态战术控制域解算方法
CN115755969A (zh) 一种基于零控拦截流型的飞行器制导方法及装置
CN115033024A (zh) 基于攻击时间和角度共同约束的多导弹三维协同制导方法
CN112504016B (zh) 一种适应协同任务规划的目标不可逃逸区可靠预测方法
CN114740883A (zh) 一种协同点侦察任务规划跨层联合优化方法
CN114035616A (zh) 一种飞行器对移动目标打击控制方法及系统
CN109543135A (zh) 一种基于信息最大化的多auv协同同时导引方法
CN116795130A (zh) 一种针对机动目标的避障制导方法
CN116227343A (zh) 拦截机动目标满足视场角约束的末角攻击制导律设计方法
CN114153143B (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN116225053A (zh) 一种面向任务分级目标运动特性的跟踪轨迹优化方法
CN115145295A (zh) 一种动态环境下无人机在线自主航迹优化控制方法
CN114662285A (zh) 一种高速飞行器火力控制模型的智能解算方法
CN112379599A (zh) 一种针对空中机动目标的三维协同中制导方法
Shin et al. Membership set-based mid course guidance: application to manoeuvreing target interception
CN113110576B (zh) 一种可连续切换的自适应固定时间收敛协同制导方法
Ni et al. Optimized Application Of Multi-UAV Electromagnetic Deception Jamming
CN113051706B (zh) 一种采用虚拟点预测的三体对抗防御制导方法
CN117707205B (zh) 制导飞行器速度控制方法及系统
CN112861314B (zh) 一种适用于水下航行器搜索策略模拟的随机目标模型
CN114200439A (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