1、自动控制原理第七章-离散控制系统教学重点l 了解线性离散系统的基本概念和基本定理,把握线性连续系统与线性离散系统的区别与联系;l 熟练掌握Z变换的方法、Z变换的性质和Z反变换;l 了解差分方程的定义,掌握差分方程的解法;l 了解脉冲传递函数的定义,熟练掌握开环与闭环系统脉冲传递函数的计算方法;l 与线性连续系统相对应,掌握线性离散系统的时域和频域分析方法和原则。2022-10-62教学难点 离散时间函数的数学表达式及采样定理,线性常系数差分方程与脉冲传递函数,采样控制系统的时域分析,采样控制系统的频域分析。2022-10-63概述:概述:近年来,随着脉冲技术、数字式元器件、数字计算机,特别是微
2、处理器的迅速发展,数字控制器在许多场合取代了模拟控制器,比如微型数字计算机在控制系统中得到了广泛的应用。离散系统理论的发展是非常迅速的。因此,深入研究离散系统理论,掌握分析与综合数字控制系统的基础理论与基本方法,从控制工程特别是从计算机控制工程角度来看,是迫切需要的。离散系统与连续系统相比,既有本质上的不同,又有分析研究方法的相似性。利用z变换法研究离散系统,可以把连续系统中的许多概念和方法,推广应用于离散系统。本章内容本章内容:主要介绍线性离散系统的分析方法。首先给出信号采样和保持的数学描述。然后介绍z变换理论与性质以及脉冲传递函数。接着研究线性离散系统稳定性、稳态误差、动态性能的分析方法,
3、并介绍最少拍系统的设计方法。最后介绍如何利用MATLAB进行线性离散系统的分析。2022-10-647.1 7.1 离散控制系统的基本概念离散控制系统的基本概念1.连续系统:如果控制系统中的所有信号都是时间变量的连续函数,也就 是说,这些信号在全部时间上都是已知的,则这样的系统称为 连续时间系统,简称连续系统。2.离散系统:如果控制系统中有一处或几处信号是脉冲序列或数码,则这样 的系统称为离散时间系统,简称离散系统。采样控制系统:系统中的离散信号是脉冲序列形式的离散系统,称为采样控制系统或脉冲控制系统。包括 数字控制系统:把数字序列形式的离散系统,称为数字控制系统或计算机控制系统。注:注:在理
4、想采样及忽略量化误差情况下,数字控制系统近似于采样控制 系统,将它们统称为离散系统。这使得采样控制系统与数字控制系统的分析与校正在理论上统一。2022-10-657.1.1 7.1.1 采样控制系统采样控制系统一般来说,采样控制系统是对传感器所采集的连续信号在某些规定的时间上取值,然后通过对这些值的比较、计算和输出,来达到控制目标的系统。采样控制系统结构构成:主要由采样器、数字控制器、保持器、执行器、被控对象和测量变送器构成,如图7-1所示。图7-1 采样控制系统方框图2022-10-661.1.信号采样信号采样()e t()e t在采样控制系统中,把连续信号转化为脉冲序列的过程称为采样。如图
5、7-2所示。图7-2 采样过程T1/sfT22/ssfT采样周期采样周期采样角频率采样角频率 采样频率采样频率2.2.信号复现信号复现在采样控制系统中,把脉冲序列转变为连续信号的过程称为信号复现。实现复现过程的装置称为保持器。()e t()he t最简单的保持器是零阶保持器,它将脉冲序列复现为阶梯信号 如图7-3所示。2022-10-67图7-3 信号复现过程7.1.2 7.1.2 数字控制系统数字控制系统数字控制系统是一种以数字计算机为控制器去控制具有连续工作状态的被控对象的闭环控制系统。其原理方框图如图7-4所示。图7-4 数字控制系统方框图2022-10-68过程分析:A/D转换器将连续
6、信号转换成数字序列,经数字控制器处理后生 成离散控制信号,再通过D/A转换器转换成连续控制信号作用于 被控对象。1.A/D1.A/D转换器转换器A/D转换器是把连续的模拟信号转换为离散数字信号的装置。A/D转换包括采样过程和量化过程。T()e t()e t采样过程 是每隔 秒对连续信号进行一次采样,得到采样信号 。量化过程 是计算机中任何数值都用二进制表示,因此,幅值上连续的离散信号()e t 信号,此过程称为量化过程。须经过编码表示成最小二进制数的整数倍,成为离散数字2022-10-692.D/A2.D/A转换器转换器D/A转换器是把离散的数字信号转换为连续模拟信号的装置。包括解码过程和复现
7、过程。解码过程就是把离散数字信号转换为离散的模拟信号。复现过程就是通过保持器,将离散模拟信号复现为连续模拟信号。2022-10-6107.2 7.2 信号的采样与保持信号的采样与保持7.2.1 7.2.1 采样过程及其数学描述采样过程及其数学描述由图7-2可写出脉冲序列 的表达式为()e t0()(0)1()1()()1()1()()1()1()()1()1()ke tette TtTtTe kTtkTtkTe kTtkTtkT式(7-2)也可写作 0()()()ke te ttkT()e t 因此,采样过程从物理意义上可以看作是脉冲调制过程。此时,采样开关相当于理想单位脉冲发生器的作用,通过
8、它将连续信号 调制成脉冲序列。2022-10-611从频率特性的角度看:假设连续信号()e t的频率特性为()()j te je t edt,该信号的频谱|()|E j如图7-5所示。图7-5 连续信号频谱离散信号()e t的拉氏变换为1()()skEsE sjkT1|()|()|skEjE jjkT()e t的傅立叶变换为 1()()skEjE jjkT即 2022-10-612如图7-6所示max2s图7-6 采样信号频谱(时)max2s图7-7 采样信号频谱(时)2022-10-6137.2.2 7.2.2 采样定理采样定理香农采样定理:()e tmax如果被采样的连续信号的频谱具有有限
9、带宽,且频谱的最高角频率为,则只要采样角频率s满足:max2s或采样频率sf满足:max2sff则通过理想滤波器,由采样得到的离散信号能够可以不失真地恢复为原连续信号。采样定理给出了采样频率下限的选取规则,对于采样频率的上限,要依据易实现性和抗干扰性来统一确定。2022-10-6147.2.3 7.2.3 信号恢复信号恢复1.1.零阶保持器零阶保持器零阶保持器是工程实践上最常用的一种保持器,它把采样时刻kT的采样值恒定不变地保持到下一个采样时刻 。(1)kT如图7-8所示。图7-8 零阶保持器的输出波形()he t()e t保持器的输出与连续输入信号之间的关系为0()()1()1(1)hke
10、kTe kTtkTtkT2022-10-615kT()tkT对于零阶保持器,在任意时刻输入单位脉冲信号,其单位脉冲响应为一个幅值为1的矩形方波。如图7-9所示。图7-9 零阶保持器的时域特性零阶保持器的频率特性为2sin12()2Tj TjhTeGjTeTjsin(0)2|()|0()2hsTTGjTTk0(0)()()2hsTGjkk 2022-10-616|()|hGj()hGj绘制幅频特性和相频特性曲线,如图7-10所示。图7-10 零阶保持器的幅频特性与相频特性2.2.一阶保持器一阶保持器其外推关系式为()(1)()()e kTe kTe kTte kTtT一阶保持器的输出波形如图7-
11、11所示。2022-10-617图7-11 一阶保持器的输出波形一阶保持器的脉冲响应如图7-12所示。图7-12 一阶保持器的脉冲响应2022-10-618 一阶保持器的幅频特性和相频特性曲线,如图7-13所示。图7-13 一阶保持器的幅频特性与相频特性2022-10-6197.3 z7.3 z变换理论变换理论Z变换是从拉氏变换引申出来的一种变换方法,是研究线性离散系统的重要数学工具,因此又称为离散拉氏变换。7.3.1 z7.3.1 z变换的定义变换的定义()e tT*()e t连续时间函数经采样周期为的采样开关后,得到采样信号 ,即*0()()()ke te kTtkT进行拉氏变换可得*0(
12、)()()kTskEsL e te kT e,引入一个新变量 ,即zTszez()E z得到以为变量的函数,即*0()()()kkE zEse kT z()E z*()e t式中,称为离散信号的z变换,记为*()e t*()()E zZ e t2022-10-6207.3.2 z7.3.2 z变换的方法变换的方法常用的求取离散函数的z变换方法有级数求和法、部分分式法和留数计算法。1.1.级数求和法级数求和法()e tT根据z变换的定义,将连续信号按周期进行采样,级数展开可得 12()(0)()(2)()nE zee T zeT ze nT z例7-1 试求单位阶跃函数()1()e tt的z变换
13、。解:因为单位阶跃函数()1()e tt在所有采样时刻上的采样值均为1,即()1,0,1,2,e kTk则*00()1()()()kke tkTtkTtkT为一理想单位脉冲序列,因此120()1()1kkkE zkT zzzz 上式为一个等比级数,若满足1|1z,则级数收敛,可写成如下闭合形式:11()11zE zzz2022-10-621()ate te(0)a 例7-2 试求衰减指数函数的z变换。解:将ate在各采样时刻上的采样值代入展开式,得1220()1akTkaTaTkaTkkE zezezezez 1|1aTez|1aTez 若,即,则可写成闭合形式:11()1aTaTzE zez
14、ze例7-3 试求函数()ke ta的z变换。解:将 在各采样时刻上的采样值代入展开式,得ka1220()1kkkkkE za zaza za z 若|1az,则可写成闭合形式:11()1zE zazza2022-10-6222.2.部分分式法部分分式法()e t()E s将连续时间函数的拉氏变换 展开成部分分式之和的形式,即12112()nniiniAAAAE sspspspsp由拉氏反变换可得原时间函数:1()inp tiie tAe直接对上式进行z变换,得1()inip TiAzE zze()e t()()aE ss sa例7-4 已知连续时间函数的拉氏变换为,试求其z变换。()E s1
15、1()()aE ss sassa解:首先将展开成部分分式的形式:然后对上式逐项求取拉氏反变换,得()1()ate tte2022-10-623根据求得的时间函数再逐项写出相应的z变换,得2(1)()1(1)aTaTaTaTzzzeE zzzezeze例7-5 利用部分分式法求取正弦函数 的z变换。解:已知正弦函数 的拉氏变换为 ,将其分解成部分分式之和的 形式,得 利用拉氏反变换求出 的原时间函数为 ,利用已知的指数函 数z变换公式可求得相应的z变换,即sin tsint22s2211sin2()2()Ltsj sjj sj 1sj()jte 2sinsin2()2()(2cos)1j Tj
16、TzzzTZtj zej zezT z 2022-10-6243.3.留数计算法留数计算法若已知连续时间函数 的拉氏变换 及其全部极点,则 的z变换 可通过留数计算求得。由拉氏反变换可得 采样后,其采样值为而 的z变换为 最后得 若满足 ,则上式可写为 由此可通过拉氏变换直接求相应的z变换函数。应用留数定理计算上式中的积 分,可得()e t()E s()e t()E z1()()2cjstcje tE s e dsj 1()(),0,1,2,2cjkTscje kTE s eds kj 0()()kkE ze kT z()e kT101()()()2cjTskcjkE zE se zdsj|T
17、sze1()()2cjTscjE s zE zdsjze 1()()inTsisszE sE zresze2022-10-625(1)若 为 的单极点,则(2)若 为 的 重极点,则例7-6 试用留数法求取拉氏变换为 的连续时间函数e(t)的z变换。解:由题意可知,的极点均为单极点,即 ,。可计 算 ,即is()E s()()()iiiTsTssssszE szE sressszezeis()E s11()1()()(1)!iiiiirrirTsTsssssizE sdzE sressszerdszeir()()aE ss sa10s 2sa,()E z00()()()()()()1TsTss
18、saTsTsssaaTazazE zresress sa zes sa zeazazssas sa zes sa zezzzze()as sa2022-10-626例7-7 已知 ,求 的z变换。解:根据题意可知,其极点为重极点,即 ,。可计算 ,即()ate tte()e t21()()E ssa1sa 12r ()E z22211()()(2 1)!()()TssaaTaTdzE zsadssazeTezze常用函数的z变换及相应的拉氏变换见表7-1。()t()as sa1ate(1)(1)()aTaTezzzekTse()tkTkz22ssin t2sin2 cos1zTzzT1s1()
19、t1zz 22sscos t2(cos)2 cos1z zTzzT11()E s()e t()E z()E s()e t()E z表7-1 z变换表2022-10-62721st2(1)Tzz 21()saatTe2()aTaTTzeze32s2t23(1)(1)T z zz22()sasinatet22sin2cosaTaTaTzeTzzeTe1saateaTzze22()sasacosatet222cos2cosaTaTaTzzeTzzeTe7.3.3 z7.3.3 z变换的性质变换的性质1、线性定理 设连续函数 、的z变换分别为 、,为常数,则有 2、时移定理 若函数 的z变换为 ,则有
20、3、初值定理 若函数 的z变换为 ,且 时,则有1()e t2()e t1()E z2()Ez,a b1212()()()()Z ae tbe taE zbE z()e t()E z()()kz e tkTzE z10()()()kknnz e tkTzE ze kT z()e t()E z0t()0e t 0(0)lim()lim()tzee tE z2022-10-6284、终值定理 若函数 的z变换为 ,且 不含有 的二重以上的极点,以及 的 极点均位于z平面的单位圆内,则有7.3.4 z7.3.4 z反变换反变换 已知z变换表达式 ,求相应离散序列 的过程,称为z反变换,记为当 时,信
21、号序列是单边的,对单边序列常用的z反变换法有部分分式法、幂级数法和反演积分法。1.1.部分分式法部分分式法 部分分式法又称查表法,根据已知的 ,通过查z变换表找出相应的 或 。()e t()E z()E z1z()E z1()lim()lim(1)()tzee tzE z()E z()e kT1()()e kTZE z0n()0e kT()E z()e t()e kT2022-10-629例7-8 已知 ,试用部分分式法求z反变换。解:首先展开成如下部分分式形式:由此可得 由表7-1查得 因此2.2.幂级数法(长除法)幂级数法(长除法)Z变换函数可以直接通过长除法得到一个无穷项幂级数的展开式。
22、根据的系 数便可以得出时间序列的值。10()(1)(2)zE zzz()101010(1)(2)12E zzzzzz1010()12zzE zzz111zZz122kzZz()10(12)ke kT 0,1,2,k 0()10(12)()kke ttkT 2022-10-630例7-9 设 ,试用幂级数法求 。解:根据题意可得 利用长除法,得 由此可得3.3.反演积分法反演积分法 反演积分法又称留数法。的幂级数展开形式为:则有反演积分公式 式中,表示函数 在极点 处的留数。留数计算方法如下:若 为单极点,则 若 为m阶重极点,则 10()(1)(3)zE zzz()e t210()43zE z
23、zz01234()01040130400E zzzzzz()10()40(2)130(3)400(4)e ttTtTtTtT()E z0()()kkE ze kT z1111()()Re ()2ikkkzzie kTE z zdzs E z zj 1Re ()ikzzs E z z1()kE z ziziz(0,1,2,)ik11Re ()lim()()iikkzzizzs E z zzz E z z11111Re ()()()(1)!iimkmkzzimz zds E z zzzE z zmdziz2022-10-631例7-10 设 ,试用反演积分法求 。解:根据前式,可得例7-11 设
24、,试用反演积分法求z反变换。解:根据题意可知,该函数有一个单极点,;有一个二重极点,得10()(1)(3)zE zzz()e kT11310()Re(1)(3)1010(1)(3)(1)(3)(1)(3)55 30,1,2,ikzzkkzzkze kTszzzzzzzzzzzk 32()(1)(5)zE zzz11z 25z 11111Re()lim(1)()16kkzzzzs E z zzE z z2111152 1212 1511Re()(5)()(1)!1(5)()(2 1)!(43)516mkmkmzzzkzkds E z zzE z zmdzdzE z zdzk2022-10-632
25、因此 相应的采样信号为总结:上述三种方法中,幂级数法最简单,但得到的z反变换是开式的,因此难以应用。而部分分式法和反演积分法得到的z反变换均为闭合形式。7.47.4离散控制系统的数学模型离散控制系统的数学模型7.4.1 7.4.1 线性常系数差分方程线性常系数差分方程 对于线性定常离散控制系统,一般可用n阶后向差分方程描述,即 同理,线性定常离散系统也可用n阶前向差分方程描述,即111(43)51(43)5()161616kkkke kT1001(43)5()()()()16()11(1)86(2)kkkke te kTtkTtkTttt11()()()nmiiijy ka y kibr kj
26、 11()()()nmiiijy kna y knibr kmj 2022-10-633工程上对于线性定常系数差分方程的求解方法通常有迭代法和z变换法。1、迭代法 若已知差分方程,并且给定输出序列的初值,则可利用递推关系,通过计算 机迭代计算出输出序列。例7-12 已知二阶后向差分方程 ,其中,输入序 列 ,初始条件为 。试用迭代法求解输出序列 。解:根据初始条件和递推关系,得()2()5(1)6(2)y kr ky ky k()1,(0)r kk(0)0,(1)1yy(),0,1,2,y kk(0)0y(1)1y(2)2(2)5(1)6(0)7yryy(3)2(3)5(2)6(1)31yry
27、y(4)2(4)5(3)6(2)115yryy2022-10-6342、z变换法 用z变换法求解线性定常系数差分方程与用拉氏变换法求解微分方程相类似。其实质是将差分方程转化为代数方程,通过代数运算及查表的方法来求出输出序列 。z变换法求解差分方程的一般步骤如下:(1)利用z变换的实数位移定理对差分方程两端取z变换,并代入相应的初始 条件,得到以z为变量的代数方程;(2)求出代数方程的解 ;(3)取z反变换,从而求得输出序列 。例7-13 试用z变换法求解下列二阶前向差分方程 其中,初始条件为 。解:对方程两端取z变换,得 即 代入初始条件,得 即 ()y k()Y z()Y z()y k(2)
28、3(1)2()0y ky ky k(0)0,(1)1yy22()(0)(1)3()3(0)2()0z Y zz yzyzY zzyY z22(32)()(0)(1)3(0)zzY zyzyyz2(32)()zzY zz2()32(1)(2)zzY zzzzz2022-10-635利用反演积分法求出z反变换,得7.4.2 7.4.2 脉冲传递函数脉冲传递函数 1.1.脉冲传递函数定义脉冲传递函数定义 在线性定常离散控制系统中,当初始条件为零时,系统离散输出信号的z变换与离散输入信号的z变换之比,称为线性定常离散控制系统的脉冲传递函数。设离散系统如图7-14所示。根据线性定常离散系统的脉冲传递函数
29、定义得:()120,1,2,ky kk 0()(12)()kky ttkT 00()()()()()nnnny nT zY zG zR zr nT z图7-14 开环采样系统2022-10-636输出是连续信号的情况,如图7-15所示。可以在系统输出端虚设一个开关,如图7-15中虚线所示。2 2、脉冲传递函数的性质、脉冲传递函数的性质(1)脉冲传递函数是复变量z的复函数(一般是有理分式)。(2)脉冲传递函数只与系统自身的结构、参数有关。(3)系统的脉冲传递函数与系统的差分方程有直接关系。(4)系统的脉冲传递函数是系统的单位脉冲响应序列的变换。(5)系统的脉冲传递函数在z平面上有对应的零、极点分
30、布。图7-15 虚设采样开关的开环采样系统2022-10-6373.3.脉冲传递函数的求法脉冲传递函数的求法 传递函数 的拉氏反变换是单位脉冲函数 ,将 离散化得到脉冲响应 序列 ,将 进行z变换即可得到 。这一变换过程可表示为例7-14 已知采样系统结构图如图7-16所示。(1)求系统的脉冲传递函数;(2)写出系统的差分方程。()G s()k t()k t()k nT()k nT()G z1()()()()()()()G sLG sk tk tk nTZ k nTG z离散化图7-16 采样系统结构图2022-10-638解:(1)首先将E(z)/z分解为部分分式形式查变换表得4.4.开环系
31、统脉冲传递函数开环系统脉冲传递函数(1 1)串联环节之间有采样开关时)串联环节之间有采样开关时 设开环离散系统如图7-17所示,在两个串联连续环节 和 之间,有理想采样开关。根据脉冲传递函数定义,有由此可得1()G s2()G s111()(0.11)10G sssss101010(1)()1(1)()TTTzzzeG zzzezze2()()()Y zG z C z1()()()C zG z R z12()()()()()Y zG zG z G zR z2022-10-639(2 2)串联环节之间无采样开关时)串联环节之间无采样开关时 设开环离散系统如图7-18所示,在两个串联连续环节G1(
32、s)和G2(s)之间没有理想采样开关。此时系统的传递函数为 上式作为一个整体进行z变换,由脉冲传递函数定义得 图7-17 串联环节之间有理想采样开关的开环采样系统12()()()G sG s G s12()()()()Y zG zGG zR z图7-18 环节之间无理想采样开关的开环采样系统2022-10-640例7-15 已知开环采样系统分别如图7-17、图7-18所示,其中 ,输入信号r(t)=1(t),试求开环脉冲传递函数和输出的z变换。解:根据z变换表,可知 对于如图7-17所示的采样系统 因此 对于如图7-18所示的采样系统1()/G sb s2()/()G sasa()1zR zz
33、1()1bbzG zZsz 2()aTaazGzZsaze212()()()(1)()aTabzG zG z G zzze32()()()(1)()aTabzY zG z R zzze12()()()abG s G ss sa12(1)()()()(1)()aTaTabbzeG zGG zZs sazze22(1)()()()(1)()aTaTbzeY zG z R zzze2022-10-6415.5.闭环系统脉冲传递函数闭环系统脉冲传递函数 在采样系统中,由于采样开关在系统中的位置不同,因此闭环离散系统结构 图的形式并不唯一。图7-19所示的是常见的系统结构图。可得该闭环离散系统脉冲传递函
34、数为闭环离散系统的误差传递函数为 ,其中 为开环离散系统的脉冲传递函数。图7-19 闭环离散系统结构图()()()()1()Y zG zW zR zGH z()1()()1()eE zW zR zGH z()1()0D zGH z 其中()GH z2022-10-642例7-16 设闭环采样系统结构图如图7-20所示,试求其闭环脉冲传递函数。解:由图可得 对D(s)离散化,有 则 又有 离散化得 即 同时,输出信号的采样拉氏变换为 进行z变换,即得图7-20 闭环采样系统结构图2()()()Y sG s Ds1()()()D sG s Es1()()()DsGs Es21()()()()Y s
35、G s Gs Es21()()()()()()()()()E sR sH s Y sR sH s G s Gs Es21()()()()()EsR sHGs Gs Es12()()1()()R sEsG s HGs212112()()()()()()()1()()Gs G s R sCsGs G s EsG s HGs1212()()()()()1()()G z GzC zW zR zGz HG z2022-10-6437.57.5离散控制系统的分析与设计离散控制系统的分析与设计7.5.1 7.5.1 离散控制系统的稳定性离散控制系统的稳定性 在线性连续系统中,判别系统稳定性是根据特征方程的根
36、在s平面的分布位 置确定的。若系统特征方程的根全部都位于s平面的左半平面,则系统稳定。1.s1.s平面与平面与z z平面之间的映射关系平面之间的映射关系 在z变换定义中,给出了s域到z域的映射关系,即 ,其中,s是复变量,可表示为 ,映射到z域则为写成极坐标形式为如图7-21所示。由此可见,可以把s平面划分为无穷多条平行于实轴的周期带,其中从 的周期带称为主带,其余的周期带称为辅带。Tszesj()jTTj Tzeee|jzz e/2/2ss图7-21 s平面与z平面之间的映射2022-10-6442.2.线性定常离散系统稳定的充要条件线性定常离散系统稳定的充要条件 离散系统稳定性的概念与连续
37、系统相同。如果一个线性定常离散系统的脉冲响应序列趋于零,则系统是稳定的,否则系统不稳定。由s域到z域的映射关系及连续系统的稳定判据,可知:(1)s左半平面映射为z平面单位圆内的区域,对应稳定区域;(2)s右半平面映射为z平面单位圆外的区域,对应不稳定区域;(3)s平面上的虚轴,映射为z平面的单位圆周,对应临界稳定情况,属不稳定。线性定常离散系统稳定的充要条件是:系统闭环脉冲传递函数的全部极点均分布在z平面上以原点为圆心的单位圆内,或者系统所有特征根的模均小于1。3.3.离散系统的稳定性判据离散系统的稳定性判据在z域中不能直接套用劳斯判据,必须引入z域到 域的线性变换,使z平面单位圆内的区域,映
38、射成s平面上的左半平面,这种新的坐标变换称为 变换。令复变量zxjyujv2022-10-645得显然由于上式的分母 始终为正,因此可得(1)等价于 ,表明平面的虚轴对应于z平面的单位圆周;(2)等价于 ,表明左半平面对应于z平面单位圆内的区域;(3)等价于 ,表明右半平面对应于z平面单位圆外的区域。z平面和平面的这种对应关系,如图7-22所示。222222()12(1)(1)xyyujvjxyxy2222()1(1)xyuxy22(1)xy0u 221xy0u 221xy0u 221xy图7-22 z平面与平面的对应关系2022-10-646例7-17 设闭环离散系统如图7-23所示,其中采
39、样周期T=0.1s,试求系统稳定时K的临界值。解:首先求系统开环传递函数的z变换 由此可得,闭环特征方程为 令 ,得 化简后,得 列出劳斯表:从劳斯表第一列系数可知,为保证系统稳定,则需 ,所以,系统稳定的临界增益值 。图7-23 闭环离散系统方框图20.632()(0.11)1.3680.368KKzG zZsszz21()(0.6321.368)0.3680G zzKz11z211(0.6321.368)0.368011K20.6321.2642.7360.6320KK2100.6322.7360.6321.26402.7360.632KKK04.33K4.33K 2022-10-6477
40、.5.2 7.5.2 离散控制系统的瞬态响应离散控制系统的瞬态响应 离散系统的动态特性,是通过在外部输入信号作用下的输出曲线来反映的。设闭环离散系统的闭环脉冲传递函数为当 时,有经部分分式展开后,可得其中,1101010101()()()()()mjmmjmnnnniizzb zb zbbM zW zN za za zaazp()1()r tt()()()()()1M zzY zW z R zN zz1(1)()()()(1)1niiic zMzY zW z R zNzzp()(1)()iiiiM pcpNp()()iizpdN zNpdz2022-10-648求z反变换,即可求得系统在采样时
41、刻的输出离散值为 式中,等号右边第一项为y(kT)的稳态分量;第二项为y(kT)的瞬态分量。1(1)()0,1,2,(1)nkiiiMy kTc pkN图7-24 闭环极点分布与瞬态分量的关系图2022-10-6491 正实轴上闭环极点当时,极点位于单位圆内的正实轴上,对应的瞬态分量指数衰减,且越靠近原点,衰减越快;当时,极点位于单位圆外的正实轴上,对应的瞬态分量指数发散,且越靠近原点,发散越快。系统不稳定;当时,极点位于单位圆上的正实轴上,对应的瞬态分量为等幅脉冲序列,系统临界稳定。2 负实轴上闭环极点当时,极点位于单位圆内的负实轴上,且当 为偶数时,为正值,当 为奇数时,为负值。01ip1
42、ip 1ip 10ip kkiic pkkiic p2022-10-650当 时,极点位于单位圆外的负实轴上,对应的瞬态分量为正、负交替的发散振荡脉冲序列;当 时,极点位于单位圆上的负实轴上,对应的瞬态分量为正、负交替的等幅振荡脉冲序列。3 Z平面上的闭环共轭复数极点当时,极点位于单位圆外的z平面上,对应的瞬态分量发散振荡,系统不稳定;当时,极点位于单位圆内的z平面上,对应的瞬态分量衰减振荡。且pi 越小,即复极点越靠近原点,衰减振荡越快;当时,极点位于单位圆上,对应的瞬态分量等幅振荡,系统处于临界稳定。1ip 1ip|1ip|1ip|1ip 2022-10-651综上所述,闭环脉冲传递函数的
43、闭环极点在单位圆内,对应的瞬态分量均为衰减,故系统稳定;当闭环极点在单位圆上或单位圆外,对应的瞬态分量等幅振荡或发散振荡,故系统不稳定。为了使离散系统具有良好的动态过程,闭环极点应尽量避免配置在单位圆的左半部,尤其不要靠近负实轴,以免产生强烈振荡。闭环极点应最好配置在单位圆的右半部,而且靠近原点的地方,这样,系统的动态过程进行较快,因而系统的快速性较好。2022-10-6527.5.3 7.5.3 离散控制系统的稳态误差离散控制系统的稳态误差1、一般方法(利用终值定理)设单位反馈离散系统结构图如图7-25所示,为系统采样信号,其z变换为()e t()()()()()()E zR zY zR z
44、G z E z图7-25 单位反馈离散系统2022-10-653系统脉冲传递函数为如果 的极点全部位于z平面上的单位圆内,即离散系统是稳定的,则可用z变换的终值定理求出系统的稳态误差。线性定常离散系统的稳态误差,与系统本身的结构和参数有关,与输入序列的形式及幅值有关,而且与采样周期的选取也有关。()1()()1()eE zW zR zG z()eW z*111(1)()lim()lim(1)()lim1()sstzzzR zee tzE zzG z2022-10-654例例7-187-18 设离散系统如图7-25所示,其中,输入连续信号分别为 和 ,试求离散系统的稳态误差。解:解:首先 的z变
45、换为系统的误差脉冲传递函数为闭环极点为 ,均位于z平面的单位圆内,因此可以应用终值定理方法求解稳态误差。当 时,()1/(1)G ss s1Ts1()tt()G s11(1)()()(1)()zeG zZ G szze21(1)(0.368)()1()0.7360.368ezzW zG zzz1,20.3680.482zj()1()r tt()/(1)R zzz21(1)(0.368)lim00.7360.368sszzzezz2022-10-655当时,2、稳态误差系数法在连续系统中,我们把开环传递函数 具有 的极点数作为划分系统型别的标准。在离散系统中,对应把开环脉冲传递函数 具有 的极点
46、数,作为划分离散系统型别的标准,类似把 中 的系统,分别称为0型、型和型离散系统。21(0.368)lim10.7360.368sszT zeTzz()r tt2()/(1)R zTzz()G s0s()G z1z()G z0,1,2n 2022-10-656(1 1)阶跃输入时的稳态误差:)阶跃输入时的稳态误差:当系统输入为阶跃函数 时,其z变换函数稳态误差为 其中 称为离散系统的稳态位置误差系数。对0型离散系统,不趋于 ,从而稳态误差 ;对型或型以上的离散系统,因而稳态误差 。()1zR zz()1()r tt11111lim1()1 lim()sszpzeG zG zK1lim()pzK
47、G zpK0sse pK 0sse 2022-10-657(2 2)斜坡输入时的稳态误差:)斜坡输入时的稳态误差:当系统输入为斜坡函数 时,其z变换函数为因而稳态误差为其中 称为离散系统的稳态速度误差系数。在斜坡输入条件下,0型系统的 ,所以,;型系统的 为有限值,存在常值速度误差;型和型以上系统 ,稳态误差为零。()r tt2()(1)TzR zz11lim(1)(1()lim(1)()sszvzTTTezG zzG zK1lim(1)()vzKzG z0vK sse vKvk 2022-10-658(3 3)加速度输入时的稳态误差:)加速度输入时的稳态误差:当系统输入为加速度函数 时,其z
48、变换函数为 因而稳态误差为 式中称为离散系统的稳态加速度误差系数。加速度输入条件下,由于0型及型系统的 ,所以 ;型系统的 为常数,则加速度误差是非零常数。2()/2r tt23(1)()2(1)T z zR zz2222211(1)lim2(1)1()lim(1)()sszazTzTTezG zzG zK21lim(1)()azKzG z0aK sse aK2022-10-659表7-2 单位反馈离散系统的稳态误差()1()r tt()r tt2()/2r tt1/pK/vT K2/aTK系统型别0型型0型00加速度误差速度误差位置误差2022-10-660例例7-197-19 设离散系统如
49、图7-25所示,其中,输入连续信号 分别为 ,和 ,试利用稳态误差系数法求该离散系统的 、以及稳态误差 。解:解:(1)输入信号 时(2)输入信号 时(3)输入信号 时()1/(1)G ss s1Ts()r t1()tt212tpKvKaKsse1111(1)lim()lim(1)()pzzzeKG zzze 101sspeK1111(1)lim(1)()lim(1)1(1)()vzzzeKzG zzzze1ssvTeK122111(1)lim(1)()lim(1)0(1)()azzzeKzG zzzze2ssaTeK()1()r tt()r tt21()2r tt2022-10-6617.5
50、.4 7.5.4 最少拍系统及其设计最少拍系统及其设计最少拍系统是指:在典型输入作用下,在各采样时刻上无稳态误差,且瞬态过程能在有限个采样周期内结束,完全跟踪控制信号的系统。在离散系统中,一个采样周期称为一拍。所谓最少拍系统,是指对于典型输入信号具有最快的响应速度,能在有限的几拍之内结束过渡过程,且在过渡过程结束后,在采样时刻上稳态误差为零的系统。最少拍系统也称为小调节时间系统或最快响应系统。2022-10-662最少拍系统的设计原则是:若系统被控对象 无延迟,且在z平面单位圆上及单位圆外无零极点,需选择闭环脉冲传递函数 ,使系统在典型输入作用下,经最少采样周期后,能使输出序列在各采样时刻的稳