CN111007867B - 一种可预设调整时间的高超声速飞行器姿态控制设计方法 - Google Patents

Abstract

本发明记载一种可预设调整时间的高超声速飞行器姿态控制设计方法，属于航空航天领域，具体方案包括以下步骤：步骤一、明确控制***的设计任务；步骤二、建立高超声速飞行器姿态***数学模型；步骤三、设计自适应滑模控制律；步骤四、闭环***分析；步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环***的性能检验。本设计方法不仅使得高超声速飞行器姿态角在有限时间内达到了精度要求，而且可以根据性能指标要求预先设定需要的调整时间。此外，该控制律还能使不确定项估计值按需增长，也减小了控制器增益，很大程度上克服了一般鲁棒控制和自适应控制设计的保守性。

Description

一种可预设调整时间的高超声速飞行器姿态控制设计方法

技术领域

本发明属于航空航天领域，具体涉及一种可预设调整时间的、可使高超声速飞行器姿态角误差在指定时间内收敛到指定精度的自适应滑模控制律设计方法。

背景技术

在高超声速飞行器的控制中，需要保证飞行器的姿态角跟踪参考指令信号，否则难以实现飞行器的飞行任务。目前对于高超声速飞行器的控制，大多只是保证姿态角在无限时间内收敛到指定参考值，少数可以保证在有限时间内收敛，但无法预先设定调整时间。为解决控制中的不确定性问题，目前应用的大多数鲁棒控制律一般需要选取较大的控制增益；而大多自适应控制律也存在估计值持续增长的问题。

发明内容

本发明的目的是提供一种可预设调整时间的高超声速飞行器姿态控制设计方法，该方法能使***姿态角误差在预设调整时间内达到指定精度要求；而且该方法能解决多数鲁棒控制律控制增益较大和多数自适应律估计值持续增长的问题，使控制器的增益按需增长。

为了实现上述目的，本发明采取的技术方案如下：

一种可预设调整时间的高超声速飞行器姿态控制设计方法，包括以下步骤：

步骤一、明确控制***的设计任务，使高超声速飞行器姿态角在预定时间内跟踪参考指令信号，并具有良好的鲁棒性和对***扰动及外界干扰的自适应能力；

步骤二、建立高超声速飞行器姿态***数学模型；

步骤三、设计自适应滑模控制律；

步骤四、闭环***分析；

步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环***的性能检验。

进一步的，步骤一中，控制***的设计任务是：给定参考信号x_d，预设的调整时间T_f和允许的跟踪误差∈，设计合适的控制律使得闭环***的状态x有界，并且当t≥T_f时，***状态的跟踪误差|x-x_d|≤∈。

进一步的，步骤二中，建立高超声速飞行器姿态***数学模型的具体步骤如下：描述高超声速飞行器三通道姿态的二阶非线性微分方程如下所示：

式中，

分别为飞行器俯仰角、偏航角和滚转角，

分别为飞行器俯仰角速度、偏航角速度和滚转角速度；

分别为飞行器俯仰角加速度、偏航角加速度和滚转角加速度，J_x,J_y,J_z分别为飞行器对弹体坐标系各轴的转动惯量，M_x,M_y,M_z分别为作用在飞行器上的合外力对质心的力矩在弹体坐标系各轴上的分量，计算公式为：

M_x＝qsl_latm_x

M_y＝qsl_latm_y

M_z＝qsl_lonm_z

式中，q,s分别表示动压和参考面积，l_lat,l_lon分别为侧向和纵向参考长度；气动力矩系数m_x,m_y,m_z为空速V，迎角α，侧滑角β和姿态角速率ω_x,ω_y,ω_z以及舵偏角δ_x,δ_y,δ_z的函数，即：

m_x＝m_x(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

m_y＝m_y(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

m_z＝m_z(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

下面定义***状态变量和控制量为：

根据式(1)-(3)即可得到飞行器姿态控制***二阶非线性模型：

式中，d_lum表示***总的不确定性，满足

d_i＞0为未知有界常数，

为其保守估计值；

G(x)＝qsLB(x)T

其中，δ_e、δ_a和δ_r分别为左右升降舵偏转角和方向舵偏转角；L＝diag(l_lat,l_lat,l_lon)为飞行器参考长度构成的矩阵；

其中，

进一步的，步骤三中，设计自适应滑模控制律的具体步骤如下：

步骤1：定义滑模面

首先，定义***实际跟踪误差轨迹与期望跟踪误差轨迹之差

z＝e-η

式中，z＝z_i＝[z₁,z₂,z₃]^T，z为时间t的函数，e为***跟踪误差；η为跟踪误差的预期轨迹，其满足以下条件：

(1)η在区间[0,∞)上有二阶连续导数；

(2)η及其一阶导数

二阶导数

均有界；

(3)η(0)＝e(0)并且

(4)当t＞T_f时，η＝0，T_f为设定的调整时间；

η依据下列公式拟合：

式中

其中，

e_i(0)(i＝1,2,3)分别为***三通道初始时刻的跟踪误差，

为e_i(0)的一阶导数，η_i(t)(i＝1,2,3)分别为***三通道预期的跟踪误差轨迹，t_f和κ为待设计参数；

定义滑模面为

其中，s＝s_i＝[s₁,s₂,s₃]^T，C＝diag(c₁,…,c_i,…,c_n)，c_i＞0为待设计参数，

为z的一阶导数，显然，由于

所以有s(0)＝0，并且滑模面s的一阶导数

式中，

为z的二阶导数，

分别为误差e的一阶导数和二阶导数，

分别为状态x的一阶导数和二阶导数，

分别为状态参考信号x_d的一阶导数和二阶导数。

容易得到，该滑模面具有下列性质：

性质1：对任意的常数

如果|s_i|＜0对

成立，那么对于

有

成立，并且如果当t→∞时s_i→0，那么当t→∞时

步骤2：设计自适应滑模控制律

根据选取的滑模面，***自适应滑模控制律设计如下：

式中，

为状态变量，k＝[k₁,k₂,k₃]^T，其中k_i＞0,i＝1,2,3为待设计参数，ε＝[ε₁,ε₂,ε₃]^T为一常数向量，SAT_ε(s)为饱和函数，

其中

ε_i满足不等式

∈_i为***期望的跟踪误差，s_i为第i个子***的滑模面。

为未知有界参数d_i的估计值，由下列自适应律给出

μ_i＞0为待设计参数；

定义估计误差为

所以有

为估计误差的导数，

为估计值

的导数。

进一步的，步骤四中，闭环***分析的具体步骤如下：

飞行器姿态控制***二阶非线性模型为：

将控制律(7)及自适应律(8)应用于上述模型，在参考指令信号x_d具有二阶连续导数且x_d、

和

均有界的条件下，有下面的结论成立：

1)

均有界

2)当t→∞时，

其中，

都为常数；

3)对于预先设定的调整时间T_f＞0和稳态时允许的跟踪误差∈_i＞0，设计参数

那么，当t≥T_f时，各通道状态实际值x_i与状态期望值x_id满足：|x_i-x_id|≤∈_i(i＝1,2,3)。

证明：将控制律(7)代入式(6)得

考虑第i个子***

定义

当|s_i|＞ε，有

下面讨论s_i(t)的变化情况，s_i(0)＝0，从t＝0开始，s_i(t)值的变化有两种情况；第一种理想情况对于

都有|s_i(t)|＜ε_i，因此，对

所以，

另一种情况为存在时间T₁＞0使得|s_i(t)|≤ε_i,t∈[0,T₁)，|s_i(T₁)|＝ε_i，并且当t＞T₁时s_i(t)超出区间[-ε_i,ε_i]。由|s_i(t)|≤ε_i,t∈[0,T₁)得，

因此

再结合条件|s_i(T₁)|＝ε_i可知

因为

所以V_i(T₁)＜ε_i ²。

由式(9)得，当s_i(t)在区间[-ε_i,ε_i]之外时，V_i(t)严格单调递减，所以，s_i(t)不可能一直保持在区间[-ε_i,ε_i]之外；假设

成立，那么有

在

时成立；所以

成立，即当

时|s_i(t)|≤ε_i成立，与假设相矛盾，所以一定存在时间T₂＞T₁使得|s_i(T₂)|＝ε_i并且当t＞T₂时s_i(t)回到区间[-ε_i,ε_i]内；又V_i(t)在[T₁,T₂]上单调递减，所以

因此当t∈[T₁,T₂]时，一定有

成立，并且V_i(T₂)＜V_i(T₁)＜ε_i ²，因此

从t＝T₂开始，s_i(t)的值又有两种变化情况，一种是对

都有|s_i(t)|≤ε_i成立，因此

所以：

另外一种情况是存在时间T₃≥T₂，使得|s_i(t)|≤ε_i,t∈[T₂,T₃)，|s_i(T₃)|＝ε_i并且当t＞T₃时s_i(t)超出区间[-ε_i,ε_i]。所以，当t∈[T₂,T₃]时，

所以

t＞T₃时的分析过程同t＞T₁时一致，这样s_i(t)不断重复上述变化过程，直到不再超出区间[-ε_i,ε_i]；所以，

对任意t＞0成立，而且对所有的t＞0，

和

均有界。

再结合性质1可得

等价于

因为当t≥T_f时η_i(t)＝0，所以有

|x_i-x_id|＜∈_i,t＞T_f

并且，

再结合x_d以及η满足的条件可知，x_id、

η_i、

均有界，因此，

有界；

有界并且

所以，存在常数

当t→∞时，

当t≥T_f时，|x_i-x_id|≤∈_i(i＝1,2,3)成立。证毕。

此外，由自适应律(8)可知，

只有在|s_i|＞ε_i时才会更新，即只有当控制增益不足以抑制不确定项和干扰项时，自适应律才会通过增大

来提供额外的增益，这样保证了控制器增益不会过大，并且估计值不会持续增长。

综上，该控制律达到了设计目标，能够使***姿态角误差在预设调整时间内达到指定精度要求；控制器增益不会过大；并且不存在估计值持续增长的问题。

本发明的优点：提供了一种可预设调整时间的高超声速飞行器姿态控制设计方法，不仅使得高超声速飞行器姿态角在有限时间内达到了精度要求，而且可以根据性能指标要求预先设定需要的调整时间。此外，该控制律还能使不确定项估计值按需增长，从而也减小了控制器增益，一定程度上克服了一般鲁棒控制和自适应控制设计的保守性。

附图说明

图1设计方法流程图；

图2控制***结构框图；

图3俯仰角及其参考信号；

图4偏航角及其参考信号；

图5滚转角及其参考信号；

图6自适应律估计值。

具体实施方式

具体实施方式一

一种可预设调整时间的高超声速飞行器姿态控制设计方法，包括以下步骤：

步骤一：明确控制***的设计任务，给出了控制***设计的目的是使高超声速飞行器姿态角信号在预定时间内跟踪参考指令信号；步骤二：建立高超声速飞行器姿态***二阶非线性数学模型；步骤三：定义了滑模面，基于滑模面设计了自适应滑模控制律；步骤四：借助计算机数值仿真工具Matlab/Simulink进行了闭环***的性能检验。经过上述步骤，设计结束。

具体步骤如下：

步骤一、明确控制***的设计任务

控制***的设计任务是：给定参考信号x_d，预设的调整时间T_f＝2s和允许的跟踪误差∈＝0.5°，设计合适的控制律使得闭环***的状态x有界，并且当t≥T_f时，***状态的跟踪误差|x-x_d|≤∈。

步骤二、建立高超声速飞行器姿态***数学模型

飞行器姿态***二阶非线性数学模型为

其中

G(x)＝qsLB(x)T

和B(x)分别为

模型中涉及到的变量说明如下：

d_i＞0为未知有界常数，取其保守估计值

飞行器的侧向和纵向参考长度l_lat,l_lon分别为24.384m和18.288m；参考面积s为334.73m²；详细的空气动力学参数计算以及气动力和力矩的公式参考文献(侯明哲,谭峰.高超声速飞行器的增益协调鲁棒参数化控制[M].北京:科学出版社,2018.)。飞行器的初始姿态设置为γ(0)＝3°，ψ(0)＝2°，

ω_x＝10°/s，ω_y＝10°/s，ω_z＝10°/s；为了验证所设计控制律的鲁棒性，假设实际气动力矩系数m_x、m_y和m_z的值都比它们的标称值增加了30％，实际转动惯量J_x、J_y和J_z的值均较其标称值增加20％。

步骤三、设计自适应滑模控制律

设计过程分为两小步：

步骤1：定义滑模面

定义滑模面为

取C＝diag(5,5,5)，此时，由

可得，0＜ε_i＜0.031，取ε_i＝0.03；又根据

可得，μ_i≥11.1，取μ_i＝15。

步骤2：设计自适应滑模控制律

根据选取的滑模面，***自适应滑模控制律设计如下：

式中，

为未知有界参数d_i的估计值，由如下自适应律给出

其中，μ_i＝15,ε_i＝0.03；另外，取k＝[1,1,1]^T，T_f＝2s,κ＝0.35,t_f＝0.3。定义估计误差为

则有

步骤四：闭环***分析及检验

借助计算机数值仿真工具Matlab/Simulink进行闭环***的性能检验。

***的控制结构框图如图2所示，***以期望的飞行器姿态角与实际姿态角的偏差作为控制输入，控制器采用设计的自适应滑模控制律，依据相应的状态偏差解算出需要的控制力，使执行机构产生相应的舵偏，调整飞行器姿态角，使得实际姿态角收敛于期望姿态角。

最终得到的控制效果如图3-5所示，(其中参考指令信号x_d由制导信号生成)，图3中，飞行器的俯仰角在预设调整时间2s内跟踪上参考指令信号，动态性能良好，跟踪误差小于0.5°，达到了设计目标；图4中飞行器的偏航角在预设调整时间2s内跟踪上参考指令信号，动态性能良好，跟踪误差小于0.5°，达到了设计目标；图5中，飞行器的滚转角在预设调整时间2s内跟踪上参考指令信号，动态性能良好，跟踪误差小于0.5°，达到了设计目标。

图6为***不确定项估计值，由图可见***不确定项估计值有界且值较小，避免了估计值持续增长的现象，达到了设计目标。

综上所述，由仿真结果可看出在***控制增益k较小的情况下，飞行器姿态角可在预先设定的2s内实现对指令信号的跟踪，跟踪误差稳定在零附近，能够达到精度要求，并且不确定项估计值有界，没有持续增长现象，很好的实现了设计要求。

Claims (3)

1.一种可预设调整时间的高超声速飞行器姿态控制设计方法，其特征在于，包括以下步骤：

步骤一、明确控制***的设计任务，使高超声速飞行器姿态角在预定时间内跟踪状态参考指令信号，并具有良好的鲁棒性和对***扰动及外界干扰的自适应能力；

步骤二、建立高超声速飞行器姿态***数学模型；定义***状态变量x和控制量u为：

式中，θ,ψ,γ分别为飞行器俯仰角、偏航角和滚转角，δ_x,δ_y,δ_z分别为控制飞行器进行滚转、偏航和俯仰运动的舵偏角；

建立高超声速飞行器姿态***数学模型

式中，d_lum表示***总的不确定性，满足

d_i＞0为未知有界常数，

为其保守估计值；

其中，

G(x)＝qsLB(x)T

式中，q,s分别表示动压和参考面积，L＝diag(l_lat,l_lat,l_lon)为飞行器参考长度构成的矩阵，l_lat,l_lon分别为侧向和纵向参考长度；δ_e、δ_a和δ_r分别为左右升降舵偏转角和方向舵偏转角；

其中，

m_x＝m_x(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

m_y＝m_y(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

m_z＝m_z(V,α,β,ω_x,ω_y,ω_z,δ_x,δ_y,δ_z)

式中，气动力矩系数m_x,m_y,m_z为空速V，迎角α，侧滑角β和姿态角速率ω_x,ω_y,ω_z以及舵偏角δ_x,δ_y,δ_z的函数；

式中，J_x,J_y,J_z分别为飞行器对弹体坐标系各轴的转动惯量；

其中，

式中，

分别为飞行器俯仰角速度、偏航角速度和滚转角速度；

步骤三、定义滑模面，基于滑模面设计自适应滑模控制律；具体步骤如下：

步骤1：定义滑模面为

其中z＝e-η为***实际跟踪误差e与预期跟踪误差轨迹η的差值，

为z的一阶导数，C＝diag(c₁,…c_i,…c_n)，c_i＞0为待设计参数；

步骤2：设计自适应滑模控制律

根据选取的滑模面，***自适应滑模控制律设计如下：

式中，x＝x_i＝[x₁,x₂,x₃]^T＝[γ,ψ,θ]^T(i＝1,2,3)，为状态变量，k＝[k₁,k₂,k₃]^T，其中k_i＞0,i＝1,2,3为待设计参数，ε＝[ε₁,ε₂,ε₃]^T为一常数向量，SAT_ε(s)为饱和函数，

其中

ε_i满足不等式

∈_i为***期望的跟踪误差；s_i代表第i个子***的滑模面；

分别为预期跟踪误差轨迹η的一阶导数和二阶导数；

为状态变量x的一阶导数；

分别为状态参考信号x_d的一阶导数和二阶导数；

为未知有界参数d_i的估计值，由下列自适应律给出

μ_i＞0为待设计参数；

步骤四、闭环***分析；

步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环***的性能检验。

2.根据权利要求1所述的一种可预设调整时间的高超声速飞行器姿态控制设计方法，其特征在于，步骤一中，所述控制***的设计任务是：给定状态参考信号x_d，预设的调整时间T_f和允许的跟踪误差∈，设计合适的控制律使得闭环***的状态变量x有界，并且当t≥T_f时，***状态的跟踪误差|x-x_d|≤∈。

3.根据权利要求1所述的一种可预设调整时间的高超声速飞行器姿态控制设计方法，其特征在于，步骤四中，闭环***分析的具体步骤如下：

飞行器姿态控制***二阶非线性模型为：

将控制律(7)及自适应律(8)应用于上述模型，在状态参考信号x_d具有二阶连续导数且x_d、

和

均有界的条件下，有下面的结论成立：

1)x,

均有界；

2)当t→∞时，

其中，

为常数；

3)对于预先设定的调整时间T_f＞0和稳态时允许的跟踪误差∈_i＞0，设计参数

那么，当t≥T_f时，各通道状态变量x_i与状态期望值x_id满足：|x_i-x_id|≤∈_i(i＝1,2,3)；

此外，由自适应律(8)可知，

只有在|s_i|＞ε_i时才会更新，即只有当控制增益不足以抑制不确定项和干扰项时，自适应律才会通过增大

来提供额外的增益，这样保证了控制器增益不会过大，并且估计值不会持续增长。

Images