CN106843265B - 有限时间收敛的三维多导弹协同制导方法与系统 - Google Patents
有限时间收敛的三维多导弹协同制导方法与系统 Download PDFInfo
- Publication number
- CN106843265B CN106843265B CN201611256552.4A CN201611256552A CN106843265B CN 106843265 B CN106843265 B CN 106843265B CN 201611256552 A CN201611256552 A CN 201611256552A CN 106843265 B CN106843265 B CN 106843265B
- Authority
- CN
- China
- Prior art keywords
- sight
- line
- missile
- acceleration
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000001133 acceleration Effects 0.000 claims abstract description 62
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 8
- 238000013461 design Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 101100491149 Caenorhabditis elegans lem-3 gene Proteins 0.000 description 1
- 244000145845 chattering Species 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 231100000225 lethality Toxicity 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/107—Simultaneous control of position or course in three dimensions specially adapted for missiles
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G3/00—Aiming or laying means
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (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
本发明公开了一种有限时间收敛的三维多导弹协同制导方法,包括:针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型;利用自适应控制、积分滑模控制及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标;通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛。本发明能够在不了解目标任何机动信息的情况下,使所有导弹以期望角度同时击中目标。同时避免抖振的发生。
Description
技术领域
本发明涉及制导技术领域,尤其涉及一种有限时间收敛的三维多导弹协同制导方法与系统。
背景技术
多导弹协同作战是通过多枚导弹相互配合、协作、共同执行作战任务来实现的。随着反导技术的快速发展,综合一体化的现代防御体系已逐步形成,单枚导弹突防愈发困难,而多枚导弹通过信息共享与功能互补,不仅极大地提高了导弹的杀伤力和突防能力,而且能够完成单枚导弹无法完成的任务,因此研究多导弹协同制导律具有非常重要的实际意义。
现有的协同制导方法针对的都是静止或低速的目标,而且需要人为地事先指定攻击时间。另外,在制导过程中各导弹之间没有动态的信息交互,不能实现真正意义上的导弹协同。
同时,由于在多导弹协同拦截过程中暗含着可能发生的碰撞,研究人员提出一种基于符号函数或滑模控制理论的有限时间一致性协议。但是,上述协议是不连续的,在实际应用中会产生抖振现象,严重影响制导性能。
发明内容
针对现有技术的上述缺陷,本发明提供一种有限时间收敛的三维多导弹协同制导方法与系统,能够在不了解目标任何机动信息的情况下,使所有导弹以期望角度同时击中目标。同时,本发明采用连续的有限时间一致性协议进行控制律设计,消除了抖振的发生。
根据本发明的一个方面,提供了一种有限时间收敛的三维多导弹协同制导方法,包括步骤:
S1.针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型;
S2.利用自适应控制、积分滑模控制、及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标;
S3.通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛。
优选地,所述视线方向的加速度制导律、及所述视线法向的加速度制导律均是连续的。
优选地,所述方法还包括:将每枚导弹的剩余时间信息在导弹间进行交互。
优选地,所述多导弹协同制导模型具体为:
其中,x1=r,x3=[qε-qεd,qβ-qβd]T,qε为视线倾角,qβ为视线偏角,aMr为导弹加速度在视线坐标系x轴的分量,aTr为目标加速度在视线坐标系x轴的分量,r为导弹与目标的相对距离,qεd为期望的视线倾角,qβd为期望的视线偏角,aMε为导弹加速度在视线坐标系y轴的分量,aMβ为导弹加速度在视线坐标系z轴的分量,aTε为目标加速度在视线坐标系y轴的分量,aTβ为目标加速度在视线坐标系y轴的分量。
优选地,在建立多导弹协同制导模型之后,设计制导律之前,所述方法还包括:将导弹的剩余时间信息引入多导弹协同制导模型,获得如下式表示的制导模型;
其中,i表示第i枚导弹,tgoi为第i枚导弹的剩余时间,
优选地,所述视线方向的加速度制导律由下式表示:
l2=2l1a+b+4a2 式20
其中,为标准一致性协议,l1、l2、l3为自适应增益,l3>0,tfi为导弹拦截目标的预测时刻,tfi(0)为tfi的初始值,w1,γ1,a,b>0,w2,γ2,l1*,l2*>0。
优选地,所述视线法向的加速度制导律由下式表示:
其中,α=diag(α1,α2),β=diag(β1,β2),α1、α2、β1、β2>0,1<λ2<2,λ1>λ2,k1=diag(k11,k12),k2=diag(k21,k22),k11,k12,k21,k22>0,0<μ<1,s为选取的非奇异快速终端滑模面,ε>0,为k的估计值,Δqi为未知正数,k3>0,k4=(k41,k42),k41,k42>0。
优选地,在建立所述视线方向的加速度制导律、及所述视线法向的加速度制导律之后,所述方法还包括:利用李雅普诺夫理论证明系统的有限时间稳定性。
根据本发明的另一方面,提供一种有限时间收敛的三维多导弹协同制导系统,包括:
模型建立单元,用于针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型;
第一制导律建立单元,用于利用自适应控制、积分滑模控制、及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标;
第二制导律建立单元,用于通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛。
优选地,所述视线方向的加速度制导律、及所述视线法向的加速度制导律均是连续的。
在本发明的技术方案中,针对三维空间中多导弹协同拦截一个机动目标的问题,提出了一种具有攻击角约束的有限时间协同制导律,该协同制导律是连续的,并且不需要知道目标机动的任何信息。具体地:首先,基于三维空间内的导弹-目标相对运动方程,建立带有攻击角约束的多导弹协同制导模型。其次,将协同制导律的设计过程分离为两个部分:一是基于自适应super-twisting算法和积分滑模控制,提出一种新的有限时间一致性协议用于设计视线方向的加速度指令,保证所有导弹同时达到目标;二是利用自适应控制和非奇异快速终端滑模控制,设计视线法向的加速度指令,保证每枚导弹与目标之间的视线角速率和视线角在有限时间收敛。
附图说明
图1是本发明的有限时间收敛的三维多导弹协同制导方法示意图;
图2是本发明的有限时间收敛的三维多导弹协同制导系统示意图;
图3是本发明实施例的三维空间导弹与目标运动关系示意图;
图4是本发明实施例的导弹通信拓扑示意图;
图5是本发明实施例的导弹与目标相对距离仿真结果示意图;
图6是本发明实施例的导弹剩余时间仿真结果示意图;
图7是本发明实施例的视线方向加速度指令仿真结果示意图;
图8是本发明实施例的视线法向加速度指令仿真结果示意图;
图9是本发明实施例的视线倾角角速率仿真结果示意图;
图10是本发明实施例的视线偏角角速率仿真结果示意图;
图11是本发明实施例的视线倾角仿真结果示意图;
图12是本发明实施例的视线偏角仿真结果示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下参照附图并举出优选实施例,对本发明进一步详细说明。然而,需要说明的是,说明书中列出的许多细节仅仅是为了使读者对本发明的一个或多个方面有一个透彻的理解,即便没有这些特定的细节也可以实现本发明的这些方面。
本发明的发明人考虑到,现有的协同制导方法针对的都是静止或低速的目标,而且需要人为地事先指定攻击时间。另外,在制导过程中各导弹之间没有动态的信息交互,不能实现真正意义上的导弹协同。并且由于采用不连续的控制律,易于产生抖振现象,严重影响制导性能。
基于上述考虑,本发明的发明人首先建立带有攻击角约束的多导弹协同制导模型。其次,设计了针对弹着时间与终端角度进行约束的制导律:基于自适应super-twisting算法和积分滑模控制,提出一种新的有限时间一致性协议用于设计视线方向的加速度指令,保证所有导弹同时达到目标;利用自适应控制和非奇异快速终端滑模控制,设计视线法向的加速度指令,保证每枚导弹与目标之间的视线角速率和视线角在有限时间收敛。
下面结合附图详细说明本发明的技术方案。
实施例一
图1示出了本发明的有限时间收敛的三维多导弹协同制导方法,如图1所示,制导方法具体如下执行:
步骤S1,针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型。
实际应用中,多导弹协同制导模型如下式表示:
其中的参数含义将在后文详细介绍。
较佳地,在设计制导律时,将每枚导弹的剩余时间信息在导弹间进行交互以实现真正地协同制导,提高制导性能。将导弹的剩余时间信息引入上述多导弹协同制导模型,可得到如下式表示的制导模型:
其中的参数含义在后文介绍。
步骤S2,利用自适应控制、积分滑模控制及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标。
具体应用中,上述视线方向的加速度制导律如下式表示:
φ1(σ)=|σ|1/2sign(σ)+l3σ
l2=2l1a+b+4a2
步骤S3,通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛。
实际应用中,上述视线法向的加速度制导律如下式所示:
特别地,上述视线方向的加速度制导律、视线法向的加速度制导律均是连续的,由此可避免抖振的发生。
较佳地,在建立上述二制导律之后,本发明还利用李雅普诺夫理论证明系统的有限时间稳定性。
这样,本发明针对带有攻击角约束的多导弹协同攻击机动目标问题,设计了上述带有攻击角约束的有限时间收敛的三维协同制导律,通过多导弹的有效协同,在不了解目标任何机动信息的情况下,控制所有导弹以期望角度同时击中目标,于此同时较好地消除了抖振现象。
实施例二
图2示出了本发明的有限时间收敛的三维多导弹协同制导系统,参见图2,多导弹协同制导系统包括模型建立单元201、第一制导律建立单元202、第二制导律建立单元203。
具体地,模型建立单元201用于针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型。第一制导律建立单元202用于利用自适应控制、积分滑模控制及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标。第二制导律建立单元203用于通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛。
特别地,上述视线方向的加速度制导律及视线法向的加速度制导律均是连续的。
实施例三
下面将详细推导、论证上述制导模型与制导律。
首先将介绍有限时间一致性与有限时间稳定性。
为了确保多导弹同时到达目标,协同拦截问题是以有限时间一致性为基础的,所以下面给出一些有关有限时间一致性的引理。
考虑如下多智能体系统,每个智能体动态描述为
其中,xi是智能体i的状态,ui表示一致性协议,是状态xi及其相邻智能体状态的函数。
多智能体之间通信网络可由图论来描述。无向图G=(v,ζ,C)代表多智能体之间的通讯拓扑关系,其中v描述图中节点组成的集合,ζ代表节点之间的连线,矩阵C=[cij]∈Rn ×n代表权系数矩阵,节点vi代表智能体i,图G中的任意边可以描述为(vi,vj),其中(vi,vj)代表无序节点对。若智能体i和智能体j之间能够进行信息交换,则否则cij=0,特别的cii=0,i∈{1,2,…,n},由于G是无向图,有cij=cji。若节点vi和节点vj之间存在一条通道,则通路上连续的节点都是相邻的。如果无向图中任意两个节点间都存在至少一条通路,则整个图是连通的。
定义1:如果在任意初始条件下,存在有限时间t*,使得任意智能体j,当t≥t*,xj(t)=x*,其中,x*为实数,则称ui为有限时间一致性协议。
引理1:考虑具有n个个体的多智能体系统1,在其通信拓扑结构图G无向且连通时,则一致性协议:
解决了多智能体有限时间一致性问题,其中,0<θi<1。
引理2:假设存在一个连续正定的函数V(t),满足其中τ>0和0<η<1。那么,系统在有限时间收敛到平衡点,其中收敛时间tf满足
引理3:假设存在一个连续正定函数V(t)满足其中α1>0,α2>0,0<δ<1。那么,系统在有限时间内快速收敛到平衡点,收敛时间tf满足
引理4:假设存在b1,b2,…,bn是正数、0<q<2,那么下面的不等式成立:
下面将展示多导弹协同制导模型的建立过程。
在三维空间,导弹与目标的相对运动关系如图3所示。图3中,T和M分别代表目标和导弹;Mxyz为惯性坐标系;Mx1y1z1为视线坐标系;r是导弹与目标之间的相对距离;qε和qβ分别为视线倾角和视线偏角。
由图1可得目标-导弹相对运动方程为:
其中[aMr,aMε,aMβ]T和[aTr,aTε,aTβ]T分别是导弹加速度和目标加速度在视线坐标系3个轴的分量。
定义qεd和qβd分别是期望的末段视线倾角和视线偏角。定义状态变量x1=r,x3=[qε-qεd,qβ-qβd]T,和结合式4、5、6得到三维制导模型如下:
其中,
引入变量tgo为导弹与目标之间末制导剩余时间,因为在实际的拦截情况中,导弹与目标之间视线方向的相对速度变化相对较小,所以末制导剩余时间tgo可以表示为:
对式8求导可以得到剩余时间的导数为:
其中,
多导弹协同拦截的目标是针对每枚导弹设计制导律aMr和aMq,使得所有的导弹在分布式的通信拓扑下同时击中机动的目标,并且以期望的视线倾角和期望的视线偏角击中目标。为了保证所有导弹能够同时击中目标,即时间的一致性,需要每枚导弹的剩余时间tgo能够通过信息交互被保持一致。
每枚导弹的状态用对应下标i(i=1,2,…,n)表示。为了控制末制导剩余时间实现多导弹同时拦截目标,将tgoi作为一个新的制导系统状态,则第i枚导弹与目标的制导系统状态方程可以表示为:
引入新的控制输入:
将式11代入状态方程10可得:
基于以上分析,由式12可知,在多导弹拦截单一目标的协同制导问题中,协同制导律的设计过程可以分成两部分:一部分是:在存在外界干扰的情况下通过设计新的控制输入使得导弹的末制导剩余时间tgoi在有限时间内达到一致,也就是说,需要设计一个有限时间一致性协议,即使存在外界干扰的情况下也能保证剩余时间tgoi是一致的,进而保证多导弹同时拦截目标。另一部分通过设计加速度指令aMqi使导弹的视线角收敛到期望的终端视线角,视线角速率收敛到零。
下面基于自适应super-twisting滑模控制方法和积分滑模控制方法,提出一种新型的有限时间一致性协议,即通过视线方向上的加速度指令保证所有导弹的剩余时间tgos在有限时间达成一致。
假设1:dri的导数满足以下不等式:
其中,Δri是未知的正常数。
引入变量tfi表示导弹i在t时刻预测的成功拦截目标的时刻,则有:
tfi=t+tgoi 式14
从14中可以看出tgoj-tgoi=tfj-tfi,所以要想使得剩余时间tgoi趋于一致,只需使得tfi趋于一致。
结合式12,可以得到:
基于标准的有限时间一致性动态方程1和2,一种新的积分换设计如下:
其中,tfi(0)是tfi的初始值,称为标准的一致性协议,表达式为:
基于上述设计的积分滑模16和自适应super-twisting算法,提出一种连续有限时间一致性协议如下,即在视线方向上的加速度指令:
其中,l3>0,l1,l2是制导律自适应增益,自适应律为:
l2=2l1a+b+4a2 式20
其中,w1,γ1,a,b>0。
定理1:针对系统15,其通信拓扑结构图G无向且连通,在假设1满足的条件下,制导律选取式18,自适应律为19和20,那么σ在有限时间收敛到零,进而使剩余时间tgoi在有限时间内达到一致。
下面给出证明:首先,结合式18,对积分滑模16进行求导,可得:
引入变量有:
取向量ζ=[|σ|1/2sign(σ)+l3σ,z]T,
根据假设1,可以得到:存在一个正常数ρi满足ρi≥2Δri和
因此,不等式24可以得到:
其中,|χi|≤ρi。
利用式25,并且注意到φ2(σ)=φ′1(σ)φ1(σ),式23能够重写为
其中,
为了证明系统15的稳定性,选择如下李雅普诺夫函数:
其中 *和γ2是正数,P是一个对称正定矩阵。
定义Vζ=ζTPζ和
运用式26,Vζ的导数为:
其中,
根据Shur补性质得到矩阵Q是正定的条件为:
l2=2l1a+b+4a2 式30
因为Vζ=ζTPζ为二次正定函数,故有:
其中,另外由式32可得:
从式28、32、33,可得:
其中,
利用式27、34,根据引理4,V1的导数为:
假设采用自适应律19和20时,l1和l2均有界,于是总是存在常数使得
因此,式35可以写成:
其中,
为了满足有限时间收敛的条件,令τ=0,从而得到增益l1,l2的自适应律为:
选取可以使得式20和式38是一致的,因为由式20可得:
式36可以写成:
根据引理2,可得V1在有限时间收敛到零,因此σ和z在有限时间收敛到零。那么从式22,也在有限时间收敛到零。
那么,意味着根据引理1可知,tfi能够在有限时间内达到一致。因为tgoj-tgoi=tfj-tfi,则剩余时间tgoi在有限时间内达到一致,而且定理1中设计的有限时间一致性协议中的可替换为:
定理证毕。
下面,本发明提出一种新的连续的自适应非奇异快速制导律,来保证视线角速率和视线角的有限时间收敛。
考虑式12中的子系统如下:
假设2:dqi满足||dqi||∞≤Δqi,其中Δqi是未知的正常数。
同一般的滑模控制设计相似,首先设计滑模面,来获得所期望的控制效果。选取如下非奇异快速终端滑模面:
其中,α=diag(α1,α2),β=diag(β1,β2),αj,βj>0(j=1,2),1<λ2<2和λ1>λ2。
对s求导可得:
为了让制导系统的状态快速地从初始状态收敛到滑模面,选择以下的快速幂次趋近律:
其中k1=diag(k11,k12),k2=diag(k21,k22),k1j,k2j>0(j=1,2),0<μ<1。
基于非奇异快速终端滑模和自适应控制,设计以下连续的自适应非奇异快速终端滑模制导律:
其中k4=(k41,k42),k4j>0(j=1,2),k3>0,ε>0,是k的估计值,且
定理2:针对系统42,在假设1和选取滑模面43的情况下,在所设计的制导律46和自适应律47的作用下,滑模面s将在有限时间内收敛到以下区域:
进一步,制导系统状态x3i和x4i将在有限时间内收敛到如下区域:
下面给出证明:令选取Lyapunov函数
对V2求导,可得:
由
和
其中,0<o<2,可得:
其中,
由式54及有界性定理,可知s和是最终一致有界的。因此,可假设其中ξ是一个正常数。
进一步,选取另一个Lyapunov函数:
对V3求导,可得:
式56能够写成以下两种形式
和
根据引理4,对于式57,如果选择和有:
其中
根据引理3,可得滑模面s将在有限时间收敛到以下区域:
针对式58,经过相似的分析,可得滑模面s也将在有限时间收敛到以下区域:
综合式60和式61,可得以下区域可在有限时间达到:
|sj|≤εj=min{ε1j,ε2j} 式62
故当系统轨迹进入区域62,有:
式63可进一步写成:
如果:
则式64也具有和非奇异快速终端滑模类似的结构,这意味着系统轨迹将一直收敛直到换句话说,x4i(j)将在有限时间内收敛到以下区域:
而且从式63和66可得x3i(j)将要在有限时间内收敛到以下区域:
定理2证毕。
下面考虑3枚导弹同时攻击一个机动目标的情况,对所设计的制导律18和46进行仿真验证。3枚导弹之间的通信网络如图4所示,该通信网络是连通的,且对应的权系数矩阵可以描述为导弹和目标之间的相对初始条件和所期望的视线角如表1所示。目标的机动形式为:aTr=50sin(0.5t)m/s2,aTβ=50cos(0.5t)m/s2。各导弹所能提供的最大加速度指令为40g,g=9.8m/s2。仿真步长取0.001s。
视线方向上的制导律18的仿真参数选取如下:
l3=0.38,w1=0.2,γ1=1,a=0.1,b=1,θi=0.3(i=1,2,3)
视线法向上的制导律46的仿真参数选取如下:
α=diag[1,1],β=diag[0.2,0.2],λ1=3,λ2=1.5,k1=diag[0.2,0.2],k2=diag[2,2],μ=0.3,k4=diag[1,3]k3=0.2,ε=2。
表1
仿真结果如图5到图12、表2所示。从图5和图6可以看到,在协同制导律的作用下,三枚导弹同时拦截目标。图7表明三枚导弹的视线方向的加速度指令在合理的范围之内、并在初始阶段有相应的变化。这是为了调整三枚导弹的末制导剩余时间快速趋于一致,随着剩余时间趋于一致,视线方向的加速度指令曲线很快也趋于一致,并随目标机动不断变化,实现对机动目标的同时拦截。图8表明三枚导弹的视线法向上的加速度指令曲线。为了能够使视线角速率收敛到零,制导初期加速度值比较大,但随着视线角速率收敛到零,加速度逐渐平缓并跟踪上目标加速度。由图9-12可知,三枚导弹相应的视线角速率收敛到零,视线角收敛到相应的期望值。由表2可见,三枚导弹都以较小的脱靶量精确命中目标,并且制导时间一致。
表2
根据本发明提供的有限时间收敛的三维多导弹协同制导方法与系统,针对三维空间中多导弹协同拦截一个机动目标的问题,提出了一种具有攻击角约束的有限时间协同制导律,可在不知道目标任何机动信息的情况下,使所有导弹以期望角度同时击中目标。同时,本发明采用连续的有限时间一致性协议进行控制律设计,消除了抖振的发生。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读取存储介质中,如:ROM/RAM、磁碟、光盘等。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种有限时间收敛的三维多导弹协同制导方法,其特征在于,包括步骤:
S1.针对未知目标,建立三维空间内的、具有攻击角约束的多导弹协同制导模型;
S2.利用自适应控制、积分滑模控制、及有限时间一致性协议建立视线方向的加速度制导律,使所有导弹可同时击中目标;
S3.通过自适应控制与非奇异快速终端滑模控制,设计视线法向的加速度制导律,使每枚导弹的视线角速率及视线角实现有限时间收敛;
其中,所述视线方向的加速度制导律、及所述视线法向的加速度制导律均是连续的;
所述方法还包括:
将每枚导弹的剩余时间信息在导弹间进行交互;
其中,所述多导弹协同制导模型具体为:
其中,qε为视线倾角,qβ为视线偏角,aMr为导弹加速度在视线坐标系x轴的分量,aTr为目标加速度在视线坐标系x轴的分量,r为导弹与目标的相对距离,qεd为期望的视线倾角,qβd为期望的视线偏角,aMε为导弹加速度在视线坐标系y轴的分量,aMβ为导弹加速度在视线坐标系z轴的分量,aTε为目标加速度在视线坐标系y轴的分量,aTβ为目标加速度在视线坐标系y轴的分量。
2.如权利要求1所述的方法,其特征在于,在建立多导弹协同制导模型之后,设计制导律之前,所述方法还包括:
将导弹的剩余时间信息引入多导弹协同制导模型,获得如下式表示的制导模型;
其中,i表示第i枚导弹,tgoi为第i枚导弹的剩余时间,
3.如权利要求2所述的方法,其特征在于,所述视线方向的加速度制导律由下式表示:
φ1(σ)=|σ|1/2sign(σ)+l3σ
l2=2l1a+b+4a2 式20
其中,ui nom为标准一致性协议,l1、l2、l3为自适应增益,l3>0,tfi为导弹拦截目标的预测时刻,tfi(0)为tfi的初始值,w1,γ1,a,b>0,w2,γ2,l1*,l2*>0。
4.如权利要求3所述的方法,所述视线法向的加速度制导律由下式表示:
其中,α=diag(α1,α2),β=diag(β1,β2),α1、α2、β1、β2>0,1<λ2<2,λ1>λ2,k1=diag(k11,k12),k2=diag(k21,k22),k11,k12,k21,k22>0,0<μ<1,s为选取的非奇异快速终端滑模面,ε>0,为k的估计值,Δqi为未知正数,k3>0,k4=(k41,k42),k41,k42>0。
5.如权利要求1-4任一所述的方法,其特征在于,在建立所述视线方向的加速度制导律、及所述视线法向的加速度制导律之后,所述方法还包括:
利用李雅普诺夫理论证明系统的有限时间稳定性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611256552.4A CN106843265B (zh) | 2016-12-30 | 2016-12-30 | 有限时间收敛的三维多导弹协同制导方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611256552.4A CN106843265B (zh) | 2016-12-30 | 2016-12-30 | 有限时间收敛的三维多导弹协同制导方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106843265A CN106843265A (zh) | 2017-06-13 |
CN106843265B true CN106843265B (zh) | 2019-12-24 |
Family
ID=59113574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611256552.4A Active CN106843265B (zh) | 2016-12-30 | 2016-12-30 | 有限时间收敛的三维多导弹协同制导方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106843265B (zh) |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107818219B (zh) * | 2017-10-31 | 2018-09-25 | 中国人民解放军国防科技大学 | 一种面向突防的多导弹协同弹道规划方法 |
CN108008736B (zh) * | 2017-11-21 | 2020-12-01 | 山东科技大学 | 飞行器协同控制方法、装置、计算机可读存储介质及终端 |
CN108052112B (zh) * | 2017-12-01 | 2020-10-02 | 哈尔滨工业大学 | 基于pn制导律辨识的多飞行器威胁度获取方法 |
CN107908895B (zh) * | 2017-12-01 | 2020-11-13 | 哈尔滨工业大学 | 一种基于imm滤波器主动防御视线策略协同制导律设计方法 |
CN108287472B (zh) * | 2017-12-20 | 2021-11-16 | 中国人民解放军海军大连舰艇学院 | 采用滚动时域框架的舰队防空决策与自动化调度系统与方法 |
CN108416098B (zh) * | 2018-02-02 | 2021-07-16 | 北京航空航天大学 | 一种拦截机动目标的攻击时间约束制导律设计方法 |
CN108362171B (zh) * | 2018-02-02 | 2019-03-29 | 北京航空航天大学 | 一种具有攻击时间和攻击角度约束的制导策略 |
CN110162084B (zh) * | 2018-10-11 | 2022-01-11 | 北京机电工程研究所 | 基于一致性理论的飞航导弹集群系统编队控制方法 |
CN109343563B (zh) * | 2018-10-15 | 2020-06-05 | 北京理工大学 | 考虑舵机失效及落角约束的飞行器制导系统及方法 |
CN109557814B (zh) * | 2018-12-04 | 2021-08-20 | 四川航天系统工程研究所 | 一种有限时间积分滑模末制导律 |
CN109977613B (zh) * | 2019-04-19 | 2021-01-01 | 哈尔滨工业大学 | 一种可预先设定调整时间的自适应滑模末制导律设计方法 |
CN110119087B (zh) * | 2019-05-05 | 2021-12-21 | 西北工业大学 | 一种有向通信下二阶多智能体系统固定时一致性跟踪方法 |
CN110083179B (zh) * | 2019-05-07 | 2021-10-15 | 西北工业大学 | 一种预定时间多智能体系统一致性跟踪控制方法 |
CN110187640B (zh) * | 2019-06-29 | 2022-04-29 | 东南大学 | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 |
CN110412874B (zh) * | 2019-07-25 | 2022-05-13 | 东南大学 | 针对机动目标和时延通信的多导弹协同制导律设计方法 |
CN110471283B (zh) * | 2019-08-16 | 2022-06-07 | 南京理工大学 | 一种带碰撞角约束的三维鲁棒制导律构建方法 |
CN110471275A (zh) * | 2019-08-30 | 2019-11-19 | 哈尔滨工业大学 | 一种非奇异终端滑模有限时间收敛角度约束制导方法 |
CN110487132B (zh) * | 2019-08-30 | 2020-05-19 | 哈尔滨工业大学 | 一种基于非奇异快速终端滑模控制的角度约束制导方法 |
CN110716585B (zh) * | 2019-10-31 | 2020-08-21 | 西安电子科技大学 | 基于分布式一致性与自组网的自主协同控制方法及系统 |
CN110687801A (zh) * | 2019-11-22 | 2020-01-14 | 深圳信息职业技术学院 | 一种基于前向制导的自适应积分滑模三维制导律设计方法 |
CN111027206B (zh) * | 2019-12-05 | 2024-02-09 | 哈尔滨工业大学 | 具有规定性能的拦截机动目标自适应滑模控制方法 |
CN111707148B (zh) * | 2020-06-28 | 2021-02-09 | 北京航空航天大学 | 一种时变时延条件下多导弹协同制导方法及系统 |
CN111930142B (zh) * | 2020-08-04 | 2022-10-11 | 西北工业大学 | 速度不可控条件下的多弹编队协同控制方法 |
CN112050692B (zh) * | 2020-08-17 | 2022-08-09 | 北京电子工程总体研究所 | 一种多导弹协同制导律设计方法 |
CN112033234B (zh) * | 2020-08-17 | 2022-08-09 | 北京电子工程总体研究所 | 一种多导弹协同制导律设计方法 |
CN112346474B (zh) * | 2020-10-20 | 2021-12-07 | 南京航空航天大学 | 一种有限时间收敛的微分对策制导律的设计方法 |
CN112344795B (zh) * | 2020-11-08 | 2022-11-11 | 西北工业大学 | 一种预定时间收敛的末制导方法 |
CN112379599B (zh) * | 2020-11-20 | 2022-06-14 | 西北工业大学 | 一种针对空中机动目标的三维协同中制导方法 |
CN112527011B (zh) * | 2020-11-20 | 2022-07-15 | 西北工业大学 | 一种编队协同中制导技术统一框架的构建方法 |
CN112577373B (zh) * | 2020-11-27 | 2022-05-27 | 南京理工大学 | 一种基于比例导引律的多导弹同时攻击静目标的制导律 |
CN112558634B (zh) * | 2020-12-08 | 2022-02-01 | 北京航空航天大学 | 一种多导弹协同制导方法及系统 |
CN112631328B (zh) * | 2020-12-12 | 2023-01-24 | 西北工业大学 | 一种多目标协同末制导律方法 |
CN112859921B (zh) * | 2021-01-21 | 2022-02-22 | 北京航空航天大学 | 一种三维同时攻击鲁棒协同制导律设计方法 |
CN112902767B (zh) * | 2021-01-28 | 2022-05-20 | 西安交通大学 | 一种多弹时间协同的导弹制导方法及系统 |
CN113110576B (zh) * | 2021-04-14 | 2022-09-13 | 西北工业大学 | 一种可连续切换的自适应固定时间收敛协同制导方法 |
CN113359813B (zh) * | 2021-05-11 | 2022-08-05 | 西北工业大学 | 一种多领弹多从弹群组协同制导方法 |
CN113625745B (zh) * | 2021-07-02 | 2022-08-23 | 西北工业大学 | 一种基于切换固定时间收敛理论的攻击时间控制制导方法 |
CN114371737A (zh) * | 2022-01-17 | 2022-04-19 | 哈尔滨工业大学 | 一种角度约束的智能弹药最优对策协同制导方法 |
CN115342815B (zh) * | 2022-08-26 | 2024-04-26 | 哈尔滨工业大学 | 反大气层内或临近空间机动目标视线角速率估计方法 |
CN115268503B (zh) * | 2022-09-28 | 2023-01-10 | 中国人民解放军国防科技大学 | 一种去奇异化的多飞行器三维协同制导方法 |
CN117193336B (zh) * | 2023-09-12 | 2024-03-22 | 中国船舶集团有限公司第七一九研究所 | 一种基于时间和角度约束的水下航行物协同制导方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245256B (zh) * | 2013-04-25 | 2015-03-11 | 北京理工大学 | 一种多导弹协同作战制导律设计方法 |
CN104266546B (zh) * | 2014-09-22 | 2015-10-07 | 哈尔滨工业大学 | 一种基于视线的有限时间收敛主动防御制导控制方法 |
-
2016
- 2016-12-30 CN CN201611256552.4A patent/CN106843265B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106843265A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106843265B (zh) | 有限时间收敛的三维多导弹协同制导方法与系统 | |
Song et al. | Three-dimensional cooperative guidance law for multiple missiles with finite-time convergence | |
Kumar et al. | Cooperative salvo guidance using finite-time consensus over directed cycles | |
Zhao et al. | Cooperative time-varying formation guidance for leader-following missiles to intercept a maneuvering target with switching topologies | |
Zhang et al. | Two-stage cooperative guidance strategy using a prescribed-time optimal consensus method | |
Dong et al. | Necessary and sufficient conditions for average formation tracking of second-order multi-agent systems with multiple leaders | |
Zhang et al. | Fault tolerance of cooperative interception using multiple flight vehicles | |
He et al. | Three-dimensional salvo attack guidance considering communication delay | |
Yu et al. | Three-dimensional adaptive fixed-time cooperative guidance law with impact time and angle constraints | |
Xu et al. | Affine formation control for heterogeneous multi-agent systems with directed interaction networks | |
CN111027206B (zh) | 具有规定性能的拦截机动目标自适应滑模控制方法 | |
Zadka et al. | Consensus-based cooperative geometrical rules for simultaneous target interception | |
Liu et al. | Three-dimensional cooperative guidance law for intercepting hypersonic targets | |
Cichella et al. | A 3D path-following approach for a multirotor UAV on SO (3) | |
Yan et al. | Limited communication consensus control of leader-following multi-UUVs in a swarm system under multi-independent switching topologies and time delay | |
Tang et al. | Cooperative circular guidance of multiple missiles: A practical prescribed-time consensus approach | |
Boukhnifer et al. | H-infinity robust control of 3-DOF helicopter | |
Wang et al. | Cooperative guidance law against highly maneuvering target with dynamic surrounding attack | |
Samiei et al. | Attitude stabilization of rigid spacecraft with minimal attitude coordinates and unknown time-varying delay | |
Wang et al. | Integrated distributed formation flight control with aerodynamic constraints on attitude and control surfaces | |
Guo et al. | Pursuit-Evasion Games of Marine Surface Vessels Using Neural Network-Based Control | |
Tao et al. | Optimal spatial-temporal cooperative guidance against a maneuvering target | |
Liu et al. | Integrated guidance and control of interceptor missile based on asymmetric barrier Lyapunov function | |
Chen et al. | Adaptive Cooperative Game Penetration Guidance for Multiple Hypersonic Vehicles with Performance Constraints | |
Zhou et al. | Distributed cooperative guidance for multi-missile system with time-varying delays and switching topologies |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Song Shenmin Inventor after: Tan Liguo Inventor after: Song Junhong Inventor after: Wang Yan Inventor before: Song Shenmin Inventor before: Song Junhong Inventor before: Wang Yan |
|
GR01 | Patent grant | ||
GR01 | Patent grant |