离散时间系统与z变换课件

第二章离散时间系统与z变换2.1取样和内插2.2离散时间信号序列2.3离散系统及其普遍关系2.4离散信号的傅氏变换2.5离散信号的z变换2.6单边z变换2.7z变换与傅氏变换的关系2.8系统的时域分析与频域分析第二章离散时间系统与z变换2.1取样和内插2.1取样和内插1.取样 将连续信号变成离散信号有各种取样方法，其中最常用的是等间隔周期取样，即每隔固定时间T取一个信号值，如图2-1所示。其中T称为取样周期，T的倒数称为取样频率或取样率。记为fS=1/T2.1取样和内插1.取样图2-1连续信号的取样图2-1连续信号的取样 取样定理——Shannon定理 任一连续信号xa(t)，设其频谱的最高频率分量为fm，则当对它进行取样时，只要选择取样率等于或大于2fm

，就可以由这个取样序列xa(nT)来唯一准确地恢复xa(t)。 设有一限带信号xa(t)。当|Ω|≤Ω

max，它的付氏变换为Xa(Ω)。将xa(t)乘一取样函数p(t)就得到xa(t)，如图2.2所示。 取样定理——Shannon定理图2-2连续信号取样的数学模型图2-2连续信号取样的数学模型离散时间系统与z变换图2-5取样过程的时域与频域关系图2-5取样过程的时域与频域关系最后需要说明一点：上述取样定理是理想取样，如果取样函数不是单位冲击函数序列，而是窄脉冲函数序列，则如图2-6所示(详细情况请参看相关资料)。最后需要说明一点：上述取样定理是理想取样，如果取样函数不是单图2-6理想取样和非理想取样的比较图2-6理想取样和非理想取样的比较2.内插 用大于奈奎斯特取样频率取样限带信号xa(t)，则被取样信号xa(t)通过理想低通滤波器，只要其截止频率Ωc满足Ωmax≤Ωc≤(Ωs－Ωmax)时，就可以恢复出原来信号，如图2-8所示。2.内插图2-8取样信号的恢复与理想低通滤波器的传输函数图2-8取样信号的恢复与理想低通滤波器的传输函数图2-9连续信号的内插表示图2-9连续信号的内插表示图2-10连续信号用三角形内插函数图2-10连续信号用三角形内插函数3.Matlab实现(1)零阶和一阶保持内插(2)三次样条内插3.Matlab实现2.2离散时间信号序列1.离散时间信号 离散时间信号是在离散的时间上取值，在两个取样间隔内数值为零的信号。2.2离散时间信号序列1.离散时间信号2.常用序列(1)单位取样序列 单位取样序列的定义为

其图形如图2-15所示。2.常用序列图2-15单位取样序列图2-15单位取样序列(2)单位阶跃序列 单位阶跃序列的定义为其图形如图2.16所示。(2)单位阶跃序列图2-16单位阶跃序列图2-16单位阶跃序列(3)矩形序列 矩形序列的定义为其图形如图2-17所示。(3)矩形序列图2-17矩形序列图2-17矩形序列(4)实指数序列 实指数序列的定义为x(n)=an其中a为不等于零的任意实数。 图2-18是0＜a＜1的一个实指数序列的图形。(4)实指数序列图2-18实指数序列图2-18实指数序列(5)正弦序列 正弦序列的定义为x(n)=sinnω0其图形如图2-19所示。(5)正弦序列图2-19正弦序列图2-19正弦序列2.3离散系统及其普遍关系1.离散系统的定义 离散系统在数学上定义为将输入序列x(n)映射成输出序列y(n)的惟一性变换或运算。亦即将一个序列变换成另一个序列的系统，记为y(n)=T［x(n)］通常将上式表示成图2-20所示的框图。2.3离散系统及其普遍关系1.离散系统的定义图2-20离散系统的模型图2-20离散系统的模型2.线性非移变系统(1)系统的线性特性 满足叠加原理的系统具有线性特性，即若对两个激励x1(n)和x2(n)有2.线性非移变系统(2)系统的非移变特性 系统的非移变是指系统的参数不随时间而变化。用数学表示为T［x(n－n0)］=y(n－n0)即不管输入信号作用的时间先后，输出信号响应的形状均相同，仅是出现的时间不同，如图2-22所示。(2)系统的非移变特性图2-22离散系统的非移变特性图2-22离散系统的非移变特性(3)线性非移变系统 线性非移变系统就是既满足迭加原理又具有非移变特性的系统，将其描绘如图2-24所示。图2-24线性非移变系统模型(3)线性非移变系统图2-24线性非移变系统模型(4)线性卷积的计算 计算线性卷积有4种方法。 ①利用两个序列的解析式直接计算式(2-34)。 ②利用两个序列的移位求和，即先把一个序列倒置。每次将它向下移一步，求出两序列重叠部分乘积之和。 ③用作图法求。 ④卷积的Matlab实现(4)线性卷积的计算3.系统的稳定性与因果性(1)稳定性 对于一个系统，当输入序列是有界时，其输出也是有界的，则称它是稳定系统。用数学描述则为如果 ｜x(n)｜＜∞对于一切n则 ｜y(n)｜＜∞对于一切n3.系统的稳定性与因果性 因为其中假设｜x(n)｜≤M。 因为2．因果性 一个系统如果其输出变化不会发生在输入变化之前，则称它是因果的。这就是说对于因果系统，如果取n0,当n<n0时，x1(n)=x2(n),则n<n0时，y1(n)=y2(n)。一个线性非移变系统当n<0时的因果充要条件是其单位取样响应等于零,即h(n)=0 n<0 这个充要条件可以从y（n）＝x(n)*h(n)的解析式中导出。2．因果性4.系统的差分方程描述(1)非递归型(FIR) 非递归型因果系统是输出的现在值仅仅取决于输入的现在值与输入的过去值的系统。非递归，即输出对输入无反馈。因此，设在n时刻输入x(n)与输出y(n)的关系为y(n)=f｛……,x(n-1),x(n),x(n+1),……｝若系统是线性非移变的，y(n)可表示为4.系统的差分方程描述(2)递归型(IIR) 递归型因果系统输出的现在值不仅取决于输入的现在值与过去值，还取决于输出的过去值。y(n)=f｛……,x(n-1),x(n),x(n+1),……｝ +g｛……,y(n-1),y(n+1),……｝(2)递归型(IIR) 同理，在系统为线性、非移变、因果时，可推得 同理，在系统为线性、非移变、因果时，可推得2.4离散信号的傅氏变换1.问题的提出 对于连续信号xa(t)与其频谱Xa(Ω)之间存在着傅氏变换关系，如图2-28所示。前边已经讨论了连续信号xa(t)的离散化，即取样的问题，已经知道取样序列的频谱是原信号频谱在Ω轴上的周期延拓，如图2-28(b)所示。2.4离散信号的傅氏变换1.问题的提出图2-28连续和离散信号的傅氏变换图2-28连续和离散信号的傅氏变换2.傅氏变换对的推导 从2.1节的讨论知道，对连续信号在时域内进行取样的结果，是频域内频谱周期的延拓，并且还得到了已取样信号xa(t)在频域内的表示，现重写为2.傅氏变换对的推导离散时间系统与z变换3.线性非移变系统的频率响应 从2.3节的讨论得到线性非移变离散系统的输入输出关系为对上式两边同时进行傅氏变换得3.线性非移变系统的频率响应离散时间系统与z变换离散时间系统与z变换4.离散信号傅氏变换的性质 首先定义两种序列，共轭对称序列与共轭反对称序列。(1)序列的卷积特性 序列的卷积特性是时域内的卷积关系映射到频域内为相乘，即4.离散信号傅氏变换的性质(2)序列傅氏变换的周期性 序列的傅氏变换是ω的周期函数，周期为2π。(3)复序列傅氏变换的对称性(2)序列傅氏变换的周期性离散时间系统与z变换离散时间系统与z变换2.5离散信号的z变换1.z变换的定义及其收敛域 对于一个序列x(n)，其z变换的定义为 其中z为复变量，也可记作Z［x(n)］=X(z)。式(2-49)的定义也称为双边z

变换；相应的还有单边z变换。2.5离散信号的z变换1.z变换的定义及其收敛域对于所有的序列或所有的z值，z变换并不总是收敛的。对于任意给定的序列，使z变换收敛的z值集合称作收敛区域：{Z：X(z)存在}＝收敛区域。其内径R-与外径R+分别取x(n)在n→∞和n→-∞时的形状。z变换收敛域的概念很重要，不同的序列可能有相同的z变换表达式，但是收敛域却不同，所以应该特别注意，只有当z变换的表达式与收敛域都相同时，才能判定两个序列相等。对于所有的序列或所有的z值，z变换并不总是收敛的。对于任意给根据以上讨论，可以概括为(1)对右边序列(n≥0存在)，｜z｜＞R-收敛，且R-是右序列的极点。(2)对左边序列(n＜0存在)，｜z｜＜R+收敛，且R+是左边序列的极点。根据以上讨论，可以概括为(3)若X(z)不只一个极点，则找与收敛域相重的那个极点，对右边序列，最外极点之外的区域为收敛域；对左边序列，最内极点之内的区域为收敛域，如图2-31所示。(4)对双边序列，若在左边序列的收敛域存在重叠部分，则这重叠部分就是它的收敛域。若不存在重迭部分，则z变换不存在。(3)若X(z)不只一个极点，则找与收敛域相重的那个极点，2.系统传递函数H(z)的频域表示 描述线性非移变系统的差分方程为对上式方程两边取z变换为2.系统传递函数H(z)的频域表示离散时间系统与z变换 在Matlab中，reqz函数计算幅度和相位响应，它有如下5种调用方式。(1)［H，w］＝freqz(b，a，N) b和a分别表示分子和分母的系数向量，与filter(b，a，x)函数中的相同。此函数在单位圆上半部上等间隔的计算N点频率响应，返回该系统的N点频率矢量w和N点复数频率响应矢量H。如果N没有说明，则缺省值为512。 在Matlab中，reqz函数计算幅度和相位响应，它有如下(2)［H，w］＝ freqz(b，a，N，’whole’) 在整个单位圆上等间隔的计算N点频率响应。(3)H＝freqz(b，a，w) 它返回矢量w指定的那些频率点上的频率响应，通常在0到π之间。(2)［H，w］＝(4)［H，F］＝freqz(b，a，N，Fs)和［H，F］＝freqz(b，a，N，'whole'，Fs) 给定取样频率Fs，单位为Hz；返回单位为Hz的频率矢量F。(5)H＝freqz(b，a，F，Fs) 给定单位为Hz的取样频率Fs，返回矢量F指定的那些频率点上的复数频率响应，单位也是Hz。(4)［H，F］＝freqz(b，a，N，Fs)和［H，F］3.z反变换 z反变换关系式可以利用柯西积分定理推导出来，柯西定理为式中c是一个逆时针方向环绕原点的围线。3.z反变换(1)幂级数法 幂级数法也就是长除法，对于给定的z变换X(z)，可以根据它的收敛域判定序列x(n)是右边序列还是左边序列，或是双边序列。(1)幂级数法(2)部分分式法 当x(z)序列为有理函数时，可将x(z)写成一个和式为因为，对于右边序列存在如下变换关系(2)部分分式法(3)留数法 我们知道式中c是X(z)收敛域内的积分围线。对于n＞0时，对应右边序列，此时极点在c内，对n＜0时，对应左边序列。此时极点在c外，根据留数定理有(3)留数法(4)Matlab

的实现 在Matlab中，函数residuez

计算有理函数的留数Rj、极点pj和直接项系数Cj，设有多项式如下：B(z)和A(z)分别是分子分母多项式，它们按z-1递增顺序排列。(4)Matlab的实现4.因果非移变系统的稳定性，收敛区与极点的关系 正如在式(2-56)中看到的，线性非移变系统的系统函数H(z)是具有实系数z的有理函数。现在来讨论系统函数的极点与系统稳定性和收敛区的关系，亦即证明系统函数的极点分布将决定系统是否稳定。

4.因果非移变系统的稳定性，收敛区与极点的关系 假设有一N阶因果系统，系统函数为H(z)，为方便起见，设H(z)只有单阶极点，这样系统的单位取样响应由式(2-67)给出 假设有一N阶因果系统，系统函数为H(z)，为方便起见，设H离散时间系统与z变换由以上的讨论清楚看到，因果稳定系统的收敛区域包括单位圆以及以外的整个z平面。因而，因果非移变的稳定系统为(1)极点都在单位圆内；(2)收敛区域为l≤z≤∞。由以上的讨论清楚看到，因果稳定系统的收敛区域包括单位圆以及以5.z变换的性质(1)线性 设X(z)与Y(z)分别是x(n)与y(n)的z变换，即5.z变换的性质(2)序列的移位 设序列x(n)的z变换为Z［x(n)］=X(z) Rx-＜|z|＜Rx+(3)乘以指数序列 如果序列x(n)乘上指数序列an(a可以是复数)，则Z[x(n)]=X(z)Rx_＜|z|＜Rx+

(2)序列的移位(4)X(z)的微分 序列x(n)之z变换的导数乘以(-z)等于x(n)经线性加权后的z变换，即(5)复序列的共轭序列(4)X(z)的微分(6)初值定理 如果n＜0时，x(n)为零，则(7)序列的卷积 如果w(n)是序列x(n)和y(n)的卷积，则w(n)的z

变换是x(n)和y(n)的z

变换的乘积，即w(n)=x(n)*y(n)则(6)初值定理(8)复卷积定理 复卷积定理与序列的卷积是对偶关系。设w(n)=x(n)y(n)则(8)复卷积定理或者式中C1是X（z/v）与Y（v）两者收敛区重叠部分的闭合围线；C2是X(v)与Y（z/v）两者收敛区重叠部分内的闭合围线。或者(9)帕斯维尔定理 设有两个序列x(n)与y(n)，则帕斯维尔定理为(9)帕斯维尔定理2.6单边z

变换1.单边z

变换的定义 单边z

变换定义为它和双边z变换的不同之处在于它只计算以序列x(n)的正向区间为系数的z-1幂级数，而不管对x(n)在n＜0时如何定义。2.6单边z变换1.单边z变换的定义2.单边z

反变换 单边z

变换也有相应的z反变换，但其解不是惟一的。3.单边z

变换的性质 在2.5节中阐明的双边z

变换的性质和定理，除了和移位特性有关的外，大都适用于单边z

变换。2.单边z反变换4.用单边z变换解线性差分方程 线性差分方程描述的因果系统，当给定适当的初值条件时，需要用单边z变换求解。例如对一阶差分方程y(n)=x(n)+ay(n-1)4.用单边z变换解线性差分方程离散时间系统与z变换2.7z变换与傅氏变换的关系 由前几节的讨论得知，一个序列x(n)的z

变换为x(n)的傅氏变换为如果令z＝ejω，则式(2-96)为2.7z变换与傅氏变换的关系 由前几节的讨论得知，一个序离散时间系统与z变换离散时间系统与z变换2.8系统的时域分析与频域分析 对一个系统来说有两种分析方法，时域分析法与频域分析法。 时城分析法就是利用y(n)=x(n)*h(n)=z-1［X(z)·H(z)］进行分析。2.8系统的时域分析与频域分析 对一个系统来说有两种分析第二章离散时间系统与z变换2.1取样和内插2.2离散时间信号序列2.3离散系统及其普遍关系2.4离散信号的傅氏变换2.5离散信号的z变换2.6单边z变换2.7z变换与傅氏变换的关系2.8系统的时域分析与频域分析第二章离散时间系统与z变换2.1取样和内插2.1取样和内插1.取样 将连续信号变成离散信号有各种取样方法，其中最常用的是等间隔周期取样，即每隔固定时间T取一个信号值，如图2-1所示。其中T称为取样周期，T的倒数称为取样频率或取样率。记为fS=1/T2.1取样和内插1.取样图2-1连续信号的取样图2-1连续信号的取样 取样定理——Shannon定理 任一连续信号xa(t)，设其频谱的最高频率分量为fm，则当对它进行取样时，只要选择取样率等于或大于2fm

，就可以由这个取样序列xa(nT)来唯一准确地恢复xa(t)。 设有一限带信号xa(t)。当|Ω|≤Ω

max，它的付氏变换为Xa(Ω)。将xa(t)乘一取样函数p(t)就得到xa(t)，如图2.2所示。 取样定理——Shannon定理图2-2连续信号取样的数学模型图2-2连续信号取样的数学模型离散时间系统与z变换图2-5取样过程的时域与频域关系图2-5取样过程的时域与频域关系最后需要说明一点：上述取样定理是理想取样，如果取样函数不是单位冲击函数序列，而是窄脉冲函数序列，则如图2-6所示(详细情况请参看相关资料)。最后需要说明一点：上述取样定理是理想取样，如果取样函数不是单图2-6理想取样和非理想取样的比较图2-6理想取样和非理想取样的比较2.内插 用大于奈奎斯特取样频率取样限带信号xa(t)，则被取样信号xa(t)通过理想低通滤波器，只要其截止频率Ωc满足Ωmax≤Ωc≤(Ωs－Ωmax)时，就可以恢复出原来信号，如图2-8所示。2.内插图2-8取样信号的恢复与理想低通滤波器的传输函数图2-8取样信号的恢复与理想低通滤波器的传输函数图2-9连续信号的内插表示图2-9连续信号的内插表示图2-10连续信号用三角形内插函数图2-10连续信号用三角形内插函数3.Matlab实现(1)零阶和一阶保持内插(2)三次样条内插3.Matlab实现2.2离散时间信号序列1.离散时间信号 离散时间信号是在离散的时间上取值，在两个取样间隔内数值为零的信号。2.2离散时间信号序列1.离散时间信号2.常用序列(1)单位取样序列 单位取样序列的定义为

其图形如图2-15所示。2.常用序列图2-15单位取样序列图2-15单位取样序列(2)单位阶跃序列 单位阶跃序列的定义为其图形如图2.16所示。(2)单位阶跃序列图2-16单位阶跃序列图2-16单位阶跃序列(3)矩形序列 矩形序列的定义为其图形如图2-17所示。(3)矩形序列图2-17矩形序列图2-17矩形序列(4)实指数序列 实指数序列的定义为x(n)=an其中a为不等于零的任意实数。 图2-18是0＜a＜1的一个实指数序列的图形。(4)实指数序列图2-18实指数序列图2-18实指数序列(5)正弦序列 正弦序列的定义为x(n)=sinnω0其图形如图2-19所示。(5)正弦序列图2-19正弦序列图2-19正弦序列2.3离散系统及其普遍关系1.离散系统的定义 离散系统在数学上定义为将输入序列x(n)映射成输出序列y(n)的惟一性变换或运算。亦即将一个序列变换成另一个序列的系统，记为y(n)=T［x(n)］通常将上式表示成图2-20所示的框图。2.3离散系统及其普遍关系1.离散系统的定义图2-20离散系统的模型图2-20离散系统的模型2.线性非移变系统(1)系统的线性特性 满足叠加原理的系统具有线性特性，即若对两个激励x1(n)和x2(n)有2.线性非移变系统(2)系统的非移变特性 系统的非移变是指系统的参数不随时间而变化。用数学表示为T［x(n－n0)］=y(n－n0)即不管输入信号作用的时间先后，输出信号响应的形状均相同，仅是出现的时间不同，如图2-22所示。(2)系统的非移变特性图2-22离散系统的非移变特性图2-22离散系统的非移变特性(3)线性非移变系统 线性非移变系统就是既满足迭加原理又具有非移变特性的系统，将其描绘如图2-24所示。图2-24线性非移变系统模型(3)线性非移变系统图2-24线性非移变系统模型(4)线性卷积的计算 计算线性卷积有4种方法。 ①利用两个序列的解析式直接计算式(2-34)。 ②利用两个序列的移位求和，即先把一个序列倒置。每次将它向下移一步，求出两序列重叠部分乘积之和。 ③用作图法求。 ④卷积的Matlab实现(4)线性卷积的计算3.系统的稳定性与因果性(1)稳定性 对于一个系统，当输入序列是有界时，其输出也是有界的，则称它是稳定系统。用数学描述则为如果 ｜x(n)｜＜∞对于一切n则 ｜y(n)｜＜∞对于一切n3.系统的稳定性与因果性 因为其中假设｜x(n)｜≤M。 因为2．因果性 一个系统如果其输出变化不会发生在输入变化之前，则称它是因果的。这就是说对于因果系统，如果取n0,当n<n0时，x1(n)=x2(n),则n<n0时，y1(n)=y2(n)。一个线性非移变系统当n<0时的因果充要条件是其单位取样响应等于零,即h(n)=0 n<0 这个充要条件可以从y（n）＝x(n)*h(n)的解析式中导出。2．因果性4.系统的差分方程描述(1)非递归型(FIR) 非递归型因果系统是输出的现在值仅仅取决于输入的现在值与输入的过去值的系统。非递归，即输出对输入无反馈。因此，设在n时刻输入x(n)与输出y(n)的关系为y(n)=f｛……,x(n-1),x(n),x(n+1),……｝若系统是线性非移变的，y(n)可表示为4.系统的差分方程描述(2)递归型(IIR) 递归型因果系统输出的现在值不仅取决于输入的现在值与过去值，还取决于输出的过去值。y(n)=f｛……,x(n-1),x(n),x(n+1),……｝ +g｛……,y(n-1),y(n+1),……｝(2)递归型(IIR) 同理，在系统为线性、非移变、因果时，可推得 同理，在系统为线性、非移变、因果时，可推得2.4离散信号的傅氏变换1.问题的提出 对于连续信号xa(t)与其频谱Xa(Ω)之间存在着傅氏变换关系，如图2-28所示。前边已经讨论了连续信号xa(t)的离散化，即取样的问题，已经知道取样序列的频谱是原信号频谱在Ω轴上的周期延拓，如图2-28(b)所示。2.4离散信号的傅氏变换1.问题的提出图2-28连续和离散信号的傅氏变换图2-28连续和离散信号的傅氏变换2.傅氏变换对的推导 从2.1节的讨论知道，对连续信号在时域内进行取样的结果，是频域内频谱周期的延拓，并且还得到了已取样信号xa(t)在频域内的表示，现重写为2.傅氏变换对的推导离散时间系统与z变换3.线性非移变系统的频率响应 从2.3节的讨论得到线性非移变离散系统的输入输出关系为对上式两边同时进行傅氏变换得3.线性非移变系统的频率响应离散时间系统与z变换离散时间系统与z变换4.离散信号傅氏变换的性质 首先定义两种序列，共轭对称序列与共轭反对称序列。(1)序列的卷积特性 序列的卷积特性是时域内的卷积关系映射到频域内为相乘，即4.离散信号傅氏变换的性质(2)序列傅氏变换的周期性 序列的傅氏变换是ω的周期函数，周期为2π。(3)复序列傅氏变换的对称性(2)序列傅氏变换的周期性离散时间系统与z变换离散时间系统与z变换2.5离散信号的z变换1.z变换的定义及其收敛域 对于一个序列x(n)，其z变换的定义为 其中z为复变量，也可记作Z［x(n)］=X(z)。式(2-49)的定义也称为双边z

变换；相应的还有单边z变换。2.5离散信号的z变换1.z变换的定义及其收敛域对于所有的序列或所有的z值，z变换并不总是收敛的。对于任意给定的序列，使z变换收敛的z值集合称作收敛区域：{Z：X(z)存在}＝收敛区域。其内径R-与外径R+分别取x(n)在n→∞和n→-∞时的形状。z变换收敛域的概念很重要，不同的序列可能有相同的z变换表达式，但是收敛域却不同，所以应该特别注意，只有当z变换的表达式与收敛域都相同时，才能判定两个序列相等。对于所有的序列或所有的z值，z变换并不总是收敛的。对于任意给根据以上讨论，可以概括为(1)对右边序列(n≥0存在)，｜z｜＞R-收敛，且R-是右序列的极点。(2)对左边序列(n＜0存在)，｜z｜＜R+收敛，且R+是左边序列的极点。根据以上讨论，可以概括为(3)若X(z)不只一个极点，则找与收敛域相重的那个极点，对右边序列，最外极点之外的区域为收敛域；对左边序列，最内极点之内的区域为收敛域，如图2-31所示。(4)对双边序列，若在左边序列的收敛域存在重叠部分，则这重叠部分就是它的收敛域。若不存在重迭部分，则z变换不存在。(3)若X(z)不只一个极点，则找与收敛域相重的那个极点，2.系统传递函数H(z)的频域表示 描述线性非移变系统的差分方程为对上式方程两边取z变换为2.系统传递函数H(z)的频域表示离散时间系统与z变换 在Matlab中，reqz函数计算幅度和相位响应，它有如下5种调用方式。(1)［H，w］＝freqz(b，a，N) b和a分别表示分子和分母的系数向量，与filter(b，a，x)函数中的相同。此函数在单位圆上半部上等间隔的计算N点频率响应，返回该系统的N点频率矢量w和N点复数频率响应矢量H。如果N没有说明，则缺省值为512。 在Matlab中，reqz函数计算幅度和相位响应，它有如下(2)［H，w］＝ freqz(b，a，N，’whole’) 在整个单位圆上等间隔的计算N点频率响应。(3)H＝freqz(b，a，w) 它返回矢量w指定的那些频率点上的频率响应，通常在0到π之间。(2)［H，w］＝(4)［H，F］＝freqz(b，a，N，Fs)和［H，F］＝freqz(b，a，N，'whole'，Fs) 给定取样频率Fs，单位为Hz；返回单位为Hz的频率矢量F。(5)H＝freqz(b，a，F，Fs) 给定单位为Hz的取样频率Fs，返回矢量F指定的那些频率点上的复数频率响应，单位也是Hz。(4)［H，F］＝freqz(b，a，N，Fs)和［H，F］3.z反变换 z反变换关系式可以利用柯西积分定理推导出来，柯西定理为式中c是一个逆时针方向环绕原点的围线。3.z反变换(1)幂级数法 幂级数法也就是长除法，对于给定的z变换X(z)，可以根据它的收敛域判定序列x(n)是右边序列还是左边序列，或是双边序列。(1)幂级数法(2)部分分式法 当x(z)序列为有理函数时，可将x(z)写成一个和式为因为，对于右边序列存在如下变换关系(2)部分分式法(3)留数法 我们知道式中c是X(z)收敛域内的积分围线。对于n＞0时，对应右边序列，此时极点在c内，对n＜0时，对应左边序列。此时极点在c外，根据留数定理有(3)留数法(4)Matlab

的实现 在Matlab中，函数residuez

计算有理函数的留数Rj、极点pj和直接项系数Cj，设有多项式如下：B(z)和A(z)分别是分子分母多项式，它们按z-1递增顺序排列。(4)Matlab的实现4.因果非移变系统的稳定性，收敛区与极点的关系 正如在式(2-56)中看到的，线性非移变系统的系统函数H(z)是具有实系数z的有理函数。现在来讨论系统函数的极点与系统稳定性和收敛区的关系，亦即证明系统函数的极点分布将决定系统是否稳定。

4.因果非移变系统的稳定性，收敛区与极点的关系 假设有一N阶因果系统，系统函数为H(z)，为方便起见，设H(z)只有单阶极点，这样系统的单位取样响应由式(2-67)给出 假设有一N阶因果系统，系统函数为H(z)，为方便起见，设H离散时间系统与z变换由以上的讨论清楚看到，因果稳定系统的收敛区域包括单位圆以及以外的整个z平面。因而，因果非移变的稳定系统为(1)极点都在单位圆内；(2)收敛区域为l≤z≤