1现代控制理论实验指导书2014.3.17

《现代控制理论》实验指导书安阳工学院电子信息与电气工程学院目录实验一控制系统的数学模型及转换..................................2实验二状态空间模型的线性变换及其标准形..........................6实验三求解系统方程..............................................9实验四系统能控性、能观性的判别.................................11实验五系统稳定性仿真实验.......................................14实验六状态反馈和状态观测器的设计...............................16实验一控制系统的数学模型及转换一．实验目的（1）熟悉线性系统的数学模型及模型转换。（2）了解MATLAB中相应的函数。二．实验条件带有MATLAB的微机一台。三．实验原理（1）由传递函数建立状态空间系统的传递函数为11101110nnnnnYsbsbsbGsUssasasa（i）系统只含单实极点时的情况。设Us可分解为：12Unssss则1niiiYscUss若令状态变量为1iiXUss其向量-矩阵形式为11122201101nnnxxxxuxx，12nnxxycccx（ii）系统含重实极点时的情况。例如Ds可分解为314Unssss则131112324111niiiYsccccUsssss若令状态变量为1iiXUss111111211213113444101001101nnnxxxxxxuxxxx1112134nycccccx（2）状态方程转化为传递函数设系统的模型如式(1－1)示。pmnRyRuRxDCxyBuAxx(1－1)其中A为n×n维系数矩阵、B为n×m维输入矩阵C为p×n维输出矩阵，D为传递阵，一般情况下为0，只有n和m维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式(1－2)示。DBASICsdensnumsG1)()()()((1－2)式(1－2)中，)(snum表示传递函数阵的分子阵，其维数是p×m；)(sden表示传递函数阵的按s降幂排列的分母。四．练习内容（1）采用MATLAB编程，求系统的传递函数阵或状态空间表达式。注意：ss2tf和tf2ss是互为逆转换的指令；（2）在MATLAB界面下调试程序，并检查是否运行正确。（3）[例1.1]已知SISO系统的状态空间表达式为(1－3)，求系统的传递函数。,631234100010321321uxxxxxx321001xxxy(1－3)程序:%首先给A、B、C阵赋值；A=[010;001;-4-3-2];B=[1;3;-6];C=[100];D=0;[num,den]=ss2tf(A,B,C,D)程序运行结果:num=01.00005.00003.0000den=1.00002.00003.00004.0000从程序运行结果得到:系统的传递函数为:43235)(232sssssSG(1－4)（4）[例1.2]从系统的传递函数(1－4)式求状态空间表达式。程序:num=[0153];den=[1234];[A,B,C,D]=tf2ss(num,den)程序运行结果:A=B=-2-3-4110000100C=D=1530五．实验内容与要求（1）在运行以上[例]程序的基础上，应用MATLAB对(1－5)系统编程，求系统的A、B、C、D阵，并运行得出结果。432352)(232ssssssSG(1－5)提示：num=[0012；0153]；（2）用两种方法验证上述结果是否正确。提示：①[num,den]=ss2tf(a,b,c,d,u)；sys=tf(num,den);②DBASICsdensnumsG1)()()(()(symss;六．讨论[例1.2]程序运行结果不等于式(1－3)中的A、B、C阵，是结果错了吗？为什么？实验二状态空间模型的线性变换及其标准形一．实验目的（1）掌握线性系统的对角线标准形、约当标准形的表示及相应变换阵的求解。（2）了解MATLAB中相应的函数。二．实验条件带有MATLAB的微机一台。三．实验原理（1）对角规范型设A阵为任意形式的方阵，且有n个互异实数特征值12,,,n，则可由非奇异线性变换化为对角阵，121nPAPP阵由A阵的实数特征向量1,2,,ipin组成：12nPppp特征向量满足iiiApp，1,2,,in程序实现：[P,D]=eig(A)%P为变换阵，D为对角阵（2）约当标准形设A阵具有m重实特征值1，其余为nm个互异实特征值，但在求解iiiApp时只有一个独立实特征向量1p，则只能使A化为约当阵J，111111010mnJPAP121mmnPppppp其中12,,,mppp是广义特征向量，满足11121211,,,,,,1mmpppAppp其中1,,mnpp是互异特征值对应的实特征向量。[P,J]=jordan(A)%P为变换阵，J为约旦阵（3）P变换若已知变换矩阵P，则APPA1,BPB1,CPC四．练习内容（1）输入状态空间模型uxx1006116100010，试做线性变换，要求变换后系统矩阵A为对角阵。A=[010;001;-6-11-6];B=[0;0;1];[P,D]=eig(A);Q=inv(P);A1=Q*A*P;B1=Q*B;（2）试将矩阵452100010A化为约旦标准型，并求出变换阵。A=[010;001;2-54];[P,J]=jordan(A);五．实验内容与要求编写程序运算以下两题，并运行得出结果。（1）输入状态空间模型01161166115A，121B，100011C，21D求A的特征多项式、特征值，A的对角或约当标准形，变换矩阵P。（2）输入状态空间模型0123A，01B，10C，0D输入变换矩阵3002P，求经P变换的模型。六．讨论对于一个系统，已知状态空间模型，如何判断变换之后的系统为对角线标准型，还是约旦标准型？能否统一用jordan语句来转换？实验三求解系统方程一．实验目的（1）掌握状态转移矩阵的求法。（2）掌握系统方程的求解方法。（3）了解MATLAB中相应的函数。二．实验条件带有MATLAB的微机一台。三．实验原理（1）状态转移矩阵的计算方法级数展开法022!1!1!21kkkkkAttAktAktAAtIe拉普拉氏变换11)(AsILeAt（2）求解系统方程线性齐次方程的解)0()(xetxAt非齐次方程的解dBuexetxttAAt)()0()(0)(三、练习内容试求矩阵4321A的状态转移矩阵Ate。symst;A=[12;34];eAt=expm(A*t)四、实验内容与要求（1）输入矩阵0123A求状态转移矩阵,并计算ｔ＝0.3时的状态转移矩阵的值。（2）已知线性系统状态方程为uxx103210，01)0(x，)(1)(ttu求系统状态方程的解。五、讨论对于0123A，如何用公式11)(AsILeAt求状态转移矩阵，写出程序并运行，得出结果。实验四系统能控性、能观性的判别一、实验目的（1）系统的能控性和能观测性的判别方法、系统的能控性和能观测性分解。（2）了解MATLAB中相应的函数。二、实验条件带有MATLAB的微机一台。三、实验原理（1）能控性判据线性定常连续系统完全能控的充分必要条件：1nrankBABABn，其中n为矩阵A的维数。（2）能观测性判据线性定常连续系统完全能观测的充分必要条件：1nCCAranknCA，其中n为矩阵A的维数。（3）系统的能控性分解不能控系统的动态方程：xAxBu，yCx变换为下列的规范表达式1111212200ccccccxxxAABPAPPBuuxxxA112ccccxxyCPCCxx其中cx为r维能控状态子向量；cx为()nr维不能控状态子向量令12yyy，则可得到子系统的动态方程，其中能控子系统动态方程11121cccxAxAxBu，11yCx不能控子系统动态方程为22ccxAx，22cyCx（4）系统的能观测性分解不能观测系统的动态方程：xAxBu，yCx变换为下列能观测分解的规范表达式111122122ˆˆ0ˆˆˆooooooxxxABTATTBuuxxxBAA11ˆ0ooooxxyCTCxx能观测子系统动态方程为111ˆˆooxAxBu，11ˆoyCxy不能观测子系统动态方程为21222ˆˆˆoooxAxAxBu，10y四．练习内容已知系统uxx0136101101600xy]100[（1）判别系统的能控性。程序如下：A=[00-6;10-11;01-6];B=[3;1;0];Qc=ctrb(A,B);%Qc为能控性矩阵n=rank(Qc);%求能控性矩阵的秩L=length(A);ifn==Lstr='系统状态完全能控'elsestr='系统状态不完全能控'end（2）求系统的能控性分解后的模型。程序如下：A=[00-6;10-11;01-6];B=[3;1;0];C=[001];[111,,,,ABCTK]=ctrbf(A,B,C)%T为变换矩阵，sum（K）%可求出能控状分量的个数五．实验内容与要求调试完所有实验内容后，输入状态空间模型001103013A，110B，012C，0D（1）判别系统的能观性。提示：Qo=obsv(A,C)（2）求系统的能观性分解后的模型。提示：[111,,,,ABCTK]=obsvf(A,B,C)六．讨论自己构造一个完全能控的系统，用编程的方法将其变换为能控标准型。实验五系统稳定性仿真实验一．实验目的（1）掌握线性系统稳定性的判别方法（2）了解MATLAB中相应的函数二．实验条件带有MATLAB的微机一台。三．实验原理（1）线性定常系统为渐进稳定的充要条件是：对给定的任一个正定对称阵Q，都存在唯一的对称正定阵P，满足Lyapunov方程：TAPPAQ。（2）线性系统的稳定性的充分必要条件是：它的微分方程的全部