一类具有短时延和数据丢包的网络化线性参数变化系统的H∞容错控制方法

摘要

本发明提供一类具有短时延和数据丢包的网络化线性参数变化系统的H∞容错控制方法，属于网络化线性参数变化系统的控制领域。本方法考虑线性参数变化系统存在执行器故障、网络短时延和数据丢包的情况，首先构造适当的Lyapunov函数得到闭环容错控制系统稳定的充分条件，然后利用李雅普诺夫稳定性定理和线性矩阵不等式分析方法得到了H∞容错控制器存在的充分条件。最后利用近似基函数和网格化技术将无限维的线性矩阵不等式组的求解问题近似为有限维线性矩阵不等式组的求解问题，并且利用Matlab LMI工具箱进行求解，给出了H∞容错控制器的增益矩阵K(ρ(k))＝YP011(ρ(k))。本发明适用于一般的网络化线性变化系统的H∞容错控制，降低了H∞容错控制器的保守性，获得了较好的性能指标γ。

说明书

技术领域

本发明属于网络化线性参数变化系统的控制领域，涉及一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法。

背景技术

线性参数变化(Linear Parameter-Varying，LPV)系统是一类参数在不断变化的系统，此类系统的状态矩阵的元素是具有时变参数的确定函数，而与函数相关的时变参数的范围是能够测量的。许多实际的系统都能用这样的模型来描述，如飞行器系统、风能转换系统等，这类复杂系统一旦发生故障，就可能造成重大的财产损失甚至是人员伤亡，因此，研究如何减少灾难性事故、环境污染和人身财产损失的发生，提高现代复杂控制系统的可靠性与安全性，具有十分重要的社会意义和经济意义。特别是由于网络的应用，控制器与执行器之间可能存在一定的距离，一旦发生故障，难以及时的进行修理。动态系统的故障诊断技术和容错控制技术为解决上述问题提供了重要的技术手段，研究网络化线性参数变化系统的容错控制具有重要意义。

针对一般的网络化系统的容错控制和线性参数变化系统的容错控制分别已有较充分的研究，本发明在前人研究的基础上主要实现了在控制器和执行器之间采用网络通信方式的网络化线性参数变化系统的容错控制，考虑了网络诱导的时延、随机丢包和随机故障等。

发明内容

针对上述现有技术中存在的问题，本发明提供了一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法。设计了相应的H_∞容错控制器，使得网络化线性参数变化系统的闭环容错控制在执行器发生随机故障时，在一定范围内仍能保持闭环系统渐进稳定，并且满足H_∞性能指标。

本发明的技术方案：

一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法，包括以下步骤：

1)建立线性参数变化系统模型

其中，x(k)∈Rⁿ为状态向量，u(k)∈R^m为控制输入向量，z(k)∈R^r为系统输出向量；ω(k)∈R^q是外部干扰输入向量，属于集合l₂[0,∞)，l₂[0,∞)是平方可积向量空间；A(ρ(k))∈R^n×n、B₀(ρ(k))∈R^n×m、B₁(ρ(k))∈R^n×m、C(ρ(k))∈R^r×n、D(ρ(k))∈R^r×q和R(ρ(k))∈R^n×q为系统矩阵，都是ρ(k)的函数，ρ(k)满足实时可测；f(k,x_k)∈Rⁿ是非线性向量函数，满足Lipschitz条件其中G₁∈Rⁿ是已知常数矩阵，ε₁是大于0的标量；

2)设计状态反馈容错控制器

u(k)＝MK(ρ(k))x(k) (2)

其中，K(ρ(k))∈R^m×n是状态反馈容错控制器的增益矩阵，M表示执行器故障矩阵：

M＝diag{m₁,m₂,…,m_n}

其中，m₁,m₂,…,m_n∈[0,1]为n个互相独立的随机变量，E()表示数学期望计算，和分别为M和m_i的期望，当m_i＝1时表示执行器正常，当m_i＝0时表示执行器完全失效，当0<m_i<1时，则表示执行器存在部分失效；

3)建立网络短时延和数据丢包的数学模型

用一步时延u(k-1)来表示网络中存在的短时延；

用两个状态的马尔可夫链来描述网络通道中数据传输的丢包过程，其状态转移矩阵为P＝[p_ij]，p_ij＝prob{σ(k+1)＝j|σ(k)＝i}，prob{·}表示事件发生概率，当σ(k)＝0时表示数据在网络通道中没有丢失，当σ(k)＝1时表示数据在网络通道中有丢失；

当σ(k)＝0，且网络诱导时延小于一个采样周期时，设增广向量在执行器发生随机故障时，具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统的模型为：

其中，F(k)＝[f^T(k,x(k))>T，0为零矩阵；

当σ(k)＝1时，在网络通道中传输的数据包丢失，此时控制输入使用前一时刻的值，即u(k)＝u(k-1)，则具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统的模型为：

其中，

I为单位阵；

结合σ(k)＝0和σ(k)＝1两种情况，针对执行器发生的随机故障，把具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统描述为马尔可夫跳变系统：

4)构造Lyapunov函数

其中P_i(ρ(k))＝diag{P_i11(ρ(k))，P_i22(ρ(k))},i＝0,1，P₀₁₁(ρ(k))∈Rⁿ，P₀₂₂(ρ(k))∈Rⁿ，P₁₁₁(ρ(k))∈Rⁿ，P₁₂₂(ρ(k))∈Rⁿ是对称正定矩阵；

5)具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统渐进稳定和存在H_∞容错控制的充分条件

当存在对称正定矩阵P₀₁₁(ρ(k))，P₀₂₂(ρ(k))，P₁₁₁(ρ(k))，P₁₂₂(ρ(k))时，实数矩阵Y(ρ(k))和标量ε₁>0，使矩阵不等式(4)和(5)成立，则在执行器发生随机故障的情况下，具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统渐进稳定，存在H_∞容错控制，状态反馈容错控制器的增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))；

其中，*代表对称位置矩阵的转置，

P₀₁₁(ρ(k))，P₀₂₂(ρ(k))，P₁₁₁(ρ(k))，P₁₂₂(ρ(k))，K(ρ(k))，Y(ρ(k))和ε₁为未知变量，其他变量都是已知的；利用Matlab>∞性能约束，H_∞容错控制器增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))，能继续进行步骤6)；当矩阵不等式(4)和(5)无解时，则闭环容错控制系统不是渐进稳定的，且不能获得H_∞容错控制器增益矩阵，不能进行步骤6)，结束；

6)最优化H_∞控制

选取基函数：f₁(ρ)＝1，f₂(ρ)＝sin(k)²，得到：

P₀₁₁(ρ(k))＝P₀₁₁₁+sin(k)²P₀₁₁₂，P₁₁₁(ρ(k))＝P₁₁₁₁+sin(k)²P₁₁₁₂，

P₀₂₂(ρ(k))＝P₀₂₂₁+sin(k)²P₀₂₂₂，P₁₂₂(ρ(k))＝P₁₂₂₁+sin(k)²P₁₂₂₂，

用网格化技术将ρ的参数空间划分为十个，根据求出闭环容错控制系统的性能指标γ，H_∞控制下最优H_∞性能指标γ_opt的条件为：

令e＝γ²，当优化问题：

min e s.t.

成立时，则闭环容错控制系统在符合H_∞容错控制条件下，系统的最优H_∞性能指标同时H_∞容错控制器的增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))；

7)实现H_∞容错控制

根据求出的H_∞容错控制器的增益矩阵K(ρ(k))，具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统实现H_∞容错控制。

本发明的有益效果：

1)本发明针对线性参数变化系统，同时考虑了外部扰动以及执行器可能出现的随机故障、网络中存在的短时延、数据丢包，通过一系列的推导、转化建立了闭环容错控制系统模型，给出了H_∞容错控制器的设计方法；

2)本发明考虑了网络信号传输中存在的短时延和数据丢包现象，更具有实际意义；

3)本发明适用于一般网络化线性参数变化系统的H_∞容错控制，提出了依赖于参数的H_∞性能约束，降低了该H_∞容错控制器设计方法的保守性。

附图说明

图1是一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法的流程图。

图2是具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制的结构图。图中x(k)∈Rⁿ为系统的状态向量，为系统的输出向量，为系统的控制输入向量，为外部干扰输入向量。

图3是时的H_∞控制状态响应图。

图4是时的H_∞控制状态响应图。

图5是时的H_∞控制状态响应图。

具体实施方式

下面结合附图对本发明的具体实施方式做进一步说明。

实施例1

如附图1所示，一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法，包括以下步骤：

步骤1：建立线性参数变化系统模型为式(1)。

步骤2：设计状态反馈控制器为式(2)。

步骤3：建立具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统为式(3)。

步骤4：构造合适的Lyapunov函数：

其中P_i(ρ(k))＝diag{P_i11(ρ(k)),P_i22(ρ(k))}，i＝0,1。

当ω_k＝0，由式(3)有：

其中，

步骤5：基于步骤2构造的Lyapunov函数，利用Lyapunov稳定性理论和线性矩阵不等式分析方法，得到具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统渐近稳定和H_∞控制器存在的充分条件。

首先判断闭环容错控制系统的稳定性，得到闭环容错控制系统渐进稳定的充分条件。

假设根据Schur补引理可得

当i＝0时，即

在式(8)的左边和右边同时乘以对角矩阵以及它的转置矩阵，可以得到式(9)。

同理可得式(10)成立。

其中，

所以，可以由式(9)和(10)得到即其中，β为的最小特征值。由此可得，因此，可以得到闭环容错控制系统(3)是渐进稳定的，并且H_∞容错控制器的增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))。

接着，根据H_∞性能约束条件，在零初始条件下定义：

其中，

可以写成如下形式：

假设根据schur补引理可得：

式(12)成立等价于式(13)成立：

同理可得式(14)成立：

可以由式(13)和(14)得到当N→∞时，k从0加到∞得到：

这也意味着闭环容错控制系统(3)对所有非零ω(k)满足H_∞性能指标γ，并且存在H_∞容错控制器，H_∞容错控制器的增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))。

步骤6：最优化H_∞控制

利用Matlab LMI工具箱求解优化问题，当有解时闭环容错控制系统在符合H_∞容错控制条件下，系统的最优H_∞性能指标同时H_∞容错控制器增益矩阵为K(ρ(k))＝YP₀₁₁(ρ(k))。

步骤7：实现H_∞容错控制

根据求出的H_∞容错控制器增益矩阵K(ρ(k))，具有短时延和数据丢包的网络化线性参数变化系统的闭环容错控制系统实现H_∞容错控制。

实施例2：

采用本发明提出的一类具有短时延和数据丢包的网络化线性参数变化系统的H_∞容错控制方法，具体实现方法如下：

某风能发电机转化为线性参数变化系统，其数学模型为式(1)，给定其系统参数为：

其中，ρ(k)＝sin(k)²为时变参数，选取扰动信号为ω_k＝1/k²，随机丢包的马尔可夫链状态转移概率矩阵为：

选取3种随机故障情形：

1)随机故障的期望为方差为

2)随机故障的期望为方差为

3)随机故障的期望为方差为

选取基函数：f₁(ρ(k))＝1，f₂(ρ(k))＝sin(k)²，可以得到：

P₀₁₁(ρ(k))＝P₀₁₁₁+sin(k)²P₀₁₁₂，P₁₁₁(ρ(k))＝P₁₁₁₁+sin(k)²P₁₁₁₂，

P₀₂₂(ρ(k))＝P₀₂₂₁+sin(k)²P₀₂₂₂，P₁₂₂(ρ(k))＝P₁₂₂₁+sin(k)²P₁₂₂₂。

用网格化技术将ρ的参数空间划分为十个，并且根据三种不同的随机故障，应用Matlab LMI工具箱求解控制器，使得H_∞性能指标γ最小化。在不同的丢包概率下相应的控制器参数以及H_∞性能指标γ在表1中给出。

给定初始状态x₀＝[0.2>T，用Matlab仿真出不同随机故障情况下的，闭环容错控制系统的状态响应，如附图3-附图5所示。

表1不同丢包情况下的控制器参数

由附图3-附图5可以看出，闭环容错控制系统的状态响应曲线经过一段时间的振荡后都收敛为零，说明按本发明方法设计的H_∞容错控制器可以很好的使闭环容错控制系统(3)渐进稳定。而且随着随机故障发生概率的增大，H_∞性能指标γ也越大，闭环系统达到稳定状态的时间也增长，说明随机故障概率对系统的性能是有重要影响的，与实际情况相符。