CN104266546A - 一种基于视线的有限时间收敛主动防御制导控制方法 - Google Patents

Abstract

一种基于视线的有限时间收敛主动防御制导控制方法，涉及一种制导控制方法，特别是涉及一种主动防御制导控制方法。为了解决防御导弹过载能力受限的问题。本发明首先对目标、防御导弹和拦截导弹相对运动建模，采用视线制导方式为防御导弹设计制导律，然后采用非奇异终端滑模控制来设计制导律，对纵向平面和侧向平面分别定义滑模变量为并对其求导，将目标、防御导弹和拦截导弹相对运动方程带入并整理后得到纵向平面和侧向平面的制导律和，并按制导律对导弹进行控制，本发明能有效降低防御导弹需用过载。本发明适用于主动防御制导控制。

Description

一种基于视线的有限时间收敛主动防御制导控制方法

技术领域

本发明涉及一种制导控制方法，特别是涉及一种主动防御制导控制方法。

背景技术

随着导弹防御系统技术的完善，弹道导弹突防面临巨大挑战。弹道导弹在大气层外的中段飞行段由于机动能力很小，弹道固定容易被导弹防御系统拦截。为了避免被导弹防御系统拦截，传统突防的弹道导弹采用了被动逃避的方式，例如电子干扰，隐身技术，释放诱饵，中段制导，多弹头攻击和机动变轨等多种方式。与此相对的是采用主动防御方法，即弹道导弹释放防御导弹，之后防御导弹与其伴飞，拦截导弹防御系统发射的拦截导弹。主动防御方法也可以结合以上提到的被动防御方法，进一步提高弹道导弹的突防概率。为了叙述简便，将突防的弹道导弹称为目标，伴飞的防御导弹称为防御导弹，导弹防御系统发射的拦截目标的导弹称为拦截导弹。整个制导场景分为拦截导弹拦截目标和防御导弹拦截拦截导弹。在导弹末制导技术中，应用最广的是基于零化视线角速率的PN制导律。但是使用PN制导律的防御导弹要成功拦截机动目标需要很大的过载能力。花文化等人在《飞行力学》，2012年第1期“主动防御协同自适应滑模制导律”提出了一种飞机防御的制导律。采用了传统滑模控制方法。飞机防御问题属于大气层内制导问题，防御导弹需用过载大，这和大气层外防御导弹机动过载能力较小不同，作者并未验证其制导律是否能直接应用到大气层外的弹道导弹防御场景中。针对以上方法存在的问题，需要设计一种能降低防御导弹需用过载的鲁棒的制导律。

发明内容

本发明为了解决防御导弹过载能力受限的问题，进而提出了一种基于视线的有限时间收敛主动防御制导控制方法。

一种基于视线的有限时间收敛主动防御制导控制方法，包括下述步骤：

步骤1：拦截导弹-目标和防御导弹-拦截导弹相对运动建模：

目标、防御导弹和拦截导弹的速度用V_t，V_d和V_m表示；假设拦截导弹和防御导弹的控制均可解耦为纵向平面和侧向平面独立制导；

拦截导弹和目标的相对运动方程建模为：

${\stackrel{\·\·}{q}}_{\mathrm{\ϵmt}}=-2\frac{{\stackrel{\·}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\·}{q}}_{\mathrm{\ϵmt}}+\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\ϵ}}}-\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\ϵ}}}---\left(1\right)$

${\stackrel{\·\·}{q}}_{\mathrm{\βmt}}=-2\frac{{\stackrel{\·}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\·}{q}}_{\mathrm{\βmt}}-\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\β}}}+\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\β}}}---\left(2\right)$

其中，R_mt为拦截导弹和目标的距离，为R_mt的一阶导数；q_εmt和q_βmt分别为拦截导弹-目标的视线高低角和方位角，相应的视线角速度为和和分别为拦截导弹-目标的视线高低角和方位角的二阶导数；和分别为目标加速度在拦截导弹-目标视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在拦截导弹-目标视线坐标系下垂直视线的分量；

防御导弹和拦截导弹相对运动建模为：

${\stackrel{\·\·}{q}}_{\mathrm{\ϵdm}}=-2\frac{{\stackrel{\·}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\·}{q}}_{\mathrm{\ϵdm}}+\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\ϵ}}}-\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\ϵ}}}---\left(3\right)$

${\stackrel{\·\·}{q}}_{\mathrm{\βdm}}=-2\frac{{\stackrel{\·}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\·}{q}}_{\mathrm{\βdm}}-\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\β}}}+\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\β}}}---\left(4\right)$

其中，R_dm为防御导弹和拦截导弹的距离，为R_dm的一阶导数；式中q_εdm和q_βdm分别为防御导弹-拦截导弹的视线高低角和方位角，相应视线角速度为和和分别为防御导弹-拦截导弹的视线高低角和方位角的二阶导数；和分别为防御导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；

步骤2：视线制导方式建模：

采用视线制导方式为防御导弹设计制导律；在视线制导方式下，防御导弹通过控制将自身的位置置于拦截导弹和目标的视线上；随着拦截导弹不断接近目标，最终将撞上防御导弹；视线制导满足如下等式

q_εdm＝-q_εmt,q_βdm＝q_βmt   (5)

${\stackrel{\·}{q}}_{\mathrm{\ϵdm}}=-{\stackrel{\·}{q}}_{\mathrm{\ϵmt}},{\stackrel{\·}{q}}_{\mathrm{\βdm}}={\stackrel{\·}{q}}_{\mathrm{\βmt}}---\left(6\right)$

在制导过程中，防御导弹-拦截导弹的视线角和视线角速度并不能精确满足公式(5)和(6)，定义防御导弹对拦截导弹的视线角和视线角速度误差为：

$x_{{\mathrm{\ϵ}}_1}=q_{\mathrm{\ϵdm}}+q_{\mathrm{\ϵmt}},x_{{\mathrm{\β}}_1}=q_{\mathrm{\βdm}}-q_{\mathrm{\βmt}}---\left(7\right)$

$x_{{\mathrm{\ϵ}}_2}={\stackrel{\·}{q}}_{\mathrm{\ϵdm}}+{\stackrel{\·}{q}}_{\mathrm{\ϵmt}},x_{{\mathrm{\β}}_2}={\stackrel{\·}{q}}_{\mathrm{\βdm}}-{\stackrel{\·}{q}}_{\mathrm{\βmt}}---\left(8\right)$

其中 $x_{{\mathrm{\ϵ}}_2}={\stackrel{\·}{x}}_{{\mathrm{\ϵ}}_1},x_{{\mathrm{\β}}_2}={\stackrel{\·}{x}}_{{\mathrm{\β}}_1};$

步骤3：定义滑模变量：

采用非奇异终端滑模控制来设计制导律；对纵向平面和侧向平面分别定义滑模变量为σ₁和σ₂，表达式如下

${\mathrm{\σ}}_1={\mathrm{\β}}_1{x}_{\mathrm{\ϵ}2}^{{\mathrm{\α}}_1}+x_{\mathrm{\ϵ}1}---\left(9\right)$

${\mathrm{\σ}}_2={\mathrm{\β}}_2{x}_{\mathrm{\β}2}^{{\mathrm{\α}}_2}+x_{\mathrm{\β}1}---\left(10\right)$

其中，α₁,α₂是两个常数，α₁和α₂均为两个正奇数的商且满足1<α₁<2，1<α₂<2；β₁,β₂是两个常数，β₁>0，β₂>0；根据非奇异终端滑模控制方法的特点，后续设计的制导律可以保证状态x_ε1和x_β1在有限时间内收敛到0；

对公式(9)和(10)求导得到滑模变量的时间导数为

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{edm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}})+x_{\mathrm{\&epsiv;}2}---\left(11\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;mt}})+x_{\mathrm{\&beta;}2}---\left(12\right)$

将相对运动方程(1)，(2)，(3)和(4)代入式(11)和(12)整理得

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{x_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+\frac{u_{m_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}-\frac{u_{d_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{mt}}})---\left(13\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}+\frac{x_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-\frac{u_{m_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{d_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&beta;}}}}{R_{\mathrm{mt}}})---\left(14\right)$

其中 $u_{m_{\mathrm{q\&epsiv;}}}=u_{m_{2\mathrm{q\&epsiv;}}}-\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&epsiv;}}}$ 和 $u_{m_{\mathrm{q\&beta;}}}=u_{m_{2\mathrm{q\&beta;}}}+\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&beta;}}};$ 在制导过程中有R_mt≈R_dm；所设计的制导律使得目标、防御导弹和拦截导弹在一条直线上，有如下关系和因此和取值很小；因为和是拦截导弹的在拦截导弹-目标视线系和防御导弹-拦截导弹视线系下的复合加速度，因此防御导弹需用过载相应也会很小，这有助于降低防御导弹的需用过载；

步骤4：设计制导律并按制导律对导弹进行控制：

和是有限的值，设和M₁、M₂为常数，具体值和拦截导弹的过载能力有关，M₁、M₂大于零，小于等于拦截导弹最大过载能力；假定目标的机动可以忽略不计；

纵向平面和侧向平面的制导律分别设计如下

$u_{d_{\mathrm{q\&epsiv;}}}=-2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+{\mathrm{\&rho;}}_1\mathrm{sign}\left({\mathrm{\&sigma;}}_1\right)---\left(15\right)$

$u_{d_{\mathrm{q\&beta;}}}=2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}-\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-{\mathrm{\&rho;}}_2\mathrm{sign}\left({\mathrm{\&sigma;}}_2\right)---\left(16\right)$

其中切换项增益ρ₁>M₁，ρ₂>M₂；sign(·)是符号函数；

按制导律(15)(16)对导弹进行控制。

本发明使用了基于视线的制导方式并采用了非奇异终端滑模控制理论，能有效降低防御导弹需用过载。PN制导律的平均脱靶量为0.24米时，过载能力为5g；而本发明制导律的平均脱靶量为0.2米时，过载能力2.5g。

附图说明

图1基于视线的有限时间收敛主动防御制导控制方法流程图；

图2目标、防御导弹和拦截导弹三者相对运动示意图；

图3目标、拦截导弹和防御导弹弹道轨迹；

图4防御导弹对拦截导弹的视线角速率；

图5目标、拦截导弹和防御导弹三点共线情况图。

具体实施方式

具体实施方式一：结合图1说明本实施方式，一种基于视线的有限时间收敛主动防御制导控制方法，包括下述步骤：

步骤1：拦截导弹-目标和防御导弹-拦截导弹相对运动建模：

目标、防御导弹和拦截导弹三者相对运动关系如图2所示。目标、防御导弹和拦截导弹的速度用V_t，V_d和V_m表示；假设拦截导弹和防御导弹的控制均可解耦为纵向平面和侧向平面独立制导；

拦截导弹和目标的相对运动方程建模为：

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\&epsiv;}}}-\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&epsiv;}}}---\left(17\right)$

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}-\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\&beta;}}}+\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&beta;}}}---\left(18\right)$

其中，R_mt为拦截导弹和目标的距离，为R_mt的一阶导数；q_εmt和q_βmt分别为拦截导弹-目标的视线高低角和方位角，相应的视线角速度为和和分别为拦截导弹-目标的视线高低角和方位角的二阶导数；和分别为目标加速度在拦截导弹-目标视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在拦截导弹-目标视线坐标系下垂直视线的分量；

防御导弹和拦截导弹相对运动建模为：

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}+\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\&epsiv;}}}-\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\&epsiv;}}}---\left(19\right)$

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\&beta;}}}+\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\&beta;}}}---\left(20\right)$

其中，R_dm为防御导弹和拦截导弹的距离，为R_dm的一阶导数；式中q_εdm和q_βdm分别为防御导弹-拦截导弹的视线高低角和方位角，相应视线角速度为和和分别为防御导弹-拦截导弹的视线高低角和方位角的二阶导数；和分别为防御导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；

步骤2：视线制导方式建模：

采用视线制导方式为防御导弹设计制导律；在视线制导方式下，防御导弹通过控制将自身的位置控制在拦截导弹和目标的视线上；随着拦截导弹不断接近目标，最终将撞上防御导弹；视线制导满足如下等式

q_εdm＝-q_εmt,q_βdm＝q_βmt   (21)

${\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}=-{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}},{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}---\left(22\right)$

在制导过程中，防御导弹和拦截导弹的视线角和视线角速度并不能精确满足公式(21)和(22)，定义防御导弹对拦截导弹的视线角和视线角速度误差为：

$x_{{\mathrm{\&epsiv;}}_1}=q_{\mathrm{\&epsiv;dm}}+q_{\mathrm{\&epsiv;mt}},x_{{\mathrm{\&beta;}}_1}=q_{\mathrm{\&beta;dm}}-q_{\mathrm{\&beta;mt}}---\left(23\right)$

$x_{{\mathrm{\&epsiv;}}_2}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}+{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}},x_{{\mathrm{\&beta;}}_2}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}---\left(24\right)$

其中 $x_{{\mathrm{\&epsiv;}}_2}={\stackrel{\&CenterDot;}{x}}_{{\mathrm{\&epsiv;}}_1},x_{{\mathrm{\&beta;}}_2}={\stackrel{\&CenterDot;}{x}}_{{\mathrm{\&beta;}}_1};$

步骤3：定义滑模变量：

采用非奇异终端滑模控制来设计制导律；对纵向平面和侧向平面分别定义滑模变量为σ₁和σ₂，表达式如下

${\mathrm{\&sigma;}}_1={\mathrm{\&beta;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1}+x_{\mathrm{\&epsiv;}1}---\left(25\right)$

${\mathrm{\&sigma;}}_2={\mathrm{\&beta;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2}+x_{\mathrm{\&beta;}1}---\left(26\right)$

其中，α₁,α₂是两个常数；β₁,β₂是两个常数；

对公式(25)和(26)求导得到滑模变量的时间导数为

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{edm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}})+x_{\mathrm{\&epsiv;}2}---\left(27\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;mt}})+x_{\mathrm{\&beta;}2}---\left(28\right)$

将相对运动方程(17)，(18)，(19)和(20)代入式(27)和(28)整理得

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{x_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+\frac{u_{m_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}-\frac{u_{d_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{mt}}})---\left(29\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}+\frac{x_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-\frac{u_{m_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{d_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&beta;}}}}{R_{\mathrm{mt}}})---\left(30\right)$

其中 $u_{m_{\mathrm{q\&epsiv;}}}=u_{m_{2\mathrm{q\&epsiv;}}}-\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&epsiv;}}}$ 和 $u_{m_{\mathrm{q\&beta;}}}=u_{m_{2\mathrm{q\&beta;}}}+\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&beta;}}};$ 在在制导中有R_mt≈R_dm；所设计的制导律使得目标、防御导弹和拦截导弹在一条直线上，有如下关系和因此和取值很小；因为和是拦截导弹的在拦截导弹-目标视线系和防御导弹-拦截导弹视线系下的复合加速度，因此防御导弹需用过载相应也会很小，这有助于降低防御导弹的需用过载；

步骤4：设计制导律并按制导律对导弹进行控制：

和是有限的值，设和M₁、M₂为常数，具体值和拦截导弹的过载能力有关，M₁、M₂大于零，小于等于拦截导弹最大过载能力；假定目标的机动可以忽略不计；

纵向平面和侧向平面的制导律分别设计如下

$u_{d_{\mathrm{q\&epsiv;}}}=-2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+{\mathrm{\&rho;}}_1\mathrm{sign}\left({\mathrm{\&sigma;}}_1\right)---\left(31\right)$

$u_{d_{\mathrm{q\&beta;}}}=2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}-\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-{\mathrm{\&rho;}}_2\mathrm{sign}\left({\mathrm{\&sigma;}}_2\right)---\left(32\right)$

其中切换项增益ρ₁>M₁，ρ₂>M₂；sign(·)是符号函数；

按制导律(31)和(32)对导弹进行控制。

具体实施方式二：本实施方式中步骤3中的

其它步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式中步骤3中的β₁＝β₂＝40。

其它步骤及参数与具体实施方式一相同。

具体实施方式四：本实施方式中，步骤4中的M₁、M₂大于零，小于等于拦截导弹最大过载能力，一般拦截导弹最大过载能力小于100m/s²，取0<M₁<100、0<M₂<100。

其它步骤及参数与具体实施方式三相同。

具体实施方式五：本实施方式中，步骤4中根据拦截导弹的过载能力，取M₁＝M₂＝30。

其它步骤及参数与具体实施方式四相同。

实施例

目标在惯性系下的初始位置为(500.0,400.0,0.0)公里，初始速度大小为7071m/s，初始弹道倾角θ_t0＝-35.3°，初始弹道偏角ψ_vt0＝0°。拦截导弹在惯性系下的初始位置为(559.371,319.188,0.746)公里，初始速度大小为3000m/s，初始弹道倾角θ_m0＝44°，初始弹道偏角ψ_vm0＝0°。防御导弹在惯性系下的初始位置为(502.903,396.206,-1.98)公里，初始速度大小为7071m/s，初始弹道倾角θ_d0＝-35°，初始弹道偏角ψ_vd0＝-2.1°。拦截导弹的最大过载为3g，其中g是重力加速度，取值为9.8m/s²。制导律参数取值如下α₁＝α₂＝9/7,β₁＝β₂＝40,M₁＝M₂＝30，ρ₁＝ρ₂＝31。

在本次仿真算例中，防御导弹的脱靶量为0.11米，达到碰撞杀伤的要求。目标、拦截导弹和防御导弹三者的弹道轨迹如图3所示。从图中可以看到目标和防御导弹的弹道轨迹很接近，实际上在整个制导期间，防御导弹一直在目标的右下方伴飞且之间的距离不小于5.7公里，因此防御导弹不会和目标碰撞。

防御导弹对拦截导弹的视线角速率变化情况如图4所示。视线高低角速度和方位角速度除制导末期外都在[-0.05,0.05]deg/s之内，很好的零化了视线角速率。

目标、拦截导弹和防御导弹三点共线情况如图5所示。角度误差q_εdm+q_εmt和q_βdm-q_βmt分别在[0,-0.4]deg和[-2.1,-1.8]deg范围内变化，三点共线情况良好。

本发明制导律称为ADNTSMG(Active Defense Nonsingular Terminal Sliding ModeGuidance)。为了验证ADNTSMG制导律能有效的降低防御导弹的需用过载，将ADNTSMG制导律和PN制导律进行仿真对比。仿真了防御导弹最大过载从2g变化到5g的情况，每种情况进行了200次蒙特卡洛仿真。拦截导弹和防御导弹的初始位置在前面仿真算例给出的初始位置基础上分别增加了2范数为3.5km和5km的随机飘移。其他初始条件与仿真算例的初始条件一样。仿真统计了平均脱靶量，脱靶量小于0.3米和小于0.5米的概率。仿真结果如表1所示。从表1中可以看出在达到相同脱靶量的情况下，ADNTSMG需要的最大过载比PN制导律要小。

表1ADNTSMG和PN制导律性能比较

Claims (4)

1.一种基于视线的有限时间收敛主动防御制导控制方法，其特征在于包括下述步骤：

步骤1：拦截导弹-目标和防御导弹-拦截导弹相对运动建模：

目标、防御导弹和拦截导弹的速度用V_t，V_d和V_m表示；假设拦截导弹和防御导弹的控制均可解耦为纵向平面和侧向平面独立制导；

拦截导弹和目标的相对运动方程建模为：

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\&epsiv;}}}-\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&epsiv;}}}---\left(1\right)$

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}-\frac{1}{R_{\mathrm{mt}}}{u}_{t_{\mathrm{q\&beta;}}}+\frac{1}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&beta;}}}---\left(2\right)$

其中，R_mt为拦截导弹和目标的距离，为R_mt的一阶导数；q_εmt和q_βmt分别为拦截导弹-目标的视线高低角和方位角，相应的视线角速度为和和分别为拦截导弹-目标的视线高低角和方位角的二阶导数；和分别为目标加速度在拦截导弹-目标视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在拦截导弹-目标视线坐标系下垂直视线的分量；

防御导弹和拦截导弹相对运动建模为：

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}+\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\&epsiv;}}}-\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\&epsiv;}}}---\left(3\right)$

${\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}=-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-\frac{1}{R_{\mathrm{dm}}}{u}_{m_{2\mathrm{q\&beta;}}}+\frac{1}{R_{\mathrm{dm}}}{u}_{d_{\mathrm{q\&beta;}}}---\left(4\right)$

其中，R_dm为防御导弹和拦截导弹的距离，为R_dm的一阶导数；式中q_εdm和q_βdm分别为防御导弹-拦截导弹的视线高低角和方位角，相应视线角速度为和和分别为防御导弹-拦截导弹的视线高低角和方位角的二阶导数；和分别为防御导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；和分别为拦截导弹加速度在防御导弹-拦截导弹视线坐标系下垂直视线的分量；

步骤2：视线制导方式建模：

采用视线制导方式为防御导弹设计制导律；视线制导满足如下等式

q_εdm＝-q_εmt,q_βdm＝q_βmt   (5)

${\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}=-{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}},{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}---\left(6\right)$

定义防御导弹对拦截导弹的视线角和视线角速度误差为：

$x_{{\mathrm{\&epsiv;}}_1}=q_{\mathrm{\&epsiv;dm}}+q_{\mathrm{\&epsiv;mt}},x_{{\mathrm{\&beta;}}_1}=q_{\mathrm{\&beta;dm}}-q_{\mathrm{\&beta;mt}}---\left(7\right)$

$x_{{\mathrm{\&epsiv;}}_2}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}+{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}},x_{{\mathrm{\&beta;}}_2}={\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}---\left(8\right)$

其中 $x_{{\mathrm{\&epsiv;}}_2}={\stackrel{\&CenterDot;}{x}}_{{\mathrm{\&epsiv;}}_1},x_{{\mathrm{\&beta;}}_2}={\stackrel{\&CenterDot;}{x}}_{{\mathrm{\&beta;}}_1};$

步骤3：定义滑模变量：

采用非奇异终端滑模控制来设计制导律；对纵向平面和侧向平面分别定义滑模变量为σ₁和σ₂，表达式如下

${\mathrm{\&sigma;}}_1={\mathrm{\&beta;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1}+x_{\mathrm{\&epsiv;}1}---\left(9\right)$

${\mathrm{\&sigma;}}_2={\mathrm{\&beta;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2}+x_{\mathrm{\&beta;}1}---\left(10\right)$

其中，α₁,α₂是两个常数，α₁和α₂均为两个正奇数的商且满足1<α₁<2，1<α₂<2；β₁,β₂是两个常数，β₁>0，β₂>0；

对公式(9)和(10)求导得到滑模变量的时间导数为

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{edm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}})+x_{\mathrm{\&epsiv;}2}---\left(11\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}({\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+{\stackrel{\&CenterDot;\&CenterDot;}{q}}_{\mathrm{\&beta;mt}})+x_{\mathrm{\&beta;}2}---\left(12\right)$

将相对运动方程(1)，(2)，(3)和(4)代入式(11)和(12)整理得

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_1={\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1{x}_{\mathrm{\&epsiv;}2}^{{\mathrm{\&alpha;}}_1-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{x_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+\frac{u_{m_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}-\frac{u_{d_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&epsiv;}}}}{R_{\mathrm{mt}}})---\left(13\right)$

${\stackrel{\&CenterDot;}{\mathrm{\&sigma;}}}_2={\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2{x}_{\mathrm{\&beta;}2}^{{\mathrm{\&alpha;}}_2-1}(-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}}{R_{\mathrm{dm}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}+2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}+\frac{x_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-\frac{u_{m_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{d_{\mathrm{q\&beta;}}}}{R_{\mathrm{dm}}}+\frac{u_{t_{\mathrm{q\&beta;}}}}{R_{\mathrm{mt}}})---\left(14\right)$

其中 $u_{m_{\mathrm{q\&epsiv;}}}=u_{m_{2\mathrm{q\&epsiv;}}}-\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&epsiv;}}}$ 和 $u_{m_{\mathrm{q\&beta;}}}=u_{m_{2\mathrm{q\&beta;}}}+\frac{R_{\mathrm{dm}}}{R_{\mathrm{mt}}}{u}_{m_{1\mathrm{q\&beta;}}};$ 在制导过程中内有R_mt≈R_dm；所设计的制导律使得目标、防御导弹和拦截导弹在一条直线上，有如下关系和

步骤4：设计制导律并按制导律对导弹进行控制：

和是有限的值，设和M₁、M₂为常数；M₁、M₂大于零，小于等于拦截导弹最大过载能力；假定目标的机动可以忽略不计；

纵向平面和侧向平面的制导律分别设计如下

$u_{d_{\mathrm{q\&epsiv;}}}=-2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&epsiv;mt}}+\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&epsiv;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_1{\mathrm{\&alpha;}}_1}+{\mathrm{\&rho;}}_1\mathrm{sign}\left({\mathrm{\&sigma;}}_1\right)---\left(15\right)$

$u_{d_{\mathrm{q\&beta;}}}=2{\stackrel{\&CenterDot;}{R}}_{\mathrm{dm}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;dm}}-2\frac{{\stackrel{\&CenterDot;}{R}}_{\mathrm{mt}}{R}_{\mathrm{dm}}}{R_{\mathrm{mt}}}{\stackrel{\&CenterDot;}{q}}_{\mathrm{\&beta;mt}}-\frac{R_{\mathrm{dm}}{x}_{\mathrm{\&beta;}2}^{2-{\mathrm{\&alpha;}}_2}}{{\mathrm{\&beta;}}_2{\mathrm{\&alpha;}}_2}-{\mathrm{\&rho;}}_2\mathrm{sign}\left({\mathrm{\&sigma;}}_2\right)---\left(16\right)$

其中切换项增益ρ₁>M₁，ρ₂>M₂；sign(·)是符号函数；

按制导律(15)(16)对导弹进行控制。

2.根据权利要求1所述的一种基于视线的有限时间收敛主动防御制导控制方法，其特征在于，步骤3中的

3.根据权利要求1或2所述的一种基于视线的有限时间收敛主动防御制导控制方法，其特征在于，步骤3中的β₁＝β₂＝40。

4.根据权利要求3所述的一种基于视线的有限时间收敛主动防御制导控制方法，其特征在于，步骤4中0<M₁<100、0<M₂<100。