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

重庆交通大学自动控制课程设计课题:ode图法控制系统设计一串联滞后校正作 者:学 号:专 业:班 级:指导教师:2013.06.28在现代科学技术的众多领域中，自动控制技术起着越来越重要的作用， 而自 动控制理论是自动控制科学的核心。自动控制理论自至今已经过了三代的发展。现代控制理论已广泛应用于制造业、农业、交通、航空及航天等众多产业部门。自动控制理论从线性近似到非线性系统的研究取得了新的成就，借助微分几何的固有非线性框架来研究非线性系统的控制，已成为目前重要研究方向之一。在控制技术需求推动下，控制理论本身也取得了显著进步。为了实现各种复杂的控制 任务首先要将被控制对象和控制装置按照一定的方式连接起来，组成一个有机整体，这就是自动控制系统。本次课程设计是利用滞后-超前校正网络来校正系统以改善系统性能，首先 应该根据原有系统和初始条件要求来确定校正系统，然后利用MATLA分析校正后的系统是否达到要求以及其性能。关键字：自动控制MATLAB滞后-超前校正 系统分析前言2一、设计任务 51、设计要求 52、设计方案分析 53、控制器的 MATLAB程序实现 6二、控制系统的模拟化设计 71、模拟控制器的离散化 7A、冲激不变法 7B、加零阶彳持器的Z变换法 8C、差分变换法 8D、双线性变换法 8E、频率预畸变双线性变换法 92、控制器的MATLAB 离散程序 93、控制器的计算机实现 9A、直接程序设计法 9B、用行程序法10C、并行程序法 11D、控制器的计算机实现流程图 11三、控制系统的MATLAB实现 121、控制系统的 MATLAB仿真 122、控制系统的 simulink仿真 13四、小结 14五、参考文献 15一、设计任务Bode图法控制系统设计-串联滞后校正、10设被控对象的传递函数为Go(s)=s(s 5)1、设计要求(1)开环增益 Kv=20(2)频率裕量丫 =70°(3)对此控制器进行离散化，并用计算机程序实现(划出流程图)(4)用Matlab对系统进行仿真，分析系统的阶跃响应2、设计方案分析系统要求使用Bode图法对控制系统进行设计，同时要求采用串联滞后校正。一、Bode图法在Bode图中的对数频率特性的低频区表征了闭环系统的稳定特性，中频区表征了系统的相对稳定性，而高频区表征了系统的抗干扰特性。在大多数实际情况中， 校正问题的实质上是一个在稳定精度和相对稳定性之间取折中的问题。为了获得比较高的开环增益及满意的相对稳定性，必须改变开环频率特性响应曲线的形状，这主要体现为：在低频区和中频区增益应足够大，且中频区的对数幅频特性的斜率应为一20d吃，并有足够的宽带，以保证适当的相角裕度；而在高频区，要使增益尽可,dec能的衰减下来，以便使高频噪声的影响达到最小。下面讨论一下基于Bode图法的串联超前校正的方法及 MATLA取现。二、串联滞后校正串联滞后校正的主要作用在不改变系统动态特性的前提下，提高系统的开环放大倍 数，使系统的稳态误差减小，并保证一定的相对稳定性。设滞后校正装置的传递函数为1 aTs.MKE1. Bode图的几何设计方法用这种方法设计滞后校正装置的步骤如下。(1) 根据稳态指标确定未校正系统Go'(s)的型别和开环增益 K,并绘制Bode图。(2) 根据动态指标要求确定滞后校正装置的参数，分两种情形加以讨论。第一种情形：给出了 0c的要求值。-, .1根据 201g |Go'(j8c)|=-20lg a,求出 a的值。a=|Go '(j -c)|为了减少滞后校正对系统的影响，通常取1/aT忠(1/5 1/10)6-并求出T =1 0 /伞c 。)第二种情形：未给出的要求值。6c若相角裕度不不足，找出满足尸的频率点作为校正系统的剪贴频率，然后按第一种情形额步骤处理。验算性能指标3、控制器的MATLAB程序实现本设计主要运用 bode图法的MATLA改计方法，主要利用伯德(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(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);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 =193.1852Transfer function:13.4 s + 1193.2 s + 1gm =Infpm =72.0162 wcg =Inf wcp =1.34174 ,13.4S 1即所得的控制器连续函数为：D S =193.2S 1、控制系统的模拟化设计1、模拟控制器的离散化方法从信号理论角度来看，模拟控制器就是模拟信号滤波器应用于反馈控制系统中作为校正装置。滤波器对控制信号中有用的信号起着保存和加强的作用，而对无用的信号起着抑制和衰减的作用。模拟控制器离散化成的数字控制器，也可以认为是数字滤波器。主要有以下几种离散化方法。A、冲激不变法冲激不变法的基本思想是：数字滤波器产生的脉冲响应序列近似等于模拟滤波器的脉冲响应函数的采样值。n设模拟控制器的传递函数为D(s) =U© ='、, aA-E(s) y s . ain在单位脉冲作用下输出响应为u(t) = L，D(s) I -Aeti 1 n ai kT u(kT) = Aei W其采样值为即数字控制器的脉冲响应序列，因此得到' n A'D(z)C 工DB、加零阶彳持器的Z变换法这种方法就是用零阶保持器与模拟控制器串联，然后再进行Z变换离散化成数字控制器,即：D(z) = Z 1 e D(s)1 s 一C、差分变换法模拟控制器若用微分方程的形式表示，其导数可用差分近似。常用的一阶差分近似方法有两种：前向差分和后向差分。对于名定D(s)U(s)一 E(s)(1).后向差分变换法=1其微分方程为： 曲 = e(t)sdt用差分代替微分，则du(t)dtu(k) -y-1)= e(k)两边取z变换得D(z)=U(z)E(z)_1(1-z )U(z)=TE(z)1一 1 -zJT可以看出，D(z)与D(s)的形式完全相同，由此可得如下等效代换关系1-z s 二T便可得到D(z),即 D(z) = D(s) s-T(2).前向差分变换法如果将微分用下面差分代替，得到du(t)dt呻 1的二 e(k)两边取Z变换得 (z-1)U(z) =TE(z) 即嗡*T由此可得如下等效代换关系可得到D(z) =D(s)前向差分变换法中稳定的D(s)不能保证变换成稳定的D(z),且不能保证有相同的脉冲响应和频率响应。D、双线性变换法双线性变换又称塔斯廷(Tustin)变换法，它是s与z关系的另一种近似式。由Z变换的定义和级数展开式可知Tsts eT Z = eTse-2因此2 1 -z s =T 1 zd Ts1 -21 一至2即 D(z) = D(s) 2(j工)s =_ Fi也一E、频率预畸变双线性变换法上述的双线性变换，将 S平面的虚轴变换到 z平面的单位圆周，因而没有混叠现象。但 是在模拟频率和离散频率3之间却存在非线性关系。2 X T1 = 一 tg T 2当coT取值0兀时，Q的值为08。这意味着，模拟滤波器的全部频率响应特性被 压缩到离散滤波器的 0< 3 T<兀的频率范围之内。这两种频率之间的非线性特性，使得由双 线性变换所得的离散频率响应产生畸变，可以采用预畸变的办法来补偿频率特性的畸变。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.06884z - 0.9995Sampling time: 0.10.06936 z -0.06884即所得离散化后控制器为：D Z =z - 0.99953、控制器的计算机实现控制器的计算机实现共有三种方法：直接程序设计法、串行程序设计法和并行程序设计 法。A、直接程序设计法D(z)=设数字控制器的一般形式为:mn U(z) c biZE(z)ajZU(z) i Oj 1mn取 Z 反变换得： u(k)=£ bie(k -i)-Z aju(k-j) i =0j 4显然，上式是物理可实现的。其中，z-1为延时环节，每计算一次 u(k)需要以前n个输出值以及当前误差值 e(k)和以前 个误差值，所以必须将这些值存储起来以备使用。这样，每计算一次 u(k)需要做n+m+1 乘法，n+m次加法，并做n+m次数据转移。可以改进其算法，以便减少延时器，即减少运算次数。mn由U(z)=£ bzE(z)£ ajz=U(z)并设 a0=1 i =0j EnmU(z)< ajZ=E(z)" bz， j =0i =0-1_定义： Q(z) E(z)"ajz"j j =0 m所以改进的直接程序设计法中：U(z) -Q(z)v bzi =0 mu(k)="biq(k -i) i=0nQ(z) = E(z) -Q(z)x ajz口 j mnq(k)=e(k) -% ajq(k - j)B、用行程序法j 12,设数字控制器的零点、极点分别为 zi, pj, i= 1, 2,，m; j= 1,D(z)U(z)E(z)K(z -