《控制系统数字仿真与cad》第8章控制系统的计算机辅助设计

1、1 本章内容(1) 利用MATLAB实现串联频率校正的三种方法；(2) 利用MATLAB实现系统状态反馈的两种方法；(3) 利用MATLAB实现系统状态观测器的两种方法；(4) 利用MATLAB实现带状态观测器的状态反馈系统；(5) 利用MATLAB实现系统的解耦；(6) 利用MATLAB实现状态反馈的线性二次型最优控制器 的设计；(7) 利用MATLAB实现输出反馈的线性二次型的最优控制。 第8章 控制系统的计算机辅助设计2 控制系统的设计，就是在系统中引入适当的环节，用以对原有系统的某些性能进行校正，使之达到理想的效果，故又称为系统的校正，下面介绍几种常用的系统校正方法的计算机辅助设计实现

2、。3l 应用频率法对系统进行校正，其目的是改变系统的频率特性形状，使校正后的系统频率特性具有合适的低频、中频和高频特性以及足够的稳定裕量，从而满足所要求的性能指标。l控制系统中常用的串联校正装置是带有单零点与单极点的滤波器，若其零点比极点更靠近原点，则称之为超前校正，否则称之为滞后校正。4l1.超前校正装置的特性l设超前校正装置的传递函数为l其频率特性为) 1(11cTsTssG）（) 1(j1j1jcTTG）（jT1T1图8-1 超前校正装置的零极点图o5l 超前校正装置的极坐标图如图8-2所示。l当0变化时，Gc(j)的相位角 0，Gc(j)的轨迹为一半圆，由图可得超前校正的最大超前相位角

3、m为l (8-3)l令l l可得对应于最大相位角m时的频率m为11sin1m0d)(dT1mj01)0(图8-2 超前校正装置的极坐标图m6l超前校正装置的对数坐标图如图8-3所示。l当 l l l l由此可见，超前校正装置是一个高通滤波器(高频通过，低频被衰减),它主要能使系统的瞬态响应得到显著改善，而稳态精度的提高则较小。越大，微分作用越强，从而超调量和过渡过程时间等也越小。0|j|lg20, 0c）（Glg20|j|lg20,c）（Glg10|j|lg20,cm）（GT1T10lg10lg20L( )dB0mm图8-3 超前校正装置的对数坐标图m7l 超前校正装置的主要作用是通过其相位超

4、前效应来改变频率响应曲线的形状，产生足够大的相位超前角,以补偿原来系统中元件造成的过大的相位滞后。因此校正时应使校正装置的最大超前相位角出现在校正后系统的开环剪切频率(幅频特性的交接频率)c处。8l(1)根据性能指标对稳态误差系数的要求,确定开环增益k；l(2)利用确定的开环增益k，画出未校正系统的Bode图，并求出其相位裕量r0和幅值裕量kg；l(3)确定为使相位裕量达到要求值，所需增加的超前相位角，即r -r0+.式中 r为要求的相位裕量，是考虑到系统增加串联超前校正装置后系统的剪切频率要向右移而附加的相位角，一般取=515；l(4)令超前校正装置的最大超前相位角=，则由下式可求出校正装置

5、的参数；ccsin1sin19l(5)若将校正装置的最大超前相位角处的频率作为校正后系统的剪切频率，则有l即l或l由此可见，未校正系统的幅频特性幅值等于-20lg时的频率即为c；0| )(j)(j|lg20coccGGlg20| )(j|lg200| )(j|lg20lg20ococGG1| )(j|coG10l(6) 根据=c，利用下式求参数Tl(7) 画出校正后系统的Bode图，检验性能指标是否已全部达到要求，若不满足要求，可增大值，从第三步起重新计算。l例例8-1 设有一单位反馈系统，其开环传递函数为l要求系统的稳态速度误差系数kv=20(1/s),相位裕量 r500,幅值裕量kg10d

6、B,试确定串联校正装置c1T)2()(0ssksG11l解解 根据l可求出k =40，即l根据串联超前校正的设计步骤，可编写以下m文件。l%Example8_1.mlnumo=40;deno=conv(1,0,1,2);lGm1,Pm1,Wcg1,Wcp1=margin(numo,deno);lr=50; r0=Pm1;w=logspace(-1,3);l mag1,phase1=bode(numo,deno,w);202)2(lim)(lim000kssksssGkssv)2(40)(0sssG12lfor epsilon=5:15lphic=( r-r0+epsilon)*pi/180;l

7、alpha=(1+sin(phic)/(1-sin(phic);li1,ii=min(abs(mag1-1/sqrt(alpha);lwc=w(ii);lT=1/(wc*sqrt(alpha);lnumc=alpha*T,1;denc=T,1;lnum,den=series(numo,deno,numc,denc);lGm,Pm,Wcg,Wcp=margin(num,den);lif(Pm= r);break;endlend 13lprintsys(numc,denc),printsys(num,den)lmag2,phase2=bode(numc,denc,w);lmag,phase=bod

8、e(num,den,w);lsubplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);lgrid;ylabel(幅值 (dB);ltitle(- Go, -. Gc, _ GoGc);lsubplot(2,1,2);semilogx(w,phase,w,phase1,-,w,phase2,-.,w,(w-180-w),:);lgrid;ylabel(相位 (度);xlabel(频率 (rad/sec)ltitle(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=

9、,num2str(Pm),);ldisp(校正前：幅值裕量=,num2str(20*log10(Gm1), dB,相位裕量=,num2str(Pm1),);ldisp(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=,num2str(Pm),);14l执行后可得如下结果及图8-4所示曲线lnum/den = l 0.22541 s + 1l -l 0.053537 s + 1lnum/den = l 9.0165 s + 40l -l 0.053537 s3 + 1.1071 s2 + 2 sl 校正前：幅值裕量=Inf dB,相位裕量=17.9642 l 校正

10、后：幅值裕量=Inf dB,相位裕量=50.719615l 图8-4 超前校正装置及校正前后系统的伯德图16l1.滞后校正装置的特性l设滞后校正装置的传递函数为l其频率特性为) 1(11cTsTssG）（) 1(j1j1(jcTTG）图8-5 滞后校正装置的零极点图jT1T1o17l 滞后校正装置的极坐标图如图8-6所示。由图可知，当=0变化时，Gc(j)的相位角0，Gc(j)的根轨迹为一半圆。l同理可求得最大滞后相位角和对应的频率分别为11sin1mT1m11)(0j0m图8-6 滞后校正装置的极坐标图18l滞后校正装置的对数坐标图如图8-7所示l由此可见，滞后校正装置是一个低通滤波器（低频

11、通过，高频被衰减），且越大，高频衰减越厉害，抗高频干扰性能越好，但使响应速度变慢，故滞后校正能使稳态得到显著提高，但瞬态响应时间却随之而增加，越大，积分作用越强，稳态误差越小。0|lg20, 0c）（jGlg20|lg20,c）（jGlg10|lg20,cm）（jGlg20lg10T1T1mmdB00)(L图8-7 滞后校正装置的极坐标图m19l 滞后校正装置的主要作用是在高频段造成幅值衰减，降低系统的剪切频率，以便能使系统获得充分的相位裕量，但应同时保证系统在新的剪切频率附近的相频特性曲线变化不大。l利用频率法设计滞后校正装置的步骤利用频率法设计滞后校正装置的步骤:l（1）根据性能指标对稳态

12、误差系数的要求，确定开环增益k；l（2）利用已确定的开环增益k，画出未校正系统的Bode图，并求出其相位裕量r0和幅值裕量kg；20l（3）如未校正系统的相位和幅值裕量不满足要求，寻找一新的剪切频率c，在c处开环传递函数的相位角应满足下式lGo(jc)= -180+ r +l式中 r为要求的相角裕量，是为补偿滞后校正装置的相位滞后而附加的相位角，一般取=512；l(4) 为使滞后校正装置对系统的相位滞后影响较小（一般限制在512），m应远离c，一般取滞后校正装置的第一个交接频率：1=1/T=(1/51/10)c（即m=r);break;end;lend24lprintsys(numc,denc

13、);printsys(num,den);lmag2,phase2=bode(numc,denc,w);mag,phase=bode(num,den,w);lsubplot(2,1,1);lsemilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);lgrid;ylabel(幅值 (dB);title(- Go, -. Gc, _ GoGc);lsubplot(2,1,2);lsemilogx(w,phase,w,phase1,-,w,phase2,-.,w,(w-180-w),:);lgrid;ylabel(相位 (度);x

14、label(频率 (rad/sec)ltitle(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=,num2str(Pm),);ldisp(校正前：幅值裕量=,num2str(20*log10(Gm1), dB,相位裕量=,num2str(Pm1),);ldisp(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=,num2str(Pm),);25l执行后可得如下结果及图8-8所示曲线。lnum/den = l 8.3842s + 1l -l 59.7135s + 1lnum/den = l 41.9208s+5l -l 14.92

15、84s4+74.8918s3+60.9635s2+sl校正前：幅值裕量= -3.8573e-015 dB，相位裕量=7.3342e-006l校正后：幅值裕量=15.8574dB，相位裕量=40.655226l图8-8 滞后校正装置及校正前后系统的伯德图27l1.滞后-超前校正装置的特性l设滞后-超前校正装置的传递函数为l上式等号右边的第一项产生超前网络的作用，而第二项产生滞后网络的作用。) 1,(1111G122211cTTsTsTsTsTs）（28l 滞后-超前校正装置的极坐标图如图8-9所示。l 由图可知，当角频率在00之间变化时, 滞后-超前校正装置起着相位滞后校正的作用；当在0之间变化

16、时,它起着超前校正的作用，对应相位角为零的频率0为2101TT010mIeR0图8-9 滞后-超前校正装置的极坐标图29l 滞后-超前校正装置的对数坐标图如图8-10所示。从图可清楚看出，当00时滞后-超前校正装置起着相位滞后校正的作用；当01)l numc1=1/w1,1;denc1=1/(w1*beta),1;l w1=w1+0.01;l numc,denc=series(numc1,denc1,numc2,denc2);l num,den=series(numo,deno,numc,denc);l mag,phase=bode(num,den,w);lendlprintsys(numc1

17、,denc1);lprintsys(numc2,denc2);lprintsys(num,den);lGm,Pm,Wcg,Wcp=margin(num,den);lmag2,phase2=bode(numc,denc,w);lmag,phase=bode(num,den,w);35lsubplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),-,w,20*log10(mag2),-.);lgrid;ylabel(幅值 (dB);ltitle(- Go, -. Gc, _ GoGc);lsubplot(2,1,2);semilogx(w,ph

18、ase,w,phase1,-, w,phase2,-., w,(w-180-w),:);lgrid;ylabel(相位 (度);xlabel(频率 (rad/sec)ltitle(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=,num2str(Pm),);ldisp(校正前：幅值裕量=,num2str(20*log10(Gm1), dB,相位裕量=,num2str(Pm1),);ldisp(校正后：幅值裕量=,num2str(20*log10(Gm), dB,相位裕量=,num2str(Pm),);36l执行后可得如下结果及图8-11所示曲线。l num/de

19、n = l 1.8817 s + 1l -l 0.18817 s + 1lnum/den = l 7.0711 s + 1l -l 70.7107 s + 1lnum/den = l 133.0595 s2 + 89.5281 s + 10l -l 6.653 s5 + 55.4084 s4 + 120.1542 s3 + 72.3989 s2 + sl校正前：幅值裕量= -10.4567 dB，相位裕量= -28.0814l校正后：幅值裕量=13.7848 dB，相位裕量=52.4219 37l图8-11 滞后超前校正装置及校正前后系统的伯德图3839l 状态反馈是将系统的状态变量乘以相应的

20、反馈系数，然后反馈到输入端与参考输入叠加形成控制作为受控系统的控制输入，采用状态反馈不但可以实现闭环系统的极点任意配置，而且也是实现解耦和构成线性最优调节器的主要手段。40l 给定控制系统的状态空间模型，则经常希望引入某种控制器，使的该系统的闭环极点移动到某个指定位置，因为在很多情况下系统的极点位置会决定系统的动态性能。l假设系统的状态空间表达式为l其中 A:nn;B:nr;C:mnl引入状态反馈，使进入该系统的信号为l u = r Kxl式中 r为系统的外部参考输入，K为rn矩阵。CxyBuAxx 41可得状态反馈闭环系统的状态空间表达式为 (8-12) 可以证明，若给定系统是完全能控的，则

21、可以通过状态反馈将该系统的闭环极点进行任意配置。假定单变量系统的n个希望极点为1,2,n，则可求出期望的闭环特征方程为 f*(s)= (s-1)( s-2) (s-n)= sn + a1sn-1 + + an这时状态反馈阵K可根据下式求得 K=0 0 1Uc-1f*(A) (8-13)式中 Uc=b Ab An-1b，f*(A)是将系统期望的闭环特征方程式中的s换成系统矩阵A后的矩阵多项式。CxyBrxBKAx)(42l例例8-4已知系统的状态方程为l采用状态反馈，将系统的极点配置到-1,-2,-3，求状态反馈阵K。l解解 MATLAB程序为l%Example8_4.mlA=-2 1 1;1

22、0 1;-1 0 1; b=1;1;1;lUc=ctrb(A,b);rc=rank(Uc);lf=conv(1,1,conv(1,2,1,3);lK=zeros(1,length(A)-1) 1*inv(Uc)*polyvalm(f,A) uxx11110110111243l执行后得lK=l -1 2 4l其实，在MATLAB的控制系统工具箱中就提供了单变量系统极点配置函数acker( )，该函数的调用格式为lK=acker(A,b,P)l式中 P为给定的极点，K为状态反馈阵。l对例8-4，采用下面命令可得同样结果lA=-2 -1 1;1 0 1;-1 0 1;b=1;1;1;rc=rank(

23、ctrb(A,b);lp=-1,-2,-3;K=acker(A,b,p)l结果显示lK=l -1 2 444l 对 于 多 变 量 系 统 的 极 点 配 置 ，MATLAB控制系统工具箱中也给出了函数place( )，其调用格式为lK=place(A,B,P)l例例8-5 已知系统的状态方程为l 求使状态反馈系统的闭环极点为-2-3,l(-1j3)/2的状态反馈阵K。uxx331134029514102033821310140045l解解 MATLAB程序为 l%Example8_5.mlA=0 0 4 1;10 13 2 8;-3 -3 0 -2;-10 -14 -5 -9;lB=-2 0

24、;4 -3;-1 1;-3 3;lP=-2;-3;(-1+sqrt(3)*j)/2;(-1-sqrt(3)*j)/2;lK=place(A,B,P)l执行后得lK=l 32.5923 65.6844 58.8332 46.6557l 55.4594 111.8348 103.6800 81.023946l 在一些特定的应用中，有时没有必要去对所有的极点进行重新配置，而只需对其中若干个极点进行配置，使得其他极点保持原来的值，例如若系统开环模型是不稳定的，则可以将那些不稳定的极点配置成稳定的值，而不去改变那些原本稳定的极点。作这样配置的前提条件是原系统没有重极点，这就能保证由系统特征向量构成的矩阵

25、是非奇异的。47l 假设xi为对应于i的特征向量，即A xi =i xi，这样可以对各个特征值构造特征向量矩阵X=x1,x2,xn，由前面的假设可知X矩阵为非奇异的，故可以得出其逆阵T=X-1，且令T的第i个行向量为Ti，且想把i配置到i的位置，则可以定义变量ri=（i-i）/bi，其中bi为向量Tb的第i个分量，这时配置全部的极点，则可以得出状态反馈阵 l 特别地，若不想对哪个极点进行重新配置，则可以将对应的项从上面的求和式子中删除就可以得出相应的状态反馈阵，它能按指定的方式进行极点配置。niiiTrK148l例例8-6 对于例8-4所示系统，实际上只有一个不稳定的极点1,若仅将此极点配置到

26、-5，试采用部分极点配置方法对其进行。l解解 MATLAB程序为l%Example8_6.mlA=-2 -1 1;1 0 1;-1 0 1;lb=1;1;1;lX,D=eig(A);lii=find(diag(D)0);T=inv(X);mu= -5;V=diag(D);lbb=T*b;lgamma=(mu-V(ii)/bb(ii);lK=real(gamma*T(ii,:)l执行后得l K=l 1.5000 -1.5000 -6.000049l1.全维状态观测器的设计极点配置是基于状态反馈，因此状态x必须可量测，当状态不能量测时，则应设计状态观测器来估计状态。l对于系统l若系统完全能观测，则

27、可构造如图8-12所示的状态观测器。CxyBuAxx 50uyBACBACL-+-图8-12 带状态观测器的系统-+xx 51l由上图可得观测器的状态方程为l即 l其特征多项式为 f（s）=|sI-(A-LC)|l由于工程上要求能比较快速的逼近x，只要调整反馈阵L，观测器的极点就可以任意配置达到要求的性能，所以，观测器的设计与状态反馈极点配置的设计类似。LyxLCBuxAxLyBuxLCAx)(52l假定单变量系统所要求的n个观测器的极点为1,2,n，则可求出期望的状态观测器的特征方程为l f*(s)=(-1)(-2)(-n)=sn + a1sn-1 + + anl这时可求得反馈阵L为l式中

28、， , f*（A）是将系统期望的观测器特征方程l中s换成系统矩阵A后的矩阵多项式。100)(10*VAfL10nCACACV53l 利用对偶原理，可使设计问题大为简化，求解过程如下：l首先构造系统式(8-14)的对偶系统l (8-15)l然后，根据下式可求得状态观测器的反馈阵L。l LT=acker(AT,CT,P)l或 LT=place(AT,CT,P)l其中 P为给定的极点，L为状态观测器的反馈阵。zBwnCAzTTT54l例例8-7 已知开环系统l其中l设计全维状态观测器，使观测器的闭环极点为l-2j23,-5。CxybuAxx 001,100,6116100010CbA55l解解 为求

29、出状态观测器的反馈阵L，先为原系统构造一对偶系统。l然后采用极点配置方法对对偶系统进行闭环极点位置的配置，得到反馈阵K，从而可由对偶原理得到原系统的状态观测器的反馈阵L。lMATLAB程序为lA=0 1 0;0 0 1;-6 -11 -6;b=0;0;1;C=1 0 0;ldisp(The Rank of Obstrabilaty Matrix)lr0=rank(obsv(A,C),A1=A;b1=C;C1=b;lP=-2+2*sqrt(3)*j -2-2*sqrt(3)*j -5;lK=acker(A1,b1,P);L=K zBwnCAzTTT56l执行后得lThe Rank of Obst

30、rabilaty Matrixlr0 =l 3lL =l 3.0000l 7.0000l -1.0000l由于rankr0=3，所以系统能观测，因此可设计全维状态观测器。57l 前面所讨论的状态观测器的维数和被控系统的维数相同，故称为全维观测器，实际上系统的输出y总是能够观测的。因此，可以利用系统的输出量y来直接产生部分状态变量，从而降低观测器的维数。假设系统是完全能观测器，若状态x为n维，输出y为m维，由于y是可量测的，因此只需对n-m个状态进行观测，也就是说用(n-m)维的状态观测器可以代替全维观测器，这样观测器的结构可以大大简化。58 状态观测器解决了受控系统的状态重构问题，为那些状态变

31、量不能直接量测得到的系统实现状态反馈创造了条件。带状态观测器的状态反馈系统由三部分组成，即原系统、观测器和控制器，图8-13是一个带有全维观测器的状态反馈系统。59uryBACBACL-+-图8-13 带状态观测器的状态反馈系统-+Kxx +-60l设能控能观测的受控系统为l （8-21）l状态反馈控制律为l （8-22）l状态观测器方程为l （8-23）l由以上三式可得闭环系统的状态空间表达式为CxyBuAxx xKruLyBuxLCAx)(CxyBrxBKLCALCxxBrxBKAxx)(61l 可以证明，由观测器构成的状态反馈闭环系统，其特征多项式等于状态反馈部分的特征多项式|sI-(A

32、-BK)|和观测器部分的特征多项式|sI-(A-LC)|的乘积，而且两者相互独立。因此，只要系统0(A,B,C)能控能观测，则系统的状态反馈阵K和观测器反馈阵L可分别根据各自的要求，独立进行配置，这种性质被称为分离特性。l同理，用降维观测器构成的反馈系统也具有分离特性62l例例8-9 已知开环系统l设计状态反馈使闭环极点为 -1.8j2.4，而且状态不可量测，因此设计状态观测器使其闭环极点为-8,-8。l解解 状态反馈和状态观测器的设计分开进行，状态观测器的设计借助于对偶原理。在设计之前，应先判别系统的能控性和能观测性，MATLAB的程序为l%Example8_9.mlA=0 1; 20.6

33、0;b=0;1;C=1 0;xyuxx011006 .201063l% Check Contrillability and Observablityldisp (The rank of Controllability Matrix)lrc=rank (ctrb(A,b)ldisp (The rank of Observability Matrix)lro=rank(obsv(A,C)l%Design RegulatorlP=-1.8+2.4*j -1.8-2.4*j;lK=acker(A,b,P)l%Design State ObserverlA1=A;b1=C;C1=b;lP1=-8 -8;

34、K1=acker(A1,b1,P1);lL=K164l执行后得lThe rank of Controllability Matrixlrc =l 2lThe rank of Observability Matrixlro =l 2lK =l 29.6000 3.6000lL =l 16.0000l 84.600065l 离散系统的极点配置和状态观测器的设计的求解过程与连续系统基本相同，在MATLAB中，可直接采用工具箱中的place( )和acker( )函数进行设计，这里不在赘述。66l 在多变量系统中，如果传递函数阵不是对角矩阵，则不同的输入与输出之间存在着耦合，即第i输入不但会对第i输出

35、有影响，而且还会影响到其他的输出，就给控制系统的设计造成了很大的麻烦，故在多变量控制系统的设计中就出现了解耦控制方法。67l假设控制系统的状态空间表达式为l (8-25) l其中 A:nn;B:nr;C:mn;D:mrl引入状态反馈l (8-26)l其中 R为r1参考输入向量，在解耦控制中实际还应要求r=m，亦即系统的输入个数等于输出个数，这时闭环系统的传递函数矩阵可以写成KxHRuHDBBKAsIDKCsRsYsG)()()()(1DuCxyBuAxx 68l若闭环系统的mr矩阵G(s)为对角的非奇异矩阵，则称该系统是动态解耦的系统，若G(0)为对角非奇异矩阵，且系统为稳定的，则称该系统是静

36、态解耦的。l在给定的控制结构下，若系统的D矩阵为0，则闭环传递函数阵G(s)可以简化成l (8-28)l由上式可见，若H矩阵为奇异矩阵，则G(s)矩阵必为奇异的，所以为使得系统可以解耦，首先应该要求H为非奇异矩阵。l其中BHBKAsICsRsYsG1)()()()(HsCBaBBKACCBsBKAsIsGnnn)(|1)(211.|011asasBKAsInnn69l首先这里将给出能解耦的条件：可以证明，若按下面方法生成的矩阵B*为非奇异的，若取H=(B*)-1，则由前面给出的控制格式得出的系统能解耦原系统。l （8-29） l式中 C,C, Cm为C矩阵的行向量，参数d,d,dm是在保证B为

37、非奇异的前提下任选区间0,n-1上的整数。若确定了di参数，则可以直接获得解耦矩阵BACBACBACBmdmdd21211111mdmdACACHK70l例例8-11 对如下系统进行解耦l解解 MATLAB程序为l%Example8_11.mlA=-1 0 0;0 -2 0;0 0 -3;B=1 0;2 3;-3 -3;C=1 0 0;1 1 1;lD=0 0;0 0;lm,n=size(C);lBB=C(1,:)*B;ld(1)=0;xyuxx11100133320130002000171lfor ii=2:ml for jj=0:n-1l BB=BB;C(ii,:)*Ajj*B;l if

38、rank(BB)=ii,d(ii)=jj;break;else BB=BB(1:ii-1,:);endl endlendlH=inv(BB)lCC=C(1,:)*A(d(1)+1);lfor ii=2:mlCC=CC;C(ii,:)*A(d(ii)+1);lendlK=H*CCln1,d1=ss2tf(A-B*K,B*H,C,D,1)ln2,d2=ss2tf(A-B*K,B*H,C,D,2)72l执行后可得lH =l 1.0000 0l -1.3333 0.3333lK =l -1.0000 0 0l 1.6667 1.3333 3.0000ln1 =l 0 1.0000 -0.0000 -0

39、.0000l 0 0.0000 0.0000 0.0000ld1 =l 1.0000 -0.0000 -0.0000 0ln2 =l 0 0 0 0l 0 0.0000 1.0000 0ld2 =l 1.0000 -0.0000 -0.0000 073l亦即系统解偶后的传递函数阵为ll 解耦控制系统的目的是将原模型变换成解耦的模型，而并不必去考虑变换之后的响应品质，因为响应品质这类问题可以在解耦之后按照单变量系统进行设计补偿，单回路的设计当然可以采用单变量系统的各种方法，例如可以采用超前滞后补偿，PI设计以及PID设计等，并能保证这样设计出来的控制器不会去影响其他回路。1001001)(223

40、ssssssG74l假设控制系统的状态空间表达式为l函数estim( )将生成下述状态和输出估计器DuCxyBuAxx xICxyxCyLxAx)( 75l在MATLAB中，函数estim( )的调用格式如下lest=estim(A,B,C,D,L)l其中 A,B,C,D为系统系数矩阵，L为状态估计增益矩阵。状态估计增益矩阵L可由极点配置函数place( )形成，或者由Kalman滤波函数kalman生成。利用以上命令可生成给定增益矩阵L下的状态空间模型A,B,C,D的输出估计器est。76l例例8-13 利用例8-7所得的状态观测器的反馈阵L，求其系统的状态估计器。l解解 MATLAB程序为

41、lA=0 1 0;0 0 1;-6 -11 -6;b=0;0;1;C=1 0 0;l L=3;7;-1lest=estim(A,b,C,0,L)l执行后得lest =l -3.0000 1.0000 0l -7.0000 0 1.0000l -5.0000 -11.0000 -6.000077l假设控制系统的状态空间表达式为l利用函数reg( )可生成下述控制器DuCxyBuAxx xKuLyxKLDBLCAx)( 78l在MATLAB中，函数reg( )的调用格式为lest=reg(A,B,C,D,K,L)l其中 A,B,C,D为系统系数矩阵，K为状态反馈增益矩阵，L为状态估计增益矩阵。利用

42、以上命令可生成给定状态反馈增益矩阵K及状态估计增益矩阵L下的状态空间模型A,B,C,D的控制器est。假定系统的所有输出可测。79l例例8-14 利用例8-7所得的状态观测器的反馈阵L，求其系统的控制器。假设状态反馈阵K=-1 2 4。l解解 MATLAB程序为lA=0 1 0;0 0 1;-6 -11 -6;b=0;0;1;C=1 0 0;lK=-1 2 4;lL=3;7;-1lest=reg(A,b,C,0,K,L)l执行后得l est =l -3 1 0l -7 0 1l -4 -13 -1080 MATLAB控制系统工具箱中也提供了很多函数用来进行系统的最优控制设计，相关函数如表8-3

43、所示。81l设线性定常系统的状态空间表达式为l (8-31)l式中 A:nn;B:nr;C:mnl并设目标函数为二次型性能指标l (8-32)l式中 Q(t)为nn半正定实对称矩阵，R(t)为rr正定实对称矩阵。一般情况下，假定这两个矩阵为定常矩阵，它们分别决定了系统暂态误差与控制能量消耗之间的相对重要性。S为对称半正定终端的加权阵，它为常数。)()()()()(tCxtytButAxtx f0d)()()()()()(21)()(21TTffTttttutRtutxtQtxtSxtxJ82l 当x(tf)值固定时，则为终端控制问题，特别是当x(tf)0 时，则为调节器问题；当t0 , tf均固定时，则为暂态过程最优控制。l 最优控制问题是为给定的线性系统式（8-31）寻找一个最优控制律u*(t)，使系统从初始状态x(t0)转移到终端状态x(tf)，且满足性能指标式（8-32）最小。它可以用变分法、极大值原理和动态规划等三种方法中的任一种求解。这里我们采用极大值原理求解u*(t)。83l MATLAB的控制系统工具箱中也提供了完整的解决线性二次型最优控制的函数，其中命令lqr( )和lqry( )可以直接求解二次型调节器问题及相关的