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

1重庆交通大学自动控制课程设计课题：三、Bode 图法控制系统设计 -串联滞后校正作 者：学 号：专 业：班 级：指导教师：2013.06.282摘 要在现代科学技术的众多领域中，自动控制技术起着越来越重要的作用，而自动控制理论是自动控制科学的核心。自动控制理论自至今已经过了三代的发展。现代控制理论已广泛应用于制造业、农业、交通、航空及航天等众多产业部门。自动控制理论从线性近似到非线性系统的研究取得了新的成就，借助微分几何的固有非线性框架来研究非线性系统的控制，已成为目前重要研究方向之一。在控制技术需求推动下，控制理论本身也取得了显著进步。为了实现各种复杂的控制任务首先要将被控制对象和控制装置按照一定的方式连接起来，组成一个有机整体，这就是自动控制系统。本次课程设计是利用滞后-超前校正网络来校正系统以改善系统性能，首先应该根据原有系统和初始条件要求来确定校正系统，然后利用 MATLAB 分析校正后的系统是否达到要求以及其性能。关键字：自动控制 MATLAB 滞后-超前校正 系统分析3目 录前言21、设计任务 51、设计要求 52、设计方案分析 53、控制器的 MATLAB 程序实现 6 二、控制系统的模拟化设计 71、模拟控制器的离散化 7A、冲激不变法 7B、加零阶保持器的 Z 变换法 8C、差分变换法 8D、双线性变换法 8E、频率预畸变双线性变换法 92、控制器的 MATLAB 离散程序 93、控制器的计算机实现 9A、直接程序设计法 9B、串行程序法 10C、并行程序法 411D、控制器的计算机实现流程图 11三、控制系统的 MATLAB 实现 121、控制系统的 MATLAB 仿真 122、控制系统的 simulink 仿真 13四、小结 14五、参考文献 155一、设计任务Bode 图法控制系统设计-串联滞后校正设被控对象的传递函数为 1、设计要求（1）开环增益 Kv=20（2）频率裕量 =70（3）对此控制器进行离散化，并用计算机程序实现（划出流程图）（4）用 Matlab 对系统进行仿真，分析系统的阶跃响应2、设计方案分析系统要求使用 Bode 图法对控制系统进行设计，同时要求采用串联滞后校正。一、Bode 图法在 Bode 图中的对数频率特性的低频区表征了闭环系统的稳定特性，中频区表征了系统的相对稳定性，而高频区表征了系统的抗干扰特性。在大多数实际情况中，校正问题的实质上是一个在稳定精度和相对稳定性之间取折中的问题。为了获得比较高的开环增益及满意的相对稳定性，必须改变开环频率特性响应曲线的形状，这主要体现为：在低频区和中频区增益应足够大，且中频区的对数幅频特性的斜率应为 ，并有足够的宽带，以保证适当的相角裕度；而在高频区，要使增益20dBec尽可能的衰减下来，以便使高频噪声的影响达到最小。下面讨论一下基于 Bode 图法的串联超前校正的方法及 MATLAB 实现。二、串联滞后校正)5(10)0sG6串联滞后校正的主要作用在不改变系统动态特性的前提下，提高系统的开环放大倍数，使系统的稳态误差减小，并保证一定的相对稳定性。设滞后校正装置的传递函数为 1()()CcaTsGsK1Bode 图的几何设计方法用这种方法设计滞后校正装置的步骤如下。（1） 根据稳态指标确定未校正系统 的型别和开环增益 K，并绘制 Bode 图。0()s（2） 根据动态指标要求确定滞后校正装置的参数，分两种情形加以讨论。第一种情形：给出了 的要求值。c根据 ，求出 a 的值。 。0lg|()|2lgcGj01|()|cGj为了减少滞后校正对系统 的影响，通常取 。并求出 1/5T。10/()cTa第二种情形：未给出的要求值。若相角裕度 不足，找出满足 的频率点作为校正系统的剪贴频率，然后按第一种 c情形额步骤处理。验算性能指标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(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);7step(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 + 1-193.2 s + 1gm =Infpm =72.0162wcg =Infwcp =1.34178即所得的控制器连续函数为： 13.492SD二、控制系统的模拟化设计1、模拟控制器的离散化方法从信号理论角度来看，模拟控制器就是模拟信号滤波器应用于反馈控制系统中作为校正装置。滤波器对控制信号中有用的信号起着保存和加强的作用，而对无用的信号起着抑制和衰减的作用。模拟控制器离散化成的数字控制器，也可以认为是数字滤波器。主要有以下几种离散化方法。A、冲激不变法 冲激不变法的基本思想是：数字滤波器产生的脉冲响应序列近似等于模拟滤波器的脉冲响应函数的采样值。设模拟控制器的传递函数为在单位脉冲作用下输出响应为 其采样值为 即数字控制器的脉冲响应序列，因此得到B、加零阶保持器的 Z 变换法 这种方法就是用零阶保持器与模拟控制器串联，然后再进行 Z 变换离散化成数字控制器，即 ：C、差分变换法 模拟控制器若用微分方程的形式表示，其导数可用差分近似。常用的一阶差分近似方法有两种：前向差分和后向差分 。（1） 后向差分变换法 对于给定 其微分方程为 ：用差分代替微分，则 两边取 Z 变换得 即 niiasAEUsD1)(itieLtunikTaieAku1)( )(1)()( 1sDzeAkTuzDniTai)(1)(sDezTsEUsD1)()(tedu()1)dutkkT)1(zzTEz)( Tzs19可以看出，D(z)与 D(s)的形式完全相同，由此可得如下等效代换关系 ：便可得到 D(z)，即 （2） 前向差分变换法如果将微分用下面差分代替，得到 ()1)()dutkuekT两边取 Z 变换得 即 由此可得如下等效代换关系 可得到前向差分变换法中稳定的 D(s)不能保证变换成稳定的 D(z)，且不能保证有相同的脉冲响应和频率响应。 D、双线性变换法 双线性变换又称塔斯廷(Tustin)变换法，它是 s 与 z 关系的另一种近似式。由 Z 变换的定义和级数展开式可知因此即 E、频率预畸变双线性变换法 上述的双线性变换，将 S 平面的虚轴变换到 Z 平面的单位圆周，因而没有混叠现象。但是在模拟频率 和离散频率 之间却存在非线性关系。当 T 取值 0 时， 的值为 0。这意味着，模拟滤波器的全部频率响应特性被压缩到离散滤波器的 0 T 的频率范围之内。这两种频率之间的非线性特性，使得由双线性变换所得的离散频率响应产生畸变，可以采用预畸变的办法来补偿频率特性的畸变。 2、控制器的 MATLAB 离散程序采用 MATLAB 编程实现连续系统的离散化，此程序基于加零阶保持器的 Z 变换法，使Tzsz1)()()1(zTEUzEzD)(Tzs1Tzsz1()212TsezsTs12Ts)1(2)(zTsDz2tg10用函数 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.9995Sampling time: 0.1即所得离散化后控制器为： 0.693 z.084Z=5D3、控制器的计算机实现控制器的计算机实现共有三种方法：直接程序设计法、串行程序设计法和并行程序设计法。A、直接程序设计法设数字控制器的一般形式为：取 Z 反变换得：显然，上式是物理可实现的。其中，z-1 为延时环节，每计算一次 u(k)需要以前 n 个输出值以及当前误差值 e(k)和以前 m个误差值，所以必须将这些值存储起来以备使用。这样，每计算一次 u(k)需要做 n+