CN113359813B - 一种多领弹多从弹群组协同制导方法 - Google Patents

一种多领弹多从弹群组协同制导方法 Download PDF

Info

Publication number
CN113359813B
CN113359813B CN202110510192.0A CN202110510192A CN113359813B CN 113359813 B CN113359813 B CN 113359813B CN 202110510192 A CN202110510192 A CN 202110510192A CN 113359813 B CN113359813 B CN 113359813B
Authority
CN
China
Prior art keywords
projectile
slave
collar
group
target
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
CN202110510192.0A
Other languages
English (en)
Other versions
CN113359813A (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.)
Northwestern Polytechnical University
Shanghai Institute of Electromechanical Engineering
Original Assignee
Northwestern Polytechnical University
Shanghai Institute of Electromechanical Engineering
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 Northwestern Polytechnical University, Shanghai Institute of Electromechanical Engineering filed Critical Northwestern Polytechnical University
Priority to CN202110510192.0A priority Critical patent/CN113359813B/zh
Publication of CN113359813A publication Critical patent/CN113359813A/zh
Application granted granted Critical
Publication of CN113359813B publication Critical patent/CN113359813B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/107Simultaneous control of position or course in three dimensions specially adapted for missiles

Abstract

本发明提供一种多领弹多从弹群组协同制导方法,包括:步骤1:建立导弹相对目标的运动关系数学模型;步骤2:预测领弹的剩余命中时间;步骤3:定义领弹的一致性误差变量;步骤4:给出领弹的协同制导律;步骤5:定义第k组第f个从弹变量为;步骤6:给出从弹的协同制导律。通过本发明设计的方法,借助弹群间领弹的通信、以及弹群内的通信,通过建立运动关系数学模型,定义协同变量,分别给出了领弹和从弹的协同制导律,使得所有的导弹最终同时命中目标。

Description

一种多领弹多从弹群组协同制导方法
技术领域
本发明涉及一种多领弹多从弹群组协同制导方法,属于飞行器制导、控制领域。具体为发明了包含多领弹多从弹的群组协同制导律,使得所有导弹最终同时命中目标。
背景技术
协同制导方法能够使得多导弹实现对目标的同时命中,具有突防能力强和命中率高的优点。根据导弹在群体中承担“角色”的不同,协同制导方法可分为无主式协同制导律和基于“领从”式架构的协同制导律。对于“领从”式协同制导,现有的方法主要集中于单领弹多从弹的布局模式,多领弹多从弹的协同制导方法较少。如何设计出多领弹多从弹的群组协同制导方法,通过这种协同作战的方式实现对目标的同时命中具有重要意义。
发明内容
本发明针对多领弹多从弹协同打击静止目标问题,提出了一种多群组协同制导方法。攻击导弹由多个弹群组成,每个弹群包含1个领弹和多个从弹。借助弹群间领弹的通信、以及弹群内领从弹的通信,在建立运动关系数学模型,定义协同变量的基础上,针对领弹和从弹设计出多群组协同制导方法,实现所有导弹对目标的同时命中。
本发明的设计思路为:为多组导弹群设计协同制导方法,保证所有导弹能够最终同时命中目标。该方法首先建立领弹相对目标、从弹相对目标的运动关系数学模型,在此基础上预测领弹的剩余命中时间,并定义剩余命中时间的协同变量,设计出领弹的协同制导律。随之,建立弹群内从弹的一致性距离协同变量,设计出可保证与领弹弹目距离一致的从弹协同制导律,从而实现所有导弹最终同时命中目标。
本发明为一种多领弹多从弹群组协同制导方法,由以下步骤组成:
步骤1:建立导弹相对目标的运动关系数学模型
本发明假设n组弹群打击静止目标,每组弹群各包含1枚领弹和mk(k=1,2,…,n)枚从弹。领弹之间,以及组内的导弹之间在满足通信条件时可以交互信息,可采用图论表示导弹之间的通信关系。领弹之间的通信关系可用矩阵A=[aij]来表示,如果第i(i=1,2,…,n)枚领弹能够和第j(j=1,2,…,n,j≠i)枚领弹建立通信则其取值aij=1,否则,aij=0。采用
Figure BDA0003059920420000021
表示第k组飞行器第f(f=1,2,…,mk)枚从弹和第h(h=1,2,…,mk,f≠h)枚从弹的通信关系,如果可以通信则
Figure BDA0003059920420000022
否则
Figure BDA0003059920420000023
一个组中的领弹只可以发送信息给从弹,但不能接收从弹的信息,采用
Figure BDA0003059920420000024
表示第k组飞行器第f枚从弹与所在组领弹之间的关系,如果可接收领弹信息则
Figure BDA0003059920420000025
否则
Figure BDA0003059920420000026
如果任意两个弹体都能找到至少一条通信路径,则该通信拓扑图为连通图。若通信链路都为双向,则为无向图,若存在单向通信链路则为有向图。
针对二维平面下多导弹对固定目标的打击,假设弹间通信拓扑图为连通图,第k枚领弹与目标的相对运动学模型可由如下的方程描述:
Figure BDA0003059920420000027
式中,rk表示第k枚领弹与目标之间的距离,Vk表示恒定的导弹飞行速度,qk为视线角,γk和φk分别为导弹的航迹角和前置角,法向加速度ak垂直于速度方向。
第k组弹群中第f个从弹相对目标的运动关系可以表示为:
Figure BDA0003059920420000031
式中,
Figure BDA0003059920420000032
表示第k组弹群中第f个从弹与目标的距离,
Figure BDA0003059920420000033
Figure BDA0003059920420000034
表示第f个从弹的视线角、航迹角、前置角和速度,与领弹不同,从弹可以分别利用切向加速度
Figure BDA0003059920420000035
和法向加速度
Figure BDA0003059920420000036
来调整速度的大小和方向。
定义
Figure BDA0003059920420000037
那么可得:
Figure BDA0003059920420000041
式中,
Figure BDA0003059920420000042
Figure BDA0003059920420000043
为第k组弹群中第f个从弹加速度沿着和垂直于视线方向的分量,
Figure BDA0003059920420000044
Figure BDA0003059920420000045
分别为从弹的速度沿着和垂直于视线方向的分量,可建立以下关系:
Figure BDA0003059920420000046
步骤2:预测领弹的剩余命中时间
第i枚领弹的剩余命中时间
Figure BDA0003059920420000047
可采用下式预测:
Figure BDA0003059920420000048
式中Ns>2表示导航比。
步骤3:定义领弹的一致性误差变量
领弹的一致性误差变量定义为:
Figure BDA0003059920420000049
步骤4:给出领弹的协同制导律
领弹的协同制导律设计为:
Figure BDA00030599204200000410
式中ki为制导增益。
步骤5:定义第k组第f个从弹变量为:
Figure BDA0003059920420000051
式中,
Figure BDA0003059920420000052
为距离一致性协同误差,
Figure BDA0003059920420000053
为虚拟控制项,
Figure BDA0003059920420000054
为虚拟控制误差。
步骤6:给出从弹的协同制导律
从弹的协同制导律设计为:
Figure BDA0003059920420000055
式中,
Figure BDA0003059920420000056
为辅助变量,wf,1,wf,2和wf,3为正实数。
本发明的有益效果为:设计了一种多领弹多从弹的群组协同制导方法。借助弹群间领弹的通信、以及弹群内的通信,通过建立运动关系数学模型,定义协同变量,分别给出了领弹和从弹的协同制导律,使得所有的导弹最终同时命中目标。
附图说明
图1弹间通信关系示意图。
图2第一组弹群运动轨迹。
图3第二组弹群运动轨迹。
图4第三组弹群运动轨迹。
图5第一组弹群弹目距离变化。
图6第二组弹群弹目距离变化。
图7第三组弹群弹目距离变化。
图8三枚领弹弹目距离变化。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,参照附图1—8对本发明做进一步说明。
所发明一种多领弹多从弹的群组协同制导方法包含以下步骤:
步骤1:建立导弹相对运动关系数学模型
第k枚领弹与目标的相对运动学模型可由如下的方程描述:
Figure BDA0003059920420000061
式中,rk表示第k枚领弹与目标之间的距离,Vk表示恒定的导弹飞行速度,qk为视线角,γk和φk分别为导弹的航迹角和前置角,法向加速度ak垂直于速度方向。
第k组弹群中第f个从弹相对目标的运动关系可以表示为:
Figure BDA0003059920420000071
式中,
Figure BDA0003059920420000072
表示第k组弹群中第f个从弹与目标的距离,
Figure BDA0003059920420000073
Figure BDA0003059920420000074
分别表示第f个从弹的视线角、航迹角、前置角和速度。与领弹不同,从弹可以分别利用切向加速度
Figure BDA0003059920420000075
和法向加速度
Figure BDA0003059920420000076
来调整速度的大小和方向。
我们定义
Figure BDA0003059920420000077
则第f个从弹相对目标的运动关系可转换为:
Figure BDA0003059920420000078
式中,
Figure BDA0003059920420000079
Figure BDA00030599204200000710
为第k组弹群中第f个从弹加速度沿着和垂直于视线方向的分量,
Figure BDA00030599204200000711
Figure BDA00030599204200000712
分别为从弹的速度沿着和垂直于视线方向的分量,可建立以下关系:
Figure BDA0003059920420000081
步骤2:预测领弹的剩余命中时间
第i枚领弹的剩余命中时间
Figure BDA0003059920420000082
可采用下式预测:
Figure BDA0003059920420000083
式中Ns>2表示导航比。
步骤3:定义领弹的一致性误差变量
领弹的一致性误差变量定义为:
Figure BDA0003059920420000084
步骤4:给出领弹的协同制导律
领弹的协同制导律设计为:
Figure BDA0003059920420000085
式中ki为制导增益。
步骤5:定义第k组第f个从弹协同变量为:
Figure BDA0003059920420000086
式中,
Figure BDA0003059920420000087
Figure BDA0003059920420000088
分别为距离一致性协同误差、虚拟控制项和虚拟控制误差。
步骤6:给出从弹的协同制导律
从弹的协同制导律设计为:
Figure BDA0003059920420000091
式中,
Figure BDA0003059920420000092
为辅助变量,wf,1,wf,2和wf,3为正实数。
借助Matlab仿真环境来验证所发明多领弹多从弹群组协同制导方法的有效性。对于实施例,以打击空中静止目标为背景,选择不同位置的3枚领弹,每枚领弹各分属三枚从弹形成一个群组,即构成三个群组。弹间通信拓扑关系如图1所示,其中Li表示第i个领弹,Fi表示第i个从弹,目标位置(9000m,12000m),三枚领弹的速度都选为200m/s,九枚从弹的速度为300m/s,初始横轴坐标位置和纵轴坐标位置分别为(1000,5000),(2000,4000)和(3000,3000),第一组从弹的位置为(3000,4000),(5000,3000)和(8000,2000),第二组从弹的位置为(1800,3800),(2000,3500)和(2200,3800),第三组从弹的位置为(2700,2700),(3000,2500)和(3300,2700),导引律参数设置为ki=1.5,wf,1=0.5,wf,2=3,wf,3=5,Ns=3,轴向过载最大限制为20g,法向过载最大限制为5g。
仿真结果见图2-图8,由平面运动轨迹和弹目距离曲线可知,在所发明的多领弹多从弹群组协同制导方法下,三枚领弹能够保证命中时间的同时性,各从弹的攻击时间也能够与所在弹群领弹的攻击时间保持一致。因此,所有的导弹最终能够实现对目标的同时命中。

Claims (1)

1.一种多领弹多从弹群组协同制导方法,其特征在于,包括如下步骤:
步骤1:建立导弹相对目标的运动关系数学模型
针对二维平面下多导弹对固定目标的打击,设弹间通信拓扑图为连通图,第k枚领弹与目标的相对运动学模型可由如下的方程描述:
Figure FDA0003683600570000011
式中,rk表示第k枚领弹与目标之间的距离,Vk表示恒定的导弹飞行速度,qk为视线角,γk和φk分别为导弹的航迹角和前置角,法向加速度ak垂直于速度方向;
第k组弹群中第f个从弹相对目标的运动关系可以表示为:
Figure FDA0003683600570000012
式中,
Figure FDA0003683600570000013
表示第k组弹群中第f个从弹与目标的距离,
Figure FDA0003683600570000014
Figure FDA0003683600570000015
表示第f个从弹的视线角、航迹角、前置角和速度,与领弹不同,从弹可以分别利用切向加速度
Figure FDA0003683600570000016
和法向加速度
Figure FDA0003683600570000017
来调整速度的大小和方向;
定义
Figure FDA0003683600570000021
那么可得:
Figure FDA0003683600570000022
式中,
Figure FDA0003683600570000023
Figure FDA0003683600570000024
为第k组弹群中第f个从弹加速度沿着和垂直于视线方向的分量,
Figure FDA0003683600570000025
Figure FDA0003683600570000026
分别为从弹的速度沿着和垂直于视线方向的分量,可建立以下关系:
Figure FDA0003683600570000027
步骤2:预测领弹的剩余命中时间
第i枚领弹的剩余命中时间
Figure FDA0003683600570000028
可采用下式预测:
Figure FDA0003683600570000029
式中Ns>2表示导航比;
步骤3:定义领弹的一致性误差变量
领弹的一致性误差变量定义为:
Figure FDA0003683600570000031
步骤4:给出领弹的协同制导律
领弹的协同制导律设计为:
Figure FDA0003683600570000032
式中ki为制导增益;
步骤5:定义第k组第f个从弹变量为:
Figure FDA0003683600570000033
Figure FDA0003683600570000034
式中,
Figure FDA0003683600570000035
为距离一致性协同误差,
Figure FDA0003683600570000036
为虚拟控制项,
Figure FDA0003683600570000037
为虚拟控制误差;
步骤6:给出从弹的协同制导律
从弹的协同制导律设计为:
Figure FDA0003683600570000038
Figure FDA0003683600570000039
Figure FDA00036836005700000310
Figure FDA00036836005700000311
式中,
Figure FDA00036836005700000312
为辅助变量,wf,1,wf,2和wf,3为正实数。
CN202110510192.0A 2021-05-11 2021-05-11 一种多领弹多从弹群组协同制导方法 Active CN113359813B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110510192.0A CN113359813B (zh) 2021-05-11 2021-05-11 一种多领弹多从弹群组协同制导方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110510192.0A CN113359813B (zh) 2021-05-11 2021-05-11 一种多领弹多从弹群组协同制导方法

Publications (2)

Publication Number Publication Date
CN113359813A CN113359813A (zh) 2021-09-07
CN113359813B true CN113359813B (zh) 2022-08-05

Family

ID=77526136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110510192.0A Active CN113359813B (zh) 2021-05-11 2021-05-11 一种多领弹多从弹群组协同制导方法

Country Status (1)

Country Link
CN (1) CN113359813B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113867388A (zh) * 2021-09-30 2021-12-31 北京机电工程研究所 多约束条件下针对机动目标的分布式协同制导方法及装置
CN114020021B (zh) * 2021-11-03 2023-07-25 西北工业大学深圳研究院 一种多导弹分布式协同制导律的设计方法及系统
CN114415722B (zh) * 2022-01-11 2023-09-26 北京航空航天大学 导弹群协同制导方法、电子设备和存储介质
CN114995518B (zh) * 2022-07-27 2022-11-18 西北工业大学 从飞行器gps目标定位失效时的主-从式协同制导方法
CN115016544B (zh) * 2022-08-04 2022-10-28 中国科学院自动化研究所 基于多个飞行器的协同制导方法、装置和电子设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017011407A1 (de) * 2017-12-11 2019-06-13 Mbda Deutschland Gmbh System und verfahren zur personenkoordinierten zielfindung eines lenkflugkörpers
CN110017729A (zh) * 2019-04-18 2019-07-16 西安交通大学 一种带碰撞角约束的多导弹时间协同制导方法
CN110187640A (zh) * 2019-06-29 2019-08-30 东南大学 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法
CN111930142A (zh) * 2020-08-04 2020-11-13 西北工业大学 速度不可控条件下的多弹编队协同控制方法
CN112129292A (zh) * 2020-09-02 2020-12-25 北京航空航天大学 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法
CN112379600A (zh) * 2020-11-20 2021-02-19 西北工业大学 一种基于通信时变延迟的分布式协同制导律构建方法
CN112558634A (zh) * 2020-12-08 2021-03-26 北京航空航天大学 一种多导弹协同制导方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106843265B (zh) * 2016-12-30 2019-12-24 哈尔滨工业大学 有限时间收敛的三维多导弹协同制导方法与系统
CN110412874B (zh) * 2019-07-25 2022-05-13 东南大学 针对机动目标和时延通信的多导弹协同制导律设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017011407A1 (de) * 2017-12-11 2019-06-13 Mbda Deutschland Gmbh System und verfahren zur personenkoordinierten zielfindung eines lenkflugkörpers
CN110017729A (zh) * 2019-04-18 2019-07-16 西安交通大学 一种带碰撞角约束的多导弹时间协同制导方法
CN110187640A (zh) * 2019-06-29 2019-08-30 东南大学 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法
CN111930142A (zh) * 2020-08-04 2020-11-13 西北工业大学 速度不可控条件下的多弹编队协同控制方法
CN112129292A (zh) * 2020-09-02 2020-12-25 北京航空航天大学 考虑攻击时间和攻击角度约束的异构飞行器协同制导方法
CN112379600A (zh) * 2020-11-20 2021-02-19 西北工业大学 一种基于通信时变延迟的分布式协同制导律构建方法
CN112558634A (zh) * 2020-12-08 2021-03-26 北京航空航天大学 一种多导弹协同制导方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Distributed Anti-Disturbance Control Strategy for Multi-Aircraft Systems With Time-Varying Topology;Li Guofei;《IEEE Access》;20190221;第7卷;第26278-26292页 *
主-从多飞行器三维分布式协同制导方法;李国飞;《航空学报》;20210104;第42卷(第11期);第1-10页 *

Also Published As

Publication number Publication date
CN113359813A (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
CN113359813B (zh) 一种多领弹多从弹群组协同制导方法
CN105910495B (zh) 基于性能指标的面向效能的导弹武器系统设计方法
CN111707148B (zh) 一种时变时延条件下多导弹协同制导方法及系统
CN111649624B (zh) 一种空间微型精确制导武器控制方法
CN112035780B (zh) 一种导弹末制导阶段杀伤效果计算方法
CN114020021B (zh) 一种多导弹分布式协同制导律的设计方法及系统
CN109376364B (zh) 基于扩展卡尔曼滤波的高速旋转弹气动参数辨识方法
CN116360489A (zh) 一种基于成型制导律的协同中制导律设计方法
CN114819055A (zh) 一种基于lstm网络的导弹拦截点预测方法
CN114153143A (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
RU2602162C2 (ru) Способ стрельбы реактивными снарядами реактивной системы залпового огня в условиях контрбатарейной борьбы
CN112464451A (zh) 基于作战仿真系统的防空导弹武器命中概率修正方法
CN111832193A (zh) 一种超视距空战中空中目标的威胁评估方法
KR101838679B1 (ko) 미사일 마스킹 영역 산출방법 및 이를 이용한 공대공 미사일 운용방법
Li et al. Simulation Research on Sec Search Algorithm for Medium Range Air-to-Air Missile Attack Area
Tao et al. Modeling study on the operational process of photoelectric-guided munitions
Dautov et al. RECONNAISSANCE-STRIKE AND RECONNAISSANCE-FIRE COMPLEXES AS A WAY TO IMPROVE THE EFFECTIVENESS OF THE ARMED FORCES
Zhao et al. Simulation design of fuze warhead system of air defense missile at very low altitude
Shi et al. The Research on Damage Effect of Target under the Action of Fragmentation Warhead
Gkritzapis et al. Computational atmospheric trajectory simulation analysis of spin-stabilised projectiles and small bullets
Jiang et al. An Analysis Method for Hit Accuracy of Gun-Launched Reconnaissance and Strike Integrated Unmanned Aerial Vehicle
Yongjie et al. Ballistic and Aerodynamic Characteristics Simulation for Trajectory Correction Projectile
Nan et al. Influencing factor analysis of interception probability and classification-regression neural network based estimation
Gao et al. Algorithm and Simulation of Airborne Fire Control Aiming Angle Correction
Zhang et al. Correction Control Time Optimized Design of Two-Dimensional Trajectory Correction Projectile Based on PPNGA

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