一种含有减速伞的导弹敏捷转弯最优控制方法

1.本发明属于导弹大角度机动转弯控制领域。涉及一种具有减速伞的导弹敏捷转弯最优控制方法。


背景技术：

2.导弹大角度机动转弯控制是导弹敏捷转弯领域内的一类重要问题，被广泛应用于控制、攻防对抗等领域。
3.随着军事科技的快速发展，导弹大角度机动转弯又叫做导弹敏捷转弯，被广泛应用于战略战术武器中，尤其应用于空空导弹、地空导弹等武器系统中。众所周知，传统导弹在进行大角度机动转弯过程中，导弹离开发射架时迅速进入初制导阶段，而此时导弹的速度指向与导弹视线之间往往具有很大的夹角。因此，在导弹的初制导的作用下，夹角会变得越来越小，直到导弹的指向与目标之间的夹角为零。所以，传统的大角度机动转弯往往仅考虑如何能完成大角度机动转弯，不会考虑如何让在保证完成大角度机动转弯的约束。因此，在此情况下本专利提出一种带有减速伞的导弹敏捷转弯最优控制方法。
4.该方法首先将导弹敏捷转弯过程分为减速段、大角度机动段以及末制导段。在进入减速段的过程中，减速伞迅速张开，在减速伞的作用下导弹速度迅速降低完成导弹极短时间减速，同时在空气舵的作用下导弹完成大角度机动转弯，大角度机动完成之后，导弹迅速进入末制导段，从而完成对目标的打击。
5.减速转弯段开始转弯时刻为减速段的起始时刻t0，减速伞迅速张开，在减速伞的作用下导弹速度迅速降低；同时在空气舵的作用下导弹完成大角度机动转弯，转弯角度范围为90
°
～180
°
。
6.因此，综合以上问题，目前缺少一种含有减速伞的导弹敏捷转弯最优控制方法。


技术实现要素：

7.有鉴于此，本发明提供了一种含有减速伞的导弹敏捷转弯最优控制方法，在进行大角度机动转弯的过程中，完成对每个阶段每个执行结构的最优控制率设计。
8.为达到上述目的，本发明的技术方案包括如下步骤：
9.步骤1：将减速伞加入导弹动力学模型中，对传统铅垂面上导弹运动学方程改进，得到铅垂面上伞-导弹动力学方程为：
[0010][0011]
其中，p发动机推力，m为导弹质量，gr导弹重力，v为导弹飞行速度、v 为导弹飞行加速度，θ为导弹弹道倾角，α为导弹攻角，为导弹俯仰角，为的一阶导，δδ为导弹舵偏，为导弹俯仰角速率，jz为导弹转动惯量，y为导弹升力，x,y为铅垂面位置，为x的一阶导，为y的一阶导，f
aero
为导弹气动阻力，f
um
为减速伞作用力，为导弹攻角产生的俯仰力矩，为导弹俯仰阻尼力矩，为导弹舵偏产生的俯仰力矩，为减速伞产生的俯仰力矩，δum 为减速伞偏角。
[0012]
其次，方程中的各个气动力与气动力矩的表达式为：
[0013][0014]
其中，s为导弹参考面积，q为导弹动压，c
x
为导弹阻力因数，c
um
为减速伞阻力因数，s
um
为减速伞面积，l
um
为减速伞重心与导弹之间的参考长度，cy为导弹升力因数，l为导弹参考长度，为导弹攻角产生的俯仰力矩因数，为导弹俯仰阻尼因数，为导弹舵偏产生的俯仰力矩因数，ρ大气密度，减速伞产生的俯仰阻尼因数。
[0015]
步骤2：定义指代向量即指代向量中指代量x1～x6分别用于指代所述铅垂面上伞-弹动力学方程转化为：
[0016][0017]
将状态方程转换到状态空间下，得到：
[0018]
其中
[0019][0020][0021]
其中f1～f6为矩阵f(x)中的元素，g
11
，g
21
，g
31
，g
41
，g
51
，g
61
，g
21
，g
22
，g
32
， g
42
，g
52
，g
62
分别为矩阵g(x)中的元素，u为控制向量：其中，u1，u2分别为减速伞控制量，舵偏控制量。
[0022]
步骤3：减速转弯段控制率设计，具体过程如下：
[0023]
步骤3.1：在进行减速伞控制率设计的过程中，采用最优控制动态规划方法，设计过程中考虑的约束条件为，导弹以最小的转弯半径r完成大角度机动转弯；建立最优性能指标函数j
*
如下：
[0024][0025]
其中t0为起始时刻，t1为末端时刻。
[0026]
步骤3.2：根据j
*
建立哈密顿函数h如下：
[0027][0028]
其中λ
t
的表达式为：
[0029]
步骤3.3：根据步骤3.1和步骤3.2得到哈密顿函数方程为：
[0030][0031]
步骤3.4：此阶段过程中，考虑在减速力与控制力不受约束时，则得到此阶段最优控制的取值条件为：
[0032][0033]
其中，fi为f(x)中的元素，i为1～6，和表达式为：
[0034][0035][0036][0037]
步骤3.5：根据步骤3.4中的方程可以得到最优控制下的表达式为：
[0038][0039]
为最优的减速伞控制率，为最优的舵偏控制量。
[0040]
步骤3.6：将步骤3.5中的方程带入步骤3.2最优性能指标函数得到如下方程：
[0041][0042]
为；
[0043]
其中得到：
[0044][0045]
由于j
*
中不包含x2,x4,x5,x6，则：
[0046]
步骤3.7：将步骤3.6中的方程进行求解得到：
[0047][0048]
进一步：
[0049][0050]
其中ξγ为两个设定比例因子。
[0051]
步骤3.8：将步骤3.7与步骤3.5结合可以得到减速转弯段段最优控制输入为：
[0052][0053]
进一步地，含有减速伞的导弹，其飞行过程中导弹的减速伞张开时间忽略不计。
[0054]
进一步地，含有减速伞的导弹，其减速伞的伞绳一直处于拉紧状态，即伞绳与导弹的连接为刚体连接。
[0055]
进一步地，含有减速伞的导弹，其伞绳脱离导弹时，伞绳与导弹的脱离时间忽略不计。
[0056]
进一步地，含有减速伞的导弹，其减速伞作用力，默认作用于导弹的x轴上。
[0057]
有益效果：
[0058]
本发明提供的一种含有减速伞的导弹敏捷转弯最优控制方法，通过将减速伞在减速方面的优点和导弹高速飞行相结合，提出的新的伞-弹控制系统；
[0059]
本发明提供的一种含有减速伞的导弹敏捷转弯最优控制方法，是一种的新形势的
导弹大角度机动转弯方法，该方法可以普遍适用于对转弯半径要求比较严格的飞行环境。
[0060]
本发明提供的一种含有减速伞的导弹敏捷转弯最优控制方法，伞-弹飞行过程中，该发明通过动态规划的设计方法，设计不同阶段下的最优控制律，根据对应的仿真结果可以看出该方法的有效性。
附图说明
[0061]
图1为本发明的一种新的大角度机动转弯过程，该过程主要包含开伞减速段以及伞控大角度机动转弯段；
[0062]
图2为本发明导弹在纵向平面的受力分析示意图。
具体实施方式
[0063]
下面结合附图并举实施例，对本发明进行详细描述。
[0064]
本发明提供了一种本发明主要利用最优控制中的动态规划方法解决伞-弹系统大角度机动敏捷转弯所遇到的最优控制问题。首先，针对整个过程中各个通道的最优控制问题，采用最优控制原理对各伞控力、气动控制力等执行机构的控制率进行求解，进而得到整个飞行过程的最优控制率。
[0065]
图1为本发明的一种新的大角度机动转弯过程，该过程主要包含开伞减速段以及伞控大角度机动转弯段。
[0066]
本发明一种具有减速伞的导弹敏捷转弯最优控制方法，具体步骤如下：
[0067]
步骤1：首先由于将减速伞加入导弹动力学模型中，考虑减速伞的作用，首先对传统的导弹动力学模型进行改进，建立新型的伞-弹敏捷转弯控制系统，在建立伞-弹敏捷转弯控制系统之前，本发明首先做如下假设：
[0068]
1)飞行过程中导弹的减速伞可以迅速张开，同时忽略张开时间；
[0069]
2)减速伞的伞绳一直处于拉紧状态，即可将伞绳与导弹的来连接视作为刚体连接；
[0070]
3)伞绳迅速脱离导弹时，伞绳与导弹的脱离时间忽略不计；
[0071]
4)减速伞作用力，默认作用于导弹的x轴上。
[0072]
图2为本发明导弹在纵向平面的受力分析示意图。
[0073]
正在以上假设的基础上，对传统铅垂面上导弹运动学方程改进，得到铅垂面上伞-导弹动力学方程为：
[0074][0075]
其中，p发动机推力，m为导弹质量，gr导弹重力，v为导弹飞行速度、为导弹飞行加速度，θ为导弹弹道倾角，α为导弹攻角，为导弹俯仰角，为的一阶导，δδ为导弹舵偏，为导弹俯仰角速率，jz为导弹转动惯量，y为导弹升力，x,y为铅垂面位置，为x的一阶导，为y的一阶导，f
aero
为导弹气动阻力，f
um
为减速伞作用力，为导弹攻角产生的俯仰力矩，为导弹俯仰阻尼力矩，为导弹舵偏产生的俯仰力矩，为减速伞产生的俯仰力矩，δum为减速伞偏角；
[0076]
其次，方程中的各个气动力与气动力矩的表达式为：
[0077][0078]
其中，s为导弹参考面积，q为导弹动压，c
x
为导弹阻力因数，c
um
为减速伞阻力因数，s
um
为减速伞面积，l
um
为减速伞重心与导弹之间的参考长度，cy为导弹升力因数，l为导弹参考长度，为导弹攻角产生的俯仰力矩因数，为导弹俯仰阻尼因数，为导弹舵偏产生的俯仰力矩因数，ρ大气密度，减速伞产生的俯仰阻尼因数；
[0079]
步骤2：定义指代向量即指代向量中指代量x1～x6分别用于指代所述铅垂面上伞-弹动力学方程转化为：
[0080][0081]
将状态方程转换到状态空间下，得到：
[0082]
其中
[0083][0084][0085]
其中f1～f6为矩阵f(x)中的元素，g
11
，g
21
，g
31
，g
41
，g
51
，g
61
，g
21
，g
22
，g
32
， g
42
，g
52
，g
62
分别为矩阵g(x)中的元素，u为控制向量：其中，u1，u2分别为减速伞控制量，舵偏控制量。
[0086]
定义指代向量即指代向量中指代量x1～x6分别用于指代所述铅垂面上伞-弹动力学方程转化为：
[0087][0088]
将状态方程转换到状态空间下，得到：
[0089]
其中
[0090][0091][0092]
其中f1～f6为矩阵f(x)中的元素，g
11
，g
21
，g
31
，g
41
，g
51
，g
61
，g
21
，g
22
，g
32
， g
42
，g
52
，g
62
分别为矩阵g(x)中的元素，u为控制向量：其中，u1，u2分别为减速伞控制量，舵偏控制量；
[0093]
步骤3：减速转弯段控制率设计，具体过程如下：
[0094]
步骤3.1：在进行减速伞控制率设计的过程中，采用最优控制动态规划方法，设计过程中考虑的约束条件为，导弹以最小的转弯半径r完成大角度机动转弯；建立最优性能指标函数j
*
如下：
[0095][0096]
其中t0为起始时刻，t1为末端时刻；
[0097]
步骤3.2：根据j
*
建立哈密顿函数h如下：
[0098][0099]
其中λ
t
的表达式为：
[0100]
步骤3.3：根据步骤3.1和步骤3.2得到哈密顿函数方程为：
[0101][0102]
步骤3.4：此阶段过程中，考虑在减速力与控制力不受约束时，则得到此阶段最优控制的取值条件为：
[0103][0104]
其中，fi为f(x)中的元素，i为1～6，和表达式为：
[0105][0106][0107][0108]
步骤3.5：根据步骤3.4中的方程可以得到最优控制下的表达式为：
[0109][0110]
为最优的减速伞控制率，为最优的舵偏控制量；
[0111]
步骤3.6：将步骤3.5中的方程带入步骤3.2最优性能指标函数得到如下方程：
[0112][0113]
为；
[0114]
其中得到：
[0115][0116]
由于j
*
中不包含x2,x4,x5,x6，则：
[0117]
步骤3.7：将步骤3.6中的方程进行求解得到：
[0118][0119]
进一步：
[0120][0121]
其中ξγ为两个设定比例因子；
[0122]
步骤3.8：将步骤3.7与步骤3.5结合可以得到减速转弯段段最优控制输入为：
[0123][0124]
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。