利用MATLAB实现极点配置、设计状态观测器(现代控制)

《利用MATLAB实现极点配置、设计状态观测器(现代控制)》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、装订线实验报告实验名称利用MATLAB实现极点配置、设计状态观测器系专业自动化班姓名学号授课老师预定时间实验时间实验台号一、目的要求1、掌握状态反馈和输出反馈的概念及性质。2、掌握利用状态反馈进行极点配置的方法。学会用MATLAB求解状态反馈矩阵。3、掌握状态观测器的设计方法。学会用MATLAB设计状态观测器。4、熟悉分离定理，学会设计带有状态观测器的状态反馈系统。二、原理简述1、状态反馈和输出反馈设线性定常系统的状态空间表达式为如果采用状态反馈控制规律u=r-Kx，其中r是参考输入，则状态反馈闭环系统的传递函数为

2、：2、极点配置如果SISO线性定常系统完全能控，则可通过适当的状态反馈,将闭环系统极点配置到任意期望的位置。MATLAB提供的函数acker()是用Ackermann公式求解状态反馈阵K。该函数的调用格式为K=acker(A,B,P)其中A和B分别为系统矩阵和输入矩阵。P是期望极点构成的向量。MATLAB提供的函数place()也可求出状态反馈阵K。该函数的调用格式为K=place(A,B,P)函数place()还适用于多变量系统极点配置，但不适用含有多重期望极点的问题。函数acker()不适用于多变量系统极点配置

3、问题，但适用于含有多重期望极点问题。三、仪器设备PC计算机，MATLAB软件13装订线四、内容步骤、数据处理题5-1　某系统状态方程如下理想闭环系统的极点为，试（1）采用直接计算法进行闭环系统极点配置；（2）采用Ackermann公式计算法进行闭环系统极点配置；（3）采用调用place函数法进行闭环系统极点配置。>>A=[010;001;-4-3-2];B=[1;3;6];C=[100];D=0;G=ss(A,B,C,D);[Q,D]=eig(A)结果：Q=0.29950.3197-0.0731i0.3197+0.

4、0731i-0.49440.0573+0.5074i0.0573-0.5074i0.8160-0.7948-0.7948D=-1.6506000-0.1747+1.5469i000-0.1747-1.5469i则矩阵A的特征根为：-1.6506，-0.1747+1.5469i，-0.1747-1.5469i13装订线程序：>>A=[010;001;-4-3-2];B=[1;3;-6];C=[100];D=0;p=[-1-2-3];k=acker(A,B,p)结果：k=1.48090.7481-0.0458验证：>>

5、A=[010;001;-4-3-2];B=[1;3;6];C=[100];D=0;k=[0.50530.70520.2299];A1=A-B*k;sys=ss(A1,B,C,D);G1=zpk(sys)结果：Zero/pole/gain:(s^2+5s+15)-------------------------(s+1)(s+1.999)(s+3.001)则其极点为-1,-2,-3(2)13装订线程序：>>A=[010;001;-4-3-2];B=[1;3;-6];C=[100];D=0;p=[-1-2-3];k=p

6、lace(A,B,p)结果：k=1.48090.7481-0.0458验证：程序：>>A=[010;001;-4-3-2];B=[1;3;-6];C=[100];D=0;p=[-1-2-3];k=place(A,B,p);A1=A-B*k;sys=ss(A1,B,C,D);G1=zpk(sys)结果：Zero/pole/gain:(s+4.303)(s+0.6972)--------------------13装订线(s+3)(s+2)(s+1)则其极点为-1，-2，-3题5-2某控制系统的状态方程描述如下：y=[

7、172424]x通过状态反馈使系统的闭环极点配置在[-30-1.2]位置上，求出状态反馈矩阵K，验证闭环极点，并绘制闭环系统的阶跃响应曲线进行对比。程序>>A=[-10-35-50-24;1000;0100;0010];B=[1;0;0;0];C=[172424];D=0;p=[-30-1.2-2.4+4i-2.4-4i];k=place(A,B,p)结果：k=26.0000172.5200801.7120759.3600验证：程序：>>A=[-10-35-50-24;1000;0100;0010];B=[1;0;

8、0;0];13装订线C=[172424];D=0;p=[-30-1.2-2.4+4i-2.4-4i];k=place(A,B,p);A1=A-B*k;sys=ss(A1,B,C,D);G1=zpk(sys)结果：Zero/pole/gain:(s+1.539)(s^2+5.461s+15.6)---------------------------------