《自动控制原理》第七章 离散控制系统复习课程.ppt

自动控制原理第七章 离散控制系统教学重点l了解线性离散系统的基本概念和基本定理，把握线性连续系统与线性离散系统的区别与联系；l熟练掌握Z变换的方法、Z变换的性质和Z反变换；l了解差分方程的定义，掌握差分方程的解法；l了解脉冲传递函数的定义，熟练掌握开环与闭环系统脉冲传递函数的计算方法；l与线性连续系统相对应，掌握线性离散系统的时域和频域分析方法和原则。7.1.1 7.1.1 采样控制系统采样控制系统一般来说，采样控制系统是对传感器所采集的连续信号在某些规定的时间上取值，然后通过对这些值的比较、计算和输出，来达到控制目标的系统。采样控制系统结构构成：主要由采样器、数字控制器、保持器、执行器、被控对象和测量变送器构成，如图7-1所示。图7-1 采样控制系统方框图1.1.信号采样信号采样在采样控制系统中，把连续信号转化为脉冲序列的过程称为采样。如图7-2所示。图7-2 采样过程采采样样周期周期采采样样角角频频率率 采采样频样频率率2.2.信号复现信号复现在采样控制系统中，把脉冲序列转变为连续信号的过程称为信号复现。实现复现过程的装置称为保持器。最简单的保持器是零阶保持器，它将脉冲序列复现为阶梯信号 如图7-3所示。图7-3 信号复现过程7.1.2 7.1.2 数字控制系统数字控制系统数字控制系统是一种以数字计算机为控制器去控制具有连续工作状态的被控对象的闭环控制系统。其原理方框图如图7-4所示。图7-4 数字控制系统方框图过程分析：A/D转换器将连续信号转换成数字序列，经数字控制器处理后生 成离散控制信号，再通过D/A转换器转换成连续控制信号作用于 被控对象。1.A/D1.A/D转换器转换器A/D转换器是把连续的模拟信号转换为离散数字信号的装置。A/D转换包括采样过程和量化过程。采样过程 是每隔 秒对连续信号进行一次采样，得到采样信号 。量化过程 是计算机中任何数值都用二进制表示，因此，幅值上连续的离散信号 信号，此过程称为量化过程。须经过编码表示成最小二进制数的整数倍，成为离散数字2.D/A2.D/A转换器转换器D/A转换器是把离散的数字信号转换为连续模拟信号的装置。包括解码过程和复现过程。解码过程就是把离散数字信号转换为离散的模拟信号。复现过程就是通过保持器，将离散模拟信号复现为连续模拟信号。7.2 7.2 信号的采样与保持信号的采样与保持7.2.1 7.2.1 采样过程及其数学描述采样过程及其数学描述由图7-2可写出脉冲序列 的表达式为式（7-2）也可写作 因此，采样过程从物理意义上可以看作是脉冲调制过程。此时，采样开关相当于理想单位脉冲发生器的作用，通过它将连续信号 调制成脉冲序列。从频率特性的角度看：假设连续信号的频率特性为，该信号的频谱如图7-5所示。图7-5 连续信号频谱离散信号的拉氏变换为的傅立叶变换为 即 如图7-6所示图7-6 采样信号频谱（时）图7-7 采样信号频谱（时）7.2.2 7.2.2 采样定理采样定理香农采样定理：如果被采样的连续信号的频谱具有有限带宽，且频谱的最高角频率为，则只要采样角频率满足：或采样频率满足：则通过理想滤波器，由采样得到的离散信号能够可以不失真地恢复为原连续信号。采样定理给出了采样频率下限的选取规则，对于采样频率的上限，要依据易实现性和抗干扰性来统一确定。7.2.3 7.2.3 信号恢复信号恢复1.1.零阶保持器零阶保持器零阶保持器是工程实践上最常用的一种保持器，它把采样时刻的采样值恒定不变地保持到下一个采样时刻 。如图7-8所示。图7-8 零阶保持器的输出波形保持器的输出与连续输入信号之间的关系为对于零阶保持器，在任意时刻输入单位脉冲信号,其单位脉冲响应为一个幅值为1的矩形方波。如图7-9所示。图7-9 零阶保持器的时域特性零阶保持器的频率特性为绘制幅频特性和相频特性曲线，如图7-10所示。图7-10 零阶保持器的幅频特性与相频特性2.2.一阶保持器一阶保持器其外推关系式为一阶保持器的输出波形如图7-11所示。图7-11 一阶保持器的输出波形一阶保持器的脉冲响应如图7-12所示。图7-12 一阶保持器的脉冲响应 一阶保持器的幅频特性和相频特性曲线，如图7-13所示。图7-13 一阶保持器的幅频特性与相频特性7.3 z7.3 z变换理论变换理论Z变换是从拉氏变换引申出来的一种变换方法，是研究线性离散系统的重要数学工具，因此又称为离散拉氏变换。7.3.1 z7.3.1 z变换的定义变换的定义连续时间函数经采样周期为的采样开关后，得到采样信号 ，即进行拉氏变换可得，引入一个新变量 ，即得到以为变量的函数，即式中，称为离散信号的z变换，记为7.3.2 z7.3.2 z变换的方法变换的方法常用的求取离散函数的z变换方法有级数求和法、部分分式法和留数计算法。1.1.级数求和法级数求和法根据z变换的定义，将连续信号按周期进行采样，级数展开可得 例7-1 试求单位阶跃函数的z变换。解：因为单位阶跃函数在所有采样时刻上的采样值均为1，即则为一理想单位脉冲序列，因此上式为一个等比级数，若满足，则级数收敛，可写成如下闭合形式：例7-2 试求衰减指数函数的z变换。解：将在各采样时刻上的采样值代入展开式，得若，即，则可写成闭合形式：例7-3 试求函数的z变换。解：将在各采样时刻上的采样值代入展开式，得若，则可写成闭合形式：2.2.部分分式法部分分式法将连续时间函数的拉氏变换 展开成部分分式之和的形式，即由拉氏反变换可得原时间函数：直接对上式进行z变换，得例7-4 已知连续时间函数的拉氏变换为，试求其z变换。解：首先将展开成部分分式的形式：然后对上式逐项求取拉氏反变换，得根据求得的时间函数再逐项写出相应的z变换，得例7-5 利用部分分式法求取正弦函数 的z变换。解：已知正弦函数 的拉氏变换为 ，将其分解成部分分式之和的 形式，得 利用拉氏反变换求出 的原时间函数为 ，利用已知的指数函 数z变换公式可求得相应的z变换，即3.3.留数计算法留数计算法若已知连续时间函数 的拉氏变换 及其全部极点，则 的z变换 可通过留数计算求得。由拉氏反变换可得 采样后，其采样值为而 的z变换为 最后得 若满足 ，则上式可写为 由此可通过拉氏变换直接求相应的z变换函数。应用留数定理计算上式中的积 分，可得（1）若 为 的单极点，则（2）若 为 的 重极点，则 例7-6 试用留数法求取拉氏变换为 的连续时间函数e(t)的z变换。解：由题意可知，的极点均为单极点，即 ，。可计 算 ，即，例7-7 已知 ，求 的z变换。解：根据题意可知，其极点为重极点，即 ，。可计算 ，即 常用函数的z变换及相应的拉氏变换见表7-1。11表7-1 z变换表7.3.3 z7.3.3 z变换的性质变换的性质1、线性定理 设连续函数 、的z变换分别为 、，为常数，则有 2、时移定理 若函数 的z变换为 ，则有3、初值定理 若函数 的z变换为 ，且 时，则有4、终值定理 若函数 的z变换为 ，且 不含有 的二重以上的极点，以及 的 极点均位于z平面的单位圆内，则有7.3.4 z7.3.4 z反变换反变换 已知z变换表达式 ，求相应离散序列 的过程，称为z反变换，记为当 时，信号序列是单边的，对单边序列常用的z反变换法有部分分式法、幂级数法和反演积分法。1.1.部分分式法部分分式法 部分分式法又称查表法，根据已知的 ，通过查z变换表找出相应的 或 。例7-8 已知 ，试用部分分式法求z反变换。解：首先展开成如下部分分式形式：由此可得 由表7-1查得 因此2.2.幂级数法（长除法）幂级数法（长除法）Z变换函数可以直接通过长除法得到一个无穷项幂级数的展开式。根据的系 数便可以得出时间序列的值。例7-9 设 ，试用幂级数法求 。解：根据题意可得 利用长除法，得 由此可得3.3.反演积分法反演积分法 反演积分法又称留数法。的幂级数展开形式为：则有反演积分公式 式中，表示函数 在极点 处的留数。留数计算方法如下：若 为单极点，则 若 为m阶重极点，则 例7-10 设 ，试用反演积分法求 。解：根据前式，可得例7-11 设 ，试用反演积分法求z反变换。解：根据题意可知，该函数有一个单极点，；有一个二重极点，得 因此 相应的采样信号为总结：上述三种方法中，幂级数法最简单，但得到的z反变换是开式的，因此难以应用。而部分分式法和反演积分法得到的z反变换均为闭合形式。7.47.4离散控制系统的数学模型离散控制系统的数学模型7.4.1 7.4.1 线性常系数差分方程线性常系数差分方程 对于线性定常离散控制系统，一般可用n阶后向差分方程描述，即 同理，线性定常离散系统也可用n阶前向差分方程描述，即工程上对于线性定常系数差分方程的求解方法通常有迭代法和z变换法。1、迭代法 若已知差分方程，并且给定输出序列的初值，则可利用递推关系，通过计算 机迭代计算出输出序列。例7-12 已知二阶后向差分方程 ,其中，输入序 列 ，初始条件为 。试用迭代法求解输出序列 。解：根据初始条件和递推关系，得 2、z变换法 用z变换法求解线性定常系数差分方程与用拉氏变换法求解微分方程相类似。其实质是将差分方程转化为代数方程，通过代数运算及查表的方法来求出输出序列 。z变换法求解差分方程的一般步骤如下：(1)利用z变换的实数位移定理对差分方程两端取z变换，并代入相应的初始 条件，得到以z为变量的代数方程；(2)求出代数方程的解 ；(3)取z反变换，从而求得输出序列 。例7-13 试用z变换法求解下列二阶前向差分方程 其中，初始条件为 。解：对方程两端取z变换，得 即 代入初始条件，得 即 利用反演积分法求出z反变换，得7.4.2 7.4.2 脉冲传递函数脉冲传递函数 1.1.脉冲传递函数定义脉冲传递函数定义 在线性定常离散控制系统中，当初始条件为零时，系统离散输出信号的z变换与离散输入信号的z变换之比，称为线性定常离散控制系统的脉冲传递函数。设离散系统如图7-14所示。根据线性定常离散系统的脉冲传递函数定义得：图7-14 开环采样系统输出是连续信号的情况，如图7-15所示。可以在系统输出端虚设一个开关，如图7-15中虚线所示。2 2、脉冲传递函数的性质、脉冲传递函数的性质（1）脉冲传递函数是复变量z的复函数（一般是有理分式）。（2）脉冲传递函数只与系统自身的结构、参数有关。（3）系统的脉冲传递函数与系统的差分方程有直接关系。（4）系统的脉冲传递函数是系统的单位脉冲响应序列的变换。（5）系统的脉冲传递函数在z平面上有对应的零、极点分布。图7-15 虚设采样开关的开环采样系统3.3.脉冲传递函数的求法脉冲传递函数的求法 传递函数 的拉氏反变换是单位脉冲函数 ，将 离散化得到脉冲响应 序列 ，将 进行z变换即可得到 。这一变换过程可表示为例7-14 已知采样系统结构图如图7-16所示。（1）求系统的脉冲传递函数；（2）写出系统的差分方程。图7-16 采样系统结构图解：（1）首先将E(z)/z分解为部分分式形式查变换表得4.4.开环系统脉冲传递函数开环系统脉冲传递函数（1 1）串联环节之间有采样开关时）串联环节之间有采样开关时 设开环离散系统如图7-17所示，在两个串联连续环节 和 之间，有理想采样开关。根据脉冲传递函数定义，有由此可得（2 2）串联环节之间无采样开关时）串联环节之间无采样开关时 设开环离散系统如图7-18所示，在两个串联连续环节G1(s)和G2(s)之间没有理想采样开关。此时系统的传递函数为 上式作为一个整体进行z变换，由脉冲传递函数定义得 图7-17 串联环节之间有理想采样开关的开环采样系统图7-18 环节之间无理想采样开关的开环采样系统例7-15 已知开环采样系统分别如图7-17、图7-18所示，其中 ，输入信号r(t)=1(t),试求开环脉冲传递函数和输出的z变换。解：根据z变换表，可知 对于如图7-17所示的采样系统 因此 对于如图7-18所示的采样系统5.5.闭环系统脉冲传递函数闭环系统脉冲传递函数 在采样系统中，由于采样开关在系统中的位置不同，因此闭环离散系统结构 图的形式并不唯一。图7-19所示的是常见的系统结构图。可得该闭环离散系统脉冲传递函数为闭环离散系统的误差传递函数为 ,其中 为开环离散系统的脉冲传递函数。图7-19 闭环离散系统结构图其中例7-16 设闭环采样系统结构图如图7-20所示，试求其闭环脉冲传递函数。解：由图可得 对D(s)离散化，有 则 又有 离散化得 即 同时，输出信号的采样拉氏变换为 进行z变换，即得图7-20 闭环采样系统结构图7.57.5离散控制系统的分析与设计离散控制系统的分析与设计7.5.1 7.5.1 离散控制系统的稳定性离散控制系统的稳定性 在线性连续系统中，判别系统稳定性是根据特征方程的根在s平面的分布位 置确定的。若系统特征方程的根全部都位于s平面的左半平面，则系统稳定。1.s1.s平面与平面与z z平面之间的映射关系平面之间的映射关系 在z变换定义中，给出了s域到z域的映射关系，即 ,其中，s是复变量，可表示为 ，映射到z域则为写成极坐标形式为如图7-21所示。由此可见，可以把s平面划分为无穷多条平行于实轴的周期带，其中从 的周期带称为主带，其余的周期带称为辅带。图7-21 s平面与z平面之间的映射2.2.线性定常离散系统稳定的充要条件线性定常离散系统稳定的充要条件 离散系统稳定性的概念与连续系统相同。如果一个线性定常离散系统的脉冲响应序列趋于零，则系统是稳定的，否则系统不稳定。由s域到z域的映射关系及连续系统的稳定判据，可知：（1）s左半平面映射为z平面单位圆内的区域，对应稳定区域；（2）s右半平面映射为z平面单位圆外的区域，对应不稳定区域；（3）s平面上的虚轴，映射为z平面的单位圆周，对应临界稳定情况，属不稳定。线性定常离散系统稳定的充要条件是：系统闭环脉冲传递函数的全部极点均分布在z平面上以原点为圆心的单位圆内，或者系统所有特征根的模均小于1。3.3.离散系统的稳定性判据离散系统的稳定性判据在z域中不能直接套用劳斯判据，必须引入z域到 域的线性变换，使z平面单位圆内的区域，映射成s平面上的左半平面，这种新的坐标变换称为 变换。令复变量得显然由于上式的分母 始终为正，因此可得（1）等价于 ，表明平面的虚轴对应于z平面的单位圆周；（2）等价于 ，表明左半平面对应于z平面单位圆内的区域；（3）等价于 ，表明右半平面对应于z平面单位圆外的区域。z平面和平面的这种对应关系，如图7-22所示。图7-22 z平面与平面的对应关系例7-17 设闭环离散系统如图7-23所示，其中采样周期T=0.1s，试求系统稳定时K的临界值。解：首先求系统开环传递函数的z变换 由此可得，闭环特征方程为 令 ，得 化简后，得 列出劳斯表：从劳斯表第一列系数可知，为保证系统稳定，则需 ，所以，系统稳定的临界增益值 。图7-23 闭环离散系统方框图7.5.2 7.5.2 离散控制系统的瞬态响应离散控制系统的瞬态响应 离散系统的动态特性，是通过在外部输入信号作用下的输出曲线来反映的。设闭环离散系统的闭环脉冲传递函数为当 时，有经部分分式展开后，可得其中，求z反变换，即可求得系统在采样时刻的输出离散值为 式中，等号右边第一项为y(kT)的稳态分量；第二项为y(kT)的瞬态分量。图7-24 闭环极点分布与瞬态分量的关系图1 正实轴上闭环极点当时，极点位于单位圆内的正实轴上，对应的瞬态分量指数衰减，且越靠近原点，衰减越快；当时，极点位于单位圆外的正实轴上，对应的瞬态分量指数发散，且越靠近原点，发散越快。系统不稳定；当时，极点位于单位圆上的正实轴上，对应的瞬态分量为等幅脉冲序列，系统临界稳定。2 负实轴上闭环极点当时，极点位于单位圆内的负实轴上，且当 为偶数时，为正值，当 为奇数时，为负值。当 时，极点位于单位圆外的负实轴上，对应的瞬态分量为正、负交替的发散振荡脉冲序列；当 时，极点位于单位圆上的负实轴上，对应的瞬态分量为正、负交替的等幅振荡脉冲序列。3 Z平面上的闭环共轭复数极点当时，极点位于单位圆外的z平面上，对应的瞬态分量发散振荡，系统不稳定；当时，极点位于单位圆内的z平面上，对应的瞬态分量衰减振荡。且pi 越小，即复极点越靠近原点，衰减振荡越快；当时，极点位于单位圆上，对应的瞬态分量等幅振荡，系统处于临界稳定。综上所述，闭环脉冲传递函数的闭环极点在单位圆内，对应的瞬态分量均为衰减，故系统稳定；当闭环极点在单位圆上或单位圆外，对应的瞬态分量等幅振荡或发散振荡，故系统不稳定。为了使离散系统具有良好的动态过程，闭环极点应尽量避免配置在单位圆的左半部，尤其不要靠近负实轴，以免产生强烈振荡。闭环极点应最好配置在单位圆的右半部，而且靠近原点的地方，这样，系统的动态过程进行较快，因而系统的快速性较好。7.5.3 7.5.3 离散控制系统的稳态误差离散控制系统的稳态误差1、一般方法（利用终值定理）设单位反馈离散系统结构图如图7-25所示，为系统采样信号，其z变换为图7-25 单位反馈离散系统系统脉冲传递函数为如果 的极点全部位于z平面上的单位圆内，即离散系统是稳定的，则可用z变换的终值定理求出系统的稳态误差。线性定常离散系统的稳态误差，与系统本身的结构和参数有关，与输入序列的形式及幅值有关，而且与采样周期的选取也有关。例例7-187-18 设离散系统如图7-25所示，其中，输入连续信号分别为 和 ，试求离散系统的稳态误差。解：解：首先 的z变换为系统的误差脉冲传递函数为闭环极点为 ，均位于z平面的单位圆内，因此可以应用终值定理方法求解稳态误差。当 时，当时，2、稳态误差系数法在连续系统中，我们把开环传递函数 具有 的极点数作为划分系统型别的标准。在离散系统中，对应把开环脉冲传递函数 具有 的极点数，作为划分离散系统型别的标准，类似把 中 的系统，分别称为0型、型和型离散系统。（1 1）阶跃输入时的稳态误差：）阶跃输入时的稳态误差：当系统输入为阶跃函数 时，其z变换函数稳态误差为 其中 称为离散系统的稳态位置误差系数。对0型离散系统，不趋于 ，从而稳态误差 ；对型或型以上的离散系统，因而稳态误差 。（2 2）斜坡输入时的稳态误差：）斜坡输入时的稳态误差：当系统输入为斜坡函数 时，其z变换函数为因而稳态误差为其中 称为离散系统的稳态速度误差系数。在斜坡输入条件下，0型系统的 ，所以，；型系统的 为有限值，存在常值速度误差；型和型以上系统 ，稳态误差为零。（3 3）加速度输入时的稳态误差：）加速度输入时的稳态误差：当系统输入为加速度函数 时，其z变换函数为 因而稳态误差为 式中称为离散系统的稳态加速度误差系数。加速度输入条件下，由于0型及型系统的 ，所以 ；型系统的 为常数，则加速度误差是非零常数。表7-2 单位反馈离散系统的稳态误差系统型别0型型0型00加速度误差速度误差位置误差例例7-197-19 设离散系统如图7-25所示，其中，输入连续信号 分别为 ，和 ，试利用稳态误差系数法求该离散系统的 、以及稳态误差 。解：解：（1）输入信号 时（2）输入信号 时（3）输入信号 时7.5.4 7.5.4 最少拍系统及其设计最少拍系统及其设计最少拍系统是指：在典型输入作用下，在各采样时刻上无稳态误差，且瞬态过程能在有限个采样周期内结束，完全跟踪控制信号的系统。在离散系统中，一个采样周期称为一拍。所谓最少拍系统，是指对于典型输入信号具有最快的响应速度，能在有限的几拍之内结束过渡过程，且在过渡过程结束后，在采样时刻上稳态误差为零的系统。最少拍系统也称为小调节时间系统或最快响应系统。最少拍系统的设计原则是：若系统被控对象 无延迟，且在z平面单位圆上及单位圆外无零极点，需选择闭环脉冲传递函数 ，使系统在典型输入作用下，经最少采样周期后，能使输出序列在各采样时刻的稳态误差为零，达到完全跟踪的目的，从而确定所需的数字控制器的脉冲传递函数 。考虑到零阶保持器的存在，广义被控对象的脉冲传递函数为系统的闭环脉冲传递函数为闭环误差脉冲传递函数为因为系统为单位反馈系统，因此有可得 或典型输入信号z变换的一般形式为 是不包含因子 的 的多项式。表7-4 典型输入信号的最少拍系统设计典型输入闭环脉冲传递函数调节时间数字控制器的脉冲传递函数a）单位阶跃信号 b）单位速度信号 c）单位加速度信号图7-26 典型输入信号的最少拍系统响应例例7-20 7-20 设单位反馈线性定常离散系统的连续部分和零阶保持器的传递函数分别为其中，采样周期T=1s。若要求系统在单位斜坡输入时实现最少拍控制，试求数字控制器脉冲传递函数D(z)。解：解：系统开环传递函数为由于故有根据 ，由表7-3可知最少拍系统应具有的闭环脉冲传递函数和误差脉冲传递函数为的零点 正好可以补偿 在单位圆中的极点 ，已包含 的传递函数延迟 。因此，上述 和 满足消去 中传递延迟 及补偿 在单位圆中极点的 约束条件，可求出最少拍系统的数字控制器脉冲传递函数为7.67.6应用应用MATLABMATLAB分析离散控制系统分析离散控制系统7.6.1 7.6.1 连续系统与离散系统模型之间的转换连续系统与离散系统模型之间的转换在MATLAB中对连续系统的离散化是通过c2dm()函数实现的，c2dm()函数的调用格式为 numd,dend=c2dm(num,den,T,zoh)其中zoh表示零阶保持器；T为采样周期，num为传递函数分子多项式系数；den为传递函数分母多项式系数。表7-5 连续系统模型与离散系统模型转换函数MATLAB函数调用格式函数说明c2dsysd=c2d(sysc,Ts,method)连续时间LTI系统模型转换成离散时间系统模型c2dmnumz,denz=c2dm(num,den,T,method)连续时间LTI传递函数模型转换成离散时间传递函数模型d2csysc=d2c(sys,method)离散时间系统模型转换成连续时间LTI系统模型d2cmA,B,C,D=d2cm(Ad,Bd,Cd,Dd,T,method)1.离散时间状态空间模型转换成连续时间状态空间模型d2dsys=d2d(sysd,T)离散时间模型转换成新的离散时间系统模型c2dtAd,Bd,Cd,Dd=c2dt(A,B,C,T,lambda)具有纯延迟输入的连续状态空间系统模型转换成离散时间状态空间系统模型其中，d表示离散系统，c表示连续系统，T表示采样周期，method表示变换方法，其具体含义如表7-5所示。表7-6 method选项及功能method选项功能zoh零阶保持器foh一阶保持器imp脉冲不变变换方法tustin双线性变换方法prewarp预先转折变换方法matched零极点匹配转换方法例例7-217-21 已知采样系统结构图如图7-27所示，其中采样周期T=1s，试利用MATLAB求开环脉冲传递函数。解：解：如用解析法求解可得图7-27 采样系统结构图利用MATLAB求解，程序如下：num=1;den=1 1 0;T=1;numz,denz=c2dm(num,den,T,zoh);printsys(numz,denz,Z)执行结果为num/den=0.36788 Z+0.26424 Z2-1.3679 Z+0.367887.6.2 7.6.2 离散系统的响应离散系统的响应在MATLAB中，求离散系统的时间响应可运用函数dstep()、dimpulse()、dslim()来实现。其分别用于求采样系统的阶跃、脉冲及任意输入时的响应。1.dstep()1.dstep()函数的调用格式函数的调用格式dstep(num,den,n)其中：num为脉冲传递函数分子多项式系数；den为脉冲传递函数分母多项式系数；n为用户指定的采样点数。例例7-227-22 已知离散系统结构图如图7-28所示，输入为单位阶跃信号，采样周期T=1s，试利用MATLAB求其输出响应。图7-28 系统结构图解：解：由可得利用MATLAB中的dstep()函数可得到输出响应，如图7-29所示。程序如下：num=0.368 0.264;den=1-1 0.632;dstep(num,den)图7-29 离散系统阶跃响应2.dimpulse()2.dimpulse()函数的调用格式函数的调用格式与dstep()函数的调用格式相同。上例中只需将最好一行程序改为dimpulse(num,den)便可得到离散系统的输出响应，结果如图7-30所示。图7-30 离散系统脉冲响应3.dslim()3.dslim()函数的调用格式函数的调用格式dslim(num,den,u)其中：num为脉冲传递函数分子多项式系数；den为脉冲传递函数分母多项式系数；u为输入信号。