当前位置:
首页 > 教材教辅 > 大学 > 鲁棒控制理论及应用pdf/doc/txt格式电子书下载

鲁棒控制理论及应用pdf/doc/txt格式电子书下载

本站仅展示书籍部分内容

如有任何咨询

请加微信10090337咨询

鲁棒控制理论及应用pdf/doc/txt格式电子书下载

书名:鲁棒控制理论及应用pdf/doc/txt格式电子书下载

推荐语:

作者:王娟、

出版社:电子工业出版社

出版时间:2011-08-01

书籍编号:30466675

ISBN:9787121142857

正文语种:中文

字数:58323

版次:1

所属分类:教材教辅-大学

全书内容:

鲁棒控制理论及应用pdf/doc/txt格式电子书下载






前言


鲁棒控制理论是处理不确定性问题的有力工具,近年来备受控制理论界关注,取得了一系列的成果和方法,并在工程领域中获得了成功的应用。该问题的研究不仅可以应用于飞行器控制、机器人控制、网络控制等其他工程领域,也可以应用于运动控制领域。


本书全面系统地介绍了不确定系统鲁棒figure_0004_0004控制和时滞系统鲁棒控制的理论基础、各种设计方法及其在工程设计中的应用等问题。主要研究内容包括:约束系统的状态及输出反馈控制方法、时滞系统的稳定性及镇定问题、鲁棒figure_0004_0005输出反馈控制器的设计方法。特别是在约束系统的鲁棒figure_0004_0006控制设计中,对工程中常见的约束控制系统的鲁棒控制设计问题,采用滚动时域策略,提出了几种figure_0004_0007鲁棒输出反馈控制器设计技术。最后给出了这些理论方法在控制工程设计中的应用实例。


本书是研究团队成员大量研究成果的总结,融入了近年来在国内外刊物上发表及国际会议上交流的多篇学术论文的精华内容。杜海英负责编写了第1章,王娟负责编写了第2、3、4章,张涛负责编写了第5章,崔艳秋负责编写了第6章和附录,徐国凯负责编写了第7章。全书由王娟统稿。


本书的出版得到了辽宁省教育厅项目和中央高校自主科研基金资助,作者对此深表谢意。


由于作者水平有限,书中的缺点错误在所难免,欢迎读者给予批评指正。


编著者

第1章 控制系统的稳定性


在控制理论中系统的稳定性是一个非常重要的概念,它是系统能否正常工作的最基本条件,因此研究系统稳定性、稳定性条件及稳定性措施是控制系统的重要内容。

1.1 稳定性概念


当系统为非线性和时变的最一般情况时,可用如下状态空间模型描述

figure_0008_0008

其中,x为n维状态向量, f(⋅)和g(⋅)分别是状态的非线性向量函数和矩阵函数,u则表示控制输入量。进而,如果系统为定常系统,则其状态方程(1.1)中将不显含t;如果系统为线性,那么方程(1.1)中 f (⋅)和g(⋅)分别为x和u的线性向量函数,此时式(1.1)化为

figure_0008_0009

研究运动稳定性问题时,我们希望系统的状态x(t)对于任意的非零初始值figure_0008_0010,都能够渐近地趋于平衡状态并且稳定在平衡点上,即

figure_0008_0011

由此可见,平衡状态即系统方程的常数解,或者系统的一种静止的运动。在大多数情况下, xe =0 ,即状态空间的原点为系统的平衡状态;除此之外,系统也可以有非零平衡状态。显然对于一个非线性系统,平衡点并非唯一。容易证明,对于孤立的平衡状态,总是可以通过移动坐标系而将其转换为空间的原点。所以为了叙述简便,下面只讨论xe =0的情况。


对于一般系统而言,要想使系统具有上述运动趋势,就必须施加控制作用。因为如果不施加驱动力,即令u=0(通常称这类系统为“自治系统”),则系统状态的动态过程将满足微分方程

figure_0008_0012

但是对于任意给定的初始条件x(t0 )≠0 ,该方程的解x(t)一般不满足式(1.3),因此需要选择适当的控制输入u使系统稳定到平衡状态。如果状态x的各元素均可检测,就可以通过状态反馈来确定控制输入,令

figure_0009_0013

其中,k为函数矩阵。这时系统的状态微分方程为

figure_0009_0014

则由初始状态所引起的运动轨迹满足上述微分方程的解,所以可以选择适当的k使得上述微分方程的解对于任意给定的初始状态都趋于零。


式(1.5)称为“状态反馈控制律”或者“控制器”,而式(1.6)称为“对应的闭环系统”,式(1.4)称为“开环系统”。显然,此时系统的控制问题就是设计控制器,使得其偏离平衡状态的受扰运动能够返回到平衡点状态,或者限制在它的有限领域内。一个实际的系统必须是稳定的,只有稳定才能付诸工程应用。虽然稳定还不是系统达到成功的工程应用的全部,但是它是系统实际应用的前提和先决条件。下面将首先讨论李雅普诺夫(Lyapunov)意义下稳定的基本概念,在此基础上研究系统的内部稳定性、外部稳定性,以及各种稳定性之间的关系和稳定性的判据;同时,还将讨论利用李雅普诺夫直接法进行系统综合的问题。


定义1.1 若个给定的状态xe满足

figure_0009_0015

则称xe是系统式(1.1)的“平衡状态”或者“平衡点”


如果没有外力作用于系统,则系统将保持平衡状态;如果系统受到外力作用,则系统能否保持这个平衡状态就是平衡状态的稳定性问题。以下给出李雅普诺夫稳定性定义。


定义1.2稳定性设figure_0009_0016为动力学系统式(1.1)的平衡态,若对任意给定的实数figure_0009_0017>0,都对应地存在一个实数figure_0009_0018,使得一切满足figure_0009_0019的系统响应x(t),在所有的时间内都满足figure_0009_0020,则称“系统的平衡状态xe是稳定的”。在此定义中,有

figure_0009_0021

figure_0009_0022选得足够小时,则由初始扰动引起的响应,在所有时间内都包括在一个超球中,即

figure_0009_0023

在二维情况下,系统稳定性的几何解释如图1-1所示。

figure_0010_0024

图1-1 系统稳定性的几何解释

如果对平衡点xe和任意给定的域figure_0010_0025,找不到满足稳定条件的相对邻域figure_0010_0026,那么系统在该平衡点是不稳定的,则称“系统是不稳定的”。


如果所取的域figure_0010_0027figure_0010_0028初始时刻t0无关,即对任何t0稳定条件不变,则称“该系统稳定状态是一致稳定的”。


定义 1.3 渐进稳定 如果系统的一个平衡状态 xe是稳定的且对靠近平衡状态 xe的任何初始点x(0)的系统解x(t)满足

figure_0010_0029

也就是说,从足够靠近xe处出发的每一个解x(t),当t→figure_0010_0030时收敛于xe,则平衡状态xe是渐进稳定的。在二维情况下,系统渐进稳定性的几何解释如图1-2所示。

figure_0010_0031

图1-2 系统渐进稳定性的几何解释

如果 xe是渐进稳定的,且系统稳定性质与初始时刻 t0无关,则系统是一致渐进稳定的。


定义1.4 全局渐进稳定 如果系统的平衡状态xe对状态空间的所有 x(0)都是稳定的且figure_0011_0032,则xe是全局渐进稳定的。


如果系统的稳定性质与初始时刻t0无关,则平衡状态xe是一致全局渐进稳定的。


从工程观点来看,感兴趣的往往是在平衡状态附近比较大范围内的稳定性问题。若在状态空间中有限的范围内满足稳定条件,则系统称为“局部稳定”。

1.2 李雅普诺夫稳定性定理


李雅普诺夫第二方法又称“李雅普诺夫直接法”,这一方法的优越性在于可在不解微分方程的条件下确定系统的稳定性。该方法建立在能量观点的基础上,若系统的某个平衡状态是渐近稳定的,则随着系统的运动其存储的能量将随时间增长而不断衰减。直至figure_0011_0033时,系统运动趋于平衡状态而能量趋于极小值。由此,李亚普诺夫创立了一个可模拟系统能量的“广义能量”函数,根据这个标量函数的性质来判断系统的稳定性。由于该方法不必求解系统的微分方程就能直接判断其稳定性,故又称为“直接法”,其最大的优点在对于任何复杂系统都适用。


对于由式(1.1)描述的系统,如果

figure_0011_0034

则系统可能的平衡状态为xe =0 ,即为坐标原点0。


为了分析系统的稳定性,李雅普诺夫引出一个虚构的能量函数,称为“李雅普诺夫函数”。分析这一函数及其一次导数的定号性而获得系统稳定性的有关信息。李雅普诺夫第二方法概念直观,方法具有一般性且物理意义清晰。下面介绍李雅普诺夫函数和李雅普诺夫稳定性定理及其应用。


李雅普诺夫函数V (x,t)是状态变量的函数,它是一个广义的“能量函数”,并应具有下列性质:


(1)V(x,t)的一次偏导数是连续的;


(2)V(x,t)是正定的,即对所有figure_0011_0035,有V(x,t)>0且V(0)=0。


对于线性系统寻找李雅普诺夫函数比较简单,选为系统状态变量x的二次型函数,即

figure_0011_0036

为了V (x,t)正定,P应是正定矩阵。


对于一个复杂的系统,要想找到一个合适的李雅普诺夫函数比较困难。设系统方程为式(1.1)并设 f (0,t)=0 ,即原点为系统的平衡点,对此系统有如下定理。


定理1.1 李雅普诺夫稳定性定理


对于式(1.1)描述的系统,如果在包含原点 O 在内的某个域 D 内存在李雅普诺夫函数V(x,t)>0,而且figure_0012_0037,则系统在原点O是稳定的。


证明:用Dε表示中心为原点O,半径为figure_0012_0038的球的内部,以figure_0012_0039表示此球的表面。设V (x,t)为正定,figure_0012_0040,V(x,t)在figure_0012_0041上的最小值为m;以figure_0012_0042表示中心为原点O,半径为figure_0012_0043的球的内部,并设在球体figure_0012_0044内各点的V(x,t)<m。设在figure_0012_0045之中有任意一点x ,下0面来研究由x0出发的轨线x(t)=x(x0,t)的变化情况。假定在某一时刻t1,x(t1)穿过figure_0012_0046,由于

figure_0012_0047

故沿运动轨线x(t),V(x,t)并不增加,因此有figure_0012_0048;另一方面,m表示在figure_0012_0049上V(x,t)的最小值,如果x(t1)超出figure_0012_0050的边界,则必有figure_0012_0051,上面的结果是矛盾的。由此可以说明随着时间的延长,x(t)不能超出figure_0012_0052的边界。因此系统在原点稳定,定理得证。


定理1.2 李雅普诺夫渐近稳定性定理


对式(1.1)描述的系统,如果在包含原点 O 在内的某个域 D 内,存在李雅普诺夫函数V (x,t)>0 ,且figure_0012_0053,则系统在原点O是渐进稳定的。


证明:设表示中心为原点O,半径为r的球的内部。若x0在Dr之内,当t>0时,设x(t,x0)不会跑出Dr之外;设figure_0012_0054为任意小的正数,由前面的定理可有figure_0012_0055。当x0figure_0012_0056之内时,x(t,x0)在figure_0012_0057内。今设x0在Dr之内,设在t>0时不落入figure_0012_0058,则x(t,x0)将落入D


δ


和Dr两球之间,即figure_0012_0059与Dr之间的区域。但在此中间层内,figure_0012_0060,因此存在某个常数k>0,使得中间层内figure_0012_0061,于是由下式

figure_0012_0062

可得

figure_0012_0063

当t不断增长时,V (x,t)渐渐减小。但V (x,t)为正定函数,不可能变负,而是渐渐趋近于0。经过一段时间之后,x(t,x0)将落入figure_0012_0064figure_0012_0065是这样选择的,当x(t,x0)落入figure_0012_0066之后不能跑出figure_0012_0067,即figure_0012_0068figure_0012_0069之内。由于figure_0012_0070是任意小的数,&nbsp;figure_0012_0071可趋近于0,因此figure_0012_0072,即系统渐进稳定,定理得证。

1.3 线性离散系统李雅普诺夫稳定性分析


前面我们以连续时间系统为对象,讨论了系统的零输入响应即自由运动的稳定性。现在我们把这一讨论推广到离散时间系统,讨论李雅普诺夫第二方法在线性离散系统稳定性分析和系统状态反馈设计中的应用。


设系统的状态方程为

figure_0013_0073

其中,X为n维状态向量,A为n × n维常数矩阵。


定理1.3 离散系统的渐进稳定 对于离散系统式(1.14),如果存在一个对于X(k)的标量函数V(X(k)) ,对于任意的X(k)满足如下条件:


(1)V(X(k))为正定;


(2)figure_0013_0074为负定。


则原点为平衡状态,即X =0为渐近稳定。


定理 1.3 给出了构造线性离散系统李雅普诺夫函数的通用方法。本章主要是基于李亚普诺夫直接法设计使线性定常系统稳定的控制器。下面应用这个定理讨论离散系统的状态反馈控制器设计问题。


设不稳定的线性定常系统的状态方程为

figure_0013_0075

其中,状态反馈律具有下述形式

figure_0013_0076

系统式(1.15)在状态反馈律如式(1.16)作用下闭环系统为

figure_0013_0077

其中, Ac =A+BK。


对于离散系统式(1.17),选下列二次型函数为系统的李雅普诺夫函数

figure_0013_0078

其中,P为n×n维对称正定矩阵,求V的增量

figure_0013_0079

由于V(X)取正定,如果要使系统渐进稳定,必须使AcTPAc+P为负定,即要求

figure_0013_0080

成立,证明完毕。


解不等式(1.20)即求解一个代数的矩阵不等式,这并不是一件容易的事情。20世纪90年代初,随着求解凸优化问题的内点法的提出,线性矩阵不等式再次受到控制界的关注并被应用到系统和控制的各个领域中。众所周知,线性矩阵不等式在控制领域得到了广泛的应用。因为它克服了黎卡提(Riccati)方程求解困难,含有不定参数,需要人为确定参数的缺点。随着求解线性矩阵不等式的内点法的提出,可以求解线性矩阵不等式的软件包的大量出现,许多控制问题都可以采用线性矩阵不等式这一工具来处理,它通过将问题转化为一个线性矩阵不等式的可行性问题或者一个具有线性矩阵不等式约束的凸优化问题。由于这里利用Matlab的线性矩阵不等式(LMI)工具箱将渐近稳定性问题等价于线性矩阵不等式的可行性问题。


式(1.19)和式(1.20)中存在两个未知矩阵变量P和K,且这两个矩阵变量以非线性的形式出现在这个矩阵不等式中,因此要直接从以上矩阵不等式中求解P和K是很困难的。可以通过矩阵的舒克尔补(Schur Complement)性质将矩阵不等式进行等价变换,然后通过一个适当的变量替换。将非线性矩阵不等式(1.20)转化为一个等价的关于新变量的线性矩阵不等式,从而可以利用Matlab的LMI(线性矩阵不等式)工具箱求得可行解。


舒克尔补性质非常有用,借助它可将某些二次非 LMI 转化为 LMI。根据矩阵的舒克尔补性质,矩阵不等式(1.20)等价于

figure_0014_0081

其中, Q2=P2−1。用diag(Q1, I )分别左乘和右乘不等式两边,并定义KQ1=Y ,得到

figure_0014_0082

矩阵不等式(1.22)是一个线性矩阵不等式,因此可以应用 LMI 工具箱中的求解器feasp来判断该线性矩阵不等式的可行性,进而可以得到使系统式(1.14)渐近稳定的反馈控制器

figure_0014_0083

1.4 约束系统的稳定性分析


控制输入饱和(执行器饱和)是控制系统最为常见的约束。引入饱和后系统具有局部渐近稳定特性,其稳定性分析的一个重要问题是估计系统的吸引域。如何获得尽可能大的吸引域估计,即在吸引域内得到尽可能大的不变集,以减少估计的保守性成为研究的热点问题。

1.4.1 状态反馈控制系统稳定性分析


1.离散时间系统


考虑约束离散系统

figure_0015_0084

其中,x∈Rn为状态变量,u∈Rm为控制输入,A和B是适当维数的常数矩阵。函数figure_0015_0085:Rm→Rm为标准饱和函数,即figure_0015_0086,其中figure_0015_0087,给定状态反馈控制器

figure_0015_0088

定义整数集M ={1, 2,…, m},并且Bi表示矩阵B的第i列, Ki表示矩阵K的第i行,则相应的闭环系统可以表示为

figure_0015_0089

定义1.5 对于x(0)=x0∈Rn为系统式(1.25)的状态轨迹,如果从某个有限的区域内出发的系统状态轨迹始终保持在这个区域内,那么这个区域被称为“不变域”。


定义1.6给定对称正定矩阵P∈Rn×n和一个正数figure_0015_0090,记椭圆figure_0015_0091=figure_0015_0092。选择V=xTPx为Lyapunov函数,若对于所有的figure_0015_0093,有figure_0015_0094,则椭圆figure_0015_0095被称为“收缩不变”。


定义1.7 由初始状态0(0)nx=x∈R出发的系统式(1.25)状态轨迹figure_0015_0096的集合

figure_0015_0097

被称为“系统的原点吸引域”。


由于获得吸引域的确切值比较困难,因此获得尽可能大的吸引域估计值是非常重要的。显然如果figure_0015_0098是收缩不变集,那么椭球figure_0015_0099位于吸引域内,即饱和系统式(1.25)的吸引域可以依靠构造这样的椭球不变集来估计。


借助椭球不变集估计系统的吸引域,将归结为以下两个问题:


(1)一个椭球集是位于吸引域内不变集的条件;


(2)最大椭球不变集的计算方法。


下面分别对这两个问题进行讨论。这一部分首先通过引入多个辅助反馈矩阵的形式对饱和项进行处理,然后利用两个松弛矩阵,结合 Lyapunov 稳定性理论给出保守性较小的闭环系统局部渐近稳定的条件。


定义 1.8[1] 给定一个整数集M ,集合N是M的所有的子集的集合,记为N={S:S⊆M}。Sc表示属于M中集合S的补集,记为Sc ={i∈M:i∉S}。


显然集合 N 中共有2m个S j ,j=1,…,2m。若S j∈N ,则S cj∈N。利用定义 1.8,可以引入多个辅助矩阵figure_0016_0100来处理饱和非线性。用figure_0016_0101表示矩阵figure_0016_0102的第i行,给出下面的饱和特性的多参数描述形式。


引理1.1 给定z∈Rmfigure_0016_0103那么

figure_0016_0104

证明:

figure_0016_0105

对于每个figure_0016_0106,可能存在如下情况:

figure_0016_0107

综合上面的分析可以得到

figure_0016_0108

为了利用引理1.1中描述的饱和特性的多参数描述形式估计闭环系统式(1.25)的吸引域,引入两个松弛变量矩阵P2和P3,容易得到

figure_0016_0109

(1.26)

通过这种方法,可以将用2m个辅助变量figure_0016_0110来描述的饱和非线性特性引入到稳定性条件figure_0016_0111中得到定理1.4。为了简化,下文中对称矩阵中的符号“*”表示沿对角线对称位置的块矩阵的转置。


定理1.4 考虑带有给定状态反馈控制式(1.24)的系统式(1.23),对于每一个Sj ∈N,j=1,…,2m ,若存在正定矩阵P1>0和矩阵P2 ,以及P3figure_0017_0112∈Rm×n,当x(k)∈figure_0017_0113(P1,1)时,使矩阵不等式(1.27)和式(1.28)成立,即

figure_0017_0114

figure_0017_0115

其中,figure_0017_0116表示figure_0017_0117的第i行。那么椭球集figure_0017_0118位于吸引域内,并且是严格不变的。


证明:考虑二次函数V=xTP1x,P1>0。若figure_0017_0119位于吸引域内,必有

figure_0017_0120

引入等式条件[2]

figure_0017_0121

其中,P2和P3是适当维的两个任意矩阵,那么

figure_0017_0122

figure_0017_0123满足条件figure_0017_0124, i∈Sj时,如果

figure_0017_0125

根据引理1.1,可以得到

figure_0017_0126

将公式(1.31)写成下面的矩阵不等式形式

figure_0017_0127

其中,figure_0017_0128figure_0017_0129。由于figure_0017_0130,故figure_0017_0131,即x(k)∈figure_0018_0132位于吸引域内。


figure_0018_0133时,figure_0018_0134成立的充要条件是[3]

figure_0018_0135

应用Schur\'s补公式,式(1.32)与不等式(1.28)等价。


如果限定P1=P3, P2=0 ,且令figure_0018_0136=H , j=1,…,2m ,那么得到推论1.1。


推论1.1考虑系统式(1.25),给定椭球figure_0018_0137,若存在矩阵H∈R并满足m×n

figure_0018_0138

figure_0018_0139

那么椭球集figure_0018_0140位于吸引域内,并且是严格不变的。


在以上讨论中给出了一种用不变椭球集描述吸引域的方法。下面的主要工作是在所有满足定理 1.4 条件的不变椭球集中寻找最大的一个。这里利用形状参考集方法度量吸引域估计值的大小。参考集X⊂Rn是先给定的有界凸集,对于集合Rfigure_0018_0141,定义

figure_0018_0142

有两种典型的XR描述,分别为椭球集figure_0018_0143,R>0和多面体集figure_0018_0144。如果figure_0018_0145,那么figure_0018_0146越大,则包含XR的集合figure_0018_0147的容积越大。因此可以通过最大化figure_0018_0148的值得到最大椭球吸引域的估计值。


定理1.4给出了椭圆集figure_0018_0149位于吸引域内的条件。用上述多面体形式的参考集,从所有满足定理1.4条件的椭圆集中选择最大的不变椭圆集figure_0018_0150,以使所估计的系统吸引域的保守性最小的问题可以描述为下面的优化问题,即

figure_0018_0151

figure_0018_0152

figure_0019_0153

使用不同的参考集估算椭球域所用的方法是相似的,本节给出用多面体集作为参考集估算系统的吸引域计算方法。给定一多面体集XR , XR =co{v1,v2,?,vl}。对每一个S j∈N ,为了将获得尽可能大的吸引域估计问题转化为带有 LMI约束条件的最优问题,设figure_0019_0154。因为XR为多面体,所以上式中的约束条件figure_0019_0155等价于

figure_0019_0156

因为 P3+P3T >P1 ,可知 P3 是非奇异的。定义 Q3=P3−1 ,并设figure_0019_0157figure_0019_0158,用QT和Q分别左乘和右乘不等式(1.27)后,得

figure_0019_0159

进一步,令figure_0019_0160,利用Schur\'s补公式得到上式的LMI描述

figure_0019_0161

不等式(1.28)等价于

figure_0019_0162

通过上面的分析,优化式(1.33)可以转化为下面的凸优化问题,即

figure_0019_0163

figure_0019_0164

figure_0020_0165

figure_0020_0166

值得注意的是,当m=1时,如果求解优化问题式(1.37),得到Q2和Q1−Q3近似等于0。


优化式(1.37)处理的是K已知时,取得最大椭圆figure_0020_0167为收缩不变的问题。这可用于设计控制反馈K,使闭环系统有一个figure_0020_0168的不变集。定义Z=KQ1,则有

figure_0020_0169

figure_0020_0170

figure_0020_0171

figure_0020_0172

求解上述的优化问题,可以在计算最大不变椭圆集figure_0020_0173的同时计算出状态反馈阵K =ZQ1−1


2.受约束连续系统的吸引域估计及控制器设计


考虑连续系统

figure_0020_0174

其中,x∈Rn,u∈Rmfigure_0020_0175是标准饱和函数同系统式(1.23),给定状态反馈控制器

figure_0020_0176

根据整数集M的定义,闭环系统可以表示为

figure_0020_0177

对于figure_0021_0178,设figure_0021_0179为系统式( 1.41 )的状态轨迹,那么原点吸引域为figure_0021_0180。对应于离散系统吸引域估计问题的分析方法,对于所有的figure_0021_0181,则椭球集figure_0021_0182被称为“收缩不变集”。


定理1.5 考虑系统式(1.41)。对于每一个S j∈N , j=1, 2,…, 2m ,若存在对称矩阵P1>0 ,以及矩阵P2、P3figure_0021_0183,当x(t)∈figure_0021_0184时,均有

figure_0021_0185

figure_0021_0186

那么椭球集figure_0021_0187位于吸引域内,并且是严格不变的。


证明:类似于离散系统的证明方法,如果式(1.42)成立,可以得到

figure_0021_0188

那么连续闭环系统在包含原点的椭球figure_0021_0189内渐近稳定。


如果限定P1=P2且P3=figure_0021_0190figure_0021_0191是一个充分小正数,且令figure_0021_0192,可以得到下面的推论。


推论1.2考虑系统式(1.41),若存在充分小的正数figure_0021_0193,对称矩阵P>0和矩阵H∈R。m×nfigure_0021_0194时,有

figure_0021_0195

figure_0021_0196

那么椭球集figure_0021_0197位于吸引域内,且是严格不变的。


因为figure_0021_0198是一个充分小的正数,如果式(1.44)成立,那么有

figure_0021_0199

这个结果等价于文献[3]的定理1.5的条件。


利用矩阵 Q 及其转置对不等式(1.42)进行相似变换,可以将估计满足定理 1.5 条件的最大椭球问题转化为LMI优化问题,即

figure_0021_0200

figure_0022_0201

figure_0022_0202

figure_0022_0203

值得注意的是,当m = 1时,如果求解上述优化问题得到Q2 =Q3,那么上述优化问题中将等价于文献[3]。


同样在用优化式(1.47)中使用变量替换Z=KQ1,这个优化问题可用于设计控制反馈K,使闭环系统有一个figure_0022_0204的不变集,方法如下

figure_0022_0205

figure_0022_0206

figure_0022_0207

figure_0022_0208

求解上述的优化问题,可以在得到最大不变椭圆集的同时计算出状态反馈阵K =ZQ1−1


3.仿真研究


例1 考虑式(1.23)所描述的一个两输入系统,即

figure_0022_0209

运用LQR方法所设计的状态反馈阵为

figure_0023_0210

figure_0023_0211

当反馈矩阵 K 和参考集 XR已知时,运用本文提出的方法求解优化式(1.37)得到figure_0023_0212=2.5652。


例2 对于单输入系统,只能定义一个辅助反馈变量H。由推论1.2可知,本文提出的算法与定理 1.5 提出的方法得到的最大稳定域应该是相同的。下面给出一个单输入系统的仿真算例。


考虑一个单输入连续系统,即m=1:

figure_0023_0213

其中,

figure_0023_0214

取v1=[−1 0.8]T ,X R =co{v1 v2},运用定理1.5提出的方法求解优化问题式(1.47)得到figure_0023_0215=4.3713。


例3 前两个算例讨论了系统的稳定性分析问题,本算例将在设计控制器的同时得到系统吸引域的估计值。为了更形象直观地比较说明本章算法的有效性,考虑一个二输入系统:

figure_0023_0216

取多面体参考集figure_0023_0217,运用本文提出的控制器设计方法求解优化问题式(1.48)得到

figure_0023_0218

选择同样的参数v1,运用文献[3]提出的方法figure_0023_0219=0.99,状态反馈为

figure_0023_0220

图1-3给出了本算例椭圆稳定域估计曲线。


文献[3]运用饱和函数的LDI描述方法得到的系统最大稳定域估计值为图中虚线表示的椭圆;本章用饱和函数的多参数描述形式估计系统的稳定域得到的最大椭圆用实线表示,显然本方法得到的吸引域估计值相对要大。


在前面提出利用多面体参考集估计吸引域,如果采用椭球参考集,那么估计吸引域的优化问题式(1.38)和式(1.48)中的约束条件figure_0024_0221应为下面的形式:

figure_0024_0222

以优化问题式(1.48)为例,选择椭球参考集为:

figure_0024_0223

则求解优化问题得到的最大椭球和采用上述多面体参考集得到的椭球稳定域相同,如图1-3最外面的实线椭圆所示。

figure_0024_0224

图1-3 例3的稳定域

图 1-4给出了所获得的系统稳定域(实线椭圆)和闭环系统的状态轨迹(虚线曲线)。


由状态轨迹可知在椭圆域内系统是渐近稳定的,且估计得到的椭圆稳定域比较接近于系统的吸引域。


这 3 个仿真算例表明对于单输入系统,本文提出的算法所得到最大稳定域与文献[3]和[4]的结果非常近似。但是对于多输入系统,本文所得到的椭圆要明显的比文献[3]和[4]的结果大,并且比较接近于系统真实的稳定域。这表明本文提出的算法减少了吸引域估计的保守性。

figure_0025_0225

图1-4 稳定域和状态轨迹

用此等式计算得到的矩阵参数构造的输出反馈控制器,使得系统式(1.55)在位于原点吸引域内的椭球域figure_0032_0300内渐近稳定。


1.4.2 输出反馈控制系统稳定性分析


上一节讨论了含控制输入饱和约束的状态反馈系统吸引域的估计问题,即假定系统状态是全息的。在实际问题中,系统的状态有时不能直接测量,难以应用状态反馈控制系统。有时即使系统的状态可以直接测量,但考虑到实施成本和系统可靠性等因素,如果可以用系统的输出反馈来达到闭环系统的性能要求,则更适合于选择输出反馈的控制方式。因此输出反馈控制器的研究更具有实际意义,动态输出反馈控制器有全阶(对象阶)和降阶的两类。在系统稳定性的研究中,全阶输出反馈控制系统的分析设计问题已经得到了较早和较多的关注;在实际控制系统中,阶次较低的控制器所需元件少,反馈增益一般也较小,容易实现且不易出故障。因此如何设计阶次较低的控制器一直是控制理论应用的一个热点问题,受到广泛关注。


1.连续时间系统


考虑状态空间描述的被控系统

figure_0025_0226

其中, x(t)∈Rn是状态变量, u(t)∈Rm是控制输入,figure_0025_0227和(n2figure_0025_0228n)是系统测量输出。饱和函数figure_0025_0229的定义相同于系统式(1.23)。As、Bs和Cs是适当维数的常数矩阵。


当通过输出y(t)只有部分的系统状态信息可以利用时,不失一般性,假设Cs是行满秩矩阵,此时存在一个非奇异的矩阵 T,使figure_0026_0230,即系统有n2个状态变量可以直接利用。定义为

figure_0026_0231

利用矩阵T可将系统式(1.49)进行转换

figure_0026_0232

figure_0026_0233

其中,



对于系统式(1.50),非线性降阶(n1阶)输出反馈控制器被定义为

figure_0026_0235

定义figure_0026_0236,这样由式(1.50)和式(1.51)组成的闭环系统可以描述为

figure_0026_0237

其中,

figure_0026_0238

针对这个系统的设计目标是构造非线性降阶控制器式(1.51),使系统式(1.49)被包含在吸引域的椭球域内,即

figure_0026_0239

渐近稳定,并且获得其中最大的椭球域,即

figure_0026_0240

由于饱和具有非线性特性,故在应用Lyapunov函数方法估计系统的最大椭球吸引域前,需要处理非线性环节。到目前为止LDI处理技术是一种比较有效的方法。运用这种技术,通过引进一个受约束的辅助矩阵,即

figure_0026_0241

并定义

figure_0027_0242

将饱和控制figure_0027_0243替代为一个多胞型线性微分包含的形式。


引理1.2[5] 对于某个给定矩阵H,如果xcl∈L(H ) ,则[S]

figure_0027_0244

其中, co表示凸组合。


基于1.4.1节中的整数集M和N的定义,下面给出几个需要用到的符号定义。


定义1.9 给定正数i∈M 和S j∈N ,如果i属于S j ,那么

figure_0027_0245

定义1.10 给定S j∈N ,那么

figure_0027_0246

利用定义1.9、定义1.10和引理1.2,对于某个给定的H

....

本站仅展示书籍部分内容

如有任何咨询

请加微信10090337咨询

本站仅展示书籍部分内容
如有任何咨询

请加微信10090337咨询

再显示