首页 理论教育多飞行器协同制导与控制的防御弹单向协同制导律

多飞行器协同制导与控制的防御弹单向协同制导律

【摘要】:针对如式所示的目标信息单向传输时的模型,类似地设计协同拦截制导律。例14-1假设目标、防御弹和攻击弹的初始位置、速度和弹道倾角同例13-1中的表13-1,3个飞行器常速飞行,攻击弹和防御弹分别采用Nm=4、Nd=4的增强比例导引律。图14-6不同防御弹制导律时qmd的变化曲线图14-7不同防

针对如式(14-5)所示的目标信息单向传输时的模型,类似地设计协同拦截制导律。令Q=,且考虑如式(14-6)所示的状态变量定义,则式(14-5)可写为

类似式(14-7),在每一时刻,将式(14-24)看作线性定常系统进行制导律设计。

选取性能指标函数为

式中,a3、b3、c3为惩罚系数,a3≥0,b3≥0和c3>0;当a3时,可使得x31在末端时刻趋于零,即qmd,实现终端拦截角约束;当b3时,可使得x32在末端时刻趋于零,即→0,实现防御弹对攻击弹的拦截;c3较大时,防御弹的控制能量较小。

由式(14-24)和式(14-25)可知,哈密顿函数为

式中,λ31和λ32为协态量,其正则方程分别为

λ31和λ32在末端时刻满足的横截条件分别为

式中,φ3=

由式(14-27)和式(14-28)可解得协态量λ31和λ32的表达式为

式中,χ=(tf-t)/rmd≈-1。

由极小值原理可知

联立式(14-26)、式(14-29)和式(14-30),可解得

式中,——防御弹的开环解。

由于λ32包含x31(tf)和x32(tf)项,不能由当前时刻的状态进行求解,因此需要求解x31(tf)和x32(tf)的表达式,将式(14-31)代入式(14-24)的第二个表达式,并从t到tf进行积分,可得

式中,Φ3(·)——状态转移矩阵

因此,x32(t)的表达式可整理为

将式(14-33)代入式(14-24)的第一个表达式,并从t到tf进行积分,可得

同样,将式(14-32)和式(14-34)整理成关于x31(tf)和x32(tf)的二元一次方程:

式中,系数A3、B3、C3、D3、E3和F3分别为

求解式(14-35),可得x31(tf)和x32(tf)为

将式(14-36)代入式(14-37),即可得到x31(tf)和x32(tf)的完整表达式;再将x31(tf)和x32(tf)代入式(14-31),即可得到防御弹加速度闭环解。

考虑防御弹以指定的角度命中攻击弹的要求,设a3和b3,同时忽略系数B3、D3和E3中的小量,则可将B3、D3和E3化简为

将化简后的系数代入式(14-37),可求得x31(tf)和x32(tf)的表达式;再将x31(tf)和x32(tf)代入式(14-31),可得最终的防御弹加速度指令wq

式中,右端第一项体现了对拦截角度的约束;第二项类似于比例导引律体现了成功拦截的要求;第三项是攻击弹在垂直于防御弹-攻击弹视线方向上的加速度。

若只考虑防御弹命中攻击弹的要求而不考虑防御弹攻击角度的约束,则可令a3=0和b3,同时将系数A3、B3、D3和E3中的小量忽略,则可将A3、B3、D3和E3化简为

将式(14-40)代入式(14-37)中的第二式,可得x32(tf)为

同时,由于a3=0,因此可将式(14-31)化简为

将式(14-41)代入式(14-42),即可得到最终的不带攻击角度约束的防御弹加速度指令wq

由式(14-43)可知,此时不带攻击角度约束的防御弹加速度指令与常用的增强比例导引律在表达形式上类似,仅导引系数不同,但两种制导律的设计思路、信息获取模式等都不一样。

例14-1 假设目标、防御弹和攻击弹的初始位置、速度和弹道倾角同例13-1中的表13-1,3个飞行器常速飞行,攻击弹和防御弹分别采用Nm=4、Nd=4的增强比例导引律。目标分别采用=15°的单向协同制导律(简称“制导律1”)、=20°的单向协同制导律(简称“制导律2”)和无拦截角约束的单向协同制导律(简称“制导律3”),三种制导律的参数设置见表14-1。

表14-1 三种制导律的参数设置

仿真结果如图14-1~图14-4所示,图例中的下标“1”“2”和“3”分别表示制导律1、制导律2和制导律3。

图14-1 不同目标制导律时的飞行轨迹

图14-2 不同目标制导律时qmd的变化曲线

图14-3 不同目标制导律时的变化曲线

图14-4 不同目标制导律时控制量的变化曲线

由图可知,当目标采用无角度约束的制导律3时,防御弹成功拦截攻击弹且拦截角度为18.64°;当目标采用具有拦截角约束的制导律1和制导律2时,防御弹-攻击弹的视线角速度趋于零,同时视线角qmd也趋于给定的理想拦截角15°和20°。这是因为,目标已知攻击弹和防御弹所采用的制导律,通过自身机动吸引攻击弹至某一合适位置,使得防御弹以期望的拦截角度对攻击弹进行拦截。由图14-4、图14-1可见,为了实现拦截角度的约束,防御弹的最大加速度较大,弹道较弯曲。

假设目标单向协同制导律中c2=105,107,同时考虑之前c2=1且拦截角约束=15°的情况,令其余参数不变,仿真结果见表14-2。

表14-2 不同c2时制导律1的仿真结果

由表14-2可知,当c2逐渐增大时,目标的加速度逐渐减小,控制能量也逐渐减小。从另一个角度讲,c2大时,目标的机动性弱,此时其无法很好地协助防御弹以指定的角度拦截攻击弹。在c2=1、c2=105和c2=107的情况下,防御弹的拦截角误差分别为0.001 1°、2.015 4°和4.180 8°,随着c2的增大而增大。

例14-2 假设目标以加速度aT=g sin(2πt/5)做蛇形机动。防御弹分别采用单向协同制导律(=30°)和式(13-48)所示的弹道成型制导律(制导阶次为0),仿真结果如图14-5~图14-8所示,图例中的下标“4”“5”分别代表防御弹单向协同制导律和弹道成型制导律。需要说明的是,在图14-5中,当防御弹采用不同制导律时,由于目标和攻击弹的运动规律是相同的,只是飞行时间略有不同,因此T4和T5的轨迹是重合的,M4和M5的轨迹也是重合的,只是D4和D5的轨迹不同。同理,在图14-8中,目标和攻击弹的控制量变化曲线重合。

图14-5 不同防御弹制导律时的飞行轨迹

由图可知,在防御弹独立飞行和考虑目标运动协同飞行两种情况下,其都按照指定的角度实现了对攻击弹的成功拦截,但是和目标协同飞行时的最大加速度要小于独立飞行时的最大加速度,弹道弯曲度更小。

图14-6 不同防御弹制导律时qmd的变化曲线

图14-7 不同防御弹制导律时的变化曲线

图14-8 不同防御弹制导律时控制量的变化曲线

假设防御弹单向协同制导律中的参数a3=1010、b3=1010=30°,c3分别设为1、103和105,则仿真结果见表14-3。

表14-3 不同c3时防御弹单向协同制导律仿真结果

由表可知,随着c3增大,防御弹的最大加速度(绝对值)逐渐减小,控制能量也逐渐减小,攻击角与期望攻击角的偏差逐渐增大。从另一个角度,如果防御弹的机动性很弱,则无法很好地协同目标飞行、以指定的角度拦截攻击弹。