一种城市排水管道的有限时间控制方法与流程

本发明属于自动化技术领域，特别是现代城市排水行业中的自动化控制，涉及一种城市排水管道的有限时间控制方法，利用输出反馈控制器对城市排水系统中的水流状态进行有限时间控制，实现对城市排水管道有效控制。

背景技术：

城市排水系统是城市地表水控制的重要基础设施，也是城市水污染防治和城市排渍防涝的重要工程，所以及时有效地收集、输送、处理和排放整个城市的污水和雨水是城市排水管道系统的重要任务。由于下水管道的蓄水能力、水泵的排水能力、污水处理厂的污水处理能力这些因素的限制，在暴雨来临时，城市排水管道很容易发生溢流，甚至发生严重的城市内涝。因此，实现对城市排水系统的实时控制,可以有效防止城市排水管道溢流和城市内涝的发生。

现有的城市排水管道系统控制大多采用状态反馈控制方法，然而现实中很多复杂情况下的水流状态变量无法测量或难以得到，因此也无法设计状态反馈控制器，这时可以采用输出反馈控制的方法。此外在夏季，雨水量很大导致现有的控制方法不能在给定的有限时间之内实现及时、有效地控制排水管道水位。因此，急需一种新的方法，对城市排水管道系统在给定的有限时间内进行及时、有效地控制。

技术实现要素：

本发明的目的是针对现有控制方法的不足，提出一种城市排水管道的控制方法，在给定的有限时间之内实现对城市排水系统及时、有效地控制。

本发明使用了基于冗余通道的数据传输方法，在一条传输通道出现水流数据丢失的情况下，采用多条冗余通道对水流状态数据进行传输，提高数据传送的成功率，并基于随机有限时间有界性分析结果设计动态输出反馈控制器，使用矩阵不等式方法求解控制器增益，从而在给定的有限时间内实现对城市排水系统的准确控制。

本发明方法的包括：

(1).建立城市排水管道系统状态空间模型；

首先，基于流体力学原理建立排水管道系统模型：

符号σ表示求和运算，π表示连乘。

x(k)表示k时刻检测到的水流状态向量，x(k)＝[x1(k),x2(k),x3(k)]^t，x1(k)、x2(k)、x3(k)分别表示k时刻所测量的水位高度值、水流速度值和水压值，上标t表示矩阵的转置；

u(k)∈r^1×1表示k时刻的控制输入量，为排水管道上游进入管道的水流量与排水管道下游流出管道的水流量之差；

y(k)∈r^1×1表示k时刻排水控制系统的测量输出，为传感器测量到的水流状态通过网络传输到远程控制中心的测量值；

z(k)∈r^1×1表示k时刻排水控制系统待估计的输出信号；城市排水控制系统利用三个传输通道将各传感器节点的水务数据信息传输到集中控制中心，在检测到主传输通道发生数据丢失的情况下，冗余通道作为备用通道进行数据传输。

a∈r³×³、b∈r³×¹、ci∈r^1×3、l∈r¹×³为已知矩阵，i＝1,2,3；其中的c1为水流状态数据从主传输通道传输到远程排水控制中心的映射矩阵，c2和c3为水流状态数据从两条冗余传输通道传输到远程排水控制中心的映射矩阵。

γj(k)是相互独立的随机序列，满足伯努利分布，j＝1,2,3；其中的γ1(k)表示主传输通道传输数据时发生数据丢包的现象，γ2(k)和γ3(k)表示两条冗余通道传输数据时发生数据丢包的现象；是γj(k)的均值，即式中e{·}表示数学期望，prob{·}表示随机事件的概率。

非线性扰动f(x(k))∈r^3×1，满足初始条件f(0)＝0和假设m∈r^3×1和n∈r^3×1都是任意向量，是已知矩阵，||·||表示向量或矩阵的欧几里得范数；

利用实测数据和计算机仿真技术，反复进行模型校验和修正。

(2).利用系统测量输出，构造动态输出反馈控制器：其中为控制器的状态向量，表示向量x(k)的估计值；f∈r^3×3、g∈r^3×1、h∈r^1×3为待求的控制器增益矩阵。

(3).动态输出反馈控制器求解；

给出闭环控制系统随机有限时间有界的充分条件，通过求解线性矩阵不等式，得到动态输出反馈控制器增益。

(3-1).引入增广向量得到增广系统：

其中，

(3-2).构建二次型函数v(η(k))＝η^t(k)pη(k)；其中，p为维数适当的正定对称矩阵。

计算二次型函数v(η(k))＝η^t(k)pη(k)沿着增广系统的差分，取期望得到：

得到：

其中，

对于任意标量μ＞0，得到有：

得到其中标量β≥1；

*表示对称矩阵中对应的对称项；

有e{v(η(k+1))}≤e{βv(η(k))}成立，通过递归计算得到：

e{v(η(k))}≤e{β^kv(η(0))}≤e{β^kη^t(0)pη(0)}。

下面考虑前述增广系统的随机有限时间有界性，即对于给定的已知正数b1、b2和正定对称加权矩阵r，若初始条件满足e{η^t(0)rη(0)}≤b1，则对于给定的正整数u有e{η^t(k)rη(k)}≤b2；

令其中是给定的正定对称加权矩阵r的平方根的逆矩阵，则由给定的初始条件确定：e{v(η(k))}≤e{β^k(λ1η^t(0)rη(0))≤β^kb1λ1；

基于构造的二次型函数得到：其中，和分别表示矩阵的最大特征值和最小特征值。

对于任意k≤u，有k∈{1,2,…,u}；

若不等式λ1b1＜λ2β^-ub2成立，则保证对于给定参数b1，b2，u和矩阵r以及初始条件e{η^t(0)rη(0)}≤b1，闭环系统是随机有限时间有界的。

由schur补引理，等价于π<0，其中，

式中i表示维数匹配的单位矩阵，其中diag{…}表示块对角矩阵；如果同时满足π<0和λ1b1＜λ2β^-ub2，则可保证闭环系统随机有限时间有界；

由schur补引理，π<0等价于

(3-3).令hs＝hw^t，gs＝sg；其中x和y为待求的正定对称矩阵，非奇异矩阵s和w满足条件sw^t＝i-xy；引入矩阵和令p＝ψ2ψ1^-1，ψ1，ψ2均可逆，且p正定对称；

令同时用和它的转置矩阵分别左乘和右乘矩阵不等式ω<0，得到矩阵不等式：

其中，

由schur补引理，ξ＜0等价于φ+ωφ^t+φω^t＜0，φ表示λ的7阶顺序主子式，即

其中，

λ61＝[lyl],

由schur补引理，φ+ωω^t+φφ^t＜0成立，当且仅当λ<0成立。

对于给定参数的b1、b2、u、r、β≥1，如果存在正定对称矩阵x和y以及矩阵fs、gs、hs、正标量μ，使得线性矩阵不等式λ<0有解，且在给定的初始条件e{η^t(0)rη(0)}≤b1下，满足有限时间有界条件λ1b1＜λ2β^-ub2，则有限时间动态输出反馈控制器的增益矩阵为：

本发明针对现有排水系统控制方法难以在给定的有限时间之内实现及时、有效地控制，提出一种城市排水管道的有限时间控制方法。由于传感器节点测量数据众多，水流数据在传输过程中容易出现丢失现象，为了提高传送数据包的成功率，本发明使用了多条冗余通道的水流数据传输方法。基于随机有限时间有界性分析结果，通过求解矩阵不等式设计动态输出反馈控制器，从而在给定的有限时间内实现对城市排水系统的准确控制。

具体实施方式

一种城市排水管道的有限时间控制方法，包括：

(1).建立城市排水管道系统状态空间模型；

首先，基于流体力学原理，建立排水管道系统模型：

其中，符号∑表示求和运算，符号π表示连乘运算；

x(k)表示k时刻检测到的水流状态向量，x(k)＝[x1(k),x2(k),x3(k)]^t，x1(k)、x2(k)、x3(k)分别表示k时刻所测量的水位高度值、水流速度值和水压值，上标t表示矩阵的转置；

u(k)∈r^1×1表示k时刻的控制输入量，为排水管道上游进入管道的水流量与排水管道下游流出管道的水流量之差；

y(k)∈r^1×1表示k时刻排水控制系统的测量输出，为传感器测量到的水流状态通过网络传输到远程控制中心的测量值；

z(k)∈r^1×1表示k时刻排水控制系统待估计的输出信号；城市排水控制系统利用三个传输通道将各传感器节点的水务数据信息传输到集中控制中心，在检测到主传输通道发生数据丢失的情况下，冗余通道作为备用通道进行数据传输；

a∈r^3×3、b∈r^3×1、ci∈r^1×3、l∈r^1×3为已知矩阵，i＝1,2,3；其中的c1为水流状态数据从主传输通道传输到远程排水控制中心的映射矩阵，c2和c3为水流状态数据从两条冗余传输通道传输到远程排水控制中心的映射矩阵；

γj(k)是相互独立的随机序列，满足伯努利分布，j＝1,2,3；其中的γ1(k)表示主传输通道传输数据时发生数据丢包的现象，γ2(k)和γ3(k)表示两条冗余通道传输数据时发生数据丢包的现象；是γj(k)的均值，即式中e{·}表示数学期望，prob{·}表示随机事件的概率，通过实验和统计方法可以得到的值。

非线性扰动f(x(k))∈r^3×1，满足初始条件f(0)＝0和假设m∈r^3×1和n∈r^3×1都是任意向量，是已知矩阵，||·||表示向量或矩阵的欧几里得范数。

最后，利用实测数据和计算机仿真技术，反复进行模型校验和修正。

(2).动态输出反馈控制器结构；

利用系统测量输出，构造动态输出反馈控制器：其中为控制器的状态向量，表示向量x(k)的估计值；f∈r^3×3、g∈r^3×1、h∈r^1×3为待求的控制器增益矩阵。

通过随机有限时间有界性分析和线性矩阵不等式方法对动态输出反馈控制器求解。

(3).动态输出反馈控制器求解：给出闭环控制系统随机有限时间有界的充分条件，通过求解线性矩阵不等式得出动态输出反馈控制器增益。

(3-1).引入增广向量得到增广系统：

其中，

(3-2).构建二次型函数v(η(k))＝η^t(k)pη(k)；其中，p为维数适当的正定对称矩阵。

计算上述二次型函数沿着增广系统的差分，并取期望得到：

由于γj(k)是相互独立的伯努利变量，得到：

其中，

根据非线性函数的假设条件，对于任意标量μ＞0，得到从而有：

进而得到，其中标量β≥1；

式中*表示对称矩阵中对应的对称项。

因此，有e{v(η(k+1))}≤e{βv(η(k))}成立，通过递归计算得到：

e{v(η(k))}≤e{β^kv(η(0))}≤e{β^kη^t(0)pη(0)}。

下面考虑前述增广系统的随机有限时间有界性，即对于给定的已知正数b1、b2和正定对称加权矩阵r，若初始条件满足e{η^t(0)rη(0)}≤b1，则对于给定的正整数u有e{η^t(k)rη(k)}≤b2。

令其中是给定的正定对称加权矩阵r的平方根的逆矩阵，则由给定的初始条件可知：e{v(η(k))}≤e{β^k(λ1η^t(0)rη(0))≤β^kb1λ1；

基于构造的二次型函数得到：其中，式中和分别表示矩阵的最大特征值和最小特征值。

从而，对于任意k≤u，有：其中，k∈{1,2,…,u}。

因此，若不等式λ1b1＜λ2β^-ub2成立，则保证对于给定参数b1，b2，u和矩阵r以及初始条件e{η^t(0)rη(0)}≤b1，在设计的输出反馈控制器作用下，闭环系统是随机有限时间有界的。

根据schur补引理可知，等价于π<0，其中，

式中i表示维数匹配的单位矩阵，其中diag{…}表示块对角矩阵。因而，如果同时满足∏<0和λ1b1＜λ2β^-ub2，则可保证闭环系统随机有限时间有界。

根据schur补引理，∏<0等价于

(3-3).令hs＝hw^t，gs＝sg；其中x和y为待求的正定对称矩阵，非奇异矩阵s和w满足条件sw^t＝i-xy。引入矩阵和令p＝ψ2ψ1^-1，ψ1，ψ2均可逆，且p正定对称,右上标-1表示矩阵的求逆运算。

令同时用和它的转置矩阵分别左乘和右乘矩阵不等式ω<0，得到矩阵不等式：

其中，

由schur补引理，ξ＜0等价于φ+ωφ^t+φω^t＜0，φ表示λ的7阶顺序主子式，即：

其中，

λ61＝[lyl],

根据基本不等式ab^t+ba^t≤aa^t+bb^t，要使φ+ωφ^t+φω^t＜0成立只需φ+ωω^t+φφ^t＜0成立。根据schur补引理可知，φ+ωω^t+φφ^t＜0成立当且仅当λ<0成立。

对于给定参数的b1、b2、u、r、β≥1，如果存在正定对称矩阵x和y以及维数适当的矩阵fs、gs、hs，以及正标量μ，使得线性矩阵不等式λ<0有解，且在给定的初始条件e{η^t(0)rη(0)}≤b1下，所求结果满足有限时间有界条件λ1b1＜λ2β^-ub2，则本发明所设计的有限时间动态输出反馈控制器的增益矩阵为：