毕业论文时滞项可微系统的时滞相关稳定性条件

1、本科毕业论文（设计）题 目：时滞项可微系统的时滞相关稳定性条件时滞项可微系统的时滞相关稳定性条件学 院： 自动化工程学院 专 业： 自动化 姓 名： xxx 指导教师： xxx 2013 年 6 月 1 日青岛大学本科生毕业论文（设计）1delay-dependent stability criteria for systems with differentiable time delays青岛大学本科生毕业论文（设计）2摘 要本文研究了带有可微时变时滞的连续系统的稳定性问题。通过使用时滞导数的信息，本文给出了时滞系统的改进的渐近稳定性。与以前的研究方法不同的是，本文考虑了时滞导数上界，即使这

2、种时滞导数的上界大于等于 1。可以证明取得的结果要比现有结论保守性更低。同时，因为涉及较少的决策变量，本文所展示稳定判据的计算复杂程度大大降低。用 matlab 证实了所得稳定条件的有效性和更低的保守性。关键字 时滞相关稳定条件 线性矩阵不等式（lmi） 时滞系统abstractthis paper studies the problem of stability for continuous-time systems with differentiable time-varying delays. by using the information of delay derivative, i

3、mproved asymptotic stability conditions for time-delay systems are presented. unlike the previous methods, the upper bound of the delay derivative is taken into consideration even if this upper bound is larger than or equal to 1. it is proved that the obtained results are less conservative than the

4、existing ones. meanwhile, the computational complexity of the presented stability criteria is reduced greatly since fewer decision variables are involved. and we use matlab illustrate the effectiveness and less conservatism of the obtained stability conditions.keywords delay-dependent stability cond

5、ition linear matrix inequality (lmi) time-delay systems青岛大学本科生毕业论文（设计）3目 录前 言.2第 1 章 绪 论.41.1 时滞系统的相关介绍.41.2 时滞系统稳定性的分析基本方法.41.3 时滞系统稳定性问题与展望.51.4 schur 补的相关知识补充.51.4.1 schur 补的定义.51.4.2 schur 引理.51.5 本文的主要研究工作.51.6 文中的符号说明.61.7 小 结.6第 2 章 lmi 工具箱介绍.72.1 线性矩阵不等式及相关术语.72.2 线性矩阵不等式的确定.92.3 线性矩阵不等式求解器.

6、16第 3 章 时滞项可微系统的时滞相关稳定性条件.213.1 主要结果.213.2 与现有结果的联系.273.3 数值算例 .343.4 结论.35结束语.36谢 辞.37参考文献.38附录 仿真程序.40青岛大学本科生毕业论文（设计）1前 言从系统理论的观点看，任何实际系统的过去状态不可避免地要对当前的状态产生影响，即系统的演化趋势不仅依赖于系统当前的状态，也依赖于过去某一时刻或若干时刻的状态，这类系统称为时滞系统。时滞产生的原因有很多，如：系统变量的测量(复杂的在线分析仪)、长管道进料或皮带传输、缓慢的化学反应过程等都会产生时滞。时滞常见于电路、光学、神经网络、生物环境及医学、建筑结构、

7、机械等领域，由于应用背景广泛，受到很多学者的关注。从理论分析的角度来看，在连续域中，时滞系统是一个无穷维的系统，特征方程是超越方程，有无穷多个特征根，而在离散域中，时滞系统的维数随时滞的增加按几何规律增长，这给系统的稳定性分析和控制器设计带来了很大的困难。因此，对于时滞系统的控制问题，无论在理论还是在工程实践方面都具有极大的挑战性。 常见的时滞系统包括奇异时滞微分系统、脉冲时滞微分系统、lurie 时滞系统、中立型时滞系统和随机时滞系统等。系统的稳定性和镇定问题是控制理论界的重要课题。若控制系统在任何足够小的初始偏差的作用下，其过渡过程（输出）随着时间的推移，逐渐衰减并趋于零，具有恢复平衡状态

8、的能力，则称该系统为稳定。镇定问题源于稳定性问题，当受控系统通过状态反馈(或者输出反馈),使的闭环系统渐近稳定,这样的问题称为镇定问题。早在 20 世纪 50 年代，就有很多学者开始研究时滞系统的镇定性问题和控制问题，其研究方法大致可分为频域方法和时域方法。在早期主要是频域方法，通过分析其特征方程根的分布以及 lyapunov 矩阵函数的解，给出时滞系统的稳定性判据和控制器设计的相应准则。频域方法在单输入输出的定常时滞系统方面已经取的了一些很好的结果。但是，对于多输入输出的时滞系统和时变时滞系统，用频域方法就很难得到结论。因此，相应的时域方法就得到发展，主要有 lyapunovkrasovsk

9、ii 泛函方法和 razumikhin 函数方法。这两种方法的基本思想都是构造 lyapunovkrasovskii 泛函或 lyapunov 函数，然后对其求导，使其导数小于零，得到稳定性判据和控制器设计基本准则。这是由krasovskii 和 razumikhin 所分别创造的，现已成为分析时滞系统镇定性和控制器设计的主要方法。尤其是在 20 世纪 90 年代，随着 riccati 方程和 matlab 中 lmi( 线性矩阵不等式)的发展，更使这两种方法得到了广泛的应用，这其中，有两类成果备受关注: 一类是时滞无关条件，一类是时滞相关条件。在 90 年代初及以前，用这两种方法所得到的条件

10、基本上都是时滞无关的，由于时滞无关条件不含时滞信息，对于小时滞系统，这类条件具有较强的保守性，于是时滞相关条件得到发展。目前对于系统的时滞相关问题的研究方法主要有: 离散 lyapunovkrasovskii 泛函方法、模型变换法、参数化模型方法、自由权矩阵法、积分不等式法。青岛大学本科生毕业论文（设计）2本文提出了消除时滞导数上界限制的新方法，同时给出了时滞系统的新的稳定条件。可以证明，新结果比现有结果具有更低的保守性。同时，得到的稳定判据有更少的决策变量，所以在数学上更简便，并且计算上更有效。在结果保守性不变的条件下，本文也给出了简化由加权矩阵和广义系统方法得到的时滞相关稳定条件的方法。青

11、岛大学本科生毕业论文（设计）3第 1 章 绪 论1.1 时滞系统的相关介绍从系统理论的观点看，任何实际系统的过去状态不可避免地要对当前的状态产生影响，即系统的演化趋势不仅依赖于系统当前的状态，也依赖于过去某一时刻或若干时刻的状态，这类系统称为时滞系统。时滞产生的原因有很多，如：系统变量的测量(复杂的在线分析仪)、长管道进料或皮带传输、缓慢的化学反应过程等都会产生时滞。时滞常见于电路、光学、神经网络、生物环境及医学、建筑结构、机械等领域，由于应用背景广泛，受到很多学者的关注。从理论分析的角度来看，在连续域中，时滞系统是一个无穷维的系统，特征方程是超越方程，有无穷多个特征根，而在离散域中，时滞系统

12、的维数随时滞的增加按几何规律增长，这给系统的稳定性分析和控制器设计带来了很大的困难。因此，对于时滞系统的控制问题，无论在理论还是在工程实践方面都具有极大的挑战性。常见的时滞系统包括奇异时滞微分系统、脉冲时滞微分系统、lurie时滞系统、中立型时滞系统和随机时滞系统等。1.2 时滞系统稳定性的分析基本方法纵观时滞系统的研究和发展，有两条主要研究途径，即时域方法和频域方法两大类。频域分析方法利用 smith 预估、变结构控制等方法设计控制器，并利用 nyquist 图等频域分析手段判断系统参数在一定范围摄动条件下闭环系统稳定性。这一类设计只适用于定常不确定系统或慢时变系统。时滞系统的时域分析方法越

13、来越成为时滞系统尤其是不确定时滞系统（包括系统矩阵的参数不确定性以及时滞本身的不确定性）稳定性分析以及控制器综合的主要方法。时域分析方法克服了频域分析不能处理时变和参数摄动的不足，而且具有方法简单、易于计算等优点，使其在实际工程应用中更加具有优势。近年来有关不确定时滞系统的结论基本上都是用时域的分析方法取得的。时域方法用得最多的是lyapunov 直接设计方法。利用 lyapunov 第二方法对时滞系统的研究主要是通过构造适当的 lyapunov 函数来求解时滞系统的无记忆反馈控制律，这是设计时变及不确定时滞系统鲁棒控制器的有效途径。基于 lyapunov 方法的无记忆反馈控制器不但设计简便，

14、在线计算量少，因而近年来受到很多学者重视。利用 lyapunov 方法对时滞系统的研究结果可以分为两大类：时滞无关结果和时滞相关结果。所谓时滞无关结果，是指所得结论都是独青岛大学本科生毕业论文（设计）4立于时滞大小的，即允许系统的滞后为无穷大，而对系统滞后的变化率一般都作了小于1 的假设。相反地，时滞依赖结果跟系统滞后的大小有关。1.3 时滞系统稳定性问题与展望1)目前有关时滞系统稳定性的分析结果很多，但是进行控制器设计时，只在个别情况下才会得到线性矩阵不等式(lmi)，多数情况下得到的是多项式矩阵不等式(pmi)或双线性矩阵不等式(bmi)。如何将多项式矩阵不等式转化为lmi，或者在无法转化

15、成lmi时，如何对其利用优化方法进行求解，是今后继续努力的方向目前发展起来的多项式优化理论有望为这一问题提供系统化方法。2)如何得到计算复杂性低，同时保守性较小的稳定性准则是未来的努力方向。其中lyapunovkrasovskii泛函的适当选取，尤其是参数依赖的lyapunov泛函的选取，将对结果的保守性产生积极影响，这方面还有大量的工作有待进行。3)基于线性矩阵不等式的稳定性准则在保守性方面难于比较，至少看起来不直观原因是线性矩阵不等式在矩阵维数、变量及变量个数方面有所不同。如何进一步寻求系统化方法进行相关分析，这方面的工作很有意义。4)近年有关时滞的讨论多数集中在线性系统，有关非线性时滞系

16、统的讨论则较少(当然也有例外)，而实际系统往往是非线性的，这也是进一步努力的方向之一。5)近年来对网络控制系统、无线通讯网络、无线传感器网络的研究蓬勃兴起，因网络中的信息必须通过通信网络分时传送，不可避免地在控制环路中引入了通讯延迟(时滞)，消除时滞对网络系统的稳定性影响是备受关注的问题，是推动时滞系统进一步研究发展的动力。1.4 schur 补的相关知识补充1.4.1 schur 补的定义 分块矩阵 m 表示为。 abcd如果 a 可逆，则 m 的 schur 补定义为； 1*dcab如果 d 可逆，则 m 的 schur 补定义为。 1*abdc1.4.2 schur 引理 分块矩阵 m

17、表示为。 abcd青岛大学本科生毕业论文（设计）5a. 如果 a 可逆，则 m0 等价为 a0 且 m 的 schur 补为正定； b. 如果 d 可逆，则 m0 等价位 d0 且 m 的 schur 补为正定。 1.5 本文的主要研究工作本文提出了消除时滞导数上界限制的新方法，同时给出了时滞系统的新的稳定条件。可以证明，新结果比现有结果具有更低的保守性。同时，得到的稳定判据有更少的决策变量，所以在数学上更简便，并且计算上更有效。在结果保守性不变的条件下 ，本文也给出了简化由加权矩阵和广义系统方法得到的时滞相关稳定条件的方法。1.6 文中的符号说明 维实向量空间；nrn 维实矩阵空间；n nr

18、nn 矩阵的转置；taa 矩阵的欧氏范数，即；aa1/2max()taa a 为正定对称矩阵；0p p 为半正定对称矩阵；0p p 具有适当维数的单位矩阵；i 矩阵中的对称部分。1.7 小 结在实际的工业生产过程和自然科学过程中，时滞现象的存在是不可避免的。特别是电力系统、机械传输系统、网络控制系统以及城市交通管理系统中，时滞现象的存在对系统造成的影响是不可忽略的。所以在分析系统的稳定性及跟踪控制问题时，考虑时间延迟对系统的影响是非常重要的。青岛大学本科生毕业论文（设计）6第 2 章 lmi 工具箱介绍线性矩阵不等式(lmi)工具箱是求解一般线性矩阵不等式问题的一个高性能软件包。由于其面向结构

19、的线性矩阵不等式表示方式，使的各种线性矩阵不等式能够以自然块矩阵的形式加以描述。一个线性矩阵不等式问题一旦确定，就可以通过调用适当的线性矩阵不等式求解器来对这个问题进行数值求解。lmi 工具箱提供了确定、处理和数值求解线性矩阵不等式的一些工具，它们主要用于： 以自然块矩阵形式来直接描述线性矩阵不等式； 获取关于现有的线性矩阵不等式系统的信息； 修改现有的线性矩阵不等式系统； 求解三个一般的线性矩阵不等式问题； 验证结果。本章将详细介绍 lmi 工具箱提供的用于解决以上各个问题的相关函数和命令。2.1 线性矩阵不等式及相关术语一个线性矩阵不等式就是具有以下一般形式的一个矩阵不等式： (2-1)0

20、)(110nnlxlxlxl其中：，是给定的对称常数矩阵，是未知变量，称为决策变量，0l1lnlnxx,1是由决策变量构成的向量，称为决策向量。ntnrxxx,1尽管表达式(1)是线性矩阵不等式的一个一般形式，但在大多数实际应用中，线性矩阵不等式常常不是以一般表示式(1)的形式出现，而是具有以下形式：),(),(11nnxxrxxl其中的和是矩阵变量的仿射函数，通过适当的代数运算，上式可以写)(l)(rnxx,1成线性矩阵不等式的一般表示式(1)的形式。例如，在系统稳定性问题中经常遇到的lyapunov 矩阵不等式 (2-2)0 xaxat青岛大学本科生毕业论文（设计）7也是一个线性矩阵不等式

21、，其中的是一个矩阵变量。我们以一个二阶矩阵x为例，将矩阵不等式(2)写成一般表示式(1)的形式。针对二阶矩阵不等式(2)，2021a对应的矩阵变量是一个二阶的对称矩阵，不等式(2)中的决策变量是矩x3221xxxxx阵中的独立元。根据对策性，矩阵变量可以写成x321,xxxx100001100001321xxxx将矩阵和上式代入矩阵不等式(2)，经整理，可得a (2-3)0400043300222321xxx这样就将矩阵不等式(2)写成了线性矩阵不等式的表示式(1)。显然，与 lyapunov 矩阵不等式(2)相比，表示式(3)缺少了许多控制中的直观意义。另外，(3)式涉及到的矩阵也比(2)式

22、中的多。如果矩阵是 n 阶的，则(3)式中的系数矩阵一般有 n(n+1)/2 个。因此，这a样的表达式在计算机中将占用更多的存储空间。由于这样的一些原因，lmi 工具箱中的函数采用线性矩阵不等式的结构表示。例如，lyapunov 矩阵不等式(2)就以矩阵变量的x不等式来表示，而不是用其一般形式(3)来表示。一般的，一个线性矩阵不等式具有块矩阵的形式，其中每一个块都是矩阵变量的仿射函数。以下通过一个例子来说明有关描述一个线性矩阵不等式的术语。考虑控制中的一个线性矩阵不等式：h0nidbdicxbxcxaxanttttt其中：、是给定的矩阵，和是问题的变量。abcdnnmtrxxr 称为外因子，块

23、矩阵n idbdicxbxcxaxaxltttt),(称为内因子。外因子可以不是一个正方矩阵，它在许多问题中常常不出现。 和是问题的矩阵变量。注意标量也可以看成是一个维的矩阵。x11青岛大学本科生毕业论文（设计）8 内因子是一个对称块矩阵。根据对称性，可以由对角线及其上方的),(xl),(xl块矩阵完全确定。 中的每一块都是矩阵变量和的仿射函数。这一函数由常数项和变量项),(xlx这两类基本项组成，其中常数项就是常数矩阵或以一些常数矩阵组成的算术表达式，例如中的 b 和 d；变量项是包含一个矩阵变量的项，例如等。),(xlixa,一个线性矩阵不等式不论多么复杂，都可以通过描述其中每一块的各项内

24、容来确定这个线性矩阵不等式。 2.2 线性矩阵不等式的确定lmi 工具可以处理具有以下一般形式的线性矩阵不等式：mxxrmnxxlnktkt),(),(11其中：是具有一定结构的矩阵变量，左、右外因子和是具有相同维数的给kxx,1nm定矩阵，左、右内因子和是具有相同块结构的对称块矩阵。)(l)(r注意在线性矩阵不等式的描述中，左边总是指不等式较小的一边，例如对线性矩阵不等式0，称为是不等式的右边，0 称为是不等式的左边，常表示成 0。xxx要确定一个线性矩阵不等式系统，需要做以下两步：1. 给出每个矩阵变量的维数和结构；kxx,12. 描述每一个线性矩阵不等式中各个项的内容。这个过程产生所描述

25、线性矩阵不等式系统的一个内部表示，它以一个单一向量的形式储存在计算机内，通常用一个名字，例如 lmisys 来表示。该内部表示 lmisys 可以在后面处理这个线性矩阵不等式时调用。下面将通过 lmi 工具箱中的一个例子来说明线性矩阵不等式系统的确定。运行lmidem 可以看到这个例子的完整描述。例例 1：考虑一个具有 4 个输入、4 个输出和 6 个状态的稳定传递函数 (2-4)basicsg1)()(和一组具有以下块对角结构的输入/输出尺度矩阵：d (2-5)5432110000000000ddddddd青岛大学本科生毕业论文（设计）9则在具有时变不确定性系统的鲁棒稳定性分析中提出了以下问

26、题：寻找一个具有结构(5)的尺度矩阵，使的。d1)(sup1djdg可以证明：这样一个问题可以转化成一个线性矩阵不等式系统的可行性问题，既寻找两个对称矩阵和，使的66rx44rddst (2-6)0sxbxbsccxaxattt (2-7)0x (2-8)is 用命令 lmivar 和 lmiterm 给出线性矩阵不等式系统(6)(8)的内部描述如下：setlmis()x=lmivar(1,6 1)s=lmivar(1,2 0;2 1)% 1st lmi lmiterm(1 1 1 x,1,a,s)lmiterm(1 1 1 s,c,c)lmiterm(1 1 2 x,1,b)lmiterm(

27、1 2 2 s,-1,1) % 2nd lmi lmiterm(-2 1 1 x,1,1)% 3rd lmilmiterm(-3 1 1 s,1,1)lmiterm(3 1 1 0,1)lmisys=getlmis 其中：函数 lmirar 定义了两个矩阵变量和，lmiterm 则描述了每一个线性矩阵不等式xs中各项的内容。getlmis 回到了这个线性矩阵不等式系统的内部表示 lmisys，lmisys 也称为是储存在机器内部的线性矩阵不等式系统的名称。以下将详细介绍这几个函数的功能和用法。青岛大学本科生毕业论文（设计）10setlmis 和和 getlmis一个线性矩阵不等式系统的描述以

28、setlmis 开始，以 getlmis 结束。当要确定一个新的系统时，输入：setlmis()如果需要将一个线性矩阵不等式添加到一个名为 lmiso 的现有的线性矩阵不等式系统中，则输入：setlmis(lmiso)当线性矩阵不等式系统被完全确定好后，输入：lmisys=getlmis该命令返回这个线性矩阵不等式系统的内部表示 lmisys。lmivar函数 lmivar 用来描述出现在线性矩阵不等式系统中的矩阵变量，每一次只能描述一个矩阵变量。矩阵变量的描述包括该矩阵变量的结构。该函数的一般表达式是：x=lmivar(type,struct)这一函数定义了一个新的矩阵变量，是该矩阵变量的变

29、量名。函数中的第一个xx输入量 type 确定了矩阵变量的类型，第二个输入量 struct 进一步根据变量的类型给xx出该变量的结构。变量的类型分成三类：type=1：对称块对角结构。这种结构对应于具有以下形式的矩阵变量：rddd00000021其中对角线上的每一个矩阵块是方阵，它可以是零矩阵、对称矩阵或者数量矩阵。这jd种结构也包含了通常意义的对称矩阵和数量矩阵（分别相当于只有一块） 。此时，struct 是一个维的矩阵。如果该矩阵的第 i 行是(m，n)，则其中的 m 表示对称矩阵块的阶2rid数，而 n 只能取 1、0 或者-1，其中 n=1 表示是一个满的对称矩阵（或者无结构的对称jd

30、矩阵） ，n=0 表示是一个数量矩阵，n=-1 表示是一个零矩阵。ididtype=2：长方形结构。这种结构对应于任意的长方矩阵。此时，srtuct=(m，n)表示矩阵的维数。青岛大学本科生毕业论文（设计）11type=3：其他结构。这种结构用来描述更加复杂的矩阵，也可以用于描述矩阵变量之间的一些关联。的每一个元或者是 0，或者是，其中是第 n 个据侧变量。相应xnxnx的，struct 是一个和变量有相同维数的矩阵，其中的每一个元取值如下：xnnxjixnxjixnjixjistruct),(,),(,0),(, 0),(如果如果如果例例 2：考虑具有三个矩阵变量、和的线性矩阵不等式系统，其

31、中1x2x3x 是一个 3 3 维的对称矩阵；1x 是一个 2 4 维 的长方矩阵；2x 其中是 5 5 维的对称矩阵，和是两个标量，2213000000ixa12 表示 2 2 维的单位矩阵。2i可以应用 lmivar 来定义这些矩阵变量。setlmis()x1=lmivar(1,3 1)x2=lmivar(2,2 4)x3=lmivar(1,5 1;1 0;2 0)lmiterm在确定了矩阵变量之后，还需要确定每一个线性矩阵不等式中各项的内容。线性矩阵不等式的项指构成这个线性矩阵不等式的块矩阵中的加项。这些项可以分成三类：1. 常数项；2. 变量项，即包含了矩阵变量的项，例如(3)式中的和

32、。一般的变量项具xatscct有形式，其中的是一个变量，和是给定的矩阵，分别称为该变量项的左系数pxqxpq和右系数；3. 外因子。青岛大学本科生毕业论文（设计）12在描述一个具有多个块的线性矩阵不等式时，lmi 工具箱提供了这样的功能，即只需要确定对角线上和对角线上方的项的内容，或者只描述对角线上和对角线下方的项的内容，其他部分项的内容可以根据线性矩阵不等式的对称性得到。用命令 lmiterm 每次可以确定线性矩阵不等式的一个项的内容。例如，对称性矩阵不等式0sxbxbsccxaxattt可以用一下一组命令来描述： lmiterm(1 1 1 x,1,a,s) lmiterm(1 1 1 s

33、,c,c) lmiterm(1 1 2 x,1,b) lmiterm(1 2 2 s,-1,1)这些命令一次描述了项、和。在每一条命令中，第 1 项是xaxatscctxbs一个四元向量，它刻画了所描述的项所在的位置和特征； 第 1 个元表示所描述的项属于哪一个线性矩阵不等式。值 m 表示第 m 个 不等式的左边，-m 表示第 m 个不等式的右边。 第 2 和 3 个元表示所描述的项所在块的位置。例如，向量1 1 2 1表示所 描述的项位于第一个线性矩阵不等式左边内因子的块(1,2)中。第 2 和第 3 个元均取零表示所描述的项在外因子中。 最后一个元表明了所描述的项是常数项还是变量项。如果是

34、变量项，则 进一步说明涉及哪一个变量。0 表示常数项，k 表示所描述的项包含第 k 个矩阵变量，-k 则表示包含矩阵变量的转置（在例 1 中，kxkxtkxx 是第 1 个变量，s 是第 2 个变量，它们按确定的先后顺序排列） 。lmiterm 的第 2 项和第 3 项包含了数据（常数项的值，外因子，变量项或者pxq中的左、右系数） 。第 4 项是可选择的，且只能是s。qpxt在描述项的内容里，有一些简化的方法。1.零块可以省略描述；2. 可以通过在命令 lmiterm 中外加一个分量s，使的可以只用一个命令 lmiterm 就能描述一个变量项与该变量项的转置的和。例如，上面的第一个命令描述了

35、。xaxat3. 可以用一个标量值来表示一个数量矩阵，即用表示数量矩阵，其中是一个i标量。如例 1 中的第 3 个不等式被描述成is 青岛大学本科生毕业论文（设计）13 lmiterm(-3 1 1 s,1,1)lmiterm(3 1 1 0,1)为了便于阅读，也可以用线性矩阵不等式和矩阵变量的名称来表示对应的线性矩阵不等式和矩阵变量。矩阵变量的变量名可以用命令 lmivar 来赋值，线性矩阵不等式的名称则可以用函数 newlmi 来确定。这些标识符可以用在命令 lmiterm 中以表示相应的线性矩阵不等式或者矩阵变量。对例 1 中的线性矩阵不等式系统，采用名称的相应描述如下：setlmis(

36、)x=lmivar(1,6 1)s=lmivar(1,2 0;2 1)brl=newlmilmiterm(brl 1 1 x,1,a,s)lmiterm(brl 1 1 s,c,c)lmiterm(brl 1 2 x,1,b)lmiterm(brl 2 2 s,-1,1)xpos=newlmilmiterm(-xpos 1 1 x,1,1) slmi=newlmilmiterm(-slmi 1 1 s,1,1)lmiterm(slmi 1 1 0,1)lmisys=getlmis其中：x 和 s 分别表示变量和，而 brl、xpos 和 slmi 则分别表示第 1、第 2 和第 3xs个线性矩

37、阵不等式。-xpos 指的是第 2 个线性矩阵不等式的右边，-x 表示变量的转置。xlmiedit青岛大学本科生毕业论文（设计）14线性矩阵不等式编辑器 lmiedit 是一个图形用户界面，它可以按符号方式直接确定线性矩阵不等式系统。输入lmiedit出现一个具有一些可编辑文本区域和各种按钮的窗户。按以下步骤来确定一个线性矩阵不等式系统：1. 在文本区域的上半部分给出每一个矩阵变量的描述（名字和结构） ，其机构是通过类型（s 表示对称块矩阵，r 表示无结构的长方矩阵，g 表示其他机构矩阵）和一个“附加”的结构矩阵（类似于 lmivar 中的 struct）来刻画的。在文本编辑区，使用一行描述一

38、个变量。2. 在文本区的下半部分，按 matlab 的表示方式给出要描述的线性矩阵不等式。例如，线性矩阵不等式0ixbxbxaxatt可以通过输入01*;*xbbxaxxa来描述。其中 x 是文本区上半部分描述矩阵变量的变量名。一个线性矩阵不等式的描x述可能需要几行，但一行中最多只能描述一个线性矩阵不等式。完成了线性矩阵不等式系统的描述后，可以通过按相应的按钮来完成以下的任务： 显示用于描述线性矩阵不等式的 lmivar/lmiterm 命令串（按钮 view commands） ；反之，通过单击按钮 describe.可以将用一串 lmivar/lmiterm 命令定义的线性矩阵不等式系统按

39、 matlab 表示式显示。 将线性矩阵不等式的符号描述存为一个 matlab 语句串（按钮 save） 。以后可以通过按钮 load 重新显示这种描述。 可以从一个文件读一串 lmivar/lmiterm 命令（按钮 read） ，然后通过单击 “describe the matrix variables”或者“describe the lmis.”显示出由这些命令确定的线性矩阵不等式系统的符号表示。 写一串用于描述一个特殊线性矩阵不等式系统的 lmivar/lmiterm 命令（按钮 write） 。 通过按钮 creat 产生线性矩阵不等式系统的内部表示，结果用一个线性矩阵不等式命名的

40、matlab 变量记录（如果线性矩阵不等式系统名字是 mylmi，则其内部表示用matlab 变量 mylmi 记录） 。内部表示 mylmi 可以被线性矩阵不等式求解器或者任何其他的线性矩阵不等式函数调用。如同命令 lmiterm 一样，可以应用简捷的方法来输入线性矩阵不等式的表示式。例如零块可以简单地输入 0，而不必定义其维数，类似地，单位矩阵只需输入字符 1 等。青岛大学本科生毕业论文（设计）15lmiedit 尽管很一般，但是它没有 lmiterm 灵活。以下是 lmiedit 的一些局限性： 在矩阵变量的两边不能使用括号。例如当 x 是一个变量名时， (a*x+b)*c+c*(a*c

41、+b)是不允许的，而(a+b)*x+x*(a+b)则是可以的。 不允许出现循环和条件语句。 当把 lmiterm 命令转换成一个线性矩阵不等式的符号描述时，如果 lmiterm 的第 1个分量不能确认就将出错。使用由 newlmi 和 lmivar 提供的线性矩阵不等式和变量标识符可以避免这样的问题。图 2.1 给出了用 lmiedit 描述例 1 中的线性矩阵不等式系统的窗口。图 2.1 lmiedit 的图形界面2.3 线性矩阵不等式求解器lmi 工具箱提供了用于求解以下三个问题的线性矩阵不等式求解器(其中表示决策x变量向量，即矩阵变量中的独立变元构成的向量)。k1,xx 可行性问题：寻找

42、一个 xrnn(或者等价的：具有给定结构的矩阵)，使满足线性矩阵的不等k1,xx 式系统青岛大学本科生毕业论文（设计）16)()(xbxa相应的求解器是 feasp。 具有线性矩阵不等式约束的一个线性目标函数的最小化问题：xctxmins.t. )()(xbxa相应的求解器是 mincx。 广义特征值的最小化问题：xmins.t.)()(xdxc)(0 xb )()(xbxa相应的求解器是 gevp。以下详细介绍 feasp 求解器的功能和使用方法。feasp求解器 feasp 的一般表达式如下：tmin,xfeas=feasp(lmisys,options,target)求解器 feasp

43、是通过求解如下的一个辅助凸优化问题mints.t.tixbxa)()(来求解该线性矩阵不等式系统 lmisys 的可行性问题。这个凸优化问题的全局最优值用 tmin 表示，作为求解器 feasp 输出的第一个分量。如果 tmin0，则系统 lmisys 是可行的。当系统 lmisys 为可行时，求解器 feasp 输出的第二个分量 xfeas 给出了线性矩阵不等式系统决策变量的一个可行解。进而，应用 dec2mat 可以得到系统 lmisys 矩阵变量的一个可行解。求解器 feasp 的输入变量 target 为 min 设置了目标值，使的只要 tmin0 表示限制决策变量在球体212rxni

44、i中，或者说向量 xfeas 的欧式范数不超过。该参数的默认值是。r910r可行域半径的设定可以避免产生具有很大数值的解 x，同时也可以加快计算过程，改进数值稳定性。options(4)：该参数用于加快迭代过程的结束，它提供了反映优化过程中迭代速度和解的精度之间的一个折中指标。当该参数取值为一个正整数时，表示在最后的次迭代jj中，如果每次迭代后 的减小幅度不超过 1%，则优化迭代过程就停止。该参数的默认值t是 10。options(5)：options(5)=1 表示不显示迭代过程中的数据，options(5)=0(默认值)则相反。将 options(i)设置为零相当于将相应的控制参数设置为默

45、认值，也可以通过忽略该输入变量来接受默认值。例例 3：求满足的矩阵，使的ip p (2-9)011 papat (2-10)022 papat (2-11)033 papat其中：,31211a7 . 23 . 15 . 18 . 02a0 . 27 . 09 . 04 . 13a为了调用 feasp，我们首先确定线性矩阵不等式系统：setlmis()p=lmivar(1,2 1)lmiterm(1 1 1 p,1,a1,s) %lmi #1lmiterm(2 1 1 p,1,a2,s) %lmi#2lmiterm(3 1 1 p,1,a3,s) %lmi#3lmiterm(-4 1 1p,1

46、,1) %lmi#4:plmiterm(4 1 1 0,1) %lmi#4:i青岛大学本科生毕业论文（设计）18lmis=getlmis然后调用 feasp 来求该现行矩阵不等式系统的一个可行决策变量：tmin,xfeas=feasp(lmis)得到 tmin=-3.1363。因此，线性矩阵不等式系统 lmis 是可行的。应用 dec2matpp=dec2mat(lmis,xfeas,p)得到问题的可行矩阵变量值：1 .1554 .1264 .1268 .270p在求解这个可行性问题的过程中，也可以附加一些约束，例如，要求矩阵的pfrobenius 范数不超过 10，且 tmin-1。也可以通

47、过调用tmin,xfeas=feasp(lmis,0,0,10,0,0,-1)来达到这些附加要求。相应的结果是 tmin=-1.1745，相应的矩阵 p 的最大特征值是。6912. 9)(maxp 如何从决策变量到矩阵变量以及从矩阵变量到决策变量如何从决策变量到矩阵变量以及从矩阵变量到决策变量当现行矩阵不等式由相应的矩阵变量描述时，线性矩阵不等式求解器涉及的是由这些矩阵变量中的独立元所组成的决策向量 x。两个函数 mat2dec 和 dec2mat 可以实现这两种变量之间的转换。考虑一个具有三个矩阵变量、的线性矩阵不等式系统。给定这些变量的1x2x3x特定值 x1、x2、x3，那么由 mat2

48、dec 可以得到相应的决策向量的值：xdec=mat2dec(lmisys,x1,x2,x3)如果 lmisys 后分量的个数和线性矩阵不等式系统 lmisys 中的矩阵变量个数不符，则系统会提示一个出错信息。这个函数在线性矩阵不等式求解器 mincx 或 gevp 的初始化中也是很有用的。例如，给定、的一个初始猜测值，mat2dec 就形成了相应决策向量的初始值 xinit。1x2x3x反之，给定决策向量的一个值 xdec，那么可以通过函数 dec2mat 给出相应的第 k 个矩阵的取值。例如，一下的表示式可以给出第 2 个矩阵变量的取值：、x2=dec2mat(lmisys,xdec,2)

49、函数 dec2mat 中的最后一个分量表明了要求的是第 2 个矩阵变量，这里也可以用lmivar 定义的相应矩阵变量的变量名。矩阵变量和决策变量的总数分别由 matnbr 和 decnbr 给出。另外，函数 decinfo 提供了决策变量和矩阵变量之间关系的一些详细信息。青岛大学本科生毕业论文（设计）19青岛大学本科生毕业论文（设计）20第 3 章 时滞项可微系统的时滞相关稳定性条件3.1 主要结果在这一节，分析时变时滞连续系统的稳定性，并利用时滞导数相关李雅普诺夫函数得到了一个充分条件。考虑下面线性系统 （3-1）( )( )( ),0dx tax ta x td tt （3-2）( )(

50、),0 x ttt 其中，( )nx tr是状态变量，daa和 是适当维数的常量矩阵，时滞项是时变连续函( )d t数并且满足 （3-3）( )d t和 ( )d t（3-4）其中, (0), 和是常数。初始条件( ) t(,0t )是连续的向量值函数。在之前的文章中，例如3和6，时滞导数的上界应该小于 1。虽然7-8中的结果可以应用到1的情况，其稳定条件与时滞导数上界无关。对于(3-1)( 3-4)所描述的时滞系统，式 （3-5）( )( )( )ttt d txs qx s ds（其中， ）常被作为李雅普诺夫函数（例如6-7,11） 。但是，如果1，0tqq则这一项就是冗余的，因为( )(

51、 )( )ttt d txs qx s ds ( )( )( )( )( )(1( )( )( )( )( )(1)( )( )ttt d tttttxs qx s dsxt qx td txtd t qx td txt qx txtd t qx td t其中(1)0。青岛大学本科生毕业论文（设计）21这说明时滞项的导数没有考虑进去，这显然是不合理的。( )d t实际上，时滞导数大于等于 1 的情况是很普遍的。例如，在网络化控制系统里，时滞项表示kti,其中(1,2,)ki k 是采样时刻。所以，这一类时滞几乎在0t 时处处满( )d t足( )1d t 。对于1的情况下，如果选择了一个正数0

52、1满足1，则有 （3-6）( )( )1d td t 且 （3-7） ( )( )( )( )( )(1)( )( )ttttt d txs qx s dsxt qx txtd t qx td t 所以此时，项的导数考虑了进去。( )d t在此事实基础上，可获得一下定理。 定理定理 1 对于给定标量(0),01，和（满足1） ，如果存在矩阵tpp，使0,0(1,2,3,4),0(1,2,3)ttiijjqqizj和z （3-8）112311132345600()*()*0000*00*0*ttdza uzzza uu 那么，(3-1)( 3-4)所描述的系统是渐近稳定的。其中，青岛大学本科生毕

53、业论文（设计）224111131121111133123141215211643331123()(1)()()()11tiidiipaa pqzzpazqzzzzqzqzqzzzzuzzz 证明证明 构造一个李雅普诺夫函数 （3-1234( )( )0120()( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )ttttttttttttt d ttd ttttttttttv xxt px txs q s dsxs q s dsxs q s dsxs q s dsxs z x s dsdxs z x s dsdxs z x s dsd 9）其中，是需要

54、确定的矩阵。0,0(1,2,3,4)0(1,2,3)ijpqizj与根据莱布尼茨-牛顿公式，以下等式对于任何维数合适的矩阵和,iiin s m y(1,2,it i 都是成立的： ,5) （3-10）( )2( ) ( )( )( )0ttt d tt n x tx td tx s ds （3-11）( )2( ) ( )()( )0t d tttt s x td tx tx s ds （3-12）( )2( ) ()( )( )0ttt d tt m x tx td tx s ds （3-13）( )2( ) ( )( )( )0tttd tt y x tx td tx s ds （3-14

55、）( )( )2( ) ( )( )( )0td ttt d tt t x td tx td tx s ds其中青岛大学本科生毕业论文（设计）23125125125 , , ,ttttttttttttnnnnssssmmmm125125 ,( )( )( )()()( )ttttttttttttttyyyytttttxt xtd txtxtxtd t且或者，以下等式正确： （3-15）( )111( )( )( )( )( )( )( )ttt d tttttt d ttxs z x s dsxs z x s dsxs z x s ds （3-16）( )221( )( )( )( )( )(

56、 )( )ttt d tttttt d ttxs z x s dsxs z x s dsxs z x s ds （3-17）33( )( )( )33( )( )( )( )( )( )( )( )( )ttttttd ttd tt d tttt d ttxs z x s dsxs z x s dsxs z x s dsxs z x s ds 对，取沿着（3-1）的轨迹的时间导数，可得0t ( )tv x4211341231( )2( )( )()()( )( )()()(1( )( )( )(1( )( )( )( )() ( )( )( )ttttiitttttttv xxt px txt

57、q x txt q x txtq x td txtd t q x td td txtd t q x td txtzzzx txs z x s ds23421134123( )( )( )( )2( )( )()()( )( )()()(1)( )( )(1)( )( )( )()tttttttttiittttxs z x s dsxs z x s dsxt px txtq x txt q x txtq x txtd t q x td txtd t q x td txtzzzx( )11232( )( )( )33( )( )( )( )( )( )( )() ( )( )( )( )( )(

58、)( )2( ) ( )( )( )2tt d tttttt d ttt d tttd ttttd tt d tttt d ttxs z x s dsxs zzzx s dsxs z x s dsxs z x s dsxs z x s dst n x tx td tx s ds( )( )( ) ( )()( )2( ) ()( )( )t d tttttt d tt s x td tx tx s dst m x tx td tx s ds青岛大学本科生毕业论文（设计）24( )( )( )1111122123132( ) ( )( )( )2( ) ( )( )( )( ) ( )tttd

59、ttd ttt d tttttttttt y x tx td tx s dst t x td tx td tx s dstnz nsz smz myz ytz ya uat (3-18)其中1131124000*(1)000=*00*0*dpaqqqq 41112()tiidqpapanysnmtmstyaaa 000通过 schur 补可得，不等式1-11112212330tttttttnznsz smz myz ytz ta ua 与下式等价 (3-19)12234+=0*t其中1111122222333333444445555541233000,ttdnsmyta unsmyta uns

60、mytnsmytnsmytdiagzzzzzu 所以，如果成立， ，则对于所有, 有成立。0 0t ( )0tv x记 ， （3-7118*t 20）青岛大学本科生毕业论文（设计）25其中2132301100001111001=000001000001100000000000000000000000000000iiiiiiiiiiiiiiii 711113172732222373332334744445555353=smyztmzytznsmzytnmytnsmyztz 7111722173274481233,nznzszszdiagzzzzz 且是非奇异的，所以当时有。10 0 相反地，如果