Bode 图法控制系统设计---串联滞后校正

1、重庆交通大学 自动控制课程设计课题：三、Bode 图法控制系统设计-串联滞后校正 作 者： 学 号： 专 业： 班 级： 指导教师： 2013.06.28摘 要在现代科学技术的众多领域中，自动控制技术起着越来越重要的作用，而自动控制理论是自动控制科学的核心。自动控制理论自至今已经过了三代的发展。现代控制理论已广泛应用于制造业、农业、交通、航空及航天等众多产业部门。自动控制理论从线性近似到非线性系统的研究取得了新的成就，借助微分几何的固有非线性框架来研究非线性系统的控制，已成为目前重要研究方向之一。在控制技术需求推动下，控制理论本身也取得了显著进步。为了实现各种复杂的控制任务首先要将被控制对象和

2、控制装置按照一定的方式连接起来，组成一个有机整体，这就是自动控制系统。本次课程设计是利用滞后-超前校正网络来校正系统以改善系统性能，首先应该根据原有系统和初始条件要求来确定校正系统，然后利用MATLAB分析校正后的系统是否达到要求以及其性能。关键字：自动控制 MATLAB 滞后-超前校正 系统分析 目 录前言21、 设计任务 5 1、设计要求 5 2、设计方案分析 5 3、控制器的MATLAB程序实现 6 二、控制系统的模拟化设计 7 1、模拟控制器的离散化 7 A、冲激不变法 7 B、加零阶保持器的Z变换法 8 C、差分变换法 8 D、双线性变换法 8 E、频率预畸变双线性变换法 9 2、控

3、制器的MATLAB离散程序 9 3、控制器的计算机实现 9 A、直接程序设计法 9 B、串行程序法 10 C、并行程序法 11 D、控制器的计算机实现流程图 11三、控制系统的MATLAB实现 12 1、控制系统的MATLAB仿真 12 2、控制系统的simulink仿真 13四、小结 14五、参考文献 15一、设计任务 Bode 图法控制系统设计-串联滞后校正 设被控对象的传递函数为 1、设计要求 （1）开环增益Kv=20 （2）频率裕量=70 （3）对此控制器进行离散化，并用计算机程序实现（划出流程图）（4） 用Matlab对系统进行仿真，分析系统的阶跃响应2、设计方案分析系统要求使用Bo

4、de 图法对控制系统进行设计，同时要求采用串联滞后校正。一、Bode图法在Bode图中的对数频率特性的低频区表征了闭环系统的稳定特性，中频区表征了系统的相对稳定性，而高频区表征了系统的抗干扰特性。在大多数实际情况中，校正问题的实质上是一个在稳定精度和相对稳定性之间取折中的问题。为了获得比较高的开环增益及满意的相对稳定性，必须改变开环频率特性响应曲线的形状，这主要体现为：在低频区和中频区增益应足够大，且中频区的对数幅频特性的斜率应为 ，并有足够的宽带，以保证适当的相角裕度；而在高频区，要使增益尽可能的衰减下来，以便使高频噪声的影响达到最小。下面讨论一下基于Bode图法的串联超前校正的方法及MAT

5、LAB实现。二、串联滞后校正 串联滞后校正的主要作用在不改变系统动态特性的前提下，提高系统的开环放大倍数，使系统的稳态误差减小，并保证一定的相对稳定性。设滞后校正装置的传递函数为1Bode图的几何设计方法用这种方法设计滞后校正装置的步骤如下。（1） 根据稳态指标确定未校正系统的型别和开环增益K，并绘制Bode图。（2） 根据动态指标要求确定滞后校正装置的参数，分两种情形加以讨论。 第一种情形：给出了的要求值。 根据，求出a的值。 为了减少滞后校正对系统的影响，通常取。并求出 。 第二种情形：未给出的要求值。若相角裕度不足，找出满足的频率点作为校正系统的剪贴频率，然后按第一种 情形额步骤处理。

6、验算性能指标3、控制器的MATLAB程序实现本设计主要运用bode图法的MATLAB设计方法，主要利用伯德（Bode）图进行系统的设计，用到的函数有：Bode伯德图作图命令；Logspace用于在某个区域中产生若干频点；Polyval求多项式的值；Ceil朝正无穷大方向取整；ng0=10;dg0=1,5,0;g0=tf(ng0,dg0);t=0:0.01:3;w=logspace(-2,2);kk=10;pm=70;mu,pu=bode(kk*ng0,dg0,w);wgc=spline(pu,w,pm+5-180),ngv=polyval(kk*ng0,j*wgc);dgv=polyval(d

7、g0,j*wgc);g=ngv/dgv;alph=abs(1/g),t=10/alph*wgc,ngc=alph*t,1;dgc=t,1;gc=tf(ngc,dgc),g0c=tf(kk*g0*gc);b1=feedback(kk*g0,1);b2=feedback(g0c,1);step(b1,t);grid on,hold on;step(b2,t),hold offfigure,bode(kk*g0,w),grid on,hold on;bode(g0c,w);hold offgm,pm,wcg,wcp=margin(g0c)wgc = 1.3397alph = 0.0694t = 19

8、3.1852 Transfer function:13.4 s + 1-193.2 s + 1 gm = Infpm = 72.0162wcg = Infwcp =1.3417即所得的控制器连续函数为：二、控制系统的模拟化设计1、模拟控制器的离散化方法 从信号理论角度来看，模拟控制器就是模拟信号滤波器应用于反馈控制系统中作为校正装置。滤波器对控制信号中有用的信号起着保存和加强的作用，而对无用的信号起着抑制和衰减的作用。模拟控制器离散化成的数字控制器，也可以认为是数字滤波器。主要有以下几种离散化方法。A、冲激不变法 冲激不变法的基本思想是：数字滤波器产生的脉冲响应序列近似等于模拟滤波器的脉冲响应

9、函数的采样值。设模拟控制器的传递函数为在单位脉冲作用下输出响应为 其采样值为 即数字控制器的脉冲响应序列，因此得到B、加零阶保持器的Z变换法 这种方法就是用零阶保持器与模拟控制器串联，然后再进行Z变换离散化成数字控制器，即 ：C、差分变换法 模拟控制器若用微分方程的形式表示，其导数可用差分近似。常用的一阶差分近似方法有两种：前向差分和后向差分 。（1）后向差分变换法 对于给定 其微分方程为 ：用差分代替微分，则 两边取Z变换得 即 可以看出，D(z)与D(s)的形式完全相同，由此可得如下等效代换关系 ：便可得到D(z)，即 （2）前向差分变换法如果将微分用下面差分代替，得到 两边取Z变换得 即

10、 由此可得如下等效代换关系 可得到前向差分变换法中稳定的D(s)不能保证变换成稳定的D(z)，且不能保证有相同的脉冲响应和频率响应。 D、双线性变换法 双线性变换又称塔斯廷(Tustin)变换法，它是s与z关系的另一种近似式。由Z变换的定义和级数展开式可知因此即 E、频率预畸变双线性变换法 上述的双线性变换，将S平面的虚轴变换到Z平面的单位圆周，因而没有混叠现象。但是在模拟频率和离散频率之间却存在非线性关系。 当T取值0时，的值为0。这意味着，模拟滤波器的全部频率响应特性被压缩到离散滤波器的0T的频率范围之内。这两种频率之间的非线性特性，使得由双线性变换所得的离散频率响应产生畸变，可以采用预畸

11、变的办法来补偿频率特性的畸变。 2、控制器的MATLAB离散程序采用MATLAB编程实现连续系统的离散化，此程序基于加零阶保持器的Z变换法，使用函数c2dm作为指令。程序和结果如下：num=13.4,1;den=193.2,1;ts=0.1;x,y=c2dm(num,den,ts);tf(x,y,ts) Transfer function:0.06936 z - 0.06884- z - 0.9995 Sampling time: 0.1即所得离散化后控制器为：3、控制器的计算机实现 控制器的计算机实现共有三种方法：直接程序设计法、串行程序设计法和并行程序设计法。A、直接程序设计法 设数字控制

12、器的一般形式为：取Z反变换得：显然，上式是物理可实现的。其中，z-1为延时环节，每计算一次u(k)需要以前n个输出值以及当前误差值e(k)和以前m个误差值，所以必须将这些值存储起来以备使用。这样，每计算一次u(k)需要做n+m+1次乘法，n+m次加法，并做n+m次数据转移。可以改进其算法，以便减少延时器，即减少运算次数。由 并设a0=1 定义： 所以改进的直接程序设计法中：B、串行程序法 设数字控制器的零点、极点分别为zi，pj，i=1，2，m；j=1，2，n；nm，则其中： 则得到：C、并行程序法 设D(z)的极点为pi，则可将D(z)表示成部分分式其中： 令 则 得到：D、控制器的计算机实

13、现流程图由此控制器形式： 可以看出，它只适合用直接程序法实现。所以，应用直接程序法 则流程图：（如下）-0.068841/z0.069360.99951/z三、控制系统的MATLAB实现1、控制系统的MATLAB仿真MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。通过MATLAB程序绘制出bode图，程序如下：ng0=10;dg0=1,5,0

14、;g0=tf(ng0,dg0);t=0:0.01:3;w=logspace(-2,2);kk=10;pm=70;mu,pu=bode(kk*ng0,dg0,w);wgc=spline(pu,w,pm+5-180),ngv=polyval(kk*ng0,j*wgc);dgv=polyval(dg0,j*wgc);g=ngv/dgv;alph=abs(1/g),t=10/alph*wgc,ngc=alph*t,1;dgc=t,1;gc=tf(ngc,dgc),g0c=tf(kk*g0*gc);b1=feedback(kk*g0,1);b2=feedback(g0c,1);step(b1,t);gr

15、id on,hold on;step(b2,t),hold offfigure,bode(kk*g0,w),grid on,hold on;bode(g0c,w);hold offgm,pm,wcg,wcp=margin(g0c)改进前后的Bode图如下其中gm = Infpm = 72.0162wcg = Infwcp =1.3417复合设计要求2、控制系统的simulink仿真 Simulink是MATLAB的重要组成部分，提供建立系统模型、选择仿真参数和数值算法、启动仿真程序对该系统进行仿真、设置不同的输出方式来观察仿真结果等功能。仿真后系统校正前后单位阶跃响应曲线如下图所示：可见，阶跃

16、响应明显改善。4、 小结通过此次课程设计学习，让我对一些重要软件的操作也有了很大的提高，在此次设计中帮到很大忙。通过它们可以进行设计仿真，检验设计的正确性。在设计的过程中，我遇到了一些难题，例如MATLAB仿真软件操作问题；通过上网查阅相关教程，锻炼了自身的自学能力，也同时与同学相互探讨，培养了自身的协作能力。要把理论转化为实践是必须的，学生必须尽力为自己寻找一些实验，动手的机会。课程设计为我们提供了这样的机会。课设过程中，大家自己独立思考，完成老师布置的题目，学习了很多东西，把自己所学用于实际，是我对自动化专业有了更深的了解，同时极大的提高了我对自动化专业的兴趣。课堂学习主要注重于理论知识，

17、而我们要将所学知识应用于实际，在此阶段，课程设计便是最好的选择了，通过课程设计，我们可以温习我们所学的理论知识，同时为将理论知识运用于实际搭建了一个很好的平台，不仅如此，通过这次的课程设计，使我知道了在当今的信息技术如此发达的世界中，我们必须运用多种渠道，去学习研究。并要很好的运用计算机和一些软件，只有这样，我们才能更好地、精确地、快速地解决问题。还有就是提高了自主解决问题的能力。这次校正系统的设计，我不但运用了以前在书本上学过的知识，还运用了新的软件MATLAB软件，在此之前我们并没有学过MATLAB软件的运用，通过查阅相关资料，终于用MATLAB软件实现了所要求的功能。通过动手实践也让我对校正方法有了更深刻的了解，这绝对是对自己自控知识的巩固与提高。MATLAB软件的强大功能确实令人兴奋，在以后的学习、生活中它将会扮演重要的角色。总之，在这次课程设计过程中，我既学习到了自动控制原