带状态观测器的控制新版系统综合设计与仿真.doc

带状态观测器控制系统综合设计与仿真一、重要技术参数：1.受控系统如图所示：U(s)X1(s)X2(s)X3(s)=Y(s)图1 受控系统方框图2.性能指标规定：（1）动态性能指标：超调量 ；超调时间 ；系统频宽 ；（2）稳态性能指标：静态位置误差（阶跃信号）静态速度误差（速度信号）二、设计思路1、按图中选定状态变量建立系统状态空间数学模型。2、对原系统在Simulink下进行仿真分析，对所得性能指标与规定性能指标进行比较。3、依照规定性能指标拟定系统综合一组盼望极点。4、假定系统状态均不可测，通过设计系统全维状态观测器进行系统状态重构。5、通过状态反馈法对系统进行极点配备，使系统满足规定动态性能指标。6、合理增长比例增益，使系统满足规定稳态性能指标。7、在Simulink下对综合后系统进行仿真分析，验证与否达到规定性能指标规定。三、实验设计环节I 、按照极点配备法拟定系统综合方案1、按图1中选定状态变量建立系统状态空间数学模型 列写每一种环节传递函数由图1有：叉乘拉式反变换得一阶微分方程组由上方程可得即拉式反变换为输出由图1可知为用向量矩阵形式表达2、对原系统在Simulink下进行仿真分析，对所得性能指标与规定性能指标进行比较原受控系统仿真图如下：图2 原受控系统仿真图原受控系统阶跃响应如下图：图3 原受控系统阶跃响应曲线很显然，原系统是不稳定。3、依照规定性能指标拟定系统综合一组盼望极点由于原系统为三阶系统，系统有3个极点，选其中一对为主导极点和，另一种为远极点，并且以为系统性能重要是由主导极点决定，远极点对系统影响很小。依照二阶系统关系式，先定出主导极点。式中，和为此二阶系统阻尼比和自振频率。可以导出：由，可得，从而有，于是选。由得由和已选得，与成果比较。这样，便定出了主导极点远极点应选取使它和原点距离远不不大于点，现取，因而拟定但愿极点为4、拟定状态反馈矩阵K由环节1所得状态空间方程知，受控系统特性多项式为而由但愿极点构成特性多项式为于是状态反馈矩阵为5、拟定放大系数L由4知，相应闭环传递函数为因此由规定跟踪阶跃信号误差，有因此对上面初步成果，再用对跟踪速度信号误差规定来验证，即显然满足规定，故。对此系统进行仿真图4 受控系统闭环系统仿真图仿真成果如下：图5 闭环系统阶跃响应曲线局部放大图：图6 闭环系统阶跃响应曲线局部放大图由仿真图得：，均满足规定。6、画出相应能控规范性闭环系统方块图已知其中，可设相应规范型状态方程为再考虑输入放大系数，最后得能控规范型闭环系统方框图如下：图7 能控规范型闭环系统方框图上述导出闭环系统方框图是相应能控规范型得到。7、拟定非奇异变换矩阵P将原受控系统传递函数方框图表达到下图形式.图8 受控系统方框图按上图选取状态变量，列状态空间方程即为依照系统能控性判据判断系统能控性则由上式知，原系统是完全能控。若做变换，那么就可建立起给定（A,B,C）和能控规范型之间关系式,。8、拟定相应于图9受控系统状态反馈矩阵K状态反馈矩阵为极点配备Matlab程序如下：A=-5 0 0;1 -10 0;0 1 0;b=1;0;0;c=0 0 1;pc=-7.07+7.07i,-7.07-7.07i,-100;K=acker(A,b,pc)运营成果为：K = 1.0e+003 *0.0991 0.4726 9.99709、画出相应于图8形式受控系统闭环方框图受控系统闭环方框图如图9示。图9 相应于图8受控系统闭环方框图仿真图形为：图10 受控系统闭环仿真图图11 闭环系统阶跃响应曲线由图可显然看出：即满足性能指标规定。II、观测器设计 假定系统状态均不可测，通过设计系统全维状态观测器进行系统状态重构1、拟定原系统能观性依照给定受控系统，求能观测性矩阵及能观测性秩则又因之前以求得系统是完全能控，因此系统即完全能控、又完全能观测。因而，系统极点可以任意配备。2、 计算观测器反馈矩阵G该设计中系统极点为取观测器极点,是观测器收敛速度是被控系统收敛速度3倍。如果仅仅对闭环极点乘以3，则阻尼比和最大超量不变，而系统上升时间和稳定期间将缩小到本来。因而，选取由所取极点，可相应闭环系统特性多项式为于是状态反馈矩阵为非奇异变换矩阵为状态反馈矩阵为因而因而观测器状态方程为3、 画出带观测器状态反馈系统闭环图带观测器状态反馈闭环系统方框图如图12所示。图12 带观测器状态反馈系统由上面计算得出带观测器状态反馈闭环系统方框图如下图13 带观测器状态反馈闭环系统方框图4、在simulink 环境下对控制系统进行仿真分析图14 带观测器状态反馈闭环系统阶跃响应曲线各状态阶跃响应曲线图15 各状态阶跃响应曲线四、参照书目1、自动控制原理 主编：李素玲 胡建 出版社：西安电子科技大学出版社2、当代控制理论 主编：王金城 出版社：化学工业出版社3、当代控制理论 主编：于长官 出版社：哈尔滨工业大学出版社4、控制系统MATLAB仿真与设计 主编：王海英 袁丽英 吴勃 出版社：高等教诲出版社5、MATLAB 7辅助控制系统设计与仿真 主编：飞思科技产品研发中心 出版社：电子工业出版社6、MATLAB控制工程工具箱技术手册 主编：魏巍 出版社：国防工业出版社7、控制系统设计与仿真 主编：赵文峰 出版社：西安电子科技大学出版社五、设计总结与心得体会不知不觉两周课程设计已经结束了，在这两周设计中，用到了所学知识涉及了自动控制原理、当代控制理论、控制系统仿真等。在设计过程中，我也懂得了必要把所学各个知识点有机结合起来，才干得出抱负成果。说实话，在最初在拿到课题时候，心里暗暗地高兴，心想这样简朴题目，几天就能解决了，谁知真正设计起来后才懂得并没有想象中简朴，每次参数选定后，按理论是可以满足设计规定，可是最后仿真分析时，不是超调量不适当就是峰值时间不能满足规定，但最后还是在通过不断地调试后选出了 合乎规定所有参数。在整个实际过程中，也不乏小小成功喜悦。在设计过程中，我以为第3步拟定盼望极点是不容易，按理论上设计极点按道理应当是满足规定，但在按所选参数完毕设计后才发现性能指标不能满足指定规定。在第四步通过状态反馈对系统进行极点配备时，按照当代控制理论上办法拟定K，但在实际设计中发现用,Matlab编程更容易实现，中间计算也就节约了大量时间。整个课程设计过程了，心中有种说不出喜悦，也许是对付出汗水承认。课程设计让我学会了学以致用，仔细想想，一学期下来，学东西还不如这两周设计。在这次设计中让我结识到做任何事情都应当认认真真，脚踏实地，积极思考，不能急于求成。附录：Matlab程序及曲线图close all clear all %The original systema=-5 0 0;1 -10 0;0 1 0;b=1;0;0;c=0 0 1;d=;cam=ctrb(a,b);rcam=rank(cam)oam=obsv(a,c);roam=rank(oam)step(a,b,c,d);hold on;grid on%The system after state feedbackpc=-100,-7.07+7.07i,-7.07-7.07i;kc=place(a,b,pc)a0=a-b*kc;k0=dcgain(a0,b,c,d);b0=b;c0=c/k0;d0=d;figure(2)step(a0,b0,c0,d0);hold on;grid on%The design of observorpo=-21,-21,-300;ko=75520 8086 327;G=koal=a-ko*c0;a2=a0 -b*kc;zeros(size(a) al;b2=b0;zeros(size(b);c2=c0 zeros(size(c);figure(3)step(a2,b2,c2,d);A=a2;onediag=eye(6);x0=1;1;1;1;1;1;K=1;B=b2;ABK=inv(A)*B*K;for t=0:0.005:1 expmat=expmdemo3(A*t); Xt=expmat*x0；%Xt=expmat*x0+(expmat-onediag)*ABK; hold on; plot(t,Xt(1),d,t,Xt(2),*,t,Xt(3),o,t,Xt(4),t,Xt(5),+,t,Xt(6),x);axis(0 1 -16 8)hold on;grid onendxlabel(times),ylabel(states vatiables)legend(y,x1(t),x2(t),x3(t),x4(t),x5(t),x6(t)附图1 原系统阶跃响应曲线附图2 经状态反馈后闭环系统阶跃响应曲线附图3 带观测器状态反馈闭环系统阶跃响应曲线附图4 带观测器状态反馈闭环系统阶跃响应曲线局部放大图