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

《自动控制原理MATLAB课程设计--滞后-超前校正.docx》由会员分享，可在线阅读，更多相关《自动控制原理MATLAB课程设计--滞后-超前校正.docx（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、自动控制原理MATLAB课程设计-滞后-超前校正 滞后-超前校正 课程设计 一、设计目的: 1. 了解控制系统设计的一般方法、步骤。 2. 掌握对系统进行稳定性的分析、稳态误差分析以及动态特性分析的方 法。 3. 掌握利用MATLAB 对控制理论内容进行分析和研究的技能。 4. 提高分析问题解决问题的能力。 二、设计内容与要求: 设计内容： 1. 阅读有关资料。 2. 对系统进行稳定性分析、稳态误差分析以及动态特性分析。 3. 绘制根轨迹图、Bode 图、Nyquist 图。 4. 设计校正系统，满足工作要求。 设计条件： 、被控制对象的传递函数是 m m 1m 2022m n sn 1n 2

2、022n b s b s b s b ()a s a a s a G S -+?+=+?（） 、参数，和，因小组而异。 设计要求： 1. 能用MATLAB 解复杂的自动控制理论题目。 2. 能用MATLAB 设计控制系统以满足具体的性能指标。 3. 能灵活应用MATLAB 的CONTROL SYSTEM 工具箱和SIMULINK 仿真软件， 分析系统的性能。 三、设计步骤： 、自学软件的基本知识，包括的基本操作命令。控制系统工具箱的用法等，并上机实验。 、基于用频率法对系统进行串联校正设计，使其满足给定的领域 性能指标。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的

3、参数，等的值。 已知开环传递函数为G(S)= 0 (2)(40)k s s s +,使用频率法设计串联滞后超前校 正装置，使系统的相角裕度大于等于40，静态速度误差系数等于20。 校正前 根据上式可化简G(S)= 0 0.0125(0.51)(0.0251) k s s s +，所以公式G(S)= 20 (0.51)(0.0251)s s s +, 所以=1，则c w = 6.1310， 相角裕度为9.3528。 校正后 串联校正滞后超前校正装置2121(1)(1) ()(1)(1) C bT S aT S G S T S T S += +, 由 111()510 C W bT =，取 221

4、1 *60.8335 bT bT =?=，又由20lg 2022lg520lg 4 a a -=?= , 121aT =,10.125T =;10.25ab b =?=;1 1.22 3.332 T bT =?= 所以校正后的函数为 (10.833)(10.5)20 ()(10.125)(1 3.33)(0.51)(0.0251)C S S G S S S s s s += + 、利用函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，为什么？ 校正前 % MATLAB PROGRAM j005.m % num=20; den= conv(1 0,conv(0.5 1,0.025 1);

5、 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.den1; 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.02

6、5 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.den1; r=roots(den); disp(r) -40.7534 -4.8304 -1.6921 + 3.4339i -1.6921 - 3.4339i -1.

7、3324 极点和特征根都在左半平面，故系统稳定 、利用作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后的动态性能指标?，e ss ,并分析其有何变化？ 脉冲响应 校正前 % 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

8、3=conv(0.5 1,conv(0.025 1,1 0);d2=conv(3.33 1,0.125 1);d1=conv(0 d3,0 d2); s1=tf(k*n1,d1); sys=feedback(s1,1); impulse(sys) 阶跃响应 校正前 % MATLAB PROGRAM j005.m % k=20;n1=1;d1=conv(1 0,conv(0.5 1,0.025 1); s1=tf(k*n1,d1); figure(2);sys=feedback(s1,1);step(sys) 由图可知：?=76%，=0.272-0.093=0.179，=0.506，=6.09，

9、e ss=1-1=0 校正后 % MATLAB PROGRAM j005.m % k=conv(0.833 1,0.4 1);n1=20; d3=conv(0.025 1,conv(0.125 1,1 0);d2=conv(0.5 1,3.33 1);d1=conv(0 d3,0 d2); s1=tf(k*n1,d1); figure(2);sys=feedback(s1,1);step(sys) 由图可知：?=36%，=0.384-0.109=0.275，=0.697，=1.27，e ss=1-1.01=-0.01 斜坡响应 校正前 % MATLAB PROGRAM j005.m % k=2

10、0;n1=1;d1=conv(1 0 0,conv(0.5 1,0.025 1); s1=tf(k*n1,d1); figure(2);sys=feedback(s1,1);step(sys) 校正后 % MATLAB PROGRAM j005.m % k=conv(0.833 1,0.4 1);n1=20; d3=conv(0.5 1,conv(0.025 1,1 0);d2=conv(3.33 1,0.125 1);d4=conv(0 d3,0 d2);d1=conv(0 d4,1 0); s1=tf(k*n1,d1); figure(2);sys=feedback(s1,1);step(

11、sys) 三条曲线关系：斜坡响应曲线的导数是阶跃，阶跃响应曲线的导数是脉冲响应曲线。 、绘制系统校正前与校正后的根轨迹图，并求其分离点，汇合点及虚轴交点的坐标和相应点的增益*值，得出系统稳定时增益K*的变化范围。绘制系统校正前与校正后的Nyquist图，判断系统的稳定性，并说明理由。 校正前 % MATLAB PROGRAM j005.m % num=20; den=conv(1 0,conv(0.5 1,0.025 1); g0=tf(num,den); rlocus(g0) 分离点 -0.987，无汇合点，虚轴交点的坐标 0.129 9.37j和相应点的增益* =0.0244 校正后 % MATLAB PROGRAM j005.m % 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); num=20; g0=tf(num,den); rlocus(g0)