[发明专利]多源干扰系统的容错抗干扰控制方法

摘要

多源干扰系统的容错抗干扰控制方法，针对含有时变故障、可建模干扰和不可建模随机干扰的多源干扰系统，设计一种容错抗干扰控制器；首先，设计故障诊断观测器估计并抵消系统中的时变故障；其次，设计干扰观测器估计并抵消多源干扰系统中的可建模干扰；再次，设计鲁棒H∞状态反馈控制器抑制多源干扰系统中的不可建模随机干扰、故障估计误差和干扰估计误差；最后，基于故障诊断观测器、干扰观测器和鲁棒H∞状态反馈控制器，设计容错抗干扰控制器；本方法具有抗干扰性强、容错性能显著、工作可靠性高等优点，可用于航空、航天领域的姿态控制子系统中。

主权项

1.一种多源干扰系统的容错抗干扰控制方法，其特征在于包括以下步骤：首先，设计故障诊断观测器估计并抵消系统中的时变故障；其次，设计干扰观测器估计并抵消多源干扰系统中的可建模干扰；再次，设计鲁棒H_∞状态反馈控制器抑制多源干扰系统中的不可建模随机干扰、故障估计误差和干扰估计误差；最后，基于故障诊断观测器、干扰观测器和鲁棒H_∞状态反馈控制器，设计容错抗干扰控制器；具体步骤如下：第一步，搭建含有多源干扰系统的动力学模型，并写成状态空间表达式针对含有时变故障、可建模干扰和不可建模随机干扰的多源干扰系统，搭建其系统动力学模型，并写成状态空间表达式如下：x·(t)+Ef(x·,t)=Ax(t)+B1(u(t)+d1(t)+F(t))+B2d2(t)]]>其中，x(t)为多源干扰系统状态变量，u(t)为控制输入，d₁(t)为可建模干扰，F(t)为时变故障，d₂(t)为不可建模随机干扰，A,E,B₁和B₂为已知维数的矩阵，为系统非线性项且满足Lipschitz条件，外部模型描述干扰d₁(t)由下列外部干扰模型∑₁表示：Σ1:d1(t)=Vw(t)w·(t)=Ww(t)+B3δ(t)]]>其中，w(t)为可建模干扰模型的状态变量，V为可建模干扰模型的输出矩阵，W表示可建模干扰模型的系统阵，δ(t)为能量有界的不可建模随机干扰，B₃为不可建模随机干扰的增益阵；第二步，设计故障诊断观测器针对多源干扰系统中的时变故障F(t)，设计故障诊断观测器对其进行实时估计，并求得估计值进而得到故障估计误差第三步，设计干扰观测器针对多源干扰系统中的可建模干扰d₁(t)，设计干扰观测器对其进行实时估计，并求得估计值进而得到干扰估计误差为w(t)的估计值；第四步，设计鲁棒H_∞状态反馈控制器针对多源干扰系统中的不可建模随机干扰d₂(t)、故障估计误差e_F(t)和干扰估计误差e_w(t)，设计鲁棒H_∞状态反馈控制器对其进行抑制，控制器结构如下：u_f(t)＝Mx(t)其中，u_f(t)为鲁棒H_∞状态反馈控制输入，M为待定状态反馈控制器增益阵；第五步，设计容错抗干扰控制器设计容错抗干扰控制器，对系统中的时变故障F(t)和可建模干扰d₁(t)进行抵消，不可建模随机干扰d₂(t)、故障估计误差e_F(t)和干扰估计误差e_w(t)进行抑制，容错抗干扰控制器结构如下：u(t)=uf(t)-F^(t)-d^1(t)]]>则多源干扰系统可表示为：x·(t)+Ef(x·,t)=(A+B1M)x(t)+B1Vew(t)+B1eF(t)+B2d2(t)]]>整理可建模干扰模型的系统估计误差方程和时变故障的系统估计误差方程如下：e·w(t)=(W+LB1V)ew(t)+LB1eF(t)+LB2d2(t)+B3δ(t)]]>e·F(t)=KB1Vew(t)+KB1eF(t)+KB2d2(t)+F·(t)]]>联列上述多源干扰系统、时变故障的系统估计误差方程和可建模干扰的系统估计误差方程得到闭环系统：x·(t)e·w(t)e·F(t)=A+B1MB1VB10W+LB1VLB10KB1VKB1x(t)ew(t)eF(t)+B2LB2KB2d2(t)+0B30δ(t)+001F·(t)-E00f(x·(t))z∞(t)=C0C1C2x(t)ew(t)eF(t)]]>其中，z_∞(t)为H_∞性能参考输出，[C₀C₁C₂]为H_∞性能可调输出矩阵；第六步，增益矩阵求解利用凸优化算法求解多源干扰系统的容错抗干扰控制器增益阵；给定初始值x(0)、e_w(0)和e_F(0)，可调输出矩阵[C₀C₁C₂]，非线性权重参数λ，干扰抑制度γ₁、γ₂和γ₃，求解以下凸优化问题：minxT(0)eT(0)P1P2xT(0)eT(0)T]]>Φ=Φ11B1G-E0B20Φ18P1C0T*Φ220P2H1R2B2P2H2λ(UB1G)TCT**-λ2I000λ(UE)T0***-γ120000****-γ220λ(UB2)T0*****-γ32I00******-I0*******-I]]>其中:e(0)＝[e_w(0)e_F(0)]^T,Φ₁₁＝(AP₁+B₁R₁)+(AP₁+B₁R₁)^T,Φ₂₂＝(P₂W₁+R₂B₁G)+(P₂W₁+R₂B₁G)^T,Φ₁₈＝λ(AP₁+B₁R₁)^T,C＝[C₁ C₂],G＝[E I],H₁＝[B₃ 0]^T,H₂＝[0 1]^T；符号*表示对称矩阵中相应部分的对称块，求解得P₁、P₂、R₁和R₂，则干扰观测器和故障诊断观测器增益阵为LK=P2-1R2,]]>状态反馈控制器增益阵为