CN109543135A - 一种基于信息最大化的多auv协同同时导引方法 - Google Patents

一种基于信息最大化的多auv协同同时导引方法 Download PDF

Info

Publication number
CN109543135A
CN109543135A CN201811333769.XA CN201811333769A CN109543135A CN 109543135 A CN109543135 A CN 109543135A CN 201811333769 A CN201811333769 A CN 201811333769A CN 109543135 A CN109543135 A CN 109543135A
Authority
CN
China
Prior art keywords
auv
information
guidance
collaboration
maximumlly
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
CN201811333769.XA
Other languages
English (en)
Other versions
CN109543135B (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
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201811333769.XA priority Critical patent/CN109543135B/zh
Publication of CN109543135A publication Critical patent/CN109543135A/zh
Application granted granted Critical
Publication of CN109543135B publication Critical patent/CN109543135B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提出一种基于信息最大化的多AUV协同同时导引方法，针对多AUV协同导引中的时间约束和观测信息最大化问题，同时考虑了AUV二阶非完整力学特性，提出了基于最优控制的协同导引律，其中性能指标函数由一步次优费希尔信息矩阵(FIM)和AUV与目标距离方差组成，实现了观测信息最大下的同时击中目标的任务，具有良好的导引性能。

Description

一种基于信息最大化的多AUV协同同时导引方法
技术领域
本发明涉及多AUV协同导引技术领域，具体为一种基于信息最大化的多AUV同时攻击同一目标的协同导引方法。
背景技术
现代反导系统的飞速发展，对主动攻击武器(导弹、鱼雷等)的突防能力提出了更高的挑战，传统的单一AUV打击已经很难满足现代战争需求。“多对一”形式的饱和攻击是一种能提高导弹或鱼雷击中目标概率的有效手段。因此，为了突破层层防御，有必要研究“多对一”形式的饱和同时攻击。针对协同导引的理论研究具有深刻的军事意义，逐渐吸引了国内外学者的研究兴趣。
由于AUV的发射位置可能不同，并且各个武器的发射时间和航行速度不可能做到绝对同步，若使用传统制导律控制AUV攻击目标会导致武器的航行时间有长有短，不能保证多AUV能够同时击中目标，因此，有必要研究一种对AUV的攻击时间进行控制的制导律。
同时，被攻击目标状态的估计精度越高，多AUV击中目标的概率越大。通常影响目标估计精度的因素不仅与所搭载传感器的量测精度有关，而且还与多AUV在空间的配置方式有关。因此，针对具有最优的观测空间配置的研究逐渐成为协同导引研究的重要方向。
另外，在实际工程中，洋流等外部随机干扰会使AUV产生侧向速度，其会对导引效果造成一定的影响。因此，针对此类二阶非完整力学特性的多AUV协同导引律的设计，是一个很有实际意义的研究问题。
综上所述，针对信息最大化的多AUV的协同同时导引，还存在以下三个问题：1、外界干扰会使AUV产生侧向速度，其会对导引效果造成一定影响；2、如何设计具有实时性、且观测信息最大的多AUV协同导引方法，具有重要的工程意义。3、由于实际水下环境存在通信能力弱等问题，很难保证多AUV间的集中式通信，因此，设计分布式多AUV协同导引具有重要的现实意义。
发明内容
发明针对多AUV协同导引中的时间约束和观测信息最大化问题，同时考虑了AUV二阶非完整力学特性，提出了一种基于信息最大化的多AUV协同同时导引方法，实现了观测信息最大下的同时击中目标的任务，具有良好的导引性能。
本发明的技术方案为：
所述一种基于信息最大化的多AUV协同同时导引方法，其特征在于：包括以下步骤：
步骤1：建立多AUV水平面运动学方程：
其中(xi,yi)为世界坐标系下第i个AUV的位置坐标，ui，vi和ψi分别为第i个AUV在体坐标系下的前向速度、侧向速度和航向角；第i个AUV与目标的相对方位角为
其中(xT,yT)为被攻击目标在世界坐标系下的位置坐标；
步骤2：对于步骤1的AUV水平面运动学方程，建立信息最大化的多AUV协同同时导引律为：
ψi＝atan2(-λyiui+λxivi,-λxiui-λyivi)
其中：
Ni为能与第i个AUV互相通信的AUV集合，为常数，其取值范围为
步骤3：利用步骤2得到的AUV航向角进行多AUV协同引导。
有益效果
本发明通过优化由AUV与目标距离方差和一步次优FIM组成的性能指标函数，来得到实时的分布式协同导引律；将导引律作用于各AUV的运动学方程，使其达到预期的导引效果。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
附图1为多AUV协同导引原理框图。
附图2为考虑侧向速度下的导引弹道图。
附图3为不考虑侧向速度下的导引弹道图。
附图4为两种导引律下的一步次优FIM变化曲线图。
附图5为两种导引律下的AUV-目标距离变化曲线。
具体实施方式
下面详细描述本发明的实施例，所述实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
本实施例中以四个AUV为研究对象，其中，各AUV的初始位置(x1(0),y1(0))＝(100,100)(x2(0),y2(0))＝(120,-80),(x3(0),y3(0))＝(-80,-140)，(x4(0),y4(0))＝(-160,100)，初始航向角均为0°，u1＝u2＝u3＝u4＝1m/s，v1＝v2＝-0.2m/s,v3＝v4＝0.2m/s。
步骤1：建立多AUV水平面运动学方程：
其中(xi,yi)为世界坐标系下第i个AUV的位置坐标，ui，vi和ψi分别为第i个AUV在体坐标系下的前向速度、侧向速度和航向角；第i个AUV与目标的相对方位角为
其中(xT,yT)为被攻击目标在世界坐标系下的位置坐标。
步骤2：针对步骤1的运动模型，采用一步最优控制方法，得到协同导引律，其中待优化的性能指标函数由FIM和距离之差组成，其表达式为：
其中，
同时，FIM的行列式为
进一步地，性能指标函数可重新写为：
其中，fi＝[uicosψi-visinψi,uisinψi+vicosψi]T，Xi＝[xi,yi]T，拉格朗日乘子向量λi＝[λxi,λyi]T。
构造哈密尔顿函数为：
Hi＝λi Tfi＝λxi(uicosψ-visinψ)+λyi(uisinψ+vicosψ)
其控制方程为：
因此，协同导引律可表示为：
ψi＝atan2(-λyiui+λxivi,-λxiui-λyivi)
为了保证拉格朗日乘子向量λi应当小于零。
由于因此，λi为常数。
其中，λi应当满足其中λi为：
其中，
以上设计的多AUV分布式协同导引律，当选择合适的参数ε时，可保证各AUV以“尽量分散”的形式同时击中同一目标。“尽量分散”的目的是使得多AUV能以一个好的空间位置以获取精确的目标估计精度。同时攻击的目的是增强AUV击中目标的概率。同时，上述考虑了AUV的二阶非完整欠驱动特性，具有一定工程意义。
针对本实施中的AUV运动模型，采用上述协同导引律进行导引，并与不考虑侧向速度的协同导引律进行了仿真对比，其导引律为ψi＝atan2(-λyi,-λxi)。从仿真图2与图3可以看出，四个AUV能同时击中同一目标，且相互之间以近似90°的角度击中同一目标，即以近似最优的观测空间配置形式攻击，仿真验证了该算法的有效性。附图4表明了考虑侧向速度的协同导引律的脱靶量精度高于不考虑侧向速度协同导引律的脱靶量精度。附图5表明了考虑侧向速度导引律具有比不考虑协同导引律更优越的观测空间配置，即前者拥有更准确的目标估计精度的可能性。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (1)

1.一种基于信息最大化的多AUV协同同时导引方法，其特征在于：包括以下步骤：
步骤1：建立多AUV水平面运动学方程：
其中(xi,yi)为世界坐标系下第i个AUV的位置坐标，ui，vi和ψi分别为第i个AUV在体坐标系下的前向速度、侧向速度和航向角；第i个AUV与目标的相对方位角为
其中(xT,yT)为被攻击目标在世界坐标系下的位置坐标；
步骤2：对于步骤1的AUV水平面运动学方程，建立信息最大化的多AUV协同同时导引律为：
ψi＝atan2(-λyiui+λxivi,-λxiui-λyivi)
其中：
Ni为能与第i个AUV互相通信的AUV集合，为常数，其取值范围为
步骤3：利用步骤2得到的AUV航向角进行多AUV协同引导。
CN201811333769.XA 2018-11-09 2018-11-09 一种基于信息最大化的多auv协同同时导引方法 Active CN109543135B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811333769.XA CN109543135B (zh) 2018-11-09 2018-11-09 一种基于信息最大化的多auv协同同时导引方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811333769.XA CN109543135B (zh) 2018-11-09 2018-11-09 一种基于信息最大化的多auv协同同时导引方法

Publications (2)

Publication Number Publication Date
CN109543135A true CN109543135A (zh) 2019-03-29
CN109543135B CN109543135B (zh) 2021-01-05

Family

ID=65846678

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811333769.XA Active CN109543135B (zh) 2018-11-09 2018-11-09 一种基于信息最大化的多auv协同同时导引方法

Country Status (1)

Country Link
CN (1) CN109543135B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991571A (zh) * 2019-04-17 2019-07-09 电子科技大学 远距离分布式无源测向定位的布阵方法
CN110764533A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 一种多水下机器人协同目标搜索方法
CN111026137A (zh) * 2019-12-26 2020-04-17 西北工业大学 一种攻击角约束下同时攻击目标的三维分布式协同导引方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8436283B1 (en) * 2008-07-11 2013-05-07 Davidson Technologies Inc. System and method for guiding and controlling a missile using high order sliding mode control
CN108490968A (zh) * 2018-03-17 2018-09-04 西北工业大学 一种基于反馈融合结构的自主水下航行器控制行为融合方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8436283B1 (en) * 2008-07-11 2013-05-07 Davidson Technologies Inc. System and method for guiding and controlling a missile using high order sliding mode control
CN108490968A (zh) * 2018-03-17 2018-09-04 西北工业大学 一种基于反馈融合结构的自主水下航行器控制行为融合方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BONG SEOK PARK: "Neural network-based output feedback control for reference tracking of underactuated surface vessels", 《AUTOMATICA》 *
张友安: "带攻击时间约束的导引律综述", 《海军航空工程学院学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991571A (zh) * 2019-04-17 2019-07-09 电子科技大学 远距离分布式无源测向定位的布阵方法
CN109991571B (zh) * 2019-04-17 2022-11-18 电子科技大学 远距离分布式无源测向定位的布阵方法
CN110764533A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 一种多水下机器人协同目标搜索方法
CN111026137A (zh) * 2019-12-26 2020-04-17 西北工业大学 一种攻击角约束下同时攻击目标的三维分布式协同导引方法
CN111026137B (zh) * 2019-12-26 2022-10-28 西北工业大学 一种攻击角约束下同时攻击目标的三维分布式协同导引方法

Also Published As

Publication number Publication date
CN109543135B (zh) 2021-01-05

Similar Documents

Publication Publication Date Title
CN108168381B (zh) 一种多枚导弹协同作战的控制方法
CN109597423B (zh) 一种基于参考视线角信号的多约束末制导律的设计方法
CN111551080B (zh) 一种攻击时间控制的制导方法
CN108362171B (zh) 一种具有攻击时间和攻击角度约束的制导策略
CN107179021B (zh) 一种驾束制导体制下多弹协同高精度制导控制方法
CN109543135A (zh) 一种基于信息最大化的多auv协同同时导引方法
CN110412874B (zh) 针对机动目标和时延通信的多导弹协同制导律设计方法
CN104881043A (zh) 一种针对多动态目标的多无人机智能协同察打方法
CN110686564B (zh) 一种红外半捷联导引头制导方法及系统
CN110345814B (zh) 一种不依赖自身导引头测量信息的末制导算法
CN114153143B (zh) 一种导弹非奇异固定时间滑模制导律的设计方法
CN115951695A (zh) 空战模拟环境中基于三方博弈的动态战术控制域解算方法
CN114415722B (zh) 导弹群协同制导方法、电子设备和存储介质
CN116663263A (zh) 一种基于滑模面的飞行器末端落速落角约束制导方法及系统
CN113064448B (zh) 一种飞行器的制导方法及系统
CN109827478A (zh) 一种带落角约束与过载约束的制导方法
CN112379599B (zh) 一种针对空中机动目标的三维协同中制导方法
CN116227343A (zh) 拦截机动目标满足视场角约束的末角攻击制导律设计方法
CN111026137B (zh) 一种攻击角约束下同时攻击目标的三维分布式协同导引方法
CN113110576B (zh) 一种可连续切换的自适应固定时间收敛协同制导方法
CN114819055A (zh) 一种基于lstm网络的导弹拦截点预测方法
CN114415724B (zh) 容许速度衰减的三维多导弹协同比例导引律设计方法
CN115289908B (zh) 一种遥控指令防空导弹引入段导引方法及装置
CN117891271B (zh) 考虑时间和角度约束的高速飞行器三维协同制导方法
CN114923365B (zh) 一种类tvm制导方法及装置

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