[0090] 下面结合具体实施例来对本发明进行进一步说明,但并不将本发明局限于这些具体实施方式。本领域技术人员应该认识到,本发明涵盖了权利要求书范围内所可能包括的所有备选方案、改进方案和等效方案。
[0091] 如图2所示,本实施例基于状态饱和的正切换系统建模一类停车场综合管理系统,并设计了一种事件触发异步滤波器用于估计各大型公共停车场的车流量,其具体步骤如下:
[0092] 步骤1、首先对停车场出入口的车流量进行数据采集,利用采集的数据建立停车场综合管理系统的状态空间模型,形式如下:
[0093] x(k+1)=sat(Aσ(k)x(k))+Bσ(k)ω(k),
[0094] y(k)=Cσ(k)x(k)+Dσ(k)ω(k),
[0095] z(k)=Eσ(k)x(k)+Fσ(k)ω(k),
[0096] 其中, 表示在第k个采样时刻停车场内的剩余车位数,n表示停车场的数量, 为驶入停车场的车辆数量,m表示停车场内的泊位数量, 表示在k时刻对y(k)的估计,即预计驶入停车场的车辆数量,s表示停车场的出入口数量,ω(k)是停车场出入口影响车流量的外部干扰,函数 是饱和函
T
数,被定义为sat(u)=[sat(u1),sat(u2),…,sat(um)] ,sat(ui)=sgn(ui)min{|ui|,1},i+
∈m,σ(k)是切换信号,其取值在一个有限集S={1,2,…,J}中,J∈Z ,
和 是已知的系统
矩阵,对于σ(k)=i,i∈S,系统矩阵满足 和
[0097] 步骤2、建立停车场出入口车辆拥堵的事件触发条件:
[0098] ||ey(k)||1>β||y(k)||1,
[0099] 其中,常数β>0,采样误差 表示采样状态,‖·‖1表示向量的1范数,即向量中所有元素的绝对值之和。
[0100] 步骤3、设计停车场综合管理系统的事件触发异步滤波器,其特征包括如下步骤:
[0101] 步骤3.1、设计事件触发异步滤波器,具体形式如下:
[0102]
[0103]
[0104] 其中,xf(k)为滤波器的状态信号,zf(k)是z(k)的估计,Afi,Bfi,Efi和Ffi是所设计的事件触发异步滤波器的增益矩阵,具体形式如下:
[0105]
[0106] 其中,1n是向量中所有元素均为1的n维向量, 表示仅第ι个元素为1、其余元素均T T T为0的n维向量,θ、h和ξ是n维向量,η和 是m维向量,θ、η、ξ和 分别表示θ、η、ξ和 的转置。
[0107] 步骤3.2、所述的基于状态饱和的停车场综合管理系统,其饱和函数满足:
[0108]
[0109] 其中, 且||H||∞≤1,||·||∞是无穷范数,表示矩阵中每行元素绝对值之和的最大值,Dl表示对角元素为0或1的n×n的对角矩阵,
I表示适当维数的单位矩阵。
[0110] 步骤3.3、定义xe(k)=xf(k)‑x(k),e(k)=zf(k)‑z(k)。根据步骤1、步骤3.1和步骤3.2,将停车场综合管理系统的状态空间模型与事件触发异步滤波器扩充为一个误差系统,具体如下:
[0111]
[0112]
[0113] 其中, 和 为停车场综合管理系统的状态空间模型与事件触发异步滤波器的系统矩阵所构成的増广形式的增益矩阵,具体形式为:
[0114]
[0115]
[0116] 其中,
[0117] 步骤3.4、具有状态饱和的停车场综合管理系统与滤波器构成的误差系统在事件触发机制下平稳运行的约束条件设计如下:
[0118] 设计常数0<μ1<1,μ2>1,λ>1,γ>0,如果存在n维向量使得下列不等式:
[0119]
[0120]
[0121]
[0122]
[0123]
[0124]
[0125]
[0126]
[0127]
[0128]
[0129]
[0130]
[0131]
[0132]
[0133]
[0134]
[0135]
[0136] 成立,那么在所设计的事件异步触发滤波器增益矩阵下,误差系统是正的,且是 增益稳定的。并且系统的切换规则满足:
[0137]
[0138]
[0139] 其中,κ‑(k0,k)表示滤波器与系统同步运行的总时间,κ+(k0,k)表示滤波器与系统异步运行的总时间,τa表示平均驻留时间,Δm表示滤波器滞后于对应子系统的最大滞后时间,Φ=I‑β1m×m,Ψ=I+β1m×m。
[0140] 进一步地,所述步骤3还包括下述步骤,用于验证所构造的误差系统在事件触发条件下的正性:
[0141] 步骤3.5、对于任意初始状态 和 步骤2中的事件触发条件在k0时刻满足:
[0142]
[0143] 其中,1m×m表示矩阵元素全为1的m×m的矩阵。
[0144] 步骤3.6、结合步骤3.3和步骤3.5,当k∈[kl,kl+Δl)时,有:
[0145]
[0146]
[0147] 当k∈[kl+Δl,kl+1)时,有:
[0148]
[0149]
[0150] 其中, 和 分别为异步时间和同步时间内下界系统的增广形式的增益矩阵,具体形式为:
[0151]
[0152]
[0153]
[0154]
[0155] 其中, i取p表示滤波器与对应的子系统同步,i取q表示滤波器与对应的子系统异步,相对应时刻的滤波器增益矩阵即为同步或异步时的矩阵。
[0156] 根据步骤3.4中的正性约束条件能够保证下界增益矩阵的正性,于是对于任意初始状态 下界系统是正的,因此误差系统是正的。
[0157] 进一步地,所述步骤3还包括下述步骤,用于验证所构造的误差系统在事件触发条件下的 增益稳定性:
[0158] 步骤3.7、结合步骤3.3和步骤3.5,当k∈[kl,kl+Δl)时,有:
[0159]
[0160]
[0161] 当k∈[kl+Δl,kl+1)时,有:
[0162]
[0163]
[0164] 定义
[0165]
[0166] 其中, 和 分别为异步时间和同步时间内上界系统的增广形式的增益矩阵,具体形式为:
[0167]
[0168]
[0169]
[0170]
[0171] 步骤3.8、设计线性余正李雅普诺夫函数:
[0172]
[0173] 其中,
[0174] 步骤3.9、结合步骤3.4、步骤3.7和步骤3.8,李雅普诺夫函数的切换满足:
[0175]
[0176] 其中,
[0177] 因此,上式通过递推可以得出:
[0178]
[0179] 步骤3.10、对于时间 由步骤3.9通过递推法可以得到如下不等式:
[0180]
[0181] 其中,Nσ(k0,T)表示切换信号σ(k)在时间[k0,T)内的切换次数,并且在平均驻留时间条件下满足
[0182] 步骤3.11、由于 令M=T‑1,于是在零初始条件下,下式成立:
[0183]
[0184] 步骤3.12、将步骤3.11中不等式两边同时乘以 可以得到:
[0185]
[0186] 步骤3.13、由步骤3.4中平均驻留时间条件和步骤3.10有:
[0187]
[0188] 因此,步骤3.12可以转化为:
[0189]
[0190] 将上式两边同时在[0,∞)求和,同时由于 可得如下不等式:
[0191]
[0192] 其中,
[0193] 因此,误差系统在性能指标γ下满足 增益性能。