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

《带状态观测器的控制系统综合设计与仿真》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、浙江大学宁波理工学院现代控制理论MATLAB仿真大作业报告题目带状态观测器的控制系统综合设计与仿真项目成员专业班级指导教师何小其学院信息科学与工程学院完成日期2015年6月18日目录摘要31主要技术参数31.1某一DC电机控制系统31.2性能指标要求42设计思路43状态空间描述53.1选定的状态变量建立系统的状态空间数学模型53.2使用Matlab得到状态空间表达式64对原系统仿真并比较性能指标65根据性能指标确定系统一组期望极点76通过状态反馈法对系统进行极点配置96.1引入状态负反馈K96.2验证状态负反馈系统的稳定性106.3使用Mat

2、lab程序求矩阵K117合理增加比例增益，使系统满足稳态指标127.1放大系数改变后系统动态性校验127.2控制系统阶跃响应指标138设计全维观测器148.1判断观测器的能观性：148.2计算观测器的反馈矩阵L158.3得到观测器的状态方程178.4对所得到的状态方程进行仿真验证178.5用Matlab求解矩阵L189在simulink下对经综合后的系统进行仿真分析1910课程设计心得体会22参考文献：23带状态观测器的控制系统综合设计与仿真摘要：状态重构器是根据系统的外部输入和输出变量的实测值，得出状态变量估计值的一类动态系统。60年代初期，

3、为了对控制系统实现状态反馈或其他需要，D.G.吕恩伯格、R.W.巴斯和J.E.贝特朗等人提出状态观测器的概念和构造方法，通过重构的途径解决了状态的不能直接量测的问题。状态观测器的出现，不但为状态反馈的技术实现提供了实际可能性，而且在控制工程的许多方面也得到了实际应用，例如复制扰动以实现对扰动的完全补偿等。关键字：系统，状态空间，matlab，稳定性，反馈，矩阵，增益，指标，仿真1主要技术参数1.1某一DC电机控制系统图1受控系统方框图1.2性能指标要求1.2.1动态性能指标超调量；超调时间秒；系统频宽；1.2.2稳态性能指标静态位置误差（阶跃

4、信号）；静态速度误差（数字信号）；2设计思路⑴按图中选定的状态变量建立系统的状态空间数学模型；⑵对原系统在simulink下进行仿真分析，对所得的性能指标与要求的性能指标进行比较；⑶根据要求的性能指标确定系统综合的一组期望极点；⑷假定系统状态均不可测,通过设计系统的全维状态观测器进行系统状态重构；⑸通过状态反馈法对系统进行极点配置，使系统满足要求的动态性能指标；⑹合理增加比例增益，使系统满足要求的稳态性能指标；⑺在simulink下对经综合后的系统进行仿真分析，验证是否达到要求的性能指标的要求。3状态空间描述3.1选定的状态变量建立系统的状态空

5、间数学模型由选定的电机控制系统要求可以写出如下关系式：由上方程可得：即拉式反变换为输出由图可知为则传递函数的状态空间表达式可写为：3.2使用Matlab得到状态空间表达式在Matlab中输入如下语句也得到状态空间表达式k=50;z=[];p=[-5-100];sys=zpk(z,p,k);G1=ss(sys)运行程序可以得到状态变量的空间数学模型G1=4对原系统仿真并比较性能指标原受控系统仿真图如图2所示：图2原受控系统仿真图原受控系统的阶跃响应如图3所示：图3原受控系统的阶跃响应曲线很显然，原系统是不稳定的。5根据性能指标确定系统一组期望极点

6、由于原系统为三阶系统，应该有三组期望极点，为了计算的方便引入两个共轭的主导极点S1、S2和一个远极点S3。由系统要求的性能指标：超调量，超调时间秒，系统频宽。可以计算求得着三个期望极点，具体过程如下。由二阶系统的各项性能指标公式式中，和为此二阶系统的阻尼比和自振频率。可以求得：⑴由，可得，从而有，于是选。⑵由，得⑶由和已选的得，与⑵的结果比较。可以确定=9.8。这样，便定出了主导极点远极点的实部应为主极点的实部的5倍以上，故选取S3=100。6通过状态反馈法对系统进行极点配置6.1引入状态负反馈K已知能控性判别矩阵为：则由上式知，因为满秩，原系

7、统是完全能控的。受控系统的特征多项式为:受控系统期望的特征多项式为：于是矩阵为:非奇异变换矩阵为：非奇异变换矩阵为：于是状态反馈矩阵为：6.2验证状态负反馈系统的稳定性在原来的开环系统中加入状态反馈可以改变系统的动态性能，状态反馈环节的添加如下图4所示：图4加入状态反馈的系统结构图根据示波器显示观察的图像如图5所示图5加状态负反馈系统输出波形显然看出系统的动态指标不能达到要求，因此还应该调整系统的放大倍数K1来达到稳态性能要求。6.3使用Matlab程序求矩阵KA=[-500;10-100;010];b=[5;0;0];c=[001];pc=[

8、-6.93+6.93i,-6.93-6.93i,-100];K=acker(A,b,pc)运行结果为K=19.77208.8690192.09967合