现代控制理论-第7章

《现代控制理论基础》第五章(讲义)1第六次课小结一、Lyapunov意义下的稳定性问题基本概念平衡状态的概念Lyapunov意义下的稳定性定义（稳定，一致稳定，渐进稳定，一致渐进稳定，大范围渐进稳定等）纯量函数的正定性，负定性，正半定性，负半定性，不定性二次型，复二次型（Hermite型）二、Lyapunov稳定性理论第一方法第二方法三、线性定常系统的Lyapunov稳定性分析应用Lyapunov方程QPAPAH来进行判别稳定性四、线性定常系统的稳定自由运动的衰减率性能估计衰减系数，一旦定出min，则可定出)(xV随时间t衰减上界。计算min的关系式五、离散时间系统的状态运动稳定性及其判据离散系统的大范围淅近稳定判据，Lyapunov稳定判据在离散系统中的应用六、线性多变量系统的综合与设计的基本问题《现代控制理论基础》第五章(讲义)2问题的提法性能指标的类型研究的主要内容七、极点配置问题问题的提出可配置条件极点配置算法《现代控制理论基础》第五章(讲义)35.2.5爱克曼公式(Ackermann’sFormula)考虑由式（5.1）给出的系统，重写为BuAxx假设该被控系统是状态完全能控的，又设期望闭环极点为nsss,,,21。利用线性状态反馈控制律Kxu将系统状态方程改写为xBKAx)((5.14)定义BKAA~则所期望的特征方程为0)())((~11121nnnnnasasassssAsIBKAsI由于凯莱-哈密尔顿定理指出A~应满足其自身的特征方程，所以0~~~)~(**11*1*IaAaAaAAnnnn(5.15)我们用式（5.15）来推导爱克曼公式。为简化推导，考虑n=3的情况。需要指出的是，对任意正整数，下面《现代控制理论基础》第五章(讲义)4的推导可方便地加以推广。考虑下列恒等式22333222~~)(~~)(~~ABKAABKBKAABKAAABKABKABKAABKAAII将上述方程分别乘以)1(,,,*0*0*1*2*3aaaaa，并相加，则可得32*1*2*3~~~AAaAaIa32*1*2*3)~()(AABKABKAaBKAaIa22~~ABKAABKBKABKAABKaABKaBKaAAaAaIa2*1*1*232*1*2*3~2~~ABKAABK（5.16）参照式（5.15）可得0)~(~~~*32*1*2*3AAAaAaIa也可得到0)(*32*1*2*3AAAaAaIa将上述两式代入式（5.16），可得BKAAABKABKaABKABKaBKaAA2*12*1*2**~~~)()~(由于0)~(*A，故《现代控制理论基础》第五章(讲义)5BKAAKKaABAKAKaKaBA2*12*1*2*)~()~~()(KAKKaAKAKaKaBAABB~~~][*12*1*22（5.17）由于系统是状态完全能控的，所以能控性矩阵][2BAABBQ的逆存在。在式（5.17）的两端均左乘能控性矩阵Q的逆，可得KAKKaAKAKaKaABAABB~~~)(][*12*1*2*12上式两端左乘[001]，可得KKAKKaAKAKaKaABAABB~~~]100[)(]][100[*12*1*2*12重写为)(][]100[*12ABAABBK从而给出了所需的状态反馈增益矩阵K。《现代控制理论基础》第五章(讲义)6对任一正整数n，有)(]][1000[*11ABAABBKn(5.18)式（5.18）称为用于确定状态反馈增益矩阵K的爱克曼方程。-------------------------------------------------[例5.1]考虑如下线性定常系统BuAxx式中100,651100010BA利用状态反馈控制Kxu，希望该系统的闭环极点为s=-2±j4和s=-10。试确定状态反馈增益矩阵K。首先需检验该系统的能控性矩阵。由于能控性矩阵为：3161610100][2BAABBQ所以得出detQ=-1，因此，rankQ=3。因而该系统是状态完全能控的，可任意配置极点。《现代控制理论基础》第五章(讲义)7下面，我们来求解这个问题，并用本章介绍的3种方法中的每一种求解。方法1：第一种方法是利用式（5.13）。该系统的特征方程为：01566511001||3221323asasasssssssAsI因此1,5,6321aaa期望的特征方程为02006014)10)(42)(42(*3*22*1323asasasssssjsjs因此200,60,14*3*2*1aaa参照式（5.13），可得《现代控制理论基础》第五章(讲义)8]855199[]6145601200[K方法2：设期望的状态反馈增益矩阵为][321kkkK并使||BKAsI和期望的特征多项式相等，可得651100010000000||sssBKAsI321[100kkk20060141)5()6(65110012312233321ssskskskskskkss因此2001,605,146123kkk从中可得8,55,199321kkk《现代控制理论基础》第五章(讲义)9或]855199[K方法3：第三种方法是利用爱克曼公式。参见式（5.18），可得)(]][100[*12ABAABBK由于IAAAA2006014)(23*11743771598855199100010001200651100010606511000101465110001023且3161610100][2BAABB《现代控制理论基础》第五章(讲义)10可得]855199[11743771598855199001016165]100[117437715988551993161610100]100[1K显然，这3种方法所得到的反馈增益矩阵K是相同的。使用状态反馈方法，正如所期望的那样，可将闭环极点配置在s=-2±j4和s=-10处。------------------------------------------------------------------------------应当注意，如果系统的阶次n等于或大于4，则推荐使用方法1和3，因为所有的矩阵计算都可由计算机实现。如果使用方法2，由于计算机不能处理含有未知参数nkkk,,,21的特征方程，因此必须进行手工计算。5.2.6注释对于一个给定的系统，矩阵K不是唯一的，而是依赖《现代控制理论基础》第五章(讲义)11于选择期望闭环极点的位置（这决定了响应速度与阻尼），这一点很重要。注意，所期望的闭环极点或所期望状态方程的选择是在误差向量的快速性和干扰、测量噪声的灵敏性之间的一种折衷。也就是说，如果加快误差响应速度，则干扰和测量噪声的影响通常也随之增大。如果系统是2阶的，那么系统的动态特性（响应特性）正好与系统期望的闭环极点和零点的位置联系起来。对于更高阶的系统，期望的闭环极点位置不能和系统的动态特性（响应特性）联系起来。因此，在决定给定系统的状态反馈增益矩阵K时，最好通过计算机仿真来检验系统在几种不同矩阵（基于几种不同的期望特征方程）下的响应特性，并且选出使系统总体性能最好的矩阵K。5.3利用MATLAB求解极点配置问题用MATLAB易于求解极点配置问题。现在我们来求解在例5.1中讨论的同样问题。系统方程为BuAxx式中《现代控制理论基础》第五章(讲义)12100651100010BA，采用状态反馈控制Kxu，希望系统的闭环极点为s=μi(i=1,2,3)，其中10,42,42321jj现求所需的状态反馈增益矩阵K。如果在设计状态反馈控制矩阵K时采用变换矩阵P，则必须求特征方程|sI-A|=0的系数1a、2a、和3a。这可通过给计算机输入语句P=poly(A)来实现。在计算机屏幕上将显示如下一组系数：则)4(3),3(2),2(1321PaaPaaPaa。为了得到变换矩阵P，首先将矩阵Q和W输入计算机，其中A=[010;001;-1-5-6]；P=poly(A)P=1.00006.00005.00001.0000《现代控制理论基础》第五章(讲义)13][2BAABBQ001011112aaaW然后可以很容易地采用MATLAB完成Q和W相乘。其次，再求期望的特征方程。可定义矩阵J，使得100004200042000000321jjJ从而可利用如下poly(J)命令来完成，即20060141)(poly];1000;0*420;00*42[QJQiiJ因此，有)4(3),3(2),2(1*3*2*1QaaaQaaaQaaa即对于*ia，可采用aai。故状态反馈增益矩阵K可由下式确定：1112233][PaaaaaaK《现代控制理论基础》第五章(讲义)14或))(inv(*]112233[PaaaaaaaaaK采用变换矩阵P求解该例题的MATLAB程序如MATLABProgram5.1所示。MATLABProgram5.1%------Poleplacement------%*****DeterminatonofstatefeedbackgainmatrixKbyuesof%transformationmatrixP*****%*****EntermatricesAandB*****A=[010;001;-1-5-6];B=[0;0;1];%*****DefinethecontrollabilitymatrixQ*****《现代控制理论基础》第五章(讲义)15Q=[BA*BA^2*B];%*****ChecktherankofmatrixQ*****rank(Q)ans=3%*****SincetherankofQis3,arbitrarypoleplacementis%possible*****%*****Obtainthecoefficientsofthecharacteristicpolynomial%|sI-A|.Thiscanbedonebyenteringstatementpoly(A)*****JA=poly(A)《现代控制理论基础》第五章(讲义)16JA=1.00006.00005.00001.0000a1=JA(2);a2=JA(3);a3=JA(4);%*****DefinematricesWandPasfollows*****W=[a2a11;a110;100];P=Q*W;%*****Obtainthedesiredchracteristicpolynomialbydefining%thefollowingmatrixJandente