CN112379599A - 一种针对空中机动目标的三维协同中制导方法 - Google Patents
一种针对空中机动目标的三维协同中制导方法 Download PDFInfo
- Publication number
- CN112379599A CN112379599A CN202011314913.2A CN202011314913A CN112379599A CN 112379599 A CN112379599 A CN 112379599A CN 202011314913 A CN202011314913 A CN 202011314913A CN 112379599 A CN112379599 A CN 112379599A
- Authority
- CN
- China
- Prior art keywords
- interceptor
- guidance
- target
- interceptors
- ith
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000001133 acceleration Effects 0.000 claims abstract description 14
- 230000033001 locomotion Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002195 synergetic effect Effects 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 description 17
- 238000005755 formation reaction Methods 0.000 description 17
- 238000004891 communication Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- PHTXVQQRWJXYPP-UHFFFAOYSA-N ethyltrifluoromethylaminoindane Chemical compound C1=C(C(F)(F)F)C=C2CC(NCC)CC2=C1 PHTXVQQRWJXYPP-UHFFFAOYSA-N 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000270295 Serpentes Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种针对空中机动目标的三维协同中制导方法,将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。本发明结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。
Description
技术领域
本发明属于制导技术领域,具体涉及一种三维协同中制导方法。
背景技术
随着反拦截器技术日趋完善,空中机动目标一方面可以通过释放诱饵、蛇形机动或者其他对抗手段对拦截器进行干扰和对抗,使得现有拦截器的拦截效率大大降低;另一方面,空中目标的性能随着人工智能、无人机技术的快速发展而大幅提升,单个拦截器的探测能力和机动能力都难以满足高精度的打击要求,需要采用协同包围攻击技术。因此,多拦截器协同突防或拦截的思路应运而生。多拦截器协同编队协同作战方式是指多个拦截器通过信息共享、功能互补,不仅能够实现对目标从多个方向进行饱和攻击,提高打击效果和命中概率,还能实现单个拦截器不易完成的协同突防等任务。针对军事价值较高的目标,多拦截器协同编队作战方式是实施突发性和高密集的精确打击的重要手段,是未来突发性地区战争的重要打击形式。
对于中远程拦截器来说,要想实现多拦截器编队进行饱和打击需要经历四个飞行阶段,分别为编队集结、编队飞行、编队中末交班和编队末制导阶段。对于中远程拦截器编队来说,要想实现协同打击目标的任务,需要从拦截器发射阶段就开始进行协同,经过协同中制导和中末制导交班,最后进行协同末制导,对目标完成协同聚能饱和打击。若单研究其中一个飞行阶段的协同制导问题是不可行的,从而无法完成协同打击目标的任务,例如:当有两枚拦截器同时指向同一目标,一枚拦截器已经进入末制导阶段,距离目标低于15KM;另一枚拦截器在中制导阶段,距离目标有30KM以上,此时,这两枚拦截器很难完成协同攻击,不符合进行协同制导的基础条件(离目标的距离误差在一定范围内)。所以,对于中远程多拦截器编队协同制导,多拦截器编队整个飞行阶段的协同都至关重要,每一阶段的协同完成效果直接影响下一阶段的协同任务完成度。因此,对于多拦截器编队协同中制导段,需要实现可以在同一时刻或者近似同一时刻到达中末制导交班区,为多拦截器编队进行协同末制导做好准备。
目前,主要的理论研究与工程方法大都集中在协同末制导段,而末制导阶段的协同与中制导阶段的协同有本质上的区别。末制导阶段的目标信息获取由拦截器本身的导引头完成,目标信息精度较高,而且多拦截器编队队形和通信拓扑结构基本不发生变化,拦截器之间距离较小,通信连接情况良好,很容易达到强连通状态。在中制导阶段,目标信息由外部传感器探测和估计所得,拦截器之间距离较大,通信连接情况较差。因此,拦截器在协同中制导段所采用的制导律与在协同末制导段采用的制导律有很大区别,主要体现在:性能指标的要求、约束条件、制导目的以及获取目标信息的方式和目标信息的质量等。
同时,因为末制导时间太短,拦截器可以调整的飞行距离和飞行速度等非常有限,必须要有非常好的协同中制导律,确保在中末制导交班时就能达到时间和空间初步的协同,即给末制导段提供非常好的协同初始条件,才能完成末制导段的协同,从而实现中远程拦截器的时间和空间协同。因此,协同中制导方法是中远程拦截器实现协同的前提和基础。
目前,对协同中制导的理论研究成果鲜有报道,已经严重影响了多拦截器编队协同控制系统的研制。因此,提出一套适合工程应用的多拦截器协同中制导方法就显得非常迫切。
发明内容
为了克服现有技术的不足,本发明提供了一种针对空中机动目标的三维协同中制导方法,将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。本发明结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:建立多拦截器和目标的运动模型:
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi:
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
步骤3:计算虚拟碰撞点;
设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
步骤4:计算三维惯性坐标系下拦截器加速度为:
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
步骤5:设置协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数,从而三维协同中制导律满足一致性要求。
本发明提出的一种针对空中机动目标的三维协同中制导方法,能够确保中远程拦截器编队几乎同时到达中末交接班区域,为拦截器编队的协同末制导阶段提供良好的初始条件。本发明方法结构简单,物理含义清晰,易于工程人员理解,同时具有协同项系数便于调节,制导律鲁棒性强等优点,易于工程实现。
附图说明
图1是本发明多个拦截器与目标在惯性系下的三维示意图。
图2是本发明实施例各拦截器与目标的运动轨迹示意图。
图3是本发明实施例各拦截器的剩余时间曲线示意图。
图4是本发明实施例各拦截器的剩余距离曲线示意图。
图5是本发明实施例拦截器1的加速度指令曲线示意图。
图6是本发明实施例拦截器2的加速度指令曲线示意图。
图7是本发明实施例拦截器3的加速度指令曲线示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
如图1所示,一种针对空中机动目标的三维协同中制导方法,包括如下步骤:
步骤1:建立多拦截器和目标的运动模型:
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi:
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
步骤3:计算虚拟碰撞点;
为了提高协同末制导对机动目标的拦截成功率,设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
步骤4:计算三维惯性坐标系下拦截器加速度为:
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
步骤5:为了使协同中制导律达到协同一致性,协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数。
以X轴为例进行分析,Y轴和Z轴证明类似,具体证明如下:
接下来分析确保不等式ηi(t+△t)<ηi(t)的具体条件:
由于
由(10)和(11)式可得:
进一步可将axi改写为axi=abi-Pxiaci,其中:
由不等式ηi(t+△t)<ηi(t),△t>0可推出:
接下来分析不等式(14)右边项是小于0的,从而Pi选择任意大于0的实数可满足ηi(t+△t)<ηi(t)。
由于一般情况下,VTx>0且xT>xi,可得:
Bi=MiVTx<0 (17)
(Ei-(xT-xi)axi)Di<2 (19)
由于abi是axi的一部分,可以近似处理为axi≈abi,因而获得:
(Ei-(xT-xmi)abi)Di<2 (20)
Wi=Bi(2-(Ei-(xT-xmi)abi)Di)<0 (21)
从(16)、(18)和(21)式可获得
因此(Pxi,Pyi,Pzi)可以选取任意大于0的实数,确保ηi(t+△t)<ηi(t)成立,从而三维协同中制导律满足一致性要求。
本发明将虚拟碰撞点作为拦截器中制导阶段的目标点,同时考虑利用在中制导阶段通过外部传感器直接获得的目标的位置信息、速度信息和加速度信息设计协同中制导方法,并通过添加协同调节项确保所有拦截器的在中制导阶段结束时,几乎同时到达中末交接班区域。
具体实施例:
本实施例通过用3个拦截器同时攻击机动目标来验证所提出的三维协同中制导导引律的有效性。实施例中3个拦截器和目标的仿真初始条件:
拦截器1:初始位置(1000m,9000m,-2300m),初始速度(1100m/s,100m/s,100m/s)
拦截器2:初始位置(500m,8500m,-500m),初始速度(1000m/s,100m/s,90m/s)
拦截器3:初始位置(2000m,7600m,3000m),初始速度(950m/s,60m/s,110m/s)
目标:初始位置(100000m,8000m,1000m),初始速度(150m/s,10m/s,15m/s)
图2为在三维平面下多拦截器与目标的轨迹,其中黑色实线为拦截器在协同中制导的轨迹,黑色点实线为拦截器为协同末制导阶段的轨迹。图3为所有拦截器的飞行剩余时间变化曲线。图4为所有拦截器与目标的剩余距离的变化曲线。图5到图7分别表示各拦截器的加速度指令变化曲线。从图2可以看出,本发明方法能确保所有拦截器几乎同时到达中末交班区域,为协同末制导提供了良好的初始条件,也为最终拦截机动目标提供的拦截环境。从图3到图4可以看出,本发明方法能确保所有拦截器的剩余时间与剩余距离在协同中制导阶段趋近收敛,达到一致的状态。从图5到图7可以看出本发明方法的加速度指令平滑,且数值较小,易于工程实现。
Claims (1)
1.一种针对空中机动目标的三维协同中制导方法,其特征在于,包括以下步骤:
步骤1:建立多拦截器和目标的运动模型:
其中,(xi,yi,zi)和(xT,yT,zT)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的位置信息,(Vxi,Vyi,Vzi)和(VTx,VTy,VTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的速度信息,(axi,ayi,azi)和(aTx,aTy,aTz)分别代表第i个拦截器和目标在惯性坐标系三个坐标方向上的加速度信息,i=1,2,...,N,N代表拦截器的个数;
步骤2:计算多拦截器到达目标的剩余时间tgoi:
其中ri代表第i个拦截器到达目标的剩余距离,计算公式如下:
步骤3:计算虚拟碰撞点;
设置虚拟碰撞点作为中制导阶段拦截器的目标点,虚拟碰撞点的计算如下:
步骤4:计算三维惯性坐标系下拦截器加速度为:
其中,tf代表中制导结束时间,(VTx(tf),VTy(tf),VTz(tf))和(Vxi(tf),Vyi(tf),Vzi(tf))分别代表目标和拦截器在惯性坐标系下三个方向的末速度,(Pxi,Pyi,Pzi)代表协同项调节系数,代表所有拦截器虚拟碰撞点坐标的平均值,计算如下:
步骤5:设置协同项调节系数(Pxi,Pyi,Pzi)为大于0的任意正实数,从而三维协同中制导律满足一致性要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011314913.2A CN112379599B (zh) | 2020-11-20 | 2020-11-20 | 一种针对空中机动目标的三维协同中制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011314913.2A CN112379599B (zh) | 2020-11-20 | 2020-11-20 | 一种针对空中机动目标的三维协同中制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112379599A true CN112379599A (zh) | 2021-02-19 |
CN112379599B CN112379599B (zh) | 2022-06-14 |
Family
ID=74589015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011314913.2A Active CN112379599B (zh) | 2020-11-20 | 2020-11-20 | 一种针对空中机动目标的三维协同中制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112379599B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113433428A (zh) * | 2021-05-10 | 2021-09-24 | 广西电网有限责任公司玉林供电局 | 一种多元混合绝缘气体的协同效应分析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106843265A (zh) * | 2016-12-30 | 2017-06-13 | 哈尔滨工业大学 | 有限时间收敛的三维多导弹协同制导方法与系统 |
CN108168381A (zh) * | 2018-01-04 | 2018-06-15 | 北京理工大学 | 一种多枚导弹协同作战的控制方法 |
CN110187640A (zh) * | 2019-06-29 | 2019-08-30 | 东南大学 | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 |
CN110412874A (zh) * | 2019-07-25 | 2019-11-05 | 东南大学 | 针对机动目标和时延通信的多导弹协同制导律设计方法 |
CN111026137A (zh) * | 2019-12-26 | 2020-04-17 | 西北工业大学 | 一种攻击角约束下同时攻击目标的三维分布式协同导引方法 |
CN111506101A (zh) * | 2019-10-21 | 2020-08-07 | 北京理工大学 | 基于通信网络拓扑结构的飞行器协同制导控制方法及系统 |
-
2020
- 2020-11-20 CN CN202011314913.2A patent/CN112379599B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106843265A (zh) * | 2016-12-30 | 2017-06-13 | 哈尔滨工业大学 | 有限时间收敛的三维多导弹协同制导方法与系统 |
CN108168381A (zh) * | 2018-01-04 | 2018-06-15 | 北京理工大学 | 一种多枚导弹协同作战的控制方法 |
CN110187640A (zh) * | 2019-06-29 | 2019-08-30 | 东南大学 | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 |
CN110412874A (zh) * | 2019-07-25 | 2019-11-05 | 东南大学 | 针对机动目标和时延通信的多导弹协同制导律设计方法 |
CN111506101A (zh) * | 2019-10-21 | 2020-08-07 | 北京理工大学 | 基于通信网络拓扑结构的飞行器协同制导控制方法及系统 |
CN111026137A (zh) * | 2019-12-26 | 2020-04-17 | 西北工业大学 | 一种攻击角约束下同时攻击目标的三维分布式协同导引方法 |
Non-Patent Citations (2)
Title |
---|
LIU ZHAOHUI 等: "Cooperative Guidance Law Design on Simultaneous Attack for Multiple Missiles under Time-Delayed Communication Topologies", 《 2019 IEEE SYMPOSIUM SERIES ON COMPUTATIONAL INTELLIGENCE (SSCI)》 * |
ZIHAO WU,等: "Three-Dimensional Cooperative Mid-Course Guidance Law Against the Maneuvering Target", 《IEEE ACCESS》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113433428A (zh) * | 2021-05-10 | 2021-09-24 | 广西电网有限责任公司玉林供电局 | 一种多元混合绝缘气体的协同效应分析方法 |
CN113433428B (zh) * | 2021-05-10 | 2023-09-01 | 广西电网有限责任公司玉林供电局 | 一种多元混合绝缘气体的协同效应分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112379599B (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111580556B (zh) | 一种时空约束下多无人机协同路径规划与制导方法 | |
CN103884237B (zh) | 基于目标概率分布信息的多对一协同制导方法 | |
CN106529073B (zh) | 基于拦截几何的高超声速目标拦截弹交接班条件分析方法 | |
CN110412874B (zh) | 针对机动目标和时延通信的多导弹协同制导律设计方法 | |
CN108534614B (zh) | 一种三维全向实时预测制导方法 | |
CN110764523B (zh) | 基于抗饱和平滑变换的比例积分前置导引攻击目标的方法 | |
CN114138000B (zh) | 考虑全捷联导引头视场约束的弹群协同制导控制一体化设计方法 | |
CN110345814B (zh) | 一种不依赖自身导引头测量信息的末制导算法 | |
CN112648886B (zh) | 一种组合制导目标拦截方法及系统 | |
CN114153143B (zh) | 一种导弹非奇异固定时间滑模制导律的设计方法 | |
CN112379600A (zh) | 一种基于通信时变延迟的分布式协同制导律构建方法 | |
CN115857538A (zh) | 三维空间下满足落角约束的多飞行器协同制导方法 | |
CN112379599B (zh) | 一种针对空中机动目标的三维协同中制导方法 | |
CN112033235A (zh) | 一种采用黄金分割搜索的空战三体对抗最优预测制导方法 | |
CN109543135B (zh) | 一种基于信息最大化的多auv协同同时导引方法 | |
CN114489101A (zh) | 一种无人机的末端制导控制方法及系统 | |
CN117932894A (zh) | 一种多模导引头信息融合目标状态估计方法 | |
CN116774714A (zh) | 一种基于事件触发机制的多约束协同制导方法 | |
CN116225048A (zh) | 基于落角和时间约束的巡飞弹末端协同打击方法 | |
CN115542939A (zh) | 空空导弹分布式协同中制导律分析方法及其导引系统 | |
Tang et al. | Research on Target State Estimation and Terminal Guidance Algorithm in the Process of Multi-UAV Cooperative Attack | |
CN115617063A (zh) | 一种带落角约束的飞行器制导控制装置及方法 | |
CN113759954A (zh) | 一种针对机动目标的复合制导方法 | |
CN113138604A (zh) | 一种反低小慢目标的无人机拦截方法 | |
CN115755955B (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 |