《自动控制原理》课程设计-基于MATLAB进行控制系统的滞后-超前校正设计.doc

武汉理工大学自动控制原理课程设计说明书 题 目: 用MATLAB进行控制系统的滞后超前校正设计 目录1 滞后-超前校正设计目的和原理11.1 滞后-超前校正设计目的11.2 滞后-超前校正设计原理12 滞后-超前校正的设计过程32.1 校正前系统的参数32.1.1 用MATLAB绘制校正前系统的伯德图32.1.2 用MATLAB求校正前系统的幅值裕量和相位裕量42.1.3 用MATLAB绘制校正前系统的根轨迹52.1.4 对校正前系统进行仿真分析52.2 滞后-超前校正设计参数计算62.2.1 选择校正后的截止频率62.2.2 确定校正参数、和62.3 滞后-超前校正后的验证72.3.1 用MATLAB求校正后系统的幅值裕量和相位裕量72.3.2 用MATLAB绘制校正后系统的伯德图82.3.3 用MATLAB绘制校正后系统的根轨迹92.3.4 用MATLAB对校正前后的系统进行仿真分析103 心得体会12参考文献13用MATLAB进行控制系统的滞后超前校正设计1 滞后-超前校正设计目的和原理1.1 滞后-超前校正设计目的所谓校正就是在系统不可变部分的基础上，加入适当的校正元部件，使系统满足给定的性能指标。校正方案主要有串联校正、并联校正、反馈校正和前馈校正。确定校正装置的结构和参数的方法主要有两类：分析法和综合法。分析法是针对被校正系统的性能和给定的性能指标，首先选择合适的校正环节的结构，然后用校正方法确定校正环节的参数。在用分析法进行串联校正时，校正环节的结构通常采用超前校正、滞后校正和滞后-超前校正这三种类型。超前校正通常可以改善控制系统的快速性和超调量，但增加了带宽，而滞后校正可以改善超调量及相对稳定度，但往往会因带宽减小而使快速性下降。滞后-超前校正兼用两者优点，并在结构设计时设法限制它们的缺点。1.2 滞后-超前校正设计原理滞后-超前校正RC网络电路图如图1所示。图1 滞后-超前校正RC网络下面推导它的传递函数：令，则其中为超前部分的参数，为滞后部分的参数。滞后-超前校正的频域设计实际是超前校正和滞后校正频域法设计的综合，基本方法是利用滞后校正将系统校正后的穿越频率调整到超前部分的最大相角处的频率。具体方法是先合理地选择截止频率，先设计滞后校正部分，再根据已经选定的设计超前部分。应用频率法确定滞后超前校正参数的步骤：1、根据稳态性能指标，绘制未校正系统的伯德图；2、选择校正后的截止频率；3、确定校正参数；4、确定滞后部分的参数；5、确定超前部分的参数；6、将滞后部分和超前部分的传递函数组合在一起，即得滞后-超前校正的传递函数；7、绘制校正后的伯德图，检验性能指标。2 滞后-超前校正的设计过程2.1 校正前系统的参数根据初始条件，调整开环传递函数：当系统的静态速度误差系数时，。则满足初始条件的最小K值时的开环传递函数为2.1.1 用MATLAB绘制校正前系统的伯德图程序：num=10;den=0.5,1.5,1,0;bode(num,den)grid得到的伯德图如图2所示。图2 校正前系统的伯德图2.1.2 用MATLAB求校正前系统的幅值裕量和相位裕量用命令margin(G)可以绘制出G的伯德图，并标出幅值裕量、相位裕量和对应的频率。用函数kg,r,wg,wc=margin(G)可以求出G的幅值裕量、相位裕量和幅值穿越频率。程序：num=10;den=0.5,1.5,1,0;G=tf(num,den);margin(G)kg,r,wg,wc=margin(G)得到的幅值裕量和相位裕量如图3所示。图3 校正前系统的幅值裕量和相位裕量运行结果： kg=0.3000 r=-28.0814 wg=1.4142 wc=2.4253即幅值裕量，相位裕量=-28.0814o。2.1.3 用MATLAB绘制校正前系统的根轨迹MATLAB中专门提供了绘制根轨迹的有关函数。p,z=pzmap(num,den)的功能是绘制连续系统的零、极点图。r,k=rlocus(num,den)的功能是绘制部分的根轨迹。程序：num=10;den=0.5,1.5,1,0;rlocus(num,den)得到校正前系统的根轨迹如图4所示。图4 校正前系统的根轨迹2.1.4 对校正前系统进行仿真分析Simulink是可以用于连续、离散以及混合的线性、非线性控制系统建模、仿真和分析的软件包，并为用户提供了用方框图进行建模的图形接口，很适合于控制系统的仿真。仿真后得到的结果如图5和图6所示。图5 校正前系统的仿真图图6 校正前系统仿真的阶跃响应曲线2.2 滞后-超前校正设计参数计算2.2.1 选择校正后的截止频率若性能指标中对系统的快速性未提明确要求时，一般对应的频率作为。从图3中得，=1.5。这样，未校正系统的相位裕量为0o，与要求值仅差+45o，这样大小的超前相角通过简单的超前校正是很容易实现的。2.2.2 确定校正参数、和由超前部分应产生超前相角而定，即。在本题中，因此取，以使滞后相角控制在-5o以内，因此，滞后部分的传递函数为。过，作20dB/dec直线，由该直线与0dB线交点坐标确定。未校正系统的伯德图在=1.5处的增益是13dB。所以过点（1.5，-13）画一条20dB/dec的直线，与0dB线的交点确定转折频率。经计算得，转折频率，另一转折频率为。所以超前部分的传递函数为。将滞后校正部分和超前校正部分的传递函数组合在一起，得滞后-超前校正的传递函数为系统校正后的传递函数为2.3 滞后-超前校正后的验证由于校正过程中，多处采用的是近似计算，可能会造成滞后-超前校正后得到的系统的传递函数不满足题目要求的性能指标。所以需要对滞后-超前校正后的系统进行验证。下面用MATLAB求已校正系统的相角裕量和幅值裕量。2.3.1 用MATLAB求校正后系统的幅值裕量和相位裕量程序：num=10,9.9,0.89;den=0.5,4.855,11.0985,6.8055,0.067,0;G=tf(num,den);margin(G)kg,r,wg,wc=margin(G)得到的校正后系统的幅值裕量和相位裕量如图7所示。图7 校正后系统的幅值裕量和相位裕量运行结果： kg=5.9195 r=47.6239 wg=3.6762 wc=1.2072即校正后系统的相位裕量，满足指标。假设验证结果不满足指标，重新选择校正后的截止频率，重复上述过程，直到满足性能指标为止。2.3.2 用MATLAB绘制校正后系统的伯德图程序：num=10,9.9,0.89;den=0.5,4.855,11.0985,6.8055,0.067,0;bode(num,den)grid得到的伯德图如图8所示。图8 校正后系统的伯德图2.3.3 用MATLAB绘制校正后系统的根轨迹程序：num=10,9.9,0.89;den=0.5,4.855,11.0985,6.8055,0.067,0;rlocus(num,den)得到的校正后系统的根轨迹如图9所示。图9 校正后系统的根轨迹2.3.4 用MATLAB对校正前后的系统进行仿真分析用Simulink对校正后的系统仿真。仿真后得到的结果如图10和图11所示。图10 校正后系统的仿真图图11 校正后系统仿真的阶跃响应曲线程序：k=10;num=conv(1,0.89,1,0.1);den=conv(conv(conv(conv(1,0,1,1),0.5,1),1,6.7),1,0.01);sys=tf(k*num,den);Lsys=feedback(sys,1,-1);y,t,x=step(Lsys);plot(t,y);ltiview得到的阶跃响应曲线如图12所示。图12 校正后阶跃响应曲线调节时间取的误差范围。由图12可知，超调量，上升时间 ，峰值时间，调节时间。对比校正前后的阶跃响应曲线可知，校正前系统是不稳定的，无法求得时域性能指标。校正后的系统是稳定的，系统的阶跃响应曲线是衰减振荡的。当调节时间取的误差范围时，调节时间。3 心得体会随着科学技术发展的日新月异，MATLAB已成为当今应用软件中空前活跃的领域，在生活中的应用可以说是无处不在，因此掌握MATLAB这个软件基本的使用方法对我们是十分有益的。MATLAB可用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。当然，MATLAB也可以用对反馈系统进行校正。此次课程设计的内容对一个单位反馈系统进行滞后-超前校正。回顾此次实践的整个过程，虽然只有短短的几天，但是真的在这个自己独立学习的过程中学到了好多东西。课程设计开始阶段比较顺利，但是做到计算校正后系统的时域性能指标这里时，遇到了不小麻烦，不会用MATLAB编程得阶跃响应曲线。后来，在同学的帮助下，终于看懂了课件了的那段程序，于是就求出了滞后-超前校正后系统的时域性能指标。这次的课程设计，不仅让我们更好的更深一步的了解MATLAB这个十分有用的软件，也能运用他对某一电路图进行仿真，与理论上相结合，从而进一步验证理论的正确性，也是理论运用于实践的很好的证明。与此同时，通过此次课程设计，加深了系统进行滞后-超前设计过程的理解，还掌握了用MATLAB编程计算系统时域性能指标和系统幅值裕量、相位裕量的方法。总而言之，这次的课程设计的确让我受益匪浅，还让我把许多新知识尽