法律状态公告日
法律状态信息
法律状态
2014-07-02
授权
授权
2012-11-28
实质审查的生效 IPC(主分类):F41G3/00 申请日:20120417
实质审查的生效
2012-10-03
公开
公开
技术领域
本发明涉及一种控制多枚导弹攻击角度和攻击时间的方法,属于制导技术 领域。
背景技术
随着科学技术的高速发展,导弹防御系统的能力不断增强,单枚导弹作战 在很多情况下已无法满足战场的要求,多枚导弹协同作战变得越来越重要。为 了提高导弹协同作战的杀伤力,通常对攻击角度有要求,而为了实现饱和攻击, 又对攻击时间有要求。如果多枚导弹能够从不同的角度同时攻击目标,将大大 增强导弹的打击能力。
国内外关于对导弹攻击角度和攻击时间同时进行控制的技术还处于初级阶 段。目前,国内外对于具有攻击角度约束的制导控制方法的研究开展得比较多, 此方面的制导控制方法主要分为两类:基于最优控制理论的方法和基于比例导 引律的方法。对具有攻击时间约束的导弹制导控制方法的研究相对来说比较少, 而同时考虑对攻击角度和攻击时间进行控制的方法就更少。在先技术[1](参见 Lee J I,Jeon I S,Tahk M J.Guidance law to control impact time and angle.IEEE Transactions on Aerospace and Electronic Systems,2007,43(1):301-310)中,基于最 优控制理论提出了可同时控制多枚导弹攻击方位角和攻击时间的制导方法,但 它只适用于导弹在二维平面攻击目标的情况。在先技术[2](参见Jeon I S,Lee J I,Tahk M J.Guidance law to controlimpact time and angle.International Conference on Control and Automation.Hungary:Budapest,2005:852-857)中,基于传统的比例导引 律设计了一种具有攻击角度和时间约束的、适用于反舰导弹的制导方法,但本 方法只适用于攻击固定目标,而且并没有对制导方法参数的取值进行说明。在 先技术[3](参见张友安,马培蓓.带有攻击角度和攻击时间控制的三维制导,航空学 报,2008,29(4):1020-1026)中,先基于李雅普诺夫稳定性理论设计了带有攻击 角度控制的三维制导律,采用此制导律将导弹导引到预定的攻击角度上,然后 通过圆弧机动飞行来对攻击时间进行控制,但此技术中需要预先设定理想攻击时 间,而理想攻击时间的合理设定一直是一个难点。
发明内容
本发明要解决的技术问题是克服上述已有技术的困难,在不需要事先指定 理想攻击时间以及不限定目标运动状态的前提下,提供一种控制多枚导弹攻击 角度和攻击时间的方法。
实现本发明的技术方案如下:
一种控制多枚导弹攻击角度和攻击时间的方法,具体步骤为:
步骤一、在协同作战的多枚导弹中,任意选取一枚作为领弹,并根据武器系 统的总体作战要求,为领弹选取任一制导律;同时以目标为球心,以领弹与目 标之间的距离为半径构造虚拟球体;
步骤二、将剩余的导弹定义为从弹,根据导弹协同作战的需求,参考领弹对 目标的攻击情况,为每一从弹设定理想攻击角度,该攻击角度都是以地面坐标 系为基准设定的;同时在虚拟球体表面为每一从弹选定虚拟点,每一虚拟点与 目标之间连线的方向为其对应从弹的理想攻击角度,其中虚拟点的位置随着领 弹与目标的逼近不断更新;
步骤三、通过控制从弹运动轨迹跟踪并逼近虚拟点轨迹,从而实现对多枚 导弹攻击角度和攻击时间的控制。
进一步地,本发明步骤三的过程为:
步骤101、在地面坐标系中,建立目标的运动模型、从弹的运动模型及虚拟 点位置模型,根据所建立的模型在地面坐标系中建立从弹与虚拟点之间的相对 运动模型;该步骤的具体过程为:
步骤I、在地面坐标系中建立目标的运动模型:
其中XT=[xT,yT,zT]表示目标的位置,VT为目标的速度,θT和ψVT为目标速度矢量 与地面坐标系坐标轴的夹角;
步骤II、由于虚拟点与目标之间的距离等于领弹与目标之间的距离,虚拟 点的运动方向为设定的从弹理想攻击角,因此,在地面坐标系中建立表征虚拟 点位置模型为:
其中:Xv=[xv,yv,zv]为虚拟点的坐标,rl为领弹与目标之间的距离,εy、εz分别 为从弹在俯仰方向和偏航方向设定的理想攻击角;
步骤III、在地面坐标系中建立从弹的运动模型为:
其中:Xm=[xm,ym,zm]表示从弹的位置,V为从弹的速度,θ和ψV为从弹的弹道倾 角和弹道偏角;
步骤IV、根据步骤I-III所建立的模型,在地面坐标系中建立从弹与虚拟点 之间的相对运动模型;
从弹和虚拟点之间的相对位置为:
对式(4)求导可得到从弹与虚拟点的相对速度差:
将式(5)写为向量形式为:
其中:Xd=[xd yd zd]T,T4=[-cosεycosεz -sinεy cosεysinεz]T,领弹与目标之 间的距离的变化率;
步骤102、建立描述从弹运动的质心动力学方程、姿态运动学方程以及姿态 动力学方程,结合以上推导的从弹与虚拟点的相对运动模型,得到对从弹攻击 角度和攻击时间进行控制的非线性控制模型,基于反馈线性化理论对所述非线 性控制模型进行精确线性化,获取线性的控制模型;
步骤103、将线性化后的控制模型看作一个终端时间无限的状态调节器,采 用线性二次型性能指标,基于最优控制理论设计从弹控制器,使控制器输出为 零,利用该设计的从弹控制器实现对虚拟点的跟踪和逼近,从而实现对从弹攻 击角度和攻击时间的控制。
进一步地,本发明当领弹按照指定的制导律飞向目标时,虚拟点轨迹为一 条直线,为了降低对从弹可用过载和控制能量的要求,按照设定的时间步长对 虚拟点位置模型进行改进设计,在地面坐标系中更新虚拟点虚拟点位置模型 为:
其中:rf为从弹与目标之间的距离,εyf、εzf分别为从弹在俯仰方向和偏航方向 的实际视线角,N为整数且其取值方法为:
其中:tgo为预估的领弹飞行时间,h为设定的时间步长;
其中:n为取遍[1,2,…,N]的整数,当领弹飞过一个时间步长时,n取1,当领弹 飞过两个时间步长时,n取2,并依次类推,直至领弹飞过N个时间步长,n取N 为止。
有益效果
本发明在控制攻击时间方面:由于每一虚拟点与领弹处于同于虚拟球上, 虚拟点跟随领弹一步步接近目标,利用从弹对其所对应的虚拟点的跟踪与逼近, 因此在无需事先对每一从弹指定理想攻击时间的情况下,实现对目标攻击时间 的控制,其可适用于攻击固定或不同运动状态的目标;在控制攻击角度方面: 由于每一虚拟点与目标之间连线的方向为其对应从弹的理想攻击角度,在从弹 跟踪与逼近虚拟点的过程中,其相对于目标的攻击角度也逐步逼近理想攻击角 度,从而实现对目标攻击角度的控制;本发明基于虚拟点的方法能够实现对从 弹攻击角度和攻击时间的有效控制,从而使从弹和领弹从预设的角度同时攻击 目标,达到良好的协同作战效果。本发明可应用于各种类型的导弹,具有很大 的应用灵活性和广阔的军事应用前景。
其次、本发明通过对虚拟点轨迹进行改进,其可满足降低从弹的过载要求。
附图说明
图1为本发明的总体示意图;
图2为三枚导弹协同攻击目标弹道图;
图3为两枚从弹攻击角度随时间变化图;
图4为两枚从弹推力随时间变化图;
图5为两枚从弹舵偏角随时间变化图;
图6为从弹1切向过载和法向过载随时间变化图;
图7为从弹2切向过载和法向过载随时间变化图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
当前,在多枚导弹协同作战研究领域中,如何控制各导弹的攻击时间和攻 击角度,来提高对目标的打击力为当前所面临的主要问题。在本发明中:
如图1所示,从多枚导弹中任意选取一枚为领弹,以目标为球心、领弹与 目标之间的距离为半径构造虚拟球体;将剩余的导弹定义为从弹,在虚拟球体 表面为每一从弹选定虚拟点,其中每一虚拟点与目标之间连线的方向为其对应 从弹的预设的理想攻击角度;通过控制从弹运动轨迹跟踪并逼近虚拟点轨迹, 从而实现对多枚导弹攻击角度和攻击时间的控制。
本发明随着领弹与目标的逼近使得虚拟球体的半径逐步缩小,通过控制从 弹运动轨迹对虚拟点的跟踪与逼近,从而达到对从弹攻击时间进行控制的目的; 同时由于各虚拟点与目标之间连线的方向为其对应从弹的预设的理想攻击角 度,在从弹的跟踪与逼近的过程中可以保证对从弹攻击角度的控制目的;因此 本发明采用构造虚拟点的方法,其很好的实现对从弹攻击角度和攻击时间的控 制。
对使用符号的注释:本发明中凡是出现符号上面带“·”的表达式,其采用 课本上的表达式,即为对该符号的求导,例如下文中为对xm对时间的求导; 凡是出现符号上带有“··”的表达式,其也采用课本上的表达式,即对该符号 进行两次求导,例如下文出现的,其表示对rl对时间进行两次求导。
本发明控制多枚导弹攻击角度和攻击时间的方法,具体的过程为:
步骤一、在协同作战的多枚导弹中,任意选取一枚作为领弹,并根据武器系 统的总体作战要求,为领弹选取任一制导律;如果强调领弹对目标的命中精度 和领弹弹道的收敛性,则可选用比例导引律;如果强调领弹的攻击角度和对目 标的毁伤效果,可选用具有角度约束的制导律;由于本发明并未对领弹采用的 制导律有限制,因此,领弹可根据需要选用任一制导律。同时以目标为球心, 以领弹与目标之间的距离为半径构造虚拟球体。
步骤二、将剩余的导弹定义从弹,根据导弹协同作战的需求,参考领弹对目 标的攻击情况,为每一从弹设定理想的攻击角度εyi和εzi,从而实现领弹和从弹 从不同的方向同时命中目标,对目标进行有效打击;其中εyi和εzi分别为从弹在 俯仰方向和偏航方向设定的理想攻击角(为了描述方便,以下省去表示从弹序 号的i);同时,在虚拟球体表面为每一从弹选定虚拟点,其中每一虚拟点与目 标之间连线的方向为其对应从弹的预设的理想攻击角度;该攻击角度都是以地 面坐标系为基准设定的。
本发明设定从弹理想的攻击角度的具体过程为:根据目标的运动情况、领弹 的初始状态、领弹采取的导引律等信息先估算出领弹的攻击角度(如领弹采用 带有攻击角度约束的制导律,则可直接得知领弹的攻击角度),再根据导弹协同 作战的需求,为每一从弹设定理想的攻击角度。例如:估算出的领弹的攻击角 为(-2°,-25°)(括号内第一项指俯仰方向攻击角,第二项指偏航方向攻击角),为 了实现对目标的立体多方向攻击,可设从弹的攻击角为(-5°,-75°)、(-10°,25°)、 (-15°,75°)等。
步骤三、通过控制从弹运动轨迹跟踪并逼近虚拟点轨迹,从而实现对多枚 导弹攻击角度和攻击时间的控制。
下面对步骤三的具体过程进行详细说明:
在地面坐标系中,建立目标的运动模型、从弹的运动模型及虚拟点位置模 型,根据所建立的模型在地面坐标系中建立从弹与虚拟点之间的相对运动模型; 该步骤的具体过程为:
步骤I、在地面坐标系中建立目标的运动模型:
其中:XT=[xT,yT,zT]表示目标的位置,VT为目标的速度,θT和ψVT为目标速度矢 量与地面坐标系坐标轴的夹角。
步骤II、由于虚拟点与目标之间的距离等于领弹与目标之间的距离,虚拟 点的运动方向为设定的从弹理想攻击角,因此,在地面坐标系中建立表征虚拟 点位置模型为:
其中:Xv=[xv,yv,zv]为虚拟点的坐标,rl为领弹与目标之间的距离,εy、εz分别 为从弹在俯仰方向和偏航方向设定的理想攻击角。
本步骤中,当领弹按照指定的制导律飞向目标时,虚拟点轨迹为一条直线, 当基于控制器对从弹进行控制,使其逼近和跟踪此虚拟点轨迹时,要求从弹具 有较大的过载和控制能量。为了降低对从弹可用过载和控制能量的要求,对虚 拟点轨迹进行改进设计,改进后的虚拟点轨迹为:
其中:rf为从弹与目标之间的距离,εyf、εzf分别为从弹的实际攻击角,N为整 数且其取值方法为:
其中:tgo为预估的领弹飞行时间,h为设定的时间步长;
其中:n为取遍[1,2,…,N]的整数,当领弹飞过一个时间步长时,n取1,当领弹 飞过两个时间步长时,n取2,并依次类推,直至领弹飞过N个时间步长,n取N 为止。
采用上述方法改进后的虚拟点轨迹先弯曲再平直,从而降低了对从弹需用 过载及控制能量的要求。
步骤III、在地面坐标系中建立从弹的运动模型为:
其中:Xm=[xm,ym,zm]表示从弹的位置,V为从弹的速度,θ和ψV为从弹的弹道倾 角和弹道偏角。
步骤IV、根据步骤I-III所建立的模型,在地面坐标系中建立从弹与虚拟点 之间的相对运动模型;
从弹和虚拟点之间的相对位置为:
对上式求导可得到从弹与虚拟点的相对速度差:
将其写为向量形式为:
其中Xd=[xd yd zd]T,T4=[-cosεycosεz -sinεy cosεysinεz]T,领弹与目标之间 的距离的变化率。
步骤102、建立描述从弹运动的质心动力学方程、姿态运动学方程以及姿态 动力学方程,结合以上推导的从弹与虚拟点的相对运动模型,得到对从弹攻击 角度和攻击时间进行控制的非线性控制模型,基于反馈线性化理论对所述非线 性控制模型进行精确线性化,获取线性控制模型。
以下对本步骤的具体过程进行说明。
(I)将描述虚拟点和从弹相对运动的方程和从弹的质心动力学方程结合起 来,进行精确线性化。
从弹的质心动力学方程为:
其中:m为从弹的质量,a为从弹的攻角,β为从弹的侧滑角,γV为从弹的倾斜 角,g为重力加速度,P为从弹发动机的推力,X、Y、Z分别是从弹的阻力、 升力和侧向力,其表达式为:
其中:cx、cy、cz分别是从弹的阻力系数、升力系数和侧向力系数,其可以根据 风洞试验获取;是从弹的动压,ρ是大气密度,S是从弹的特征面积。
将推力P作为控制量,然后将式(9)变形为:
为了表述方便,令
则,式(11)结合第(6)式,得到
接下来基于反馈线性化理论对式(12)进行精确线性化。设输出为控制模 型中从弹与虚拟点的相对距离,即y=Xd=[xd yd zd]T,对控制模型中质心动力 学方程进行微分同胚坐标变换,有:
其中:表示求0阶李导数,表示求1阶李导数,x=[xd,yd,zd,V,θ,ψV]T,
其中:
由式(15)可看出,对Z1求二阶导数可显现控制量P。在式(15)中,α、 β、γV为根据导弹的姿态角和几何关系方程求出,在此将其看作由式(15)表示 的系统的输入。由式(15)可知,输出Z1到P、α、β、γV的向量关系度为 {2,2,2,2},则系统关系度r=2+2+2+2=8。而质心动力学系统的状态变量 有六个,即系统状态维数n=6。系统维数不等于系统关系 度即n≠r,为了实现精确线性化,可将P进行二阶动态扩展,从而实现关系 度。令:
其中,ζ1、ζ2是动态扩展后引入的辅助状态,v1是引入的新的控制量。将式(15) 与式(17)合并可得:
(II)将表征从弹姿态运动学方程和姿态动力学方程进行反馈线性化。
导弹的姿态运动学方程为:
其中:ψ、γ分别为从弹的俯仰角、偏航角和倾斜角,ωx、ωy、ωz为建立 在从弹上的弹体坐标系相对地面坐标系的转动角速度在弹体坐标系各轴上的分 量。
将式(19)变形为:
导弹的姿态动力学方程为
其中:Jx、Jy、Jz分别为从弹对于弹体坐标系(即惯性主轴系)各轴的转动惯量。 Mx、My、Mz为作用在从弹上的所有外力对质心的力矩在弹体坐标系各轴上的 分量。
气动力矩的计算公式为
式(23)中,L为特征长度,mx、my、mz分别是滚转力矩系数、偏航力矩 系数和俯仰力矩系数,其可以根据风洞试验获取。
且有
式(24)中,δx、δy和δz分别为副翼、方向舵和升降舵的偏转角,也是整个 控制模型的三个输入。mx0、mz0为由于从弹外形不对称引起的滚动力矩系数和俯 仰力矩系数。为气动力矩系数导 数,例如是mx相对于β的导数。
将式(23)和(24)带入式(22),并变换形式为:
对式(20)和(21)求导并将式(25)、(26)带进去,得到姿态角加速度 方程为:
式(27)和(28)中,
接下来对导弹姿态运动学方程和导弹姿态动力学方程进行线性化,为了实 现俯仰和偏航运动的解耦控制,将滚转通道的运动独立出来进行分析。表征导 弹俯仰和偏航通道运动的方程为式(20)和式(25)即:
为了表述方便,令
接下来对由式(29)和(30)表征的导弹姿态运动非线性方程组进行线性化, 其输出为
对系统进行微分同胚坐标变换,并引入新的状态变量,结果为:
其中:为了简便将简写为fp=[F3 F4′]T,对式(32)进 行求导得:
表征导弹滚转通道的运动动力学方程组为式(21)和(26),即
为了表述方便,令
设滚转通道的输出为γ,类似于俯仰通道和偏航通道,对系统进行微分同 胚坐标变换并引入新的状态,变换结果为:
其中:fγ=[F5 F6′]T。
为了实现系统的关系度,将俯仰和偏航方向的姿态运动与进行了二阶动态 扩展的辅助动态进行合并,结合式(33)和(18),合并后的新的状态为:
此时Z3和Z4分别包含3个状态变量,对其求导可得
其中:
(III)将分步线性化后的从弹与虚拟点相对运动模型、从弹质心动力学方 程、从弹姿态运动学方程、从弹姿态动力学方程进行合并,得到:
此时,在式(41)表示的线性方程组中的新状态为:
从弹的输入为u=[P,δz,δy,δx],输出为y=[xd yd zd γ]T。引入等效输入:
且令:
则式(41)化为更直观的形式:
此时可知,系统的输出为yo=[Z1 Z5]T=[xd yd zd γ]T,系统线性化前的控 制输入为u=[P δx δy δz]T,系统线性化后的等效输入为o=[o1 o2]T,系统的 状态为Z=[Z1 Z2 Z3 Z4 Z5 Z6]T,共14个状态变量即系统的维数n=14。对系 统输出Z1求四阶导数,可显现控制量v1、δy、δz。对系统输出Z5求二阶导数, 可显现控制量δx。所以输出Z1到v1、δy、δz的向量关系度为{4,4,4},输出Z5到 δx的关系度为2。从而系统的关系度为r=4+4+4+2=14,则有系统的维数等于 系统的关系度即n=r=14,因此,考虑从弹非线性动力学特性的控制模型可进 行精确线性化。
(IV)对式(43)表示的综合考虑了从弹与虚拟点相对运动模型、从弹的 质心动力学方程、从弹姿态运动学方程和从弹姿态动力学方程的控制模型进行 整体精确线性化。针对式(43)表示的系统,其状态为fz,由于滚转通道独立于 俯仰和偏航通道求解,因此,将输出分解为
y1=[xd yd zd]T (44)
y2=γ
然后,针对两部分输出分别进行微分同胚坐标变换,可得到新的状态变量 为L=[L1,L2,L3,L4,L5,L6]T,其中:
令式(45)所示系统的等效输入为q=[q1 q2]T 其中
则线性化后的方程组即式(45)可写成状态方程为
其中,系统的状态矩阵为:
系统的控制输入矩阵为:
系统的输出矩阵为:
步骤103、将线性化后的控制模型看作一个终端时间无限的状态调节器,采 用线性二次型性能指标,基于最优控制理论设计从弹控制器,使控制器输出为 零,利用该设计的从弹控制器实现对虚拟点的跟踪和逼近,从而实现对从弹攻 击角度和攻击时间的控制。
将线性化后的控制模型看作一个终端时间无限的状态调节器,此状态调节 器问题的线性二次型性能指标为
由于控制终点所以最优控制性能指标为:
在求解最优控制时,矩阵Q取常对角半正定阵,R取常对角正定阵。针对此 状态调节器问题,基于最优控制理论,求得使式(50)所示的二次型性能指标 的最小的解为:
其中K为最优控制黎卡提方程的解。求得等效控制输入q=[q1 q2]T后,再带入 式(46)和(47),求出o1、o2,再根据式(42)可求得[v1 δz δy δx]T,对v1进 行两次积分,即得到推力P,至此,得到了对从弹的控制指令,从而实现对从 弹攻击角度和攻击时间的控制。
以下对基于虚拟点轨迹设计的控制方法的验证。
以三枚导弹(包括一枚领弹和两枚从弹)协同攻击一个高价值指挥所为例, 考虑到三枚导弹协同作战时,既能保证对目标的命中精度,又能从不同方向同 时到达目标,设领弹采用传统的比例导引律,从弹1的理想攻击角度 εy1=-20°、εz1=40°,从弹2的理想攻击角度为εy2=-20°、εz2=-85°。考虑到导弹的 自身的发动机能力限制、操纵机构限制以及过载限制,设控制量满足约束: P≤10000N、|δz|≤15°、|δy|≤15°,切向过载和法向过载均不超过10。假设各枚导弹发 射后先爬升,然后转入平飞,领弹在距目标10km时发现了目标,三枚导弹一 齐转入末制导对目标进行协同攻击。由于领弹采用比例导引律攻击目标,从弹 采用本发明涉及的控制方法进行控制,因此,主要关注从弹的运动参数变化情 况。末制导阶段,领弹和从弹均在48.7秒时命中目标(弹目距离均小于1.2米), 领弹和从弹能够保持倾斜稳定,领弹和从弹的弹道以及从弹的各特征参数变化 如图2-7所示。
机译: 具有拒绝服务攻击和分布式拒绝服务攻击功能的网络接口卡,以及一种通过使用网络接口卡来拒绝拒绝服务攻击和分布式服务攻击的方法,特别是与负载有关的一个数字电脑系统
机译: 接触相检测方法,攻击角度测量方法,接触相位检测装置和攻击角度测量装置
机译: 可以通过预测攻击消息转送时间来针对安全控制器区域网络的消息洪水攻击的缓解方法