留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

多点源诱偏系统对抗反辐射导弹效能分析

赵娟娟 马东立 郑江安

刘畅, 杨锁昌, 汪连栋, 等 . 基于快速自适应超螺旋算法的制导律[J]. 北京航空航天大学学报, 2019, 45(7): 1388-1397. doi: 10.13700/j.bh.1001-5965.2018.0654
引用本文: 赵娟娟, 马东立, 郑江安等 . 多点源诱偏系统对抗反辐射导弹效能分析[J]. 北京航空航天大学学报, 2009, 35(9): 1139-1143.
LIU Chang, YANG Suochang, WANG Liandong, et al. Guidance law based on fast adaptive super-twisting algorithm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2019, 45(7): 1388-1397. doi: 10.13700/j.bh.1001-5965.2018.0654(in Chinese)
Citation: Zhao Juanjuan, Ma Dongli, Zheng Jiang-anet al. Effectiveness analysis of multi-source decoying system against ARM[J]. Journal of Beijing University of Aeronautics and Astronautics, 2009, 35(9): 1139-1143. (in Chinese)

多点源诱偏系统对抗反辐射导弹效能分析

详细信息
    作者简介:

    赵娟娟(1982-),女,陕西商洛人,博士生,juanjuan0410@ase.buaa.edu.cn.

  • 中图分类号: TJ 762.2;TN 974

Effectiveness analysis of multi-source decoying system against ARM

  • 摘要: 基于有源诱偏系统对抗反辐射导弹(ARM,Anti-Radiation Missile)的理论,导出诱偏误差公式,以3种典型布局方案(三点源布局1、三点源布局2及四点源布局)为例,对误差公式进行简化处理,通过数值计算得到诱偏误差曲线,分析不同布局方案下辐射源功率比对诱偏误差的影响;同时分析ARM的实际攻击过程,引入ARM机动修正,通过仿真分析得到ARM落点图和雷达保护概率曲线,对比不同布局方案下的诱偏效果,结果表明四点源布局为最优方案,对研究多点源诱偏系统对抗ARM的作战效能具有指导意义.

     

  • 随着现代战场的日益复杂,传统比例导引已经远远不能满足制导系统的发展需求。人们将先进的控制理论应用于导弹的制导系统中,如最优控制、逆系统控制、H控制、微分几何控制、随机系统最优控制、滑模变结构控制等,以解决强对抗条件下的精确制导问题。最优制导律虽然在理论上可以实现零脱靶量,但形式复杂,需要信息多,且对信息误差相当敏感。逆系统控制、H控制、微分几何控制、随机系统最优控制等方法虽然具有一定的鲁棒性,但均存在形式复杂、需要信息多的缺点。

    滑模变结构以其抗干扰特性和设计简单而成功的应用于导弹制导律设计,但受到滑模控制固有的抖振和渐进稳定特性影响,滑模变结构制导律普遍存在抖振严重、收敛速度慢等局限[1]。克服抖振的方法之一是以饱和函数或边界层代替符号函数,但这使得系统轨迹稳定在滑模面附近而不是滑模面上,进而丧失了滑模控制的鲁棒性[2]。另一种抑制抖振的方法是采取高阶滑模控制,高阶滑模在抑制抖振的同时还能保持对干扰的鲁棒性,消除了相对阶的限制, 提高了控制精度[3]

    二阶滑模是目前应用最为广泛的高阶滑模控制方法,因为它的控制器结构简单且需要信息少。螺旋算法、超螺旋(Super-Twisting, ST)算法、次优算法和给定收敛律算法是二阶滑模中常用的4种算法[4],与其他二阶滑模算法相比,ST算法仅需要滑模变量而不需要滑模变量的一阶导数,因此广泛应用于制导律设计。ST算法的有限时间稳定条件需要已知系统不确定性上界,实际应用中这个界的精确值很难确定,为得到系统的稳定控制,往往选取尽可能大的参数,结果带来系统损坏和剧烈抖振[5]。为此多种自适应ST算法相继被提出[6-8],但这些算法大多假设系统不确定性满足一定形式的假设,且在系统状态距离平衡点较远时收敛速度较慢。

    针对上述问题,本文在标准ST算法的基础上引入线性项[9],设计了一种新的自适应律,提出了一种新的快速自适应超螺旋(Fast Adaptive Super-Twisting, FAST)算法,该算法不需要已知系统不确定性的边界且收敛速度较快。首先,利用类二次型Lyapunov函数[10-12]证明了系统的有限时间稳定性,并给出了收敛时间公式。随后,将FAST算法应用于末制导问题,提出了一种新的二阶滑模制导律。最后,通过数字仿真将所提制导律与自适应滑模制导律、ST制导律和光滑二阶滑模制导律进行了对比。仿真结果表明,本文算法具有更高的命中精度、更快的收敛速度和更强的鲁棒性。

    为了研究导弹拦截目标过程中的制导律,首先建立导弹与目标的相对运动模型。对于非滚转导弹,在末制导过程中,姿态控制系统可以控制导弹不滚转,因而弹目相对运动可以解耦成纵向平面和横向平面的运动。以纵向平面为例,弹目相对运动如图 1所示,横向平面的推导与之类似。

    图  1  导弹和目标相对运动示意图
    Figure  1.  Schematic diagram of relative motion of missile and target

    图 1中,MT分别代表导弹和目标所处位置;r为平面内两者的相对距离;q为弹目视线角;vmvt分别为导弹和目标的速度;θmθt分别为弹道倾角和目标航向角;amat分别为导弹和目标的法向加速度。由图 1可以得出弹目相对运动方程为[13]

    (1)

    (2)

    式中:分别为rq对于时间的导数。

    为便于推导,令导弹与目标的相对径向速度,相对法向速度,将其代入式(1)和式(2)后,对式(1)和式(2)相对于时间求一阶导数,得

    (3)

    (4)

    (5)

    其中:atramr分别为目标加速度和导弹加速度在视线方向上的分量;atqamq分别为目标加速度和导弹加速度在视线法向方向上的分量。

    整理式(1)~式(5),得

    (6)

    (7)

    代入式(7),得到系统方程:

    (8)

    在制导律设计时,将amqatq分别视为系统的控制量和干扰量。根据准平行接近原理,设计制导律的关键在于如何通过amq控制视线角速率,令其趋近于0。由式(5)可知:

    (9)

    在末制导过程中,由于受到过载能力的限制,导弹和目标实际所能提供的最大侧向加速度是有限的。同时受到导引头角跟踪系统的功率、接收机过载等因素的限制,导引头存在最小作用距离r0,当弹目相对距离小于或等于r0时,制导回路断开。记末制导开始时刻为0,不失一般性制导过程满足如下假设。

    假设1   存在常数Am>0, At>0, A1>0, A2>0,使得

    (10)

    假设2   系统(8)中的时变参数r(t)满足:

    (11)

    引理1[14]   对于非线性系统,有

    (12)

    假设存在连续可微函数V(x):UR满足:

    1) V(x)为正定函数。

    2) 存在正实数ζ1>0和α∈(0, 1),以及包含原点的开邻域U0U,使得下式成立:

    (13)

    则系统(12)有限时间稳定;若U=U0=Rn,则系统(12)全局有限时间稳定。收敛时间treach满足:

    (14)

    引理2[15]   对于非线性系统(12),假设存在连续可微函数V(x):UR满足:

    1) V(x)为正定函数。

    2) 存在正实数ζ1>0,ζ2>0和α∈(0, 1),以及包含原点的开邻域U0U,使得下式成立:

    (15)

    则系统(12)有限时间稳定;若U=U0=Rn,则系统(12)全局有限时间稳定。收敛时间treach满足:

    (16)

    考虑一阶系统:

    (17)

    式中:xR为系统状态(同时也是滑模变量);uRER分别为控制输入和不确定项。定义E=E1(x, t)+E2(x, t),E1(x, t)表示不可微的不确定性,E2(x, t)表示可微的不确定性。

    系统(17)的ST算法可表示为[16]

    (18)

    式中:sgn(x)为符号函数;k1k2为待设计的参数。若E1(x, t)=0,δ>0,则k1k2的取值满足[4]

    (19)

    显然,k1k2需要根据不确定性的上界确定。

    将式(18)代入式(17),令x1=xx2=u1,化简后的控制系统为

    (20)

    ρ1=E1(x1, t),,式(20)的等价形式为

    (21)

    式(21)的有限时间稳定性证明及收敛时间估计见文献[4]。

    假设3   E1(x, t)和E2(x, t)满足:

    (22)

    式中:g1g2为未知正数; ϕ1(x)和ϕ2(x)为滑模变量的函数,对不同的控制律有不同的形式。

    针对系统(17),设计FAST算法如下:

    (23)

    式中:

    (24)

    自适应参数控制器为

    (25)

    式中:abc1c2为任意正数。

    将式(23)代入式(17),得到控制系统为

    (26)

    式中:

    (27)

    自适应参数控制器为

    (28)

    ξT=[ϕ1(x1), x2],由

    (29)

    可得

    (30)

    式中:

    由假设3可知:

    (31)

    式(31)的等效形式为

    (32)

    式(32)可改写为

    (33)

    式中:D=[1,0];

    定理1   当系统(26)满足假设3时,存在k1*k2*使得x1在有限时间treach内从任意初始位置收敛到0。

    证明

    (34)

    (35)

    (36)

    式中:ε1为足够小的正常数。若

    (37)

    Q11 < 0,Q22 < 0,Q12=Q21=0,于是Q为半负定矩阵。令

    (38)

    由分块矩阵的性质可以得到Φ为半负定矩阵。

    选取类二次型Lyapunov函数V=ξT,由式(29)和式(30)得

    (39)

    结合

    (40)

    可得

    (41)

    ϕ1(x1)=|x1|1/2sgn(x1)+x1可得

    (42)

    (43)

    因此

    (44)

    (45)

    由引理2可知,若k1*k2*满足式(37),则ξ=[ϕ1(x1), x2]T在有限时间内收敛到0,即x1在有限时间内收敛到0,且收敛时间满足:

    (46)

    证毕

    由定理1可知,对于未知常数g1g2,存在满足式(37)约束的k1*k2*使得x1在有限时间内收敛到0。

    定理2   当系统(26)满足假设3,控制器参数k1k2满足式(28),g1g2为未知数时,x1在有限时间treach内收敛到0。

    证明

    1) 若k1k1*, k2k2*,且,由定理1可知,定理2显然成立。

    2) 若k1 < k1*, k2 < k2*,选取类二次型Lyapunov函数:

    (47)

    (48)

    式中:

    (49)

    显然

    (50)

    由定理1的证明过程可知:

    (51)

    ① 若|x1|>0,则

    (52)

    (53)

    结合式(51),得到

    (54)

    式中:

    (55)

    于是有

    (56)

    式中:

    (57)

    (58)

    由引理1可知,系统(26)有限时间稳定,且收敛时间满足:

    (59)

    ② 若|x1|=0,则系统状态到达滑模面,且,因此k1=k1*k2=k2*,由引理1可知,系统有限时间收敛。证毕

    由于s很难严格收敛至0,而是在0的一个极小邻域内波动,使得k1k2过大,进而引起系统不稳定。因此,将自适应参数控制器改进为

    (60)

    式中:ε2为很小的正数。

    由准平行接近原理可知,希望在制导过程中趋近于0,因此选取滑模面为

    (61)

    式中:c0=0.1。

    结合式(7)可得

    (62)

    式中:

    (63)

    由定理3及式(9)得到FAST制导律为

    (64)

    式中:

    (65)

    其中:k1k2的值满足式(60)。

    FAST制导律对外界干扰具有鲁棒性,且能够在有限时间内收敛。由制导律的形式可以看出,参数k1k2随着s的变化实时改变,且不需要已知外部干扰的上界。

    为验证所提制导律的有效性,将FAST制导律(Fast Adaptive Super-Twisting Guidance, FASTG)与自适应滑模制导律(Adaptive Sliding Mode Guidance,ASMG)、光滑二阶滑模制导律(Smooth Second Order Sliding Mode Guidance,SSOSMG)、ST制导律(Super-Twisting Guidance,STG)进行对比。其中,ASMG为一阶滑模制导律,其他为二阶滑模制导律,且4种制导律均采取式(61)作为滑模面。

    ASMG为[17]

    (66)

    式中:k1=k2=2。

    SSOSMG为[18]

    (67)

    式中:N=4,k1=3,k2=1,atq的估计值并假设其已知。

    STG为[19]

    (68)

    式中:k1=4,k2=3。

    FASTG形式如式(64)所示,参数满足:a=0.5, b=1, c1=c2=1, ε2=0.1。

    制导律所需的r等信息均由导弹导引头测量得到,且不考虑噪声和干扰。导弹自动驾驶仪视为理想环节,即不考虑延迟。仿真终止时间为弹目相对距离最近的时间点,仿真步长为0.01 s。

    为了全面分析4种制导律的特点,分别在目标无机动和目标有机动2种情形下对制导性能进行比较。

    仿真条件:导弹初始位置为(0, 0) km,初始速度为500 m/s,初速度方向为30°,初始过载约束为15g,切向过载为0,法向过载由制导律得到。目标初始位置为(10, 5) km,初始速度为300 m/s,初速度方向为180°,切向过载和法向过载均为0。

    MATLAB仿真仿真结果如表 1图 2~图 4所示。表 1为情形1的仿真结果,其中Δ为脱靶量,treach分别为滑模变量实际和理论收敛时间,tf为制导时间。图 2为导弹弹道曲线,黑色虚线为目标运动曲线,绿色、黑色、红色、蓝色实线分别为ASMG、SSOSMG、STG、FASTG四种制导律下的导弹运动曲线。图 3(a)图 3(b)分别为滑模变量s(t)及滑模变量的一阶导数变化曲线。图 4为导弹法向过载变化曲线。

    表  1  仿真实验结果(情形1)
    Table  1.  Simulation experimental results (Case 1)
    制导律 Δ/m treach/s tf/s
    ASMG 0.462 2 7.36 10.53 14.94
    SSOSMG 0.518 6 12.2 21.83 14.91
    STG 1.175 6 6.65 12.25 14.92
    FASTG 0.616 0 5.30 14.64 14.90
    下载: 导出CSV 
    | 显示表格
    图  2  导弹弹道曲线(情形1)
    Figure  2.  Missile ballistic curves (Case 1)
    图  3  滑模变量及其一阶导数变化曲线(情形1)
    Figure  3.  Variation curves of sliding-mode variable and its first-order derivative (Case 1)
    图  4  导弹法向过载变化曲线(情形1)
    Figure  4.  Variation curves of missile normal overload (Case 1)

    图 2表 1中的Δtf可知,4种制导律的弹道平直且差异很小,制导时间接近,脱靶量STG最大,ASMG最小,SSOSMG和FASTG略高于ASMG,由于脱靶量远小于导弹杀伤半径,因此可视为直接命中目标。由图 3表 1中的treach可知,4种制导律均能够在有限时间内使得s(t)及收敛至0,其中,FASTG收敛速度最快,STG和ASMG略慢,SSOSMG最慢且时间远大于其他方法。由此可以看出,FASTG在收敛速度方面较其他方法有一定的提高,但优势并不明显。由图 4可知,在导弹发射的前5 s内,导弹过载a由大变小并逐渐稳定到0,其中FASTG的过载远大于其他制导律且变化最为剧烈,最大过载高于其他制导律4倍以上。

    总之,在目标无机动的情形下,由于弹道接近直线,与其他方法相比,FASTG的脱靶量和制导时间相近,收敛速度最快,过载最大,其优势并不明显。

    仿真条件:导弹初始位置为(0, 0) km,初始速度为500 m/s,初速度方向为45°,最大过载为30g,切向过载为0,法向过载由制导律得到。目标初始位置为(2, 2) km,初始速度为300 m/s,初速度方向为0°,切向过载为0,法向过载nt=(5cos t)g

    MATLAB仿真结果如表 2图 5~图 7所示,其中曲线和变量含义与情形1相同。

    表  2  仿真实验结果(情形2)
    Table  2.  Simulation experimental results (Case 2)
    制导律 Δ/m treach/s tf/s
    ASMG 0.549 4 20.31 12.57
    SSOSMG 0.566 2 10.91 21.30 11.92
    STG 0.319 7 20.43 11.97
    FASTG 0.875 3 1.79 18.98 11.74
    下载: 导出CSV 
    | 显示表格
    图  5  导弹弹道曲线(情形2)
    Figure  5.  Missile ballistic curve (Case 2)
    图  6  滑模变量及其一阶导数变化曲线(情形2)
    Figure  6.  Variation curves of sliding-mode variable and its first-order derivative (Case 2)
    图  7  导弹法向过载变化曲线(情形2)
    Figure  7.  Variation curves of missile normal overload (Case 2)

    图 5图 2相比,导弹弹道弯曲,其中ASMG最高,STG和SSOSMG次之,FASTG最低。表 2中的脱靶量STG最小,ASMG和SSOSMG次之,FASTG最大,由于均远小于导弹杀伤半径,均可视为直接命中目标。制导时间FASTG最短,ASMG最长,STG和SSOSMG介于两者之间。由图 6可知,FASTG能够使得s(t)及在很短的时间内收敛至0,SSOSMG的收敛速度很慢,收敛时间接近制导时间,而ASMG和STG不能收敛,且s(t)大幅度波动。由此可见,FASTG在收敛速度方面明显优于其他制导律。由图 7可知,由于s(t)的波动,过载在0附近上下波动。初始阶段,FASTG由于s(t)变化最为剧烈,因此过载最大。中间阶段4种制导律均小于5g。制导末段,STG和ASMG的过载迅速增加并达到最大过载,FASTG的过载在10g左右,SSOSMG的过载为5g左右。

    总之,目标机动情形下,弹道弯曲,FASTG制导时间最短,脱靶量相近,收敛特性明显优于其他方法,仅在初始阶段的过载较大。因此,FASTG在此情形下的表现更为优异。

    1) 在标准ST算法的基础上,增加了自适应参数控制器和线性项,提出了FAST算法。在系统不确定性上界未知的前提下,一方面控制器参数能够自适应调节,避免参数过大造成系统不稳定;另一方面系统在远离平衡点时具有更快的收敛速度,提升了标准ST算法的收敛特性。

    2) 利用二次型Lyapunov函数证明了FAST算法的有限时间稳定性,与其他证明方法相比,该方法计算较为简单,且能够得到收敛时间的估计公式。

    3) 将FAST算法成功地应用于制导律设计。仿真结果表明,本文算法在保留标准ST算法有效抑制抖振、鲁棒性强等优点的同时,具有快速收敛特性且不需要已知不确定性的上界,使得制导系统拥有更高的命中精度和稳定性。

  • [1] McWhorter L V. Response of various monopulse seekers to a multi-source environment Southeastlon 1981 Conference Proceedings. San Diego,California: IEEE, 1981:698-709 [2] 顾尔顺. 有源诱偏反辐射导弹的理论[J]. 现代防御技术,1993(3):25-34 Gu Ershun. The theory of active decoying system against anti-radiation missile [J]. Modern Defence Technology,1993(3):25-34(in Chinese) [3] 顾尔顺. 对反辐射导弹导引头进行有源诱偏的原理[J]. 现代防御技术,1992(3):40-57 Gu Ershun. The theory of active decoying system against anti-radiation missile seeker [J]. Modern Defence Technology, 1992(3):40-57(in Chinese) [4] 关成斌,陈铁柱,王国宏,等. 用不同高度有源诱偏系统对抗ARM [J]. 雷达科学与技术,2004,2(4):193-197 Guan Chengbin, Chen Tiezhu, Wang Guohong, et al. ARM countering of active decoying system with different height baits [J].Radar Science and Technology,2004,2(4):193-197(in Chinese) [5] 谭帆,苏东林,姜铁华. 诱偏反辐射导弹四点源诱饵系统布站仿真[J]. 信息技术,2005(12):96-97 Tan Fan, Su Donglin, Jiang Tiehua. Simulation of four-source decoy system distribution for decoy anti-radiation missile [J]. Information Technology,2005(12): 96-97(in Chinese) [6] 闵涛,姜文利,周一宇,等. 诱偏反辐射导弹三点源诱饵系统布站研究[J]. 航天电子对抗,2000(3):13-15 Min Tao, Jiang Wenli, Zhou Yiyu, et al. The study on the topology of decoying system with three resources against anti-radiation missile [J]. Aerospace Electronic Warfare,2000(3):13-15(in Chinese) [7] 姜铁华,王国玉,苏东林,等. 反辐射导弹技术参数对诱偏系统性能的影响研究[J]. 航空学报,2007,28(增刊):S146-S150 Jiang Tiehua, Wang Guoyu, Su Donglin, et al. Study on effects of ARM technical property on decoying performances of active decoying system [J]. Acta Aeronautica et Astronautica Sinica, 2007, 28(Supplement):S146-S150 (in Chinese)
  • 期刊类型引用(13)

    1. 贾晓涵,付丽霞. 四旋翼飞行器多模式稳定控制. 小型微型计算机系统. 2024(03): 570-575 . 百度学术
    2. 李雪,王青. 时间约束三维超螺旋滑模制导控制一体化方法. 宇航学报. 2024(02): 283-294 . 百度学术
    3. 周蒙,钱惟贤,任侃. 多约束超螺旋滑模变结构制导律. 兵工学报. 2023(03): 799-805 . 百度学术
    4. 赵玮. 考虑外部干扰的超螺旋ZNN多机械臂协调运动控制. 机械设计与制造. 2023(08): 139-147 . 百度学术
    5. 贺玉晓,王丽梅. 永磁直线同步电动机迭代超螺旋滑模控制. 电气技术. 2022(06): 24-29 . 百度学术
    6. 李中奇,张俊豪. 基于快速自适应超螺旋算法的高速列车最优黏着控制. 铁道科学与工程学报. 2022(08): 2143-2150 . 百度学术
    7. 金爱娟,王硕勋,李少龙,刘付鑫,吴铭毅. 基于改进超螺旋算法的永磁同步电动机控制. 包装工程. 2022(19): 198-207 . 百度学术
    8. 齐照辉,张今,王远卓,王嘉,徐梦荣,张程. 基于自适应超螺旋滑模弹道跟踪控制方法. 西北工业大学学报. 2022(05): 1109-1115 . 百度学术
    9. 王永尚,高国琴,方志明. 变质心喷砂除锈移动平台自适应全局鲁棒超螺旋滑模控制. 软件导刊. 2022(10): 156-164 . 百度学术
    10. 李军,廖宇新,李珺. 三维自适应有限时间超螺旋滑模制导律. 系统工程与电子技术. 2021(03): 779-788 . 百度学术
    11. 华艺欣,邹泉,田海铭. 软式自主空中加油控制策略仿真. 北京航空航天大学学报. 2021(02): 262-270 . 本站查看
    12. 余伶俐,况宗旭,王正久,周开军. 智能车辆队列横纵向有限时间滑模控制. 控制理论与应用. 2021(08): 1299-1312 . 百度学术
    13. 文成馀,江驹,余朝军,朱平. 高超声速飞行器快速平滑自适应二阶滑模控制. 哈尔滨工业大学学报. 2020(12): 35-41 . 百度学术

    其他类型引用(15)

  • 加载中
计量
  • 文章访问数:  3621
  • HTML全文浏览量:  112
  • PDF下载量:  969
  • 被引次数: 28
出版历程
  • 收稿日期:  2008-09-17
  • 网络出版日期:  2009-09-30

目录

/

返回文章
返回
常见问答