[发明专利]基于奇异摄动分析的弹性飞行器鲁棒有限时间控制方法

摘要

本发明涉及一种基于奇异摄动分析的弹性飞行器鲁棒有限时间控制方法，针对表征刚性体的慢变子系统，采用动态面控制设计舵偏角控制器，利用神经网络估计系统不确定信息，通过构造一阶滤波器与辅助信号将学习评价信息引入神经网络权重更新律并实现学习误差有限时间收敛，同时在控制律中加入误差分数阶项以保证跟踪误差有限时间收敛；针对表征系统弹性模态的快变子系统设计滑模控制算法进行弹性模态抑制。针对速度子系统设计PID控制器实现速度跟踪。

主权项

1.一种基于奇异摄动分析的弹性飞行器鲁棒有限时间控制方法，其特征在于步骤如下：/n步骤1：考虑弹性飞行器纵向通道动力学模型：/n /n /n /n /n /n /n所述的运动学模型由七个状态量和两个控制输入U＝[δ_e,Φ]^T组成；其中，V表示速度，h表示高度，γ表示航迹倾角，α表示攻角，q表示俯仰角速度，η和表示弹性模态，δ_e表示舵偏角，Φ表示节流阀开度；m、I_yy和g分别表示质量、俯仰轴的转动惯量和重力引起的加速度；ζ、ω和N分别表示弹性模态的阻尼比、自然振动频率和广义力；/n力与力矩以及各系数的表达式为：/nT＝A₁+B₁η,/nD＝A₂+B₂η,/nL＝A₃+B₃η,/nM_yy＝A₄+B₄η,/n /n其中，表示动压，表示平均气动弦长，z_T表示推力矩臂长，S表示气动参考面积，和均为气动参数，N^α、N⁰为表征弹性体动力学的相关系数；/n步骤2：定义高度跟踪误差为e_h＝h-h_d，设计航迹角指令γ_d为：/n /n其中，h_d表示高度参考指令，表示高度参考指令的一阶微分，k_h>0和k_i>0为设计参数；/n根据时标分离，将速度看作慢动态，设计航迹角指令的一阶微分为：/n /n其中，表示高度参考指令的二阶微分；/n步骤3：取x₁＝γ，x₂＝θ，x₃＝q，其中θ＝α+γ表示俯仰角；姿态子系统(3)-(6)写为以下形式：/n /n其中，/n定义ρσ＝η，ρB₂＝β₁；姿态子系统(9)写为以下形式：/n /n设置ρ＝0，姿态子系统(10)写为以下慢变子系统形式：/n /n /n /n /n其中，‘s’表示慢变子系统，δ_es表示慢变子系统的控制输入；/n将式(14)代入式(10)中，慢变子系统(11)-(14)写为以下形式：/n /n慢变子系统(15)可进一步写为以下严格反馈形式：/n /n其中，f_i，i＝1,3是由式(15)得到的未知平滑非线性函数，g_i,i＝1,3由式(15)得到的已知非线性函数；/n步骤4：定义ψ₁＝σ-σ_s，式(6)写为以下形式：/n /n其中，δ_ef＝δ_e-δ_es表示快变子系统的控制输入；/n将式(14)代入式(17)中，快变子系统(17)写为以下形式：/n /n将式(18)进一步写为如下矩阵形式：/n /n其中，ψ＝[ψ₁,ψ₂]^T，/n步骤5：第1步：定义航迹角跟踪误差为：/ne₁＝x_1s-γ_d (20)/n设计俯仰角虚拟控制量为：/n /n其中，表示神经网络最优权重向量的估计值，θ₁表示神经网络基函数向量，k₁>0、k_f1>0、0<η₁<1；/n设计一阶滤波器为：/n /n其中，表示通过式(22)所表达的滤波器后获得的信号，为滤波后得到的信号的微分信号，α₂>0；/n定义设计一阶滤波器如下：/n /n式中x_f1、θ_f1为x_1s、θ₁通过式(23)所示的滤波器后获得的信号，设计辅助矩阵P₁、辅助向量Q₁、W₁如下：/n /n式中，l₁>0，神经网络权重通过如下自适应律获得：/n /n其中，γ₁>0，Γ₁为正定常数对角矩阵；/n第2步：定义俯仰角跟踪误差为：/n /n设计俯仰角速率虚拟控制量为：/n /n其中，k₂>0、k_f2>0、0<η₂<1；/n设计一阶滤波器为：/n /n其中，表示通过式(28)所表达的滤波器后获得的信号，为滤波后得到的信号的微分信号，α₃>0为设计参数；/n第3步：定义俯仰角速率跟踪误差为：/n /n设计慢变子系统舵偏角为：/n /n其中，表示神经网络最优权重向量的估计值，θ₃表示神经网络基函数向量，k₃>0、k_f3>0、0<η₃<1；/n令设计一阶滤波器如下：/n /n式中x_f3、θ_f3为x_3s、θ₃通过式(31)所示的滤波器后获得的信号，设计辅助矩阵P₃、辅助向量Q₃、W₃如下：/n /n式中，l₃>0，神经网络权重通过如下自适应律获得：/n /n其中，γ₃>0，Γ₃为正定常数对角矩阵；/n步骤6：定义滑模切换函数为：/nc＝Gψ (34)/n其中，G∈R^2×2为设计的矩阵；/n设计快变子系统舵偏角为：/nδ_ef＝(GQ_f)⁺[-G(P_fψ)-K_fsign(c)] (35)/n其中，‘+’表示矩阵的摩尔彭罗斯逆，K_f为设计的正定矩阵；/n步骤7：定义速度跟踪误差为：/n /n其中，V_d为速度参考指令；/n设计节流阀开度Φ为：/n /n其中，k_pV>0、k_iV>0和k_dV>0为设计参数；/n步骤8：根据得到的慢变子系统的舵偏角δ_es和快变子系统的舵偏角δ_ef，得到姿态子系统的舵偏角δ_e＝δ_es+δ_ef，结合速度子系统的节流阀开度Φ，返回到飞行器动力学模型(1)-(6)，对高度和速度进行跟踪控制。/n