扩展结构大系统的鲁棒控制研究Word文档下载推荐.doc
- 文档编号:827170
- 上传时间:2023-04-29
- 格式:DOC
- 页数:37
- 大小:892.56KB
扩展结构大系统的鲁棒控制研究Word文档下载推荐.doc
《扩展结构大系统的鲁棒控制研究Word文档下载推荐.doc》由会员分享,可在线阅读,更多相关《扩展结构大系统的鲁棒控制研究Word文档下载推荐.doc(37页珍藏版)》请在冰点文库上搜索。
Abstract II
1绪论 1
1.1引言 1
1.2选择课题的目的及意义 4
2线性矩阵不等式(LMI) 5
2.1线性矩阵不等式一般形式 5
2.2LMI工具箱简介 7
3LMI算法应用 11
3.1基于LMI算法的有机结构控制 11
3.2LMI编程中的拆分矩阵法 15
4扩展结构系统的建模与设计 20
4.1一类扩展结构大系统的数学描述 20
4.2一类扩展结构大系统的鲁棒分散关联镇定 22
5鲁棒关联镇定算例及仿真 26
结论 29
致谢 30
参考文献 31
附录程序清单 33
`
1绪论
近年来随着科学技术的发展,大系统或称关联系统、组合系统越来越多地出现在各种应用工程中,如电力网系统,交通系统,经济系统,生物工程系统,社会系统,生态系统等。
本文的研究基于LMI鲁棒分散控制在多区域互联大系统中的应用。
1.1引言
众所周知,随着现代科技的迅速发展,大型工程技术的需求,人们提出了大系统的模型,它最初面临的问题是克服与其相应的数学模型维数日益增大和复杂带来的困难。
人们最初对大系统概念的认识,是把一个系统分解成相互联结的子系统,若能由子系统的性质组合得到整个系统的性质时,就把这个系统视为大系统。
然而什么叫大系统,至今还没有一个公认的严格定义。
但是大系统具有一些比较明显的特性:
规模庞大、结构复杂、功能综合、因素众多等。
对于大系统,集中控制将会使得整个控制系统信息变换异常复杂,通讯费用也将十分昂贵,因而使得集中控制大系统变得非常不切实际。
正因为如此,分散控制便于上世纪70年代应运而生,它把一个大系统分解为若干个子系统,对某一子系统的输入输出回路来讲,将其他子系统与该子系统的互联项视为可测干扰对此回路的影响,这样,在充分考虑各个子系统之间互联部分的影响后,对每个子系统单独设计控制器。
该控制器只能使用来自子系统自身的状态和信息,而不使用来自其他子系统的状态和信息。
这样做是由于大系统维数高,要分析和控制这样的系统所需的计算量将随之剧增,因此将它分解为各个子系统进行控制有利于计算的简化;
大系统分布面广,使系统内信息传递受到地理条件的限制,这种在信息上的结构约束也使得集中控制难以实现;
分散控制对于互联部分的各种结构化和非结构化的扰动本身具有一定的鲁棒性;
分散控制不需在各个子系统间传递信息,使得子系统具有相对的自治性,所以分散控制策略在处理大系统的控制问题时是非常有效的,它不但解决了集中控制难以解决的问题,还可以提高系统的可靠性和容错性,减少系统的投资。
对于互联系统来讲,将一个高维大系统分解为若干个维数较低的子系统是分散控制的先决条件。
如果按照实际子系统的物理边界来划分,那么这种分解可以得到子系统之间相互作用的重要信息,但同时也带来了计算上的低效率。
为了计算方便,人们已经研究了许多不受子系统物理边界约束的分解方法。
如嵌套的e分解,系统的包含原理及平衡BBD分解等。
近年来,包含原理引起了很多大系统领域的研究者的关注。
它已成功地应用到大系统的模型降阶、分散自动发电控制、互联电力系统的特殊重叠结构分解中以及多区域重叠互联电力系统的重叠分散控制中。
系统的结构分解解决了复杂系统的分解问题,这只是问题的一个方面。
更重要的是对其实施控制,使系统能够稳定、可靠的运行。
因此控制的稳定性就成了首要问题。
我们知道分散控制对子系统互联部分的各种扰动来说,本身就具有鲁棒性,也就是说分散控制的子系统对其局部的输入输出来说是稳定的。
但是这种局部稳定的子系统并不能保证总的互联系统稳定性,因此当沿着子系统边界运行时,系统可能会出现不稳定。
这主要是由于互联部分的不确定性造成的。
这就提到了分散控制问题的另一个性能指标,系统的鲁棒性。
所谓的鲁棒性是指系统具有较大的稳定裕度和很好抗干扰性能。
我们在设计控制率时要综合考虑系统的稳定性和鲁棒性等各指标的要求。
近年来,鲁棒控制方法已引起了学者们的广泛关注。
针对系统参数的不确定性,文[1]同时考虑了小参数不确定性,给出了一种适合于所有可允许的参数不确定的控制器设计方法,即把鲁棒控制方法和自适应控制方法结合起来。
鲁棒控制方法用来处理小参数的不确定性,而自适应控制方法则用来处理大参数的不确定性。
文[2,3]给出了负载频率控制(LoadFreqencycontrol,即LFC)的鲁棒设计方法,它们分别采用了奇异值分解SVD和Riccati稳定理论结合的技术,匹配条件和Lyapunov稳定性理论结合的技术,并对两区域互联电力系统进行了仿真研究。
文[4]把鲁棒控制问题转换成最优控制问题。
如果满足一定的匹配条件,最优控制问题的解就是鲁棒控制问题的解。
文[5]基于Q-参数理论提出了可用于电力系统自动发电控制的鲁棒控制器设计方法。
在Q-参数化方法中,电力系统的所有镇定控制器都带有Q-参数,通过选择这个自由参数来满足鲁棒稳定性和其它设计要求。
针对系统中不确定的结构扰动,Siljak提出了互联大系统关联稳定的概念,并对该问题进行了深入的研究,且又在此基础上提出了有机结构控制方法[6-8]。
有机结构控制方法是一种以LMI算法为工具的鲁棒控制方法。
它是将子系统间的互联项当作结构扰动处理,根据Lyapunov理论和Schur补定理将鲁棒控制问题转化成几个线性矩阵不等式,通过求解这些线性矩阵不等式就可以得到控制器的反馈阵和允许的最大扰动。
这种方法的主要特点是计算简单,而且在控制器的设计中很容易加入其它设计要求,诸如增益矩阵的结构和大小,来保证控制器的可实现性。
通过此方法设计的控制器可以保证系统在结构扰动下(某些子系统脱离了大系统或又重新连接上)仍是稳定,即保证系统的关联稳定性。
文[9]将此控制方法应用到电力系统当中,并对系统进行了仿真研究。
考虑到系统中有些状态可能不可测,文[10]又给出了带有观测器的自主分散控制方法。
以H技术和LMI方法作为工具,应用于互联系统的鲁棒控制研究是有意义的。
在这方面,已取得了一定的研究成果。
H控制问题可以表述为:
寻找内稳定的控制器使由扰动到误差信号的传递函数的无穷范数趋于极小。
H控制有两种控制规律:
输出反馈和状态反馈。
它的求解方法有很多,包括Riccati方法和LMI方法等。
前一种方法需要对可变参数进行选择,只有选取了合适的参数才能得到方程的解,而参数的选择又没有好的办法,因此这种方法有一定的盲目性。
而LMI方法无需对参数进行调整,因而具有数值易解性,是目前求解H控制器的一种比较好的方法。
文[11]将最优H控制思想应用于电力系统。
文[12]借用LMI针对一类互联双线性随机时滞、时变和参数不确定的电力系统,重点分析了系统鲁棒稳定性和分散鲁棒镇定问题。
文[13]介绍一种适合于多区域互联电力系统的AGC设计方法,即分散的H输出反馈控制器的设计方法。
它的基本思路为将多区域互联电力系统分解为两两互联子系统,利用包含原理扩展系统,在扩展空间中根据LMI算法[14]设计分散鲁棒的H控制器,然后收缩并应用于原系统。
文[15]中作者首先利用包含原理的约束条件对系统进行重叠结构分解,而后利用基于LMI算法的H控制,实现了以两两区域互联子系统控制为基础的多区域系统分散控制器的设计。
结果表明系统的鲁棒性明显提高,并具有令人满意的动态性能。
文[16]基于随机包含原理,给出了几种特殊的重叠分解方法,基于此研究了两区域互联电力系统的分散LQG控制问题。
它的研究主要致力于重叠结构分解方法中分解因子的选择与分散LQG控制性能之间的关系。
文[17]在给出分散的H控制器设计方法基础上,提出了一种比较新颖的鲁棒设计方法,这种设计方法更加实用。
它的主要特点是应用遗传算法对受到LMI约束的比例积分控制器的参数进行寻优。
此外LMI方法还应用在线性不确定时滞系统的鲁棒控制、不确定性系统滤波、变结构控制和模糊控制等方面。
由于现代计算机技术、人工智能微电子学等科学的高速发展,人工智能技术逐渐成为人们的研究热点。
智能控制主要包括模糊控制、人工神经网以及专家系统等。
相对于经典控制理论和现代控制理论来说,此种控制方法无须知道被控对象(或生产过程)的数学模型。
模糊控制是根据经验控制行为,遵循反馈及反馈控制思想,总结成一系列条件语句,即控制规则,运用微机的程序来实现这些控制规则。
而电力系统本身是一个含有众多约束条件的非线性系统,难以用精确的某个数学模型来描述,并且在已进行的AGC控制中积累了大量的数据和经验,将这些数据和经验量化为控制手段时,模糊控制是一个比较好的选择。
在这方面已经取得了一定的研究成果。
文[18]利用神经网络辨识电力系统的动态模型,并将模型逻辑技术与神经网络相结合、通过动态寻优确定最优性能指标下的PID控制器参数,具体设计了一种两区域负荷频率的神经模糊自适应PID控制器。
使两区域负荷频率控制既有非线性控制作用和自学习自适应能力,又有PID控制的广泛适用性。
针对电网大的负荷变化对自动发电控制的影响,文[19]利用人工神经网络研究互联电力系统的自动发电控制,将BP算法用作ANN学习规则。
这种方法在人工智能技术迅速发展的今天,也是重叠电力系统分散控制的可行方法之一。
综上所述,在互联大系统的设计中,鲁棒分散控制占有重要地位。
如何利用LMI算法,结合分散控制思想,为系统设计鲁棒分散控制器,使整个系统镇定,同时在个别子系统断开或加入大系统后,如何保持整个系统的控制性能不受影响,这些对大系统的分散控制都是十分关键的问题。
目前,虽然对于此类问题已经有一些研究,但是仍然不够系统和深入。
同时,由于互联系统的广泛存在,使得这些研究具有重要的理论价值及实际指导意义。
基于此,在这篇文章中我们以互联大系统为例对这些问题进行了仿真研究。
1.2选择课题的目的及意义
本文主要研究扩展结构大系统的鲁棒控制设计,主要考虑了扩展结构系统的鲁棒分散关联镇定问题,利用Lyapunov稳定性理论和线性矩阵不等式(LMI)方法,给出了扩展结构大系统的鲁棒关联镇定的条件,在不改变原系统结构分散控制律的基础上,来设计新加入子系统的分散控制律,使新加入的子系统以及整个扩展后的大系统都能被镇定。
具体内容如下:
第二章主要阐述了线性矩阵不等式的一般形式和LMI工具箱的介绍,第三章介绍了LMI算法在有机结构控制中的应用和LMI中的拆分矩阵法,第四章研究了扩展结构系统的建模与设计,第五章用MATLAB仿真软件对扩展结构系统的鲁棒关联镇定问题进行了仿真研究。
本研究课题具有重要的实际意义,可以间接的带来明显的经济效益和社会效益,这方面的研究成果可以广泛的用于城市交通控制、大型电力系统以及国民经济计划管理等大系统领域。
2线性矩阵不等式(LMI)
随着控制技术的发展,当今在反馈控制系统的设计中,常常需要考虑系统的不确定性,还要考虑系统的鲁棒性。
在处理不确定性系统的鲁棒控制问题及其控制理论中引起的其它控制问题时,都可转化成一种称为线性矩阵不等式LMI(linearmatrixinequation)或带有线性矩阵不等式限制条件的最优化问题。
LMI方法以其高效的求解算法和能获得全局最优解的特点,已引起控制界的关注,为鲁棒控制分析与设计的重要方法。
下面主要介绍线性矩阵不等式的一些基本概念和MATLAB中的LMI工具箱。
2.1线性矩阵不等式一般形式
随着控制技术的发展,系统鲁棒性能的分析和综合问题受到了人们的广泛重视和研究,并取得了一系列的研究成果和方法。
特别地,随着线性矩阵不等式及求解凸优化问题的内点法的提出,为许多控制问题的分析和求解提供了有效工具。
MATLAB软件中线性矩阵不等式工具箱的推出使得各种线性矩阵不等式问题的求解更加方便、直接,从而,进一步推动了线性矩阵不等式处理方法在系统和控制领域中的应用,成为鲁棒控制分析与设计的重要方法。
下面首先就线性矩阵不等式问题作简单的介绍。
在控制工程中,许多控制问题尤其是鲁棒控制问题,都可转化成一种称为线性矩阵不等式或带有线性矩阵不等式限制条件的最优化问题的求解。
线性矩阵不等式一般形式如下:
(2.1)
其中是变量,,是已知的实对称阵。
实际应用时,通常遇到的LMI并不呈现式(2.1)的形式,其中变量不是向量而是一个(或多个)实矩阵,但它可以等价地转化成式(2.1)的形式。
例如,对于以下矩阵不等式
(2.2)
其中,矩阵变量,A,B,C为适当维数的已知实矩阵,B为对称阵。
将X按矩阵空间的简单基展开,即,为X的一个元素,则式(2.2)可以写成LMI的标准形如式(2.1)所表示的,即
(2.3)
在具体问题上常常直接采用LMI的紧凑形式,如式(2.2)。
在控制理论中,经常遇到的两种矩阵不等式为
a、李亚普诺夫(Lyapunov)不等式
,(2.4)
b、黎卡提(Riccati)不等式
(2.5)
显然,式(2.4)是线性矩阵不等式,式(2.5)由于含二次项,故此式是二次矩阵不等式而不是线性矩阵不等式,但利用下面schur补引理,可容易将其变成线性矩阵不等式,即
(2.6)
LMI的基本变换引理如下:
schur补引理[20]:
对给定的对称矩阵,其中是维的。
以下三个条件是等价的:
(ⅰ)(2.7)
(ⅱ)(2.8)
(ⅲ)(2.9)
利用该引理可以将一些非线性矩阵不等式转化为线性矩阵不等式。
S—过程(S—procedure)[21]:
设是实对称阵,要求它们之间满足以下关系:
对于
(2.10)
有
(2.11)
通常很难确定以上关系成立的解析条件。
但如果引入以下充分条件,则以上关系成立,即若存在标量,使以下LMI
(2.12)
成立。
显然式(2.12)是一个LMI可行解问题。
在控制系统的鲁棒分析和鲁棒综合中,我们常常要用S—procedure来将一些不是凸约束的问题转化成线性矩阵不等式约束。
对线性矩阵不等式的求解一般可以归纳为以下三类问题。
1.可行性问题
寻找一个(或等价的:
具有给定结构的矩阵),使得满足线性矩阵不等式系统
(2.13)
2.具有线性矩阵不等式约束的一个目标函数的最小化问题
(2.14)
3.具有线性矩阵不等式限制条件的广义特征值最小化问题
满足于(2.15)
在控制理论中,大多数控制问题都可以转化成上述三种矩阵不等式问题中的一种,从而使问题得到解决。
2.2LMI工具箱简介
在60年代,已经提出了线性矩阵不等式,但由于求解形如式(2.13)~(2.15)所描述的线性矩阵不等式的算法还不够成熟。
再加上求解量大,因而线性矩阵不等式在实际中未得到广泛应用。
近几年来,由于线性矩阵不等式的理论不断完善,求解算法也不断成熟,加上计算机的广泛应用,线性矩阵不等式的求解变得很方便,因此线性矩阵不等式在实际工程中尤其在控制工程理论中得到广泛的应用。
由于用线性矩阵不等式求解控制理论中的问题是当今控制理论发展的一个重要方向,因此出现了许多计算机应用软件,其中以美国MathsWorks.Inc公司用C语言开发的MATLAB软件最为流行;
到目前为止,已相继推出了几个版本,其中在MATLAB5.1、MATLAB5.3、MATLAB6.5等版本中,增加了用于求解线性矩阵不等式的线性矩阵不等式控制工具箱。
线性矩阵不等式工具箱提供了在鲁棒控制设计中所遇到的凸最优化问题的解,同时给出了一个用于求解线性矩阵不等式的集成环境。
由于这个工具箱功能强大和友好的用户界面,因此可以开发自己的应用程序。
这里我们只介绍工具箱中几个重要函数。
1.setlmis([]):
初始化新的LMI系统。
2.lmivar(type,struct):
增加新的矩阵变量X到当前的LMI系统中。
其中,type(类型):
根据变量X的不同类型设置(1~3),1表示矩阵变量X为对称块对角阵,2表示矩阵变量X为满秩阵,3表示矩阵变量X为其它;
struct(结构):
若type=1,则struct的第i行描述X的第i个块对角阵,其中struct(i,1)代表块的大小,struct(i,2)代表块的性质,如果是尺度块t*I,则struct(i,2)取0,如果是满块,则取1,如果是0块,则取-1。
若type=2,假如X是M×
N矩阵,则struct=[M,N]。
若type=3,则struct是一个与X同维的矩阵,其中,struct(i,j)取值为:
当X(i,j)=0时,struct(i,j)=0,当X(i,j)为第n个待求变量时,struct(i,j)=+n,当X(i,j)为第n个待求变量乘上(-1)时,struct(i,j)=-n。
3.lmiterm(termID,A,B,flag):
给当前描述的LMI系统中的某个LMI增加一项。
其中,termID为4输入向量,用来指定项的位置和性质。
对于termID
(1):
若该项位于第n个LMI的左边,则termID
(1)=+n,若该项位于第n个LMI的右边,则termID
(1)=-n。
对于termID(2:
3):
若该项属于LMI的第(i,j)块,则termID(2:
3)=[i,j],若该项属于外部因子,则termID(2:
3)=[00]。
对于termID(4):
若该项属于常数项,则termID(4)=0,若该项属于变量项A*X*B,则termID(4)=m,若该项属于变量项A*XT*B:
termID(4)=-m,其中,m为由函数lmivar返回的变量X的标识。
A可以是外部因子,常数项或者变量项A*X*B或A*XT*B的左系数,B是变量项A*X*B或A*XT*B的右系数。
Flag:
设置flag=’s’,在一
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 扩展 结构 系统 鲁棒控制 研究