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

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

2、SO线性定常系统完全能控， 则可通过适当的状态反馈，将闭环系统极点配置到任意期望的位 置。MATLAB提供的函数acker()是用Ackermann公式求解状态反馈阵 K。该函数的调用格式为K=acker(A,B,P)其中A和B分别为系统矩阵和输入矩阵。P是期望极点构成的向量。MATLAB提供的函数place()也可求出状态反馈阵K。该函数的调用格式为K=place(A,B,P)函数place()还适用于多变量系统极点配置，但不适用含有多重期望极点的问题。函数acker()不适用于多变量系统极点配置问题，但适用于含有多重期望极点问题。二、仪器设备PC计算机，MATLAB软件四、内容步骤、数据处

3、理题5-1某系统状态方程如下13 u6理想闭环系统的极点为123，试(1) 采用直接计算法进行闭环系统极点配置；(2) 采用Ackermann公式计算法进行闭环系统极点配置;(3) 采用调用place函数法进行闭环系统极点配置。 A=0 1 0;0 0 1;-4 -3 -2;B=1；3；6;C=1 0 0;D=0;G=ss(A,B,C,D);Q,D=eig(A)结果：Q =- + -0 + 00 0 -则矩阵A的特征根为：，+ ,- 程序： A=0 1 0;0 0 1;-4 -3 -2;B=1；3； -6;C=1 0 0;D=0;p=-1 -2 -3;k=acker(A,B,p)结果：k =验

4、证： A=0 1 0;0 0 1;-4 -3 -2;B=1；3；6;C=1 0 0;D=0;k=;A仁 A-B*k;sys=ss(A1,B,C,D);G仁 zpk(sys)结果：Zero/pole/ga in:(sA2 + 5s + 15)(s+1) (s+ (s+则其极点为-1 ,-2 ,-3程序: A=0 1 0;0 0 1;-4 -3 -2;B=1；3； -6;C=1 0 0;D=0; p=-1 -2 -3; k=place(A,B,p) 结果：k = 验证： 程序： A=0 1 0;0 0 1;-4 -3 -2;B=1;3; -6;C=1 0 0;D=0;p=-1 -2 -3; k=p

5、lace(A,B,p); A仁 A-B*k; sys=ss(A1,B,C,D);G仁 zpk(sys)结果：Zero/pole/ga in:(s+ (s+(s+3) (s+2) (s+1)则其极点为-1 , -2 , -3题5-2某控制系统的状态方程描述如下:10355024110000xxu0100000100y=1 7 24 24x通过状态反馈使系统的闭环极点配置在-302.4 4j 位置上，求出状态反馈矩阵 K,验证闭环极点，并绘制闭环系统的阶跃响应曲线进行对比。程序 A=-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0;B=1；0；0;0;C=1 7 24

6、 24;D=0;p=-30 +4i ;k=place(A,B,p)结果：k =验证：程序： A=-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0;B=1;0;0;0;C=1 7 24 24;D=0;p=-30 +4i ;k=place(A,B,p);A仁 A-B*k;sys=ss(A1,B,C,D);G仁 zpk(sys)结果：Zero/pole/ga in:(s+ (sA2 + +(s+30) (s+ (sA2 + + A=-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0;B=1；0；0；0；C=1 7 24 24;D=0;G=ss

7、(A,B,C,D);p=-30 +4i ;k=place(A,B,p);A1=A-B*k;G仁 ss(A1,B,C,D);t=0:20;u=on es(size(t);y2=lsim(G1,u,t);y1=lsim(G,u,t);plot(t,y1,:,t,y2,-)Q9a.eazOLE0.50.4.3J.10蓝色为配置前，绿色为配置后题5-3某系统状态空间描述如下0 10 1X 001 x 3 u4326y 1 0 0 x设计全维状态观测器，要求状态观测器的极点为123。程序 A=0 1 0;0 0 1;-4 -3 -2;B=1；3；-6：C=1 0 0;D=0;P=-1 -2 -3;L=(

8、acker(A,C,p)结果：L = 40-10题5-4已知系统0100x001 x0 u12y=6 6 0x1671(1)求系统的零点，极点和传递函数，并判断系统的能控性和能观测性。 分别选取K=0 3 0.K=1 3 2,K=0 3 1为状态反馈矩阵，求解闭环系统的零点,极点和传递函数，并判断系统的能控性和能观测性。它们是否发生改变为什么(1)程序 A=0 1 0;0 0 1;-12 -16 -7;B=0；0；1；C=6 6 0;D=0;sys=ss(A,B,C,D);G=zpk(sys)结果：Zero/pole/ga in:6 (s+1)(s+2)A2 (s+3)则系统的零点为：-1则系

9、统的极点为：-2，-2，-3程序 A=0 1 0;0 0 1;-12 -16 -7;B=0;0;1;D=0;Uc=ctrb(A,B);ra nk(Uc)结果：ans = 3则系统能控程序： A=0 1 0;0 0 1;-12 -16 -7;B=0；0；1；C=6 6 0;D=0;Vo=obsv(A,C);ran k(Vo)结果：ans =3则系统能观(2) 当k=0 3 0 时：程序 A=0 1 0;0 0 1;-12 -16 -7;B=0；0；1；C=6 6 0;D=0;k=0 3 0;A仁 A-B*k;sys=ss(A1,B,C,D);G=zpk(sys)结果：Zero/pole/ga i

10、n:(s+ 理2 + +则零点为-1 pole(G) ans =+ Uc=ctrb(A,B);ra nk(Uc) ans =3 Vo=obsv(A,C);ran k(Vo) ans =3 当k=1 3 2 时：程序 A=0 1 0;0 0 1;-12 -16 -7;B=0；0；1；D=0;k=1 3 2;A仁 A-B*k; sys=ss(A1,B,C,D);G=zpk(sys) 结果：Zero/pole/ga in:6 (s+1)(s+ (sA2 + +则闭环系统的零点为：-1 pole(G)ans =+ Uc=ctrb(A,B);ra nk(Uc)ans =3则系统具有能控性 Vo=obsv

11、(A,C);ran k(Vo)ans =3则系统具有能观测性 当k=0 1 3 时： A=0 1 0;0 0 1;-12 -16 -7;B=0；0；1；C=6 6 0;D=0;k=0 1 3;A仁 A-B*k; sys=ss(A1,B,C,D);G=zpk(sys)Zero/pole/ga in:6 (s+1)(s+ (sA2 + + Uc=ctrb(A,B);ra nk(Uc)ans =则系统具有能控性 Vo=obsv(A,C);ran k(Vo)ans =3则系统具有能观测性分析：系统完全能控则可以任意配置极点，配置极点不改变系统的能控性，但可能改变系统的能观测性。不存在零极相消的情况下，则不改变系统的能观测性。五、分析讨论通过本次试验，掌握了状态反馈和输出反馈的概念及性质。掌握了利用状态反馈进行极点配置的方法。学会了用