CN112129292B - 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法 - Google Patents

考虑攻击时间和攻击角度约束的异构飞行器协同制导方法 Download PDF

Info

Publication number
CN112129292B
CN112129292B CN202010909551.5A CN202010909551A CN112129292B CN 112129292 B CN112129292 B CN 112129292B CN 202010909551 A CN202010909551 A CN 202010909551A CN 112129292 B CN112129292 B CN 112129292B
Authority
CN
China
Prior art keywords
aircraft
target
slave
attack
master
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
CN202010909551.5A
Other languages
English (en)
Other versions
CN112129292A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN202010909551.5A priority Critical patent/CN112129292B/zh
Publication of CN112129292A publication Critical patent/CN112129292A/zh
Application granted granted Critical
Publication of CN112129292B publication Critical patent/CN112129292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Traffic Control Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供考虑攻击时间和攻击角度约束的异构飞行器协同制导方法,包括如下步骤:步骤1:建立飞行器数学模型以及相对运动关系;步骤2:为主飞行器设计指定时间命中制导律;步骤3:从飞行器制导律的设计。其中,飞行器群由单个主飞行器和多个从飞行器构成,仅有主飞行器可接收攻击时间指令。该方法首先建立了主飞行器和从飞行器的数学模型,给出了主飞行器与从飞行器相对目标的运动关系,在此基础上为主飞行器设计了指定时刻命中目标的导引律;为了确保从飞行器与主飞行器命中目标的同时性,设计协同制导律使得从飞行器与目标的距离和主飞行器与目标的距离协同一致,并且从飞行器相对目标的视线角趋于预设的期望角。

Description

考虑攻击时间和攻击角度约束的异构飞行器协同制导方法
技术领域
本发明涉及一种考虑攻击时间和攻击角度约束的异构飞行器协同制导方法,属于飞行器制导、控制领域,具体为设计协同制导策略实现整个异构飞行器群在指定时刻同时命中目标。
背景技术
在军事作战中,高价值的目标通常是有防御系统保护的,因此,采用单枚飞行器独立制导的方式打击此类目标易被防御系统拦截,作战成功率低。多飞行器协同制导被视作突破防御系统打击目标的有效方法之一,是实现由“以平台为中心”的单个飞行器作战方式向“以网络为中心”的多飞行器作战方式转变的重要方向。通过协同作战形成各飞行器相互配合的打击效果,使得整个飞行器群大系统的作战能力超过各飞行器作战能力的简单和。因此,协同制导的目的是利用多飞行器协同攻击提高导弹突防能力和综合作战效能。
然而,现有的协同制导方法大都集中于同构飞行器群的设计,对于异构飞行器的协同制导方法研究尚少。此外,对于多个飞行器指定时间打击目标,已有的制导方法是给每个飞行器分配相同的攻击时间指令,要求每个飞行器都具备接收攻击指令的能力,且飞行器之间没有信息交互,制导律独立设计,因此信息共享和协同配合能力弱。因此,如何在只有个别飞行器具备接收攻击指令的条件下,面向异构飞行器群设计协同制导方法实现指定时刻同时命中目标具有重要意义。
发明内容
本发明针对异构飞行器的指定时间打击协同制导问题,提出了一种攻击时间控制的“主-从”式协同制导方法,其中主飞行器推力发动机关闭,只能通过法向过载设计导引律。从飞行器装备有低功率推力发动机,可同时利用切向过载和法向过载。主飞行器可以给满足通信条件的从飞行器发送信息,但不接收从飞行器信息,满足通信条件的从飞行器之间可互相传送信息,部分在主飞行器通信范围内的从飞行器可接收主飞行器与目标的距离信息。因此,主飞行器的制导律独立设计,不依赖于从飞行器的状态信息。从飞行器的制导律除了自身状态信息之外,还与其它可通信交互的飞行器状态有关。
本发明的技术构思为:为主飞行器独立设计制导律,保证其能在指定时刻命中目标,为从飞行器设计协同制导律使得其与目标的距离和主飞行器与目标的距离协同一致,并与目标保持特定的视线角,实现主飞行器与从飞行器同时命中目标。
本发明基于固定时间收敛控制理论和一致性协同控制理论,为主飞行器设计攻击时间控制制导方法,实现其对目标的指定时间打击。设计协同制导律使得从飞行器与目标的距离和主飞行器与目标的距离保持协同一致,并为每个从飞行器分配不同的视线角,从而实现主飞行器和从飞行器对目标的指定时间打击。
本发明为一种考虑攻击时间和攻击角度约束的异构飞行器协同制导方法,包括以下步骤:
步骤1:建立飞行器数学模型以及相对运动关系
本发明考虑n+1个飞行器打击静止目标,飞行器由1个推力发动机关闭的主飞行器和n个配置小功率推力发动机的从飞行器组成。从飞行器之间的通信关系可以采用图论中的邻接矩阵A=[aij]来表示,如果第i(i=1,2,…,n)个从飞行器能够和第j个从飞行器建立通信关系则aij=1,否则,aij=0。为了方便表示,与第i个从飞行器相关的变量采用下标i表示,与主飞行器相关的变量用下标0来表示。当通信拓扑中的任何两个飞行器都能找到至少一条通信路径,则该通信拓扑图为连通图。如果通信拓扑图中的链路都是双向的,即aij=aji,则称该图为无向图,如果存在单向的链路则称之有向图。主飞行器与从飞行器之间的通信关系采用mi来表示,如果第i个从飞行器能够接收主飞行器信息则mi=1,否则mi=0。飞行器与目标的作战几何结构如图1所示。
主飞行器相对目标的运动关系可以表示为:
Figure BDA0002662763720000031
式中,r0为主飞行器与目标之间的距离,V0表示恒定的主飞行器速度,θ0为主飞行器相对目标的视线角,γ0为主飞行器航迹角,φ0为主飞行器前置角,法向角速度a0垂直于速度方向,
Figure BDA0002662763720000032
Figure BDA0002662763720000033
表示r0,θ0和γ0的导数,此表示方法以下类同。
从飞行器相对目标的运动关系可以表示为:
Figure BDA0002662763720000041
式中,ri表示第i个从飞行器与目标的距离,θi表示第i个从飞行器相对目标的视线角,γi为从飞行器航迹角,φi为从飞行器前置角,Vi为从飞行器速度,与主飞行器不同,配置推力发动机的从飞行器可以分别利用切向加速度at,i和法向加速度an,i来控制速度的幅值和方向。
我们定义
Figure BDA0002662763720000042
则第i个从飞行器相对目标的运动关系还可表示为:
Figure BDA0002662763720000043
式中,ar,i和aθ,i为第i个从飞行器加速度沿着和垂直于相对目标视线方向的分量,Vr,i和Vθ,i分别为从飞行器与目标的相对速度沿着和垂直于视线方向的分量,其满足以下关系:
Figure BDA0002662763720000051
步骤2:为主飞行器设计指定时间命中制导律
主飞行器的制导律独立设计,不受从飞行器的影响,主飞行器的剩余命中时间可采用下式进行预测:
Figure BDA0002662763720000052
式中,
Figure BDA0002662763720000053
表示预测的剩余命中时间,Ns>2表示导航比。
Figure BDA0002662763720000054
的动态满足关系:
Figure BDA0002662763720000055
由于φ0通常较小,可等效建立
Figure BDA0002662763720000056
和φ0=sinφ0,因此有
Figure BDA0002662763720000057
Figure BDA0002662763720000058
定义攻击时间误差为:
Figure BDA0002662763720000059
式中,t表示当前时刻,Td为期望的攻击时刻。
主飞行器的指定时间打击制导律设计为:
Figure BDA00026627637200000510
其中,0<μ1<1为正实数,k为增益。
步骤3:从飞行器制导律的设计
定义:
Figure BDA0002662763720000061
式中,ξi,1为距离一致性协同误差,αi为虚拟控制项,ξi,2为虚拟控制误差,si为定义的滑模变量,li为正实数,θf,i为期望的视线角。
设计从飞行器制导律为:
Figure BDA0002662763720000062
式中,0<bj<1(j=1,2,3)为正实数,ki,1,ki,2,…,ki,6为正实数;
本发明的有益效果为:虽然仅有主飞行器才能接收到攻击时间指令,但在所设计的攻击时间控制协同制导律下,整个异构飞行器群都能在指定时刻同时命中目标。并且,从飞行器相对目标能保持期望的视线角。
附图说明
图1是飞行器与目标的作战几何结构图。
图2是飞行器的飞行轨迹曲线。
图3是飞行器与目标之间的距离变化图。
图4是从飞行器相对目标的视线角响应图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,参照附图1—4对本发明做进一步说明。
所发明攻击时间控制的异构飞行器协同制导方法包含以下步骤:
步骤1:建立飞行器数学模型以及相对运动关系
主飞行器相对目标的运动关系可以表示为:
Figure BDA0002662763720000071
式中,r0为主飞行器与目标之间的距离,V0表示恒定的主飞行器速度,θ0为主飞行器相对目标的视线角,γ0为主飞行器航迹角,φ0为主飞行器前置角,法向角速度a0垂直于速度方向。
从飞行器相对目标的运动关系可以表示为:
Figure BDA0002662763720000072
式中,ri表示第i个从飞行器与目标的距离,γi为从飞行器航迹角,φi为从飞行器前置角,Vi为从飞行器速度,θi表示第i个飞行器相对目标的视线角,与主飞行器不同,配置推力发动机的从飞行器可以分别利用切向加速度at,i和法向加速度an,i来控制速度的幅值和方向。
我们定义
Figure BDA0002662763720000081
则第i个从飞行器相对目标的运动关系还可表示为:
Figure BDA0002662763720000082
Figure BDA0002662763720000083
Figure BDA0002662763720000084
Figure BDA0002662763720000085
式中,ar,i和aθ,i为第i个从飞行器加速度沿着和垂直于相对目标视线方向的分量,其满足以下关系:
Figure BDA0002662763720000091
步骤2:为主飞行器设计指定时间命中制导律
主飞行器的制导律独立设计,不受从飞行器的影响,主飞行器的剩余命中时间可采用下式进行预测:
Figure BDA0002662763720000092
式中,
Figure BDA0002662763720000093
表示预测的剩余命中时间,Ns>2表示导航比。
Figure BDA0002662763720000094
的动态满足关系:
Figure BDA0002662763720000095
由于φ0通常较小,可等效建立
Figure BDA0002662763720000096
和φ0=sinφ0,因此有
Figure BDA0002662763720000097
Figure BDA0002662763720000098
定义攻击时间误差为:
Figure BDA0002662763720000099
式中,t表示当前时刻,Td为期望的攻击时刻。
主飞行器的指定时间打击制导律设计为:
Figure BDA00026627637200000910
其中,0<μ1<1为正实数,δ>0为连续切换边界,k为正常数。
步骤3:从飞行器制导律设计为:
Figure BDA0002662763720000101
ξi,2=Vr,ii
Figure BDA0002662763720000102
式中,αi为虚拟控制项,li为正实数,θf,i为期望的视线角。
设计从飞行器制导律为:
Figure BDA0002662763720000103
式中,0<bj<1(j=1,2,3)为正实数,ki,1,ki,2,…,ki,6为正实数。
为验证所发明考虑碰撞规避的异构飞行器协同制导方法的有效性,利用Matlab仿真平台执行所设计算法。对于实施例,本发明所设计协同制导方法的参数选取为:Td=65s,k=6,μ1=0.5,Ns=4,θf,1=π/3,θf,2=π/4,θf,3=π/6,li=0.1,ki,1=2,ki,2=0.01,ki,3=0.5,ki,4=0.001,ki,5=ki,6=2,bi=0.5。
所有飞行器的飞行轨迹和与目标的距离如图2和图3所示,可看出主飞行器和从飞行器都在指定时刻Td=65s击中目标。图4给出了各从飞行器相对主飞行器的视线角变化曲线,由结果可看出,θi在所设计的制导律下能够收敛至θf,i,仿真结果与理论分析一致。

Claims (2)

1.一种考虑攻击时间和攻击角度约束的异构飞行器协同制导方法,其特征在于,包括以下步骤:
步骤1:建立飞行器数学模型以及相对运动关系
考虑n+1个飞行器打击静止目标,飞行器由1个推力发动机关闭的主飞行器和n个配置小功率推力发动机的从飞行器组成;从飞行器之间的通信关系采用图论中的邻接矩阵A=[aij]来表示,如果第i个从飞行器能够和第j个从飞行器建立通信关系则aij=1,i=1,2,…,n;否则,aij=0;
主飞行器相对目标的运动关系表示为:
Figure FDA0002662763710000011
式中,r0为主飞行器与目标之间的距离,V0表示恒定的主飞行器速度,θ0为主飞行器相对目标的视线角,γ0为主飞行器航迹角,φ0为主飞行器前置角,法向角速度a0垂直于速度方向,
Figure FDA0002662763710000012
Figure FDA0002662763710000013
表示r0,θ0和γ0的导数,此表示方法以下类同;
从飞行器相对目标的运动关系表示为:
Figure FDA0002662763710000021
式中,ri表示第i个从飞行器与目标的距离,θi表示第i个从飞行器相对目标的视线角,γi为从飞行器航迹角,φi为从飞行器前置角,Vi为从飞行器速度,与主飞行器不同,配置推力发动机的从飞行器分别利用切向加速度at,i和法向加速度an,i来控制速度的幅值和方向;
定义
Figure FDA0002662763710000022
则第i个从飞行器相对目标的运动关系还表示为:
Figure FDA0002662763710000031
Figure FDA0002662763710000032
Figure FDA0002662763710000033
Figure FDA0002662763710000034
式中,ar,i和aθ,i为第i个从飞行器加速度沿着和垂直于相对目标视线方向的分量,Vr,i和Vθ,i分别为从飞行器与目标的相对速度沿着和垂直于视线方向的分量,其满足以下关系:
Figure FDA0002662763710000035
步骤2:为主飞行器设计指定时间命中制导律
主飞行器的制导律独立设计,不受从飞行器的影响,主飞行器的剩余命中时间采用下式进行预测:
Figure FDA0002662763710000036
式中,
Figure FDA0002662763710000037
表示预测的剩余命中时间,Ns>2表示导航比;
Figure FDA0002662763710000038
的动态满足关系:
Figure FDA0002662763710000039
由于φ0较小,等效建立
Figure FDA00026627637100000310
和φ0=sinφ0,因此有
Figure FDA0002662763710000041
Figure FDA0002662763710000042
定义攻击时间误差为:
Figure FDA0002662763710000043
式中,t表示当前时刻,Td为期望的攻击时刻;
主飞行器的指定时间打击制导律设计为:
Figure FDA0002662763710000044
其中,0<μ1<1为正实数,k为增益;
步骤3:从飞行器制导律的设计
定义:
Figure FDA0002662763710000045
ξi,2=Vr,ii
Figure FDA0002662763710000046
式中,ξi,1为距离一致性协同误差,αi为虚拟控制项,ξi,2为虚拟控制误差,si为定义的滑模变量,li为正实数,θf,i为期望的视线角;
设计从飞行器制导律为:
Figure FDA0002662763710000051
式中,0<bj<1(j=1,2,3)为正实数,ki,1,ki,2,…,ki,6为正实数。
2.根据权利要求1所述的考虑攻击时间和攻击角度约束的异构飞行器协同制导方法,其特征在于:当通信拓扑中的任何两个飞行器都能找到至少一条通信路径,则该通信拓扑图为连通图;如果通信拓扑图中的链路都是双向的,即aij=aji,则称该图为无向图,如果存在单向的链路则称之有向图;主飞行器与从飞行器之间的通信关系采用mi来表示,如果第i个从飞行器能够接收主飞行器信息则mi=1,否则mi=0。
CN202010909551.5A 2020-09-02 2020-09-02 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法 Active CN112129292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010909551.5A CN112129292B (zh) 2020-09-02 2020-09-02 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010909551.5A CN112129292B (zh) 2020-09-02 2020-09-02 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法

Publications (2)

Publication Number Publication Date
CN112129292A CN112129292A (zh) 2020-12-25
CN112129292B true CN112129292B (zh) 2021-06-04

Family

ID=73847755

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010909551.5A Active CN112129292B (zh) 2020-09-02 2020-09-02 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法

Country Status (1)

Country Link
CN (1) CN112129292B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113110576B (zh) * 2021-04-14 2022-09-13 西北工业大学 一种可连续切换的自适应固定时间收敛协同制导方法
CN113359813B (zh) * 2021-05-11 2022-08-05 西北工业大学 一种多领弹多从弹群组协同制导方法
CN114489109B (zh) * 2021-09-13 2023-06-30 西北工业大学 一种三维主从式多飞行器攻击时间控制协同制导方法
CN113867388A (zh) * 2021-09-30 2021-12-31 北京机电工程研究所 多约束条件下针对机动目标的分布式协同制导方法及装置
CN113741548A (zh) * 2021-10-20 2021-12-03 北京机电工程研究所 一种无人机编队非线性协同导引方法、装置及存储介质
CN114638298B (zh) * 2022-03-15 2023-06-09 天津大学 飞行器攻击行为预测方法及电子设备
CN114384808B (zh) * 2022-03-24 2022-07-15 西北工业大学 一种基于自适应神经网络的到达时间控制三维制导方法
CN115016544B (zh) * 2022-08-04 2022-10-28 中国科学院自动化研究所 基于多个飞行器的协同制导方法、装置和电子设备
CN117193336B (zh) * 2023-09-12 2024-03-22 中国船舶集团有限公司第七一九研究所 一种基于时间和角度约束的水下航行物协同制导方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108718463A (zh) * 2018-05-09 2018-10-30 南京邮电大学 一种h-cran下基于多时间尺度协同优化的资源分配方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706217B (zh) * 2012-04-17 2014-07-02 北京理工大学 一种控制多枚导弹攻击角度和攻击时间的方法
CN103245256B (zh) * 2013-04-25 2015-03-11 北京理工大学 一种多导弹协同作战制导律设计方法
CN108168381B (zh) * 2018-01-04 2019-10-08 北京理工大学 一种多枚导弹协同作战的控制方法
CN109084641B (zh) * 2018-08-15 2019-07-26 北京航空航天大学 导弹制导方法及装置
CN111274740B (zh) * 2020-01-10 2021-02-12 中国人民解放军国防科技大学 一种多飞行器协同突防轨迹优化设计方法
CN111551080B (zh) * 2020-04-03 2021-02-05 北京航空航天大学 一种攻击时间控制的制导方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108718463A (zh) * 2018-05-09 2018-10-30 南京邮电大学 一种h-cran下基于多时间尺度协同优化的资源分配方法

Also Published As

Publication number Publication date
CN112129292A (zh) 2020-12-25

Similar Documents

Publication Publication Date Title
CN112129292B (zh) 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法
Yang et al. UAV air combat autonomous maneuver decision based on DDPG algorithm
Zhang et al. Three-dimensional path planning for uninhabited combat aerial vehicle based on predator-prey pigeon-inspired optimization in dynamic environment
Xu et al. Virtual motion camouflage based phantom track generation through cooperative electronic combat air vehicles
Li et al. A Multi-UCAV cooperative occupation method based on weapon engagement zones for beyond-visual-range air combat
CN114489109B (zh) 一种三维主从式多飞行器攻击时间控制协同制导方法
CN113741500B (zh) 仿哈里斯鹰智能捕食优化的无人机空战机动决策方法
CN114063644A (zh) 基于鸽群反向对抗学习的无人作战飞机空战自主决策方法
Kung et al. Missile guidance algorithm design using particle swarm optimization
Xiang et al. Research on UAV swarm confrontation task based on MADDPG algorithm
CN115951695A (zh) 空战模拟环境中基于三方博弈的动态战术控制域解算方法
CN113671825A (zh) 一种基于强化学习的机动智能决策规避导弹方法
CN113741186A (zh) 一种基于近端策略优化的双机空战决策方法
You et al. Distributed synergetic guidance law for multiple missiles with angle-of-attack constraint
CN115857340A (zh) 基于递进打靶法的双航天器博弈路径规划方法
Sun et al. The study on intelligent advanced fighter air combat decision support system
Wei et al. UCAV formation online collaborative trajectory planning using hp adaptive pseudospectral method
CN113110576B (zh) 一种可连续切换的自适应固定时间收敛协同制导方法
Xu et al. The application of situation function in differential game problem of the air combat
Fu et al. Survey of manned/unmanned air combat decision technology
Gao et al. Real-time calculation of tactical control range in beyond visual range air combat
McManus et al. Situation assessment in the Paladin tactical decision generation system
CN116679755B (zh) 基于分布式模型预测控制的多智能体协同突防打击方法
CN116679742B (zh) 一种多六自由度飞行器协同作战决策方法
Xiaoyu et al. Intelligent Air Combat Maneuvering Decision Based on TD3 Algorithm

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
CB03 Change of inventor or designer information

Inventor after: Lv Jinhu

Inventor after: Li Guofei

Inventor after: Liu Kexin

Inventor after: Zhu Guoliang

Inventor after: Wang Xiaomo

Inventor after: Wang Chengcai

Inventor before: Lv Jinhu

Inventor before: Li Guofei

Inventor before: Liu Kexin

Inventor before: Zhu Guoliang

CB03 Change of inventor or designer information