线性系统理论的matlab应用PPT课件

《线性系统理论的matlab应用PPT课件》由会员分享，可在线阅读，更多相关《线性系统理论的matlab应用PPT课件（50页珍藏版）》请在装配图网上搜索。

1、MATLAB的应用的应用曾水平控制系统的分析方法 早期的控制系统分析过程复杂而耗时，如想得到一个系统的冲激响应曲线，首先需要编写一个求解微分方程的子程序，然后将已经获得的系统模型输入计算机，通过计算机的运算获得冲激响应的响应数据，然后再编写一个绘图程序，将数据绘制成可供工程分析的响应曲线。MATLAB控制系统工具箱和SIMULINK辅助环境的出现，给控制系统分析带来了福音。控制系统的分析包括系统的稳定性分析、时域分析、频域分析及根轨迹分析。本课件是关于状态空间的时域分析一一 利用利用MATLAB进行模型转换进行模型转换1-1 传递函数与状态空间表达式之间的转换传递函数与状态空间表达式之间的转换

2、1）.连续系统状态空间表达式连续系统状态空间表达式 MATLAB是当今世界上最优秀的科技应用软件之一，它以强大是当今世界上最优秀的科技应用软件之一，它以强大的科学计算能力和可视化功能，简单易用的编程语言以及开放式的的科学计算能力和可视化功能，简单易用的编程语言以及开放式的编程环境等一些显著的优点，使得它在当今许许多多科学技术领域编程环境等一些显著的优点，使得它在当今许许多多科学技术领域中成为计算机辅助分析和设计、算法研究和应用开发的基本工具和中成为计算机辅助分析和设计、算法研究和应用开发的基本工具和首选平台。在本书中，用它作为系统分析和设计的软件平台，更显首选平台。在本书中，用它作为系统分析和

3、设计的软件平台，更显示出独特的优势。示出独特的优势。本节利用本节利用MATLAB实现数学模型的转换。实现数学模型的转换。可以用可以用ss命令来建立状态空间模型。对于连续系统，其格命令来建立状态空间模型。对于连续系统，其格式为式为 sys=ss(A,B,C,D)，其中，其中A，B，C，D为描述线性连续系为描述线性连续系统的矩阵。统的矩阵。当当sys1是一个用传递函数表示的线性定常系统时，可以用是一个用传递函数表示的线性定常系统时，可以用命令命令sys=ss(sys1)，将其转换成为状态空间形式。也可以用命，将其转换成为状态空间形式。也可以用命令令sys=ss(sys1,min)计算出系统计算出系

4、统sys的最小实现。的最小实现。例例1 1 控制系统微分方程为控制系统微分方程为uuuuyyyyy2424724503510)4(求其状态空间表达式。求其状态空间表达式。解解可以先将其转换成传递函数可以先将其转换成传递函数 2450351024247)()()(23423ssssssssusysG输入下列命令输入下列命令语句执行结果为语句执行结果为这个结果表示，该系统的状态空间表达式为这个结果表示，该系统的状态空间表达式为uyu01875.0375.04375.01000102000040000161875.07813.0188.210 xxx 注意，在输入命令中，注意，在输入命令中，sys=

5、ss(G)也可以改用也可以改用A,B,C,D=tf2ss(num,den)，在本例中其作用和，在本例中其作用和sys=ss(G)近近似，也可以计算出矩阵似，也可以计算出矩阵A、B、C、D。2）.离散系统的状态空间表达式（不要求）离散系统的状态空间表达式（不要求）离散系统的状态空间表达式为离散系统的状态空间表达式为)()()()()()1(kdukCxkykHukGxkx 和连续系统状态空间表达式的输入方法相类似，如果要输入和连续系统状态空间表达式的输入方法相类似，如果要输入离散系统的状态空间表达式，首先需要输入矩阵离散系统的状态空间表达式，首先需要输入矩阵G、H、C、d，然后输入语句然后输入语

6、句 ，即可将其输入到，即可将其输入到MATLAB的的workspace中，并且用变量名来表示这个离散系统，其中中，并且用变量名来表示这个离散系统，其中T为采为采样时间。如果样时间。如果Gyu表示一个以脉冲传递函数描述的离散系统，也表示一个以脉冲传递函数描述的离散系统，也可以用可以用ss(Gyu)命令，将脉冲传递函数模型转换成状态空间表达命令，将脉冲传递函数模型转换成状态空间表达式。式。),(TdCHGsssys 例例2 2 假设某离散系统的脉冲传递函数为假设某离散系统的脉冲传递函数为47.022.298.323.389.038.057.031.0)(23423zzzzzzzzGyu采样周期为采

7、样周期为 ，将其输入到，将其输入到MATLAB的的workspace中，并且中，并且绘制零、极点分布图。并且将该离散系统脉冲传递函数模型转换绘制零、极点分布图。并且将该离散系统脉冲传递函数模型转换成状态空间表达式。成状态空间表达式。sT1.0 解解 输入下列语句输入下列语句语句执行的结果为语句执行的结果为再输入语句再输入语句 ，绘制出零、极点分布图如下，绘制出零、极点分布图如下在执行完上述语句后，在执行完上述语句后，Gyu已经存在于已经存在于MATLAB的的workspace中，这时再执行语句中，这时再执行语句执行结果为执行结果为 结果表示，离散系统的状态空间表达式为结果表示，离散系统的状态空

8、间表达式为)(0001)(05.00000200002235.0555.099.123.3)1(kukkxx)(445.0095.0285.0310)(k.kyx1-2 求传递函数矩阵求传递函数矩阵 在已知线性定常系统中的在已知线性定常系统中的A、B、C和和D矩阵之后，则该矩阵之后，则该系统的传递函数矩阵可以按下式求出系统的传递函数矩阵可以按下式求出DBAsICsusysG1)()()(例例3 3 已知系统状态方程为已知系统状态方程为21212111013210uuxxxx2121321100003121112uuxxyyy输入以下语句输入以下语句 解解 其中其中inv()函数是求矩阵的逆矩阵

9、，而函数是求矩阵的逆矩阵，而simple()函数是函数是对符号运算结果进行简化。对符号运算结果进行简化。执行结果如下执行结果如下这表示这表示ssssAsI213231211)1/(1()1/(3)2/(1)2/(2)1/(13)1/(3()(sssssssG1-3.线性变换线性变换1）化为对角矩阵化为对角矩阵 函数函数eig()可以计算出矩阵可以计算出矩阵A的特征值以及将的特征值以及将A阵转换成对角阵的阵转换成对角阵的线性变换矩阵。其语句格式为线性变换矩阵。其语句格式为Q,D=eig(A)，则，则D为对角阵并且为对角阵并且对角线上各元素为矩阵对角线上各元素为矩阵A的特征值，满足的特征值，满足

10、，因为，因为 即：即：。DAQQ11 PQDPAP1例例4 4 线性控制系统的状态方程为线性控制系统的状态方程为 u1006116100010 xx 试作线性变换试作线性变换 ,要求变换后系统矩阵要求变换后系统矩阵A为对角阵。为对角阵。xQxPx1解解先求出系统矩阵的特征值，先求出系统矩阵的特征值，Q阵可以选择为由特征值构成阵可以选择为由特征值构成的范德蒙特矩阵。的范德蒙特矩阵。输入语句输入语句可以求出可以求出A阵的特征值为阵的特征值为 1、2和和3。因此因此 9413211111PQ输入以下语句输入以下语句执行结果如下执行结果如下 由以上计算数据可得系统经过线性变换后的方程为由以上计算数据可

11、得系统经过线性变换后的方程为 u5.015.0300020001xx也可以输入语句也可以输入语句运行结果为运行结果为 再计算线性变换矩阵再计算线性变换矩阵P，并且验证结果如下，并且验证结果如下 可见，两种线性变换虽然不同，却都可以将可见，两种线性变换虽然不同，却都可以将A阵转换为对角阵阵转换为对角阵2）化为约当矩阵化为约当矩阵 在在MATLAB中用函数命令中用函数命令jordan()来求矩阵的约当标准形。其来求矩阵的约当标准形。其命令格式为：命令格式为：Q,J=jordan(A)。输入参量输入参量A是系数矩阵，输出参是系数矩阵，输出参量量J是矩阵是矩阵A 的约当标准形矩阵，而的约当标准形矩阵，

12、而 就是线性变换矩阵，就是线性变换矩阵，满足满足 。1 QP11PAPQAQJ例例5 5 将将 化为标准形矩阵。化为标准形矩阵。452100010A解解首先输入语句首先输入语句 运行结果为运行结果为可见，不满秩，即矩阵可见，不满秩，即矩阵A的特征值中有重特征值的特征值中有重特征值,并且并且A的独立特的独立特征向量的个数小于征向量的个数小于n。因此输入语句因此输入语句语句执行结果为语句执行结果为计算结果表明，矩阵计算结果表明，矩阵A的约当阵为的约当阵为 。100110002J我们验证如下我们验证如下 执行结果为执行结果为所计算出的结果表明，满足所计算出的结果表明，满足 1PAPJ二二 用用MAT

13、LAB求解系统方程求解系统方程2.1 线性齐次状态方程的解线性齐次状态方程的解 使用使用MATLAB可以方便地求出状态方程的解。我们通过例可以方便地求出状态方程的解。我们通过例子来说明。子来说明。例例 已知线性系统齐次状态方程为已知线性系统齐次状态方程为 xx3210初始条件初始条件01)0(x求系统状态方程的解。求系统状态方程的解。解解用以下用以下MATLAB程序计算齐次状态方程的解，其中程序计算齐次状态方程的解，其中collect()函数的作用是合并同类项，而函数的作用是合并同类项，而ilaplace()函数的作用是求取拉函数的作用是求取拉普拉斯逆变换，函数普拉斯逆变换，函数det()的作

14、用是求方阵的行列式。的作用是求方阵的行列式。程序执行结果程序执行结果这表示这表示ttttttttt2222e2ee2e2eeee2)(ttttt22e2e2ee2)(x2.2 线性非齐次状态方程的解线性非齐次状态方程的解通过以下例子说明。通过以下例子说明。例例2 2 已知系统状态方程为已知系统状态方程为u103210 xx 01)0(x)(1)(ttu解解 用以下用以下MATLAB程序求系统方程的解。其中，语句程序求系统方程的解。其中，语句phi=subs(phi0,t,(t-tao)表示将符号变量表示将符号变量phi0中的自变量中的自变量t用用(t-tao)代换就构成了符号变量代换就构成了符

15、号变量phi，而语句，而语句x2=int(F,tao,0,t)表示符号变量表示符号变量F对对tao在在0到到t的积分区间上求积分，运算结果返回到的积分区间上求积分，运算结果返回到x2。程序执行结果为程序执行结果为 ttttt22eee5.0e5.0)(x这表示这表示2.3 连续系统状态方程的离散化（不要求）连续系统状态方程的离散化（不要求）在在MATLAB中，函数中，函数c2d（）的功能就是将连续时间的系统模型转（）的功能就是将连续时间的系统模型转换成离散时间的系统模型。其调用格式为：换成离散时间的系统模型。其调用格式为：sysd=c2d(sysc,T,method)。其中，输入参量。其中，输

16、入参量sysc为连续时间的系统为连续时间的系统模型；模型；T为采样周期（秒）；为采样周期（秒）；method用来指定离散化采用的方法用来指定离散化采用的方法。zoh采用零阶保持器；采用零阶保持器；foh采用一阶保持器；采用一阶保持器；tustin采用双线性逼近方法；采用双线性逼近方法；prewarm采用改进的采用改进的tustin方法；方法；matched采用采用SISO系统的零极点匹配方法；系统的零极点匹配方法；当当method为缺省时（即：调用格式为为缺省时（即：调用格式为sysd=c2d(sysc,T)时），默认时），默认的方法是采用零阶保持器。的方法是采用零阶保持器。例例2-32-3

17、某线性连续系统的状态方程为某线性连续系统的状态方程为BuAxxDuCxy其中其中 6116100010A201201B112011C0000D采用零阶保持器将其离散化，设采样周期为采用零阶保持器将其离散化，设采样周期为0.1秒。求离散化的状态秒。求离散化的状态方程模型。方程模型。解解 输入以下语句，其中输入以下语句，其中D=zeros(2)表示，将表示，将D赋值为赋值为22维的全零维的全零矩阵。矩阵。语句执行的结果为语句执行的结果为计算结果表示系统离散化后的计算结果表示系统离散化后的状态方程为状态方程为)(1936.01164.00902.01959.00047.01099.0)(5112.0

18、8366.04429.00738.09541.00246.00041.00984.09991.0)1(kkkuxx)(0000)(112011)(kkkuxy三三 MATLAB的应用的应用3-1 判断线性系统的能控性和能观测性判断线性系统的能控性和能观测性 用用MATLAB可以很方便地求出线性控制系统的能控性矩阵和能可以很方便地求出线性控制系统的能控性矩阵和能观测性矩阵，并且求出它们的秩。从而判断系统的能控性和能观测观测性矩阵，并且求出它们的秩。从而判断系统的能控性和能观测性。函数性。函数ctrb()和和obsv()分别计算系统的能控性矩阵和能观测性矩分别计算系统的能控性矩阵和能观测性矩阵。格

19、式为：阵。格式为：Qc=ctrb(A,B),Qo=obsv(A,C)。例例 判断下面的线性系统是否能控？是否能观测？判断下面的线性系统是否能控？是否能观测？CxBAxxyu其中其中 103021101A201201B010001C解解 先分别计算系统的能控性矩阵和能观测性矩阵。然后，再用先分别计算系统的能控性矩阵和能观测性矩阵。然后，再用rank()函数计算这两个矩阵的秩。函数计算这两个矩阵的秩。输入以下语句输入以下语句这些语句的执行结果为这些语句的执行结果为 从计算结果可以看出，系统能控性矩阵和能观测性矩阵的秩都从计算结果可以看出，系统能控性矩阵和能观测性矩阵的秩都是是3，为满秩，因此该系统

20、是能控的，也是能观测的。，为满秩，因此该系统是能控的，也是能观测的。注：当系统的模型用注：当系统的模型用sys=ss(A,B,C,D)输入以后，也就是当系统模输入以后，也就是当系统模型用状态空间的形式表示时，我们也可以用型用状态空间的形式表示时，我们也可以用Qc=ctrb(sys)，Qo=obsv(sys)的形式求出该系统的能控性矩阵和能观测性矩阵。的形式求出该系统的能控性矩阵和能观测性矩阵。3-2 线性系统按能控性或者能观测性分解线性系统按能控性或者能观测性分解在用在用MATLAB进行结构分解时，不能控（不能观）的系统，其结构进行结构分解时，不能控（不能观）的系统，其结构分解的系统方程形式有

21、点不同。分解的系统方程形式有点不同。)ctrbf(CB,A,KT,C,B,An)rank(CQ当系统能控性矩阵的秩当系统能控性矩阵的秩 时，我们可以使用函数命令时，我们可以使用函数命令ctrbf()可以对线性系统进行能控性分解。其调用格式为可以对线性系统进行能控性分解。其调用格式为 。其中，。其中，T为相似变换矩阵。为相似变换矩阵。CCAAAA210CBB0CCCCC 输出为一个向量，输出为一个向量，sum(K)可以求出能控的状态分量的个数。可以求出能控的状态分量的个数。类似地，当系统能观测性矩阵的秩类似地，当系统能观测性矩阵的秩 时，我们可以使用时，我们可以使用函数命令函数命令obsvf()

22、可以对线性系统进行能观测性分解。其调用格式可以对线性系统进行能观测性分解。其调用格式为为 。n)rank(OQ)obsvf(CB,A,KT,C,B,A其中，其中，T 为相似变换矩阵。为相似变换矩阵。O12OAAAA0输出为一个向量，输出为一个向量，sum(K)可以求出能观测的状态分量的个数。可以求出能观测的状态分量的个数。OOBBBOCC0例例 系统方程为系统方程为 CxBAxxyu其中其中 6101101600A013B100C试按能控性进行结构分解。试按能控性进行结构分解。解解 输入下列语句输入下列语句 语句执行结果为语句执行结果为从输出的向量可以看出有两个状态分量是能控的。可以验证从输出

23、的向量可以看出有两个状态分量是能控的。可以验证 ，输入语句，输入语句TTATA得到的结果为得到的结果为 可见，可见，A1=Abar，所得到的结果是正确的。，所得到的结果是正确的。3-3线性系统转换成能控标准形和能观标准形线性系统转换成能控标准形和能观标准形 下面通过两个例子来说明将系统变换成能控标准形和能观下面通过两个例子来说明将系统变换成能控标准形和能观标准形的方法。标准形的方法。例例 系统方程为系统方程为 CxBAxxyu其中其中 231120121A110B101C求线性变换，将其变换成能控标准形。求线性变换，将其变换成能控标准形。解解 1）判断系统是否能控，并且求出）判断系统是否能控，

24、并且求出A 阵的特征多项式阵的特征多项式输入下面语句输入下面语句运行结果为运行结果为表明系统为能控，因此可以变换成能控标准形。而且求出表明系统为能控，因此可以变换成能控标准形。而且求出A 的特的特征多项式为征多项式为11125det23AI（即：（即：，）110a121a52a2）计算变换矩阵）计算变换矩阵0010151512001011221CCQQQaaa1 QP输入以下语句输入以下语句计算结果为计算结果为3）计算出能控标准形）计算出能控标准形 输入以下语句输入以下语句计算结果为计算结果为表明经过变换以后的系统方程为表明经过变换以后的系统方程为u10051211100010 xxx1510

25、y例例 系统方程为系统方程为 CxBAxxyu其中其中 101325103A201B112C求线性变换，将其变换成能观标准形。求线性变换，将其变换成能观标准形。解解 1）判断系统是否为能观测，并且求出）判断系统是否为能观测，并且求出A阵的特征多项式阵的特征多项式 输入下面语句输入下面语句运行结果为运行结果为 表明系统为能观测，因此可以变换成能观标准形。而且求出的表明系统为能观测，因此可以变换成能观标准形。而且求出的特征多项式为特征多项式为4106det23AI（即：（即：，）40a101a62a2）计算变换矩阵）计算变换矩阵 OOQQP0010161610001011221aaa输入以下语句输

26、入以下语句计算结果为计算结果为3）计算出能观标准形）计算出能观标准形 输入以下语句输入以下语句计算结果为计算结果为表明经过变换以后的系统方程为表明经过变换以后的系统方程为u4046101001400 xxx100y四四 MATLAB的应用的应用4.1 极点配置极点配置 线性系统是状态能控时，可以通过状态反馈来任意配置系统的线性系统是状态能控时，可以通过状态反馈来任意配置系统的极点。把极点配置到极点。把极点配置到S左半平面所希望的位置上，则可以获得满意左半平面所希望的位置上，则可以获得满意的控制特性。的控制特性。状态反馈的系统方程为状态反馈的系统方程为 BvxBKAx)(Cxy 在在MATLAB

27、中，用函数命令中，用函数命令place()可以方便地求出状态反馈可以方便地求出状态反馈矩阵矩阵K；该命令的调用格式为：；该命令的调用格式为：K=place(A,b,P)。P为一个行向量，其各分量为所希望配置的各极为一个行向量，其各分量为所希望配置的各极点。即：该命令计算出状态反馈阵点。即：该命令计算出状态反馈阵K，使得（，使得（A-bK）的特征值为向）的特征值为向量量P的各个分量。使用函数命令的各个分量。使用函数命令acker()也可以计算出状态矩阵也可以计算出状态矩阵K，其作用和调用格式与其作用和调用格式与place()相同，只是算法有些差异。相同，只是算法有些差异。例例 线性控制系统的状态

28、方程为线性控制系统的状态方程为uBAxx其中其中 0100016116A001B要求确定状态反馈矩阵，使状态反馈系统极点配置为要求确定状态反馈矩阵，使状态反馈系统极点配置为 101s112s123s解解 首先判断系统的能控性，输入以下语句首先判断系统的能控性，输入以下语句 语句执行结果为语句执行结果为 这说明系统能控性矩阵满秩，系统能控，可以应用状态反馈，这说明系统能控性矩阵满秩，系统能控，可以应用状态反馈，任意配置极点。任意配置极点。输入以下语句输入以下语句 语句执行结果为语句执行结果为 计算结果表明，状态反馈阵为计算结果表明，状态反馈阵为 131435127K注注：如果将输入语句中的：如果

29、将输入语句中的 K=place(A,B,P)改为改为 K=acker(A,B,P)，可以得到同样的结果。，可以得到同样的结果。4.2 状态观测器设计状态观测器设计 在在MATLAB中，可以使用函数命令中，可以使用函数命令acker()计算出状态观测器计算出状态观测器矩阵矩阵。调用格式。调用格式 ，其中，其中AT 和和 CT 分别是分别是A 和和B 矩阵的转置。矩阵的转置。P为一个行向量，其各分量为所希望的状态观测器为一个行向量，其各分量为所希望的状态观测器的各极点。的各极点。GT为所求的状态观测器矩阵为所求的状态观测器矩阵G 的转置。的转置。),acker(PCAGTTT例例 线性控制系统的状

30、态方程为线性控制系统的状态方程为uBAxxCxy其中其中 200120001A101B011C要求设计系统状态观测器，其特征值为要求设计系统状态观测器，其特征值为:3,4,5。解解 首先判断系统的能观测性，输入以下语句首先判断系统的能观测性，输入以下语句 语句运行结果为语句运行结果为说明系统能观测，可以设计状态观测器说明系统能观测，可以设计状态观测器输入以下语句输入以下语句 语句运行结果为语句运行结果为 计算结果表明，状态观测器矩阵为计算结果表明，状态观测器矩阵为210103120G状态观测器的方程为状态观测器的方程为uyBGxGCAx)(uy1012101031202210210110510

31、30120119x五五 单级倒立摆系统单级倒立摆系统MATLAB仿真1 1 建立单级倒立摆系统的状态空间表达式。建立单级倒立摆系统的状态空间表达式。单级倒立摆系统是许多重要的宇宙空间应用的一个简单模型。单级倒立摆系统是许多重要的宇宙空间应用的一个简单模型。在水平方向，应用牛顿第二定律：在水平方向，应用牛顿第二定律：ulytmtyM)sin(dddd2222在垂直于摆杆方向，应用牛顿第二定律：在垂直于摆杆方向，应用牛顿第二定律：sin)sin(dd22mglytm而有：而有：)(cos)(sinddt cos)sin()(sindd222t)sin()(cosddt)sin()cos()(cos

32、dd222t1cos线性化：当线性化：当 和和 较小时较小时，有，有sin02化简后，得化简后，得umlymM )(mgmlym 求解得：求解得：uMMmgy1 uMlMlgmM1)(选择状态变量选择状态变量 ，为系统输入，为系统输入，为系统输出为系统输出yx 1yxx 123x 34xxuy;0100010000000010114321)(4321uxxxxxxxxMlMMlgmMMmg43210001xxxxy状态图为状态图为2 单级倒立摆系统的极点配置与状态观测器设计单级倒立摆系统的极点配置与状态观测器设计1.状态反馈系统的极点配置及其状态反馈系统的极点配置及其MATLAB/Simuli

33、nk仿真仿真选取适当参数，单级倒立摆系统的状态方程为选取适当参数，单级倒立摆系统的状态方程为uxxxxxxxx1010011001000010000104321432143210001xxxxy 首先，使用首先，使用MATLAB，判断系统的能控性矩阵是否为满秩。，判断系统的能控性矩阵是否为满秩。输入以下程序输入以下程序 计算结计算结果为果为 根据判别系统能控性的定理，该系统的能控性矩阵满秩，所以根据判别系统能控性的定理，该系统的能控性矩阵满秩，所以该系统是能控的。因为系统是能控的，所以，可以通过状态反馈来该系统是能控的。因为系统是能控的，所以，可以通过状态反馈来任意配置极点。任意配置极点。不失

34、一般性，不妨将极点配置在不失一般性，不妨将极点配置在 61s5.62s73s5.74s在在MATLAB中输入命令中输入命令得到计算结果为得到计算结果为因此，求出状态反馈矩阵为因此，求出状态反馈矩阵为 175.1495.488175.12275.204K 采用采用MATLAB/Simulink构造单级倒立摆状态反馈控制系统的构造单级倒立摆状态反馈控制系统的仿真模型，如下图所示。仿真模型，如下图所示。首先，在首先，在MATLAB的的Command Window中输入各个矩阵的值，中输入各个矩阵的值，并且在模型中的积分器中设置非零初值。然后运行仿真程序。并且在模型中的积分器中设置非零初值。然后运行仿

35、真程序。得到的仿真曲线得到的仿真曲线从仿真结果可以看出，可以将倒立摆的杆子与竖直方向的偏角控制从仿真结果可以看出，可以将倒立摆的杆子与竖直方向的偏角控制在在 （即小球和杆子被控制保持在竖直倒立状态）。（即小球和杆子被控制保持在竖直倒立状态）。03.状态观测器实现状态反馈极点配置及其仿真状态观测器实现状态反馈极点配置及其仿真首先，使用首先，使用MATLAB，判断系统的能观性矩阵是否为满秩。输，判断系统的能观性矩阵是否为满秩。输入以下程序入以下程序 计算结果为计算结果为 因为该系统的能观测性矩阵满秩，所以该系统是能观测的。因因为该系统的能观测性矩阵满秩，所以该系统是能观测的。因为系统是能观测的，所

36、以，可以设计状态观测器。而系统又是能控为系统是能观测的，所以，可以设计状态观测器。而系统又是能控的，因此可以通过状态观测器实现状态反馈。的，因此可以通过状态观测器实现状态反馈。234s223s212s 设计状态观测器矩阵，使的特征值的实部均为负，且其绝对值设计状态观测器矩阵，使的特征值的实部均为负，且其绝对值要大于状态反馈所配置极点的绝对值。通过仿真发现，这样才能保要大于状态反馈所配置极点的绝对值。通过仿真发现，这样才能保证状态观测器有足够快的收敛速度，才能够保证使用状态观测器所证状态观测器有足够快的收敛速度，才能够保证使用状态观测器所观测到的状态与原系统的状态充分接近。不妨取状态观测器的特征

37、观测到的状态与原系统的状态充分接近。不妨取状态观测器的特征值为：值为：201s输入以下命令输入以下命令计算结果为计算结果为求出状态观测器矩阵为求出状态观测器矩阵为 T24312040590278090G 如果采用如果采用MATLAB/Simulink构造具有状态观测器的单级倒立构造具有状态观测器的单级倒立摆状态反馈控制系统的仿真模型，如下图所示。摆状态反馈控制系统的仿真模型，如下图所示。首先，在首先，在MATLAB的的Command Window中输入各个矩阵中输入各个矩阵的值，并且在模型中的积分器中设置非零初值。然后运行仿的值，并且在模型中的积分器中设置非零初值。然后运行仿真程序。得到的仿真曲线。真程序。得到的仿真曲线。比较两个仿真结果，具有状态观测器的单级倒立摆状态反馈比较两个仿真结果，具有状态观测器的单级倒立摆状态反馈系统的控制效果和没有状态观测器的控制系统的控制效果。系统的控制效果和没有状态观测器的控制系统的控制效果。