现代控制工程大作业

1、目 录题目一：1一、系统设计及仿真分析11、确定被控系统状态空间表达式12、系统能控、能观性判别13、系统极点配置14、确定状态反馈增益矩阵25、确定输入变换线性放大器K26、验证跟踪单位斜坡输入信号的稳态误差27、利用SIMULINK建立控制系统的动态仿真模型3（1）跟踪单位阶跃信号的动态仿真分析3（2）跟踪单位斜坡信号的动态仿真分析3二、采用全维状态观测器的状态反馈系统41、配置闭环系统状态观测器极点42、确定观测器偏差反馈增益矩阵G43、利用SIMULINK建立系统动态仿真模型5（1）对单位阶跃输入信号仿真分析5（2）对单位斜坡输入信号仿真分析64、与采用直接状态反馈的控制系统仿真结果比

2、较7（1）且时的仿真结果比较7（2）且时的仿真结果比较75、单位阶跃输入作用下状态估值误差收敛性分析8（1）收敛过程分析8（2）收敛速度与状态观测器极点的关系及对系统性能的影响9三、实现状态反馈的降维观测器设计101、通过线性变换使状态向量按输出可检测性分解102、确定降维观测器的反馈增益矩阵113、列写降维观测器状态方程并重构原系统状态114、利用SIMULINK建立系统动态仿真模型12（1）跟踪单位阶跃信号的动态仿真分析12（2）跟踪单位斜坡信号的动态仿真13题目二：15一、线性二次型最优控制律设计151、判断系统是否具有最优控制律152、非零给定点的最优控制设计及仿真分析153、权矩阵的

3、各权值对动态性能的影响分析17（1）固定 ，使 分别取1，10,100,100017（2）固定 ，使 分别取1，10,100,100018（3）固定 ，使 分别取1，10,100,1000194、使系统具有良好动、静态性能的线性二次型最优控制20题目一：设被控系统的传递函数为，现在采用带有输入变换的状态反馈控制（如教材图5-15）改善系统的性能，要求控制系统满足如下性能指标：超调量；调节时间秒；跟踪阶跃输入信号的稳态误差；跟踪单位斜坡输入信号的稳态误差。(1)根据控制系统综合指标设计状态反馈增益矩阵及输入变换线性放大器，并在SIMULINK中建立控制系统的动态仿真模型，对系统的动态和静态性能进

4、行仿真分析。(2)根据控制系统的响应速度，合理确定状态观测器系统矩阵期望特征值的分布，设计对状态向量进行估计的全维状态观测器以实现状态反馈，并保证控制系统的性能仍满足要求；在SIMULINK中建立采用状态观测器实现状态反馈的控制系统的动态仿真模型，对系统的动态和静态性能进行仿真分析，并与采用直接状态反馈的控制系统仿真结果比较；观察单位阶跃输入作用下状态估计值与状态真实值的误差的收敛过程，并研究状态估计误差收敛速度与状态观测器极点的关系及其对系统性能的影响。(3)设系统输出量可准确测量，试设计实现状态反馈的降维观测器，在SIMULINK中建立采用降维状态观测器实现状态反馈的控制系统的动态仿真模型

5、，对系统的动态和静态性能进行仿真分析。一、系统设计及仿真分析1、确定被控系统状态空间表达式被控系统传递函数为则其能观型状态空间表达式为2、系统能控、能观性判别系统能控性判别矩阵则，所以系统完全能控；又因为其状态空间表达式为能观标准型，所以系统完全能观。综上所述，此系统既能控又能观。3、系统极点配置根据系统性能指标要求：，由，可推出，取；由，得，取。根据经典控制理论，二阶系统闭环传递函数标准形式为，将，代入式，得期望的闭环主导极点为；选择一个期望的闭环非主导极点离虚轴为主导极点的5倍以上，即。则期望的闭环特征多项式为4、确定状态反馈增益矩阵由上已知，系统完全能控，所以可以通过状态反馈任意配置系统

6、闭环极点。应用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)求得状态反馈增益矩阵5、确定输入变换线性放大器K根据题目要求，跟踪单位阶跃参考输入信号的稳态误差，求得，则带有输入变换的状态反馈系统传递函数6、验证跟踪单位斜坡输入信号的稳态误差将带有输入变换的状态反馈系统等效为具有单位负反馈的闭环系统，系统方框图如下所示：图1-1 等效单位负反馈闭环系统方框图则系统传递函数求得，对于单位斜坡输入信号，；则所以，所设计系统能满足要

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

8、单位斜坡响应局部放大图由图t=1s的时刻知，跟踪单位斜坡信号的稳态误差，故系统满足设计要求。二、采用全维状态观测器的状态反馈系统基于复合系统特征值的分离特性，只要被控系统能控能观，则用状态观测器估值形成状态反馈时，可对系统的状态反馈控制器及状态观测器分别按各自的要求进行独立设计，由（一）中已求得。1、配置闭环系统状态观测器极点系统完全能观，所以闭环系统状态观测器极点可任意配置。按通常选择观测器的响应速度比所考虑的状态反馈闭环系统快25倍这一选择原则，取观测器期望极点为,2、确定观测器偏差反馈增益矩阵G应用MATLAB极点配置函数求解观测器偏差反馈增益矩阵G，代码如下：A=0 0 0;1 0 -

9、72;0 1 -18,C=0 0 1;P=-19;-19;-116;Gt=acker(A,C,P);G=Gt求得偏差反馈增益矩阵3、利用SIMULINK建立系统动态仿真模型观测器的状态方程为：利用SIMULINK建立系统动态仿真模型如下：图1-7 采用全维状态观测器的状态反馈系统（注：Scope1可观测误差的收敛过程）（1）对单位阶跃输入信号仿真分析在t=0.5s时刻输入单位阶跃信号，对于和的情况下仿真结果如下图所示：图1-8 且时单位阶跃响应图1-9 且时单位阶跃响应局部放大图静态性能分析：由图1-8可知，系统满足跟踪阶跃输入信号的稳态误差。动态性能分析：从局部放大图1-9中易知，超调量，调

10、节时间秒（注意：为便于观察，本实验中单位阶跃信号在t=0.5s时刻输入），系统满足设计要求。（2）对单位斜坡输入信号仿真分析在t=0时刻输入单位斜坡信号，对于各状态变量初值且情况下的仿真结果如下图所示：图1-10且时单位斜坡信号响应图1-11且时单位斜坡信号响应局部放大图由图t=1s的时刻知，跟踪单位斜坡信号的稳态误差，故系统满足设计要求。4、与采用直接状态反馈的控制系统仿真结果比较（1）且时的仿真结果比较（a）直接状态反馈的单位阶跃响应 （b）采用状态观测器的单位阶跃响应（c）直接状态反馈的单位斜坡响应 （d）采用状态观测器的单位斜坡响应图1-12 且的仿真结果比较从以上仿真结果图可以看出，

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

12、两种情况下的系统进行收敛性仿真分析，结果如下所示：图1-14 时误差收敛情况 图1-15 ，时误差收敛情况由图可知，当初始值相同，即时，状态估值没有误差；当初始值不同，即，时，开始一段时期存在估值误差，但是，随着时间的推移，即估值误差收敛到0。从理论上讲，状态观测器之所以出现估值误差是因为观测器与系统初始状态在实际工程中难以做到完全一致所造成的。估值误差为根据上式，当初始时刻，时，显然根据上式；当，时，由于我们通过选择输出偏差反馈增益矩阵使的所有特征值均位于复平面的左半开平面，尽管初始时刻时，存在差异，但观测器的状态仍将以一定精度和速度渐进逼近系统的实际状态。因此，仿真结果与理论分析相符。（2

13、）收敛速度与状态观测器极点的关系及对系统性能的影响现重新选取全维观测器期望极点为,应用MATLAB极点配置函数求解观测器偏差反馈增益矩阵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求得偏差反馈增益矩阵将代入采用全维状态观测器的状态反馈系统中，并在SIMULINK中进行动态仿真，当， 时，各状态变量的收敛情况如下图所示：图1-16 时的各状态变量的收敛情况 由（1）中已知情况下的误差收敛情况如下图所示： 图1-17 时的各状态变量的收敛情况从图1-16和图1-17中各状态变量的收敛情况可以看出

14、观测器极点取不同值时，的误差收敛情况。由图可知，全维观测器闭环极点离虚轴的距离越远，状态估值误差收敛速度越快。但是，同时也使得反馈增益矩阵中的元素需要很大，从而导致物理实现困难。理论上讲，估值误差收敛速度越快，系统性能越好，但是，在实际工程应用中，观测器的响应速度太快会出现放大量测噪声，使系统无法正常工作。因此，观测器期望极点的选择应兼顾快速性、抗干扰性等折中考虑，通常选择观测器的响应速度比所考虑的状态反馈闭环系统快25倍。三、实现状态反馈的降维观测器设计系统为能观标准型，所以完全能观，又，故可构造维状态观测器。1、通过线性变换使状态向量按输出可检测性分解构造非奇异变换矩阵则作线性变换，则将变

15、换成，即由于，故只需设计二维观测器重构。将分块，得，2、确定降维观测器的反馈增益矩阵配置观测器极点为则期望多项式比较和各项系数，联立方程并解得3、列写降维观测器状态方程并重构原系统状态则所对应的状态向量的估值为将变换为原系统状态空间，得到原系统的状态重构为4、利用SIMULINK建立系统动态仿真模型图1-18 采用降维状态观测器的状态反馈系统（1）跟踪单位阶跃信号的动态仿真分析在t=0.5s时刻输入单位阶跃信号，对于原系统状态变量初始值，降维观测器状态初始值的仿真结果如下：图1-19 采用降维观测器的单位阶跃响应图1-20 采用降维观测器的单位阶跃响应局部放大图静态性能分析：由图可知，系统满足

16、跟踪阶跃输入信号的稳态误差。动态性能分析：从局部放大图1-20中易知，超调量，调节时间秒（注意：为便于对于观察，本实验中单位阶跃信号在t=0.5s时刻输入），系统满足设计要求。（2）跟踪单位斜坡信号的动态仿真在t=0时刻输入单位斜坡信号，对于原系统状态变量初始值，降维观测器状态初始值的仿真结果如下：图1-21 采用降维观测器的单位斜坡响应图1-22 采用降维观测器的单位斜坡响应局部放大图由图t=1s的时刻知，跟踪单位斜坡信号的稳态误差，故系统满足设计要求。可通过观察整个系统的状态变量随时间的变化规律，进一步得出两种输入信号作用下，系统的估值误差都具有很好的收敛性，从而保证了由降维观测器构成的状

17、态反馈系统具有很好的稳定性。题目二：设被控系统状态空间表达式为：设计使系统阶跃响应具有良好的动、静态特性的线性二次型最优控制律，对闭环系统的阶跃响应进行仿真，并研究系统二次型性能指标泛函中权矩阵的不同选取对动态性能的影响。一、线性二次型最优控制律设计1、判断系统是否具有最优控制律要使系统阶跃响应具有良好的动、静态特性，则可按非零给定点的最优控制律设计，即，由于的维数等于的维数，所以。系统能控性判别矩阵则，所以系统完全能控，最优控制存在。选取性能泛函中的对称正定矩阵，。2、非零给定点的最优控制设计及仿真分析利用SIMULINK建立原控制系统的动态仿真模型如下：图2-1 原控制系统在t=0.5s时

18、刻，输入单位阶跃信号，仿真分析结果如下图所示：图2-2 原系统单位阶跃响应由图2-2可知，系统输出响应发散,引入最优控制：选取，应用MATLAB程序求解，代码如下：A=0 1 0;0 0 1;0 -2 -5;B=0;0;1;C=1 0 0;Q=100 0 0;0 1 0;0 0 1;R=1;K=lqr(A,B,Q,R);W=inv(-C*inv(A-B*K)*B) %与题中wc互为倒数求得，利用SIMULINK建立非零给定点最优控制系统的动态仿真模型如下：图2-3 非零给定点最优调节系统在t=0.5s时刻，输入单位阶跃信号，仿真分析结果如下图所示：图2-4 非零给定点最优调节系统单位阶跃响应由图可知，通过引入最优控制，系统获得了较好的动静态性能。