自动控制原理MATLAB课程设计--滞后-超前校正.doc

PAGE 5 滞后-超前校正 ——课程设计 一、设计目的: 了解控制系统设计的一般方法、步骤。 掌握对系统进行稳定性的分析、稳态误差分析以及动态特性分析的方法。 掌握利用MATLAB对控制理论内容进行分析和研究的技能。 提高分析问题解决问题的能力。 二、设计内容与要求: 设计内容： 阅读有关资料。 对系统进行稳定性分析、稳态误差分析以及动态特性分析。 绘制根轨迹图、Bode图、Nyquist图。 设计校正系统，满足工作要求。 设计条件： １、被控制对象的传递函数是 （ｎｍ） 　　　２、参数ａ０，ａ１，ａ２，．．．ａｎ和ｂ０，ｂ１，ｂ２，．．．ｂｍ因小组而异。 设计要求： 能用MATLAB解复杂的自动控制理论题目。 能用MATLAB设计控制系统以满足具体的性能指标。 能灵活应用MATLAB的CONTROL SYSTEM工具箱和SIMULINK仿真软件，分析系统的性能。 三、设计步骤： １、自学ＭＡＴＬＡＢ软件的基本知识，包括ＭＡＴＬＡＢ的基本操作命令。控制系统工具箱的用法等，并上机实验。 ２、基于ＭＡＬＡＢ用频率法对系统进行串联校正设计，使其满足给定的领域性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数Ｔ，α等的值。 已知开环传递函数为G(S)= ,使用频率法设计串联滞后—超前校正装置，使系统的相角裕度大于等于40°，静态速度误差系数等于20。 校正前 根据上式可化简G(S)= ，所以公式G(S)= ,所以=1，则= 6.1310，相角裕度为9.3528。 校正后 串联校正滞后—超前校正装置, 由，取，又由 , ,;; 所以校正后的函数为 ３、利用ＭＡＴＬＡＢ函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，为什么？ 校正前 % MATLAB PROGRAM j005.m % num=[20]; den= conv([1 0],conv([0.5 1],[0.025 1])); g=tf(num,den); sys=feedback(g,1); pzmap(g); den= conv([1 0],conv([0.5 1],[0.025 1])); t=tf(num,den); pzmap(t); [p,z]=pzmap(g); den=sys.den{1}; r=roots(den); disp(r) -41.0006 -0.4997 + 6.2269i -0.4997 - 6.2269i 系统没有零极点在右边，所以系统开环稳定。 校正后 % MATLAB PROGRAM j005.m % num=conv([0.833 1],[6.6 20]); d3=conv([0.5 1],conv([0.025 1],[1 0]));d2=conv([0.125 1],[3.33 1]);den=conv([0 d3],[0 d2]); g=tf(num,den); sys=feedback(g,1); pzmap(g); conv([0.5 1],conv([0.025 1],[1 0]));d2=conv([0.125 1],[4 1]);den=conv([0 d3],[0 d2]); t=tf(num,den); pzmap(t); [p,z]=pzmap(g); den=sys.den{1}; r=roots(den); disp(r) -40.7534 -4.8304 -1.6921 + 3.4339i -1.6921 - 3.4339i -1.3324 极点和特征根都在左半平面，故系统稳定 ４、利用ＭＡＴＬＡＢ作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后的动态性能指标?％，ｔｒ，ｔｐ，ｔｓ，ess,并分析其有何变化？ 脉冲响应 校正前 % MATLAB PROGRAM j005.m % k=20;n1=1;d1=conv([1 0],conv([0.5 1],[0.025 1])); s1=tf(k*n1,d1); sys=feedback(s1,1); impulse(sys) 校正后 % MATLAB PROGRAM j005.m % k=conv([0.833 1],[0.4 1]);n1=20; d