多采样率MIMO网络时滞控制系统的建模与稳定性研究.docx

1、多采样率MIMO网络时滞控制系统的建模与稳定性研究多采样率MIMO网络时滞控制系统的建模与稳定性研究摘 要：首先建立一个MIMO网络控制系统，然后利用提升技术使系统由线性时变的多采样率系统转化为高维的线性时不变的单采样率系统，推导出了该类网络控制系统的时滞离散时间数学模型。将整个MIMO网络控制系统分成若干个子系统进行稳定性分析，若每一个子系统都是稳定的，则整个系统是稳定的。最后给出了系统稳定的条件，并导出了使系统稳定的最大时延。通过MATLAB进行系统仿真，实验结果证明本文所设计的多采样率MIMO网络控制系统是稳定的。关键词：多采样率；MIMO网络控制系统；时滞；稳定性 0 引言 网络控制系

2、统（NCS， Networked Control System)是指以网络作为信息传输的通道从而将被控对象、传感器、控制器和执行器连接起来而形成的反馈控制系统。目前，对MIMO(Multi Input & Multi Output, MIMO)网络控制系统的研究也有了一些研究成果。Zhang W针对传感器采用多包传输的网络控制系统，在不考虑网络时延情况下建立了切换系统模型，分析了传感器数据封装为两个数据包时系统的稳定性；李静等针对MIMO网络控制系统中同时存在时延、数据包丢失以及多包传输问题进行了研究，根据数据寄存和静态调度方法建立系统模型，并利用Lyapunov函数和线性矩阵不等式方法(LM

3、Is)导出了系统稳定的充分条件；张俊等针对多延时的MIMO网络控制系统进行了建模与稳定性分析，并用Lyapunov函数和Razumikhin定理得出了系统稳定的时延参数和稳定性条件。 本文将针对一类多采样率MIMO网络时滞控制系统，利用提升技术建立准确合理的多采样周期 NCSs 数学模型， 并对所建立的模型进行了稳定性分析，并利用Lyapunov函数和线性矩阵不等式方法(LMIs)导出了系统稳定的充分条件。通过MATLAB仿真实验，证明系统在多采样率的情况下，不仅具有较强的分散控制能力并仍可以保持较高的稳定性。 1 多采样率MIMO网络控制系统的数学模型 MIMO网络控制系统结构图如图1所示。

4、在建立数学模型时先做如下假设：假设同一回路采用相同的周期采样和保持，且各回路状态参数之间相互独立；通过网络传输数据时，同一回路的输入或输出数据用单包传输，不考虑数据包错序与丢失；假设传感器、控制器和执行器都为时钟驱动。 图1 MIMO网络控制系统结构 被控对象为线性时不变连续控制系统，其状态方程描述如下： p(t)=A px p(t)+B pu p(t) y p(t)=C px p(t) (1) 其中， x p(t)R n p ，u p(t)R m，y p(t)R r分别是被控对象的状态向量、输入向量和输出向量。设第i个输入通道的输入采样周期为T i，T i=q iT,其中T为基本采样周期，整

5、个系统的循环周期为T 0=qT，q=LCM(q 1,q 2,q r)。以T 0 为采样周期，并设图1中所有的保持器均是零阶保持器，可以得到连续时间被控对象 G p 在采样各采样点时刻的离散时间状态方程为： x p(k+1)T 0=Ax p(kT 0)+Bu p(kT 0) y p(kT 0)=Cx p(kT 0) (2) 被控对象状态方程中的系数矩阵都是以 T 0 为周期的时变矩阵，这就给闭环系统的分析和控制系统的分析带来了很大的不便，为了克服这一困难，可以利用提升技术定义被控对象状态方程中的扩展输入与输出向量： u p(kT 0)=u 1 p(kT 0) u 1 p(kT 0+T 1) u

6、1 pkT 0+(M 11)T 1 u m p(kT 0) u m p(kT 0+T m) u m pkT 0+(M m1)T m ， y p(kT 0)=y 1 p(kT 0) y 1 p(kT 0+T 1) y 1 pkT 0+(M 11)T 1 y m p(kT 0) y m p(kT 0+T m) y m pkT 0+(M m1)T m 同时定义相应的矩阵： A p=A； B p=B 1,B m为n prj=1M j矩阵， C p为mi=1M in p 矩阵。 由此我们得到多采样率MIMO网络控制系统被控对象的状态空间描述为： x p(k+1)T 0= A px p(kT 0)+ B

7、p u p (kT 0) y p (kT 0)= C p x p(kT 0) (3) 若考虑控制器的处理时延，即采样 c 来表示，多采样率MIMO网络控制系统控制单元采用状态反馈方程控制方式，其状态反馈控制方程可以描述为： x c(k+1)T 0=A cx c(kT 0)+B c u c (kT 0) y c (kT 0)= K c u c(kT 0 c) (4) 其中， x cR n c ， u c, y cR mi=1M i 分别是控制器的状态、输入和输出向量。当考虑网络时延且其值小于一个采样周期时，根据图1的信号关系，并利用(3)和(4)式，则： u c(kT 0 c(kT 0)=y 1

8、 p(kT 0 11 sc 11 c) y 1 p(kT 0 12 sc 12 c) y 1 p(kT 0 1M 1 sc 1M 1 c) y m p(kT 0 m1 sc m1 c) y m p(kT 0 m2 sc m2 c) y m p(kT 0 mM m sc mM m c)= ri=1E i ix p(kT 0 11 sc 11 c) x p(kT 0 12 sc 12 c) x p(kT 0 1M 1 sc 1M 1 c) x p(kT 0 m1 sc m1 c) x p(kT 0 m2 sc m2 c) x p(kT 0 mM m sc mM m c) (5) u p(kT 0)

9、= y c(kT 0 ca (kT 0)=m=1ri=1F E i x p(kT 0 11 sc 11 ca 11 c) x p(kT 0 12 sc 12 ca 12 c) x p(kT 0 1M 1 sc 1M 1 ca 1M 1 c) x p(kT 0 m1 sc m1 ca m1 c) x p(kT 0 m2 sc m2 ca m2 c) x pkT 0 mM m sc mM m ca mM m c) (6) 将上面 u p(kT 0)与x p 的关系式代入式被控对象Gp和控制单元Gc状态方程中，被控对象Gp状态方程变为： x p(k+1)T 0= A px p(kT 0)+ B Pm

10、=1ri=1F E i x p(kT 0 11 sc 11 ca 11 c) x p(kT 0 12 sc 12 ca 12 c) x p(kT 0 1M 1 sc 1M 1 ca 1M 1 c) x p(kT 0 m1 sc m1 ca m1 c) x p(kT 0 m2 sc m2 ca m2 c) x p(kT 0 mM m sc mM m ca mM m c) (7) 控制单元Gc的状态方程变为： x c(k+1)T 0=A cx c(kT 0)+B cri=1E i x p(kT 0 11 sc ) x p(kT 0 12 sc ) x p(kT 0 1M 1 sc ) x p(kT

11、 0 m1 sc ) x p(kT 0 m2 sc ) x pkT 0 mM m sc (kT 0) (8) 取 x(kT)=x p(kT 0) x c(kT 0)R n p+n c 作为网络控制系统的广义状态空间向量，则有： x(k+1)T 0= A p0 0A cx p(kT 0) x c(kT 0)+00 B cri=1E i0 x p(kT 0 11 sc ) x p(kT 0 mM m sc ) x c(kT 0 11 sc ) x c(kT 0 mM m sc )+ B Pm=1ri=1F E i0 00 x p(kT 0 11 sc 11 ca 11 c) x p(kT 0 mM

12、 m sc mM m ca mM m c) x c(kT 0 11 sc 11 ca 11 c) x c(kT 0 mM m sc mM m ca mM m c) (9) 其中，令 A 0= A p0 0A c ， A 1=00 B cri=1E i0 ， A 2= B Pm=1ri=1F E i0 00 。 2 网络控制系统的稳定性分析 多采样率多输入多输出网络控制系统是由多个采样和保持周期不同的控制回路组成，若将每个回路看成一个子系统，若其中的任意一个子系统稳定则整个系统是稳定的。 从系统(9)中取其中的第 j 条子系统网络控制回路来分析，其状态方程描述为： x(k+1)T 0=A j 0

13、x(kT 0)+A j 1x(kT 0 j sc )+ A j 2x(kT 0 j sc j ca j c) (10) 其中， A j 0= A p0 0A c，A j 1=B cE i0 0，A j 2= B pF E i0 00 。 为便于讨论令 kT 0=k ，则式（10）可以表示成： x(k+1)=A j 0x(k)+2i=1A j ix(k i) (11) 网络中由于受网络带宽的限制，总会存在一定大小的时延，由上述式（11）为一个时延不一致的时滞系统。 子系统（11）稳定的条件是： (A j 0)r(1) (12) 网络控制系统极点落在以坐标原点为圆心，以 r为半径的圆形区域内，则称

14、该子系统为稳定的。其中=1r i+1 2i=1Q 1 0A j iQ 0，0r1，A j 0可对角化，且Q 1 0A j iQ 0=diag (A j 0)， (•)表示矩阵第 个特征值。 证明：由线性控制系统的相关理论知识可知，系统稳定的充要条件是系统极点必须分布在以坐标原点为圆心，r为半径的D(0,r)，若极点值落在圆外，则系统不稳定。 |z|r det (zIA j 02i=1A j iz i )0 i0, (13) 若令 z=r ，则可得 1 并将之代入上述式(13)，得： 1r(A j 0)+1r2i=1A j i(r) i (14) 显然，只要不等式： 1r(A j 0

15、)+1r2i=1A j i(r) i 1 (15) 成立，子系统极点便分布在圆 D(0,r) 内。 由文献中引理和定理和上述知 1，所以0r|r| 。则式(16)可变为： 可以得到： 1r(A j 0)1r(A j 0)+1r2i=1A j i(r) i 1r2i=1Q 1 0A j iQ 0(r) i = (16) 由上述说假设的条件 (A j 0)r(1) 可知则： 1r(A j 0)1 (17) 由式（16）变化可得： a1r(A j 0)1r(A j 0)+1r2i=1A j i(r) i 1r(A j 0)+1r2i=1A j i(r) i 1r(A j 0) (18) 由式(17)

16、和式(18)可知： 1r(A j 0)+1r2i=1A j i(r) i +1r(A j 0) (1+)=1 (19) 即由上述式（19）可知，该子系统极点分布在以原点为圆心的单位圆内，所以可证上述整个多采样率MIMO网络时滞系统是稳定的。 3 仿真实例 假设多采样率控制系统的线性时不变的状态空间模型为： x(k+1)T 0=A rx(kT 0)+B ru r(kT 0) y r(kT 0)=C rx(kT 0)+D ru r(kT 0) (21) 其中系统中各常数矩阵为： A r=1.37670.42630.0413 0.80470.99170.1925 1.99470.71541.0603

17、； B r=0.15070.06540.13250.11760.1054 0.12480.03390.09830.07630.058 0.25090.08460.18330.12470.0735 C r=101 1.47830.21831.0053 2.02220.47681.0225 2.64710.78201.0538； D r=00000 0.17880000 0.242300.178800 0.315800.24230.17880 选择所要配置的特征向量集合为：0.1227,0.5219,0.6524，它所对应的特征向量构成的矩阵 P 为： p=0.76280.35370.4453 0

18、.19250.02750.8935 0.61760.93500.0576 该系统中 C r为列满秩，B r为行满秩，根据 2 = defB + r(pp 1 A r)和K(C r+D r 2)= 2两式可可计算出闭环系统的反馈增益K 为： K=4.932000 0.0207000 04.229800 1.12310.72144.65150 010.8761.84622.0211 再结合上述稳定性分析可知： =1,1.8326,0.6235 ； Q 0=1.0000.02730.0256 00.08460.0813 00.95380.8724 ； A 1=0.12031.23160.3216 0

19、.23645.6211.832 0.33616.12422.1347 。 根据上述所得结果与式（12）可知，当取半径 r=0. 9 •，使系统稳定时的最大时延的取值范围为： max 0.035s。 本文在MATLAB6.5环境下，得出了多输入多输出网络控制系统的状态曲线，如图2所示，实验结果表明多采样率在MIMO网络控制系统的应用能够使系统的具有较好的稳定性。 图2 多输入多输出NCSs 的状态曲线图 4 结束语 本文对一类多采样率MIMO网络时滞控制系统进行了研究，通过采用提升技术将多采样率系统的时变特性转化为高维的时不变系统，建立基于状态反馈的离散时间数学模型。从此数学模型中取

20、其中任意一个子系统来分析，根据线性控制系统稳定性理论的相关知识，得到了单个子系统的极点分布在以原点为坐标原点的单位圆内，从而证明了整个系统的稳定性，并求出了使系统稳定的最大时延。最后用MATLAB进行系统仿真，证明了本文所设计的多采样率MIMO网络控制系统是稳定的。 参考文献： 1 LIAN FENGLI,JAMES MOYNE,DAWN TIBURY.Analysis and modeling of networked control system:MIMO case with multiple time delayC.Proceeding of the American Conteol C

21、onference. Arlington,VA,USA,2001. 2 MA X H,XIE J Y,WEI Z.Research on network delays of networked control systemsJ.Journal of shanghai Jiaotong University (Science),2005(1). 3 ZHANG W,BRANICKY M S,PHILLIPS S M.Stability of networked control systemsJ.IEEE Control Systems Magazine,2001(1). 4 刘鲁源,吕伟杰,陈玉柱.MIMO网络控制系统的稳定性分析J.信息与控制,2006(3). 5 李静,左斌,胡云安,等.MIMO网络控制系统的建模与稳定性分析J.华中科技大学学报（自然科学版）,2008(10). 6 张俊,罗大庸.一类多延时多输入多变量网络控制系统的建模及稳定性 J.控制与决策,2011(8). 7 肖建.多采样率数字控制系统M.北京,科学出版社,2003.