现代控制工程大作业

1、题目一：1.一、系统设计及仿真分析1.1、确定被控系统状态空间表达式 12、系统能控、能观性判别 13、系统极点配置 14、确定状态反馈增益矩阵 F 25、确定输入变换线性放大器 K 26、验证跟踪单位斜坡输入信号的稳态误差ev 27、利用SIMULINK建立控制系统的动态仿真模型 3(1)跟踪单位阶跃信号的动态仿真分析 3(2)跟踪单位斜坡信号的动态仿真分析 3二、采用全维状态观测器的状态反馈系统 41、配置闭环系统状态观测器极点 42、确定观测器偏差反馈增益矩阵G 43、利用 SIMULINK 建立系统动态仿真模型 5(1)对单位阶跃输入信号仿真分析 5(2)对单位斜坡输入信号仿真分析 6

2、4、与采用直接状态反馈的控制系统仿真结果比较 7(1) x(0) 0且？(0) 0时的仿真结果比较 7(2) x(0) 0且洌0) 2 2 2时的仿真结果比较 75、单位阶跃输入作用下状态估值误差收敛性分析 8(1)收敛过程分析 8(2)收敛速度与状态观测器极点的关系及对系统性能的影响 9三、实现状态反馈的降维观测器设计 101、通过线性变换使状态向量按输出可检测性分解 102、确定降维观测器的反馈增益矩阵 103、列写降维观测器状态方程并重构原系统状态 114、利用 SIMULINK 建立系统动态仿真模型 11(1)跟踪单位阶跃信号的动态仿真分析 12(2)跟踪单位斜坡信号的动态仿真 12题

3、目二：14.一、线性二次型最优控制律设计 141、判断系统是否具有最优控制律 142、非零给定点的最优控制设计及仿真分析 143、权矩阵的各权值对动态性能的影响分析 16(1)固定q22q331,使q11 分别取1,10,100,100016(2)固定q11q331,使q22分别取1,10,100,100017(3)固定qnq221,使q33分别取1,10,100,1000 184、使系统具有良好动、静态性能的线性二次型最优控制 19题目一:设被控系统的传递函数为Go(s) 1 ,现在采用带有输入变换的状态反馈控s(s 6)(s 12)制(如教材图 5-15)改善系统的性能，要求控制系统满足如

4、下性能指标：超调量5% ；调节时间ts 0.5秒；跟踪阶跃输入信号的稳态误差ep 0 ；跟踪单位斜坡输入信号的稳态误差 ev 0.2。(1)根据控制系统综合指标设计状态反馈增益矩阵F及输入变换线性放大器K ,并在SIMULINK中建立控制系统的动态仿真模型，对系统的动态和静态性能进行仿真分析。(2)根据控制系统的响应速度，合理确定状态观测器系统矩阵期望特征值的分布，设计 对状态向量x进行估计的全维状态观测器以实现状态反馈，并保证控制系统的性能仍满足要求；在SIMULINK中建立采用状态观测器实现状态反馈的控制系统的动态仿真模型，对系 统的动态和静态性能进行仿真分析，并与采用直接状态反馈的控制系

5、统仿真结果比较；观察单位阶跃输入作用下状态估计值攵与状态真实值x的误差的收敛过程，并研究状态估计误差收敛速度与状态观测器极点的关系及其对系统性能的影响。(3)设系统输出量可准确测量，试设计实现状态反馈的降维观测器，在SIMULINK中建立采用降维状态观测器实现状态反馈的控制系统的动态仿真模型，对系统的动态和静态性能进行仿真分析。一、系统设计及仿真分析1、确定被控系统状态空间表达式被控系统传递函数为Go(s)1s(s 6)(s 12)则其能观型状态空间表达式为x 1 072 x 0 u0 1180y 0 0 1x2、系统能控、能观性判别1 0 0系统能控性判别矩阵 Qc B AB A2B 0 1

6、 00 0 1则rankQc 3 n ,所以系统完全能控; 统完全能观。又因为其状态空间表达式为能观标准型,所以系综上所述，此系统既能控又能观。3、系统极点配置根据系统性能指标要求：5% , ts 0.5由 e 丫1 25%,可推出 0.690107,取 0.7；由 ts 2 0.5 ,n 10,取 n 11。2根据经典控制理论，二阶系统闭环传递函数标准形式为 pn2 ,将S 2 nS n0.7, n 11代入式s2 2 nS n2 0 ,得期望的闭环主导极点为 1,2 7.7 j7.8556; 选择一个期望的闭环非主导极点离虚轴为主导极点的5倍以上，即3 7.7 6 46.2。则期望的闭环特

7、征多项式为*32p (s) (s 1 )(s2 )(s3 ) S3 61.6s2 832.48s 5590.24、确定状态反馈增益矩阵F由上已知，系统完全能控，所以可以通过状态反馈任意配置系统闭环极点。应用 MATLAB极点配置函数求解状态反馈增益矩阵，代码如下：A=0 0 0;1 0 -72;0 1 -18;B=1;0;0;P=-7.7+j*7.8556;-7.7-j*7.8556;-46.2;F=place(A,B,P)求得状态反馈增益矩阵 F f1 f2 f343.624.3 2888.85、确定输入变换线性放大器K根据题目要求，跟踪单位阶跃参考输入信号的稳态误差ep 1 C( A BF

8、) 1BK 0,求p得K 55902 ,则带有输入变换的状态反馈系统传递函数c ,、 K5590.2GF(s) *32p (s) s 61.6s832.48s 5590.26、验证跟踪单位斜坡输入信号的稳态误差ev将带有输入变换的状态反馈系统等效为具有单位负反馈的闭环系统，系统方框图如下所示：Rs)*图1-1等效单位负反馈闭环系统方框图则系统传递函数(s) GGF(s)1 G(s) F5590.23_2s3 61.6s2 832.48s 5590.2求得，G(s)5590.2332s 61.6s832.48s对于单位斜坡输入信号,R(s) 1/s2；则 ev lim sE(s) lim sR(

9、s) 0.149 0.2s 0s 01 G(s)所以，所设计系统能满足要求系统性能指标。7、利用SIMULINK建立控制系统的动态仿真模型图1-2带有输入变换的状态反馈系统(1)跟踪单位阶跃信号的动态仿真分析,.二一门=割二7-S1为便于输出图形数据的观察比较，在t=0.5s处输入单位阶跃响应，仿真结果如下图所示:Ini图1-3 单位阶跃响应图1-4 单位阶跃响应局部放大图静态性能分析：由图1-3可知，系统满足跟踪阶跃输入信号的稳态误差ep 0。动态性能分析：从局部放大图1-4中易知，超调量 1.05 1 5% ,调节时间ts 0.5秒1(注意：为便于观察，本实验中单位阶跃信号在t=0.5s时

10、刻输入)，系统满足设计要求。(2)跟踪单位斜坡信号的动态仿真分析在t=0处输入单位阶跃响应，仿真结果如下图所示:口 511.522时间1 ( ?)图1-5单位斜坡响应1.31.21.10190 6 070.6 920.94D.9C。，典11021.011.Q6时间t (小图1-6单位斜坡响应局部放大图由图t=1s的时刻知，跟踪单位斜坡信号的稳态误差ev1 0.80.2,故系统满足设计要求。二、采用全维状态观测器的状态反馈系统基于复合系统特征值的分离特性，只要被控系统能控能观，则用状态观测器估值形成状态反馈时，可对系统的状态反馈控制器及状态观测器分别按各自的要求进行独立设计，由(一)中已求得 F

11、 43.624.3 28888。1、配置闭环系统状态观测器极点系统完全能观，所以闭环系统状态观测器极点可任意配置。按通常选择观测器的响应速度比所考虑的状态反馈闭环系统快2 5倍这一选择原则，取观测器期望极点为1,2 2.5 ( 7.7)19 , 3 2.5 ( 46.2)1162、确定观测器偏差反馈增益矩阵G应用MATLA瞅点配置函数求解观测器偏差反馈增益矩阵G,代码如下：A=0 0 0;1 0 -72;0 1 -18,C=0 0 1;P=-19;-19;-116;Gt=acker(A',C',P);G=Gt'求得偏差反馈增益矩阵 G g1 g2 g341876 469

12、7 1363、利用SIMULINK建立系统动态仿真模型观测器的状态方程为：0 0X (A GC)? Gy Bu 1 00 1418764187614769 X4697y0u1541360O利用SIMULINK 建立系统动态仿真模型如下:*入信母1 B.C4 。 U0.OL.小 二)J=以密温宦slaf图1-7采用全维状态观测器的状态反馈系统(注：Scopel可观测误差的收敛过程)(1)对单位阶跃输入信号仿真分析在t=0.5s时刻输入单位阶跃信号，对于x(0) 0和？(0) 0的情况下仿真结果如下图所示：14图1-8 x(0) 0且X(0) 0时单位阶跃响应0.7AS03I1.1121.3叼同t

13、G图1-9 x(0) 0且？(0) 0时单位阶跃响应局部放大图静态性能分析：由图1-8可知，系统满足跟踪阶跃输入信号的稳态误差ep 0。动态性能分析：从局部放大图1-9中易知，超调量 1.0： 1 5% ,调节时间ts 0.5秒(注意：为便于观察，本实验中单位阶跃信号在t=0.5s时刻输入)，系统满足设计要求。(2)对单位斜坡输入信号仿真分析在t=0时刻输入单位斜坡信号，对于各状态变量初值x(0) 0且？(0) 0情况下的仿真结果如下图所示：00311.5225时间I <«>图1-10 x(0) 0且然0) 0时单位斜坡信号响应I.LbQ.酒 03 口.各 0 H09。第

14、1tC&1,1时埠C)二V F.一像衿萧乌今图1-11 x(0) 0且？(0) 0时单位斜坡信号响应局部放大图由图t=1s的时刻知，跟踪单位斜坡信号的稳态误差ev1 0.80.2,故系统满足设计要求。4、与采用直接状态反馈的控制系统仿真结果比较(1) x(0) 0且？(0)0时的仿真结果比较(a)直接状态反馈的单位阶跃响应(b)采用状态观测器的单位阶跃响应(a)直接状态反馈的单位阶跃响应L(b)采用状态观测器的单位阶跃响应U1仃：IR%】聘I(c)直接状态反馈的单位斜坡响应(d)采用状态观测器的单位斜坡响应图1-12x(0) 0且X(0) 0的仿真结果比较从以上仿真结果图可以看出，在各

15、状态变量x(0) 0且X(0) 0的情况下，直接状态反馈与采用状态观测的反馈系统仿真结果几乎一样。(2)x(0) 0且？(0) 2 2 2时的仿真结果比较(c)直接状态反馈的单位斜坡响应(d)采用状态观测器的单位斜坡响应图1-13 x(0) 0且*0)2 2 2的仿真结果比较由图1-13可知，尽管初始时刻x(0) 0和？(0) 2 2 2存在差异，但随着时间的推移, 初始值不同情况下的响应也将趋于同样的稳定状态。这其实与后面讨论的估值误差的收敛性质有关。对于复杂的实际工程应用，一般不可能采用直接状态反馈，而是采用观测器进行状态反馈以获得控制精度并且降低成本。5、单位阶跃输入作用下状态估值误差收

16、敛性分析(1)收敛过程分析对x(0)洌0) 0和x(0) 0,父(0) 2 2 2两种情况下的系统进行收敛性仿真分析, 结果如下所示：QHi1”2时间ocoon Boom 2 i-1I-2图1-14x(0) ?(0) 0时误差收敛情况图1-15 x(0) 0,就0) 2 2 2时误差收敛情况由图可知，当初始值相同，即 x(0)父(0) 0时，状态估值没有误差；当初始值不同，即x(0) 0,X(0) 2 2 2时，开始一段时期存在估值误差，但是，随着时间的推移x(t) ?(t),即估值误差收敛到0。从理论上讲，状态观测器之所以出现估值误差是因为观测器与系统初始状态在实际工程 中难以做到完全一致所

17、造成的。估值误差为x(t) eAGCt?(0) x(0)根据上式，当初始时刻 t0 0, X(0) ?(0)时，显然根据上式 x(t) 0；当X(0) 0, ?(0) 2 2 2时，由于我们通过选择输出偏差反馈增益矩阵G使A GC的所有特征值均位于复平面的左半开平面，尽管初始时刻t0 0时，x(0) 0, ?(0) 2 2 2存在差异，但观测器的状态？(t)仍将以一定精度和速度渐进逼近系统的实际状态Mt)。因此，仿真结果与理论分析相符。(2)收敛速度与状态观测器极点的关系及对系统性能的影响现重新选取全维观测器期望极点为1,2 4.5 ( 7.7)35, 3 4.5 ( 46.2)208应用MA

18、TLA瞅点配置函数求解观测器偏差反馈增益矩阵G,代码如下：A=0 0 0;1 0 -72;0 1 -18,C=0 0 1;P=-35;-35;-208;Gt=acker(A',C',P);G=Gt'求得偏差反馈增益矩阵 Gig1 g2 g3254800 15713 260将Gi代入采用全维状态观测器的状态反馈系统中，并在 SIMULINK中进行动态仿真， 当x(0) 0 ,x(0)2 2 2时，各状态变量的收敛情况如下图所示：I.W项2L00-2E04IJU名co"WTO-12OT1.5时间tG图1-161 235, 3208时的各状态变量的收敛情况由(1)中

19、已知1,219,116情况下的误差收敛情况如下图所示:型等用3M欲降e图1-171,219, 3116时的各状态变量的收敛情况从图1-16和图1-17中各状态变量的收敛情况可以看出观测器极点取不同值时，贷t）的误差收敛情况。由图可知，全维观测器闭环极点离虚轴的距离越远，状态估值误差收敛速度越快。但是，同时也使得反馈增益矩阵中的元素需要很大，从而导致物理实现困难。理论上 讲，估值误差收敛速度越快，系统性能越好，但是，在实际工程应用中，观测器的响应速度 太快会出现放大量测噪声，使系统无法正常工作。 因此，观测器期望极点的选择应兼顾快速性、抗干扰性等折中考虑，通常选择观测器的响应速度比所考虑的状态反

20、馈闭环系统快25倍。三、实现状态反馈的降维观测器设计系统为能观标准型，所以完全能观，又m rankC 1 ,故可构造n m 2维状态观测器。1、通过线性变换使状态向量按输出可检测性分解构造n n非奇异变换矩阵tTn mC作线性变换x T又，则将(A,B,C)变换成(A,B,C),即x TAT 1 xy CT 1x0 0TBu 1 00 10 0 ix0172 x 0 u180由于xx3y ,故只需设计二维观测器重构xXi一。将A,B分块，得X20 0A11, A121 00,A210 1 ,A2272118,场 ，B2002、确定降维观测器的反馈增益矩阵G1 g1 g2g12 一 f( ) d

21、et I (A11 G1A21)detg2 g11 g2配置观测器极点为1,2 2.5 ( 7.7)19则期望多项式f ( ) (19)22 38361比较f()和f ()各项系数，联立方程并解得G1361 38T3、列写降维观测器状态方程并重构原系统状态w (An GA21)w (B1 G1B2)u(A11 G1A2DG1A12G1A22 y0361w11u138w20?X1wG1yw1X2w27220y471361Wi361 y38 yw238y则 (A,B,C)所对应的状态向量X的估值为?X3yW1 361yw2 38yy将夕变换为原系统状态空间，得到原系统(A,B,C)的状态重构为1

22、0 0 w1 361y 1 oG T M 0 1 0 w2 38y 0 0 1 yw1 361 yw2 38 yy4、利用SIMULINK建立系统动态仿真模型图1-18采用降维状态观测器的状态反馈系统(1)跟踪单位阶跃信号的动态仿真分析在t=0.5s时刻输入单位阶跃信号，对于原系统状态变量初始值x(0) 0,降维观测器状态初始值w(0) 2 2的仿真结果如下：图1-19 采用降维观测器的单位阶跃响应粕A僚日 粕上不应-0 353Suati061 05110.9口珏时间B二；IU同国之0H图1-20采用降维观测器的单位阶跃响应局部放大图静态性能分析:由图可知，系统满足跟踪阶跃输入信号的稳态误差e

23、p 0。动态性能分析:从局部放大图1-20中易知，超调量秒(注意：为便于对于观察，本实验中单位阶跃信号在 求。空5% ,调节时间ts 0.51t=0.5s时刻输入)，系统满足设计要(2)跟踪单位斜坡信号的动态仿真在t=0时刻输入单位斜坡信号，对于原系统状态变量初始值 初始值w(0) 2 2的仿真结果如下：x(0) 0,降维观测器状态05I1522.5时间t (s)图1-21采用降维观测器的单位斜坡响应11 II 05 1 匚茹09 口价as a?507O.S? m 136 。的 I 11)2 I K I 061C6时间tG图1-22采用降维观测器的单位斜坡响应局部放大图由图t=1s的时刻知，跟

24、踪单位斜坡信号的稳态误差ev1 0.80.2,故系统满足设计要求。可通过观察整个系统的状态变量随时间的变化规律，进一步得出两种输入信号作用下，系统的估值误差都具有很好的收敛性，从而保证了由降维观测器构成的状态反馈系统具有很好的稳定性。题目二:设被控系统状态空间表达式为：0 100x 0 01 x 0 u0251y 1 0 0 x设计使系统阶跃响应具有良好的动、静态特性的线性二次型最优控制律，对闭环系统的阶跃响应进行仿真，并研究系统二次型性能指标泛函中权矩阵的不同选取对动态性能的影 响。一、线性二次型最优控制律设计1、判断系统是否具有最优控制律要使系统阶跃响应具有良好的动、静态特性，则可按非零给

25、定点的最优控制律设计，即1u (t) Kx(t) Ug ,由于y的维数等于u的维数，所以u (t) Kx(t)也(0)yw。001系统能控性判别矩阵 Qc B AB A2B0 15c 15 23则rankQc 3 n ,所以系统完全能控，最优控制u (t)存在。q1100选取性能泛函中的对称正定矩阵Q 0 q22 0 , R 1。00 q332、非零给定点的最优控制设计及仿真分析利用SIMULINK建立原控制系统的动态仿真模型如下：图2-3非零给定点最优调节系统在t=0.5s时刻，输入单位阶跃信号，仿真分析结果如下图所示:图2-4非零给定点最优调节系统单位阶跃响应由图可知，通过引入最优控制，系

26、统获得了较好的动静态性能。3、权矩阵的各权值对动态性能的影响分析按2中方法，在t=0.5s时刻输入单位阶跃信号，改变变换权矩阵Q的值，分析其对系统动态性能的影响。(1)固定 q22 q33 1 ,使 qii分别取 1, 10,100,1000。仿真结果分别如下所示：(a)qii1eLh "rrqbll 年(b) qn 10(c) qn 100- - - - - = «_-:1 ' - E(d)q11 1000图2-5 q22 q33 1, q11取不同值时系统单位阶跃响应由图2-5可知，随着q11的增大，调节时间变短；到 qn 100以后，继续增大q11,由 于超调量的增大，调节时间反而增大。(2)固定 q