资源描述
现代控制理论第五次上机实验报告实验六利用MATLAB计状态观测器实验目的:1、学习观测器设计算法;2、通过编程、上机调试,掌握基于观测器的输出反应控制系统设计方法。实验步骤1、基于观测器的输出反应控制系统的设计,采用 MATLAB的m-文件编程;2、在MATLAB界面下调试程序,并检查是否运行正确。实验要求1 .在运行以上例程序的根底上,考虑图6.3所示的调节器系统,试针对被控对象设计基于全阶观测器和降阶观测器的输出反应控制器。设极点配置局部希望的闭环极点是尢2 = 2 j 2寸3,希望的观测器极点是(a) 对于全阶观测器, 比=-8和 巴=-8;(b) 对于降阶观测器,k = -8o比拟系统对以下指定初始条件的响应:(a)对于全阶观测器:为(0) =1, X2(0) =0,e(0) =1, 62(0) =0(b)对于降阶观测器:为(0) =1x(0) =0,e(0) =1进一步比拟两个系统的带宽。i(3+2)图6.3调节器系统设计闭环极点: a=0 1;0 -2;b=0;1; c=4 0;v1=-2+j*2*sqrt(3) -2-j*2*sqrt(3);K=acker(a,b,v1)16.00002.0000全阶状态观测器: v2=-8 -8;G=(acker(a-b*K),c,v2)G =30降阶状态观测器: T1 =0 1;4 0; T =0 0.25;1 0; a1 =T1*a*Tb1 =T1*b;c1 =c*T;Aaa=-2;Aab=0;Aba=4;Abb=0;Ba=1;Bb=0;v3=-8;l=(acker(Aaa,Aba,v3)Ahat=Abb-l*AabBhat=Ahat*l+Aba-l*AaaFhat=Bb-l*Baa1 =-2040l =1.5000Ahat =0Bhat =Fhat =-1.5000全阶观测器响应: sys=ss(a-b*K,0;0,eye(2),0);t=0:0.01:4;x=initial(sys,1;0,t);x1=1 0 *x;x2=0 1 *x;subplot(4,2,1);plot(t,x1),gridtitle(x1)ylabel(x1)subplot(4,2,2);plot(t,x2),grid title(x2)ylabel(x2)sys1=ss(a-G*c,0;0,eye(2),0);t=0:0.01:4;e=initial(sys1,1;0,t);e1=1 0 *x+x1;e2=0 1 *x+x2;subplot(4,2,3);plot(t,e1),grid title(全阶状态观测x1) ylabel(x1guji)subplot(4,2,4);plot(t,e2),grid title(全阶状态观测x2) ylabel(x2guji)降阶观测器响应: sys2=ss(Aaa-l*Aba,1,eye(1),0);t=0:0.01:4;e3=initial(sys2,1,t);e4=e3+x1;e5=x2;subplot(4,2,5);plot(t,e4),gridtitle(降阶状态观测x1)ylabel(x1guji)subplot(4,2,6);plot(t,e5),gridtitle(降阶状态观测x2)ylabel(x2guji)由上述两图比照可得,降阶观测器衰减速度比全阶观测器快,带宽比全阶观测器窄2.假设SISO受控系统的开环传递函数为1 G(s) s(1) 假设根据系统的性能指标要求,希望配置的系统极点为W=-3,/,3=-2-j2求受控系统的状态反应矩阵。(2) 设计观测器反应系数矩阵L ,使全维状态观测器的希望极点均为-3.状态反应矩阵: A=0 1 0;0 0 1;0 0 -1;B=0;0;1;C=1 0 0;J= -3 -2+j*2 -2-j*2;K=acker(A,B,J)24206全维观测器: v=-3 -3 -3;L=acker(A,B,v)L =27278实验心得:matlab编程仿真更这次实验做的很艰难,但一点点弄出来之后对观测器的理解也表的更加深刻,通过 加形象地了解了全阶观测器和降阶观测器的特点。
展开阅读全文