数字信处理第三版课件

《数字信处理第三版课件》由会员分享，可在线阅读，更多相关《数字信处理第三版课件（71页珍藏版）》请在装配图网上搜索。

1、第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版第第2章章 时域离散信号和系统的频域分析时域离散信号和系统的频域分析 2.1 引言引言 2.2 序列的傅里叶变换的定义及性质序列的傅里叶变换的定义及性质 2.3 周期序列的离散傅里叶级数及傅里叶变换表示式周期序列的离散傅里叶级数及傅里叶变换表示式2.4 时域离散信号的傅里叶变换与模拟时域离散信号的傅里叶变换与模拟 信号傅里叶变换之间的关系信号傅里叶变换之间的关系 2.5 序列的序列的Z变换变换 2.6 利用利用Z变换分析信号和系统的频域特性变换分析信号和系统的频域特性 第第1章章 时域离散信号和时域离散系统时域

2、离散信号和时域离散系统 数字信号处理第三版2.2 序列的傅里叶变换的定义及性质序列的傅里叶变换的定义及性质 2.2.1 序列傅里叶变换的定义序列傅里叶变换的定义 定义定义()( )jj nnX ex n e(2.2.1) 为序列为序列x(n)的傅里叶变换，用的傅里叶变换，用FT(Fourier Transform)表表示。示。 FT成立的成立的充分必要条件充分必要条件是序列是序列x(n)绝对可和，即满绝对可和，即满足下式：足下式： ( )nx n (2.2.2) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 FT反变换定义为：反变换定义为：(2.2.4)

3、 (2.2.1)和和(2.2.4)式组成式组成一对傅里叶变换公式一对傅里叶变换公式。一些绝对不可和的序列（如周期序列），其傅里叶一些绝对不可和的序列（如周期序列），其傅里叶变换可用冲激函数的形式表示出来。变换可用冲激函数的形式表示出来。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版【例例2.2.1】设x(n)=RN(n)，求x(n)的傅里叶变换。解解 当N=4时，其幅度与相位随频率的变化曲线如图2.2.1所示。 )2/sin()2/sin(e )ee (e)ee (ee1e1 ee )()e (2)1( j2/j2/j2/j2/j2/j2/jjj10jjj

4、NnRxNNNNNnNnnN第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版图2.2.1R4(n)的幅度与相位曲线 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2.2.2 序列傅里叶变换的性质序列傅里叶变换的性质1. FT的周期性的周期性 在定义式在定义式(2.2.1)中，中，n取整数，下式成立：取整数，下式成立： (2)()( ),jjM nnX ex n eM为整数为整数 (2.2.6) 序列的傅里叶变换是频率序列的傅里叶变换是频率的周期函数，周期是的周期函数，周期是2。 这样这样X(ej)可以展成傅里叶级数，可以

5、展成傅里叶级数，(2.2.1)式就是傅里叶级式就是傅里叶级数的形式，数的形式，x(n)是其系数。是其系数。 ()( )jj nnX ex n e由于由于FT的周期性，一般的周期性，一般只分析只分析或或02之间的之间的FT第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2. 线性线性 11221212()( ),()( ),( )( )()()jjjjX eFT x nXeFT x nFT ax nbx naX ebXe则则 设设 式中式中a, b为常数。为常数。 3. 时移与频移时移与频移 设设X(e j) = FTx(n)，则：，则：(2.2.7)0000

6、( ()()( )()j njjnjFT x nneX eFT ex nX e (2.2.8) (2.2.9) ）第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版4. 对称性对称性 先了解先了解共轭对称共轭对称与与共轭反对称共轭反对称以及它们的性质：以及它们的性质： 定义：定义：设序列设序列xe(n)满足满足 xe(n)=x*e(-n) 则称则称xe(n)为共轭为共轭对称序列对称序列。共轭对称序列的共轭对称序列的性质性质： 将将xe(n)用其实部与虚部表示：用其实部与虚部表示： xe(n) = xer(n)+jxei(n) 两边两边 n 用用 n 代替，并取

7、共轭，得：代替，并取共轭，得： x*e(-n)=xer(-n)-jxei(-n) xer(n) = xer(-n) (2.2.11) xei(n) = -xei(-n) (2.2.12)v 共轭对称序列其实部是偶函数，而虚部是奇函数。共轭对称序列其实部是偶函数，而虚部是奇函数。对比两式，对比两式，得：得：第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版共轭反对称序列的共轭反对称序列的性质性质： 将将x0(n) 用实部与虚部表示：用实部与虚部表示： xo(n) = xor(n)+jxoi(n) 得：得： xor(n) = -xor(-n) (2.2.14) x

8、oi(n) = xoi(-n) (2.2.15)v 共轭反对称序列的实部是奇函数，共轭反对称序列的实部是奇函数， 而虚部是偶函数而虚部是偶函数。定义：定义：满足下式的序列称满足下式的序列称共轭反对称序列：共轭反对称序列： xo(n) = -x*o(-n) (2.2.13)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版v 一般序列可用共轭对称与共轭反对称序列之和表示一般序列可用共轭对称与共轭反对称序列之和表示， 即：即： x(n) = xe(n)+xo(n) (2.2.16) xe(n), xo(n)可以分别用原序列可以分别用原序列x(n)求出求出 将将(2

9、.2.16)式中的式中的n用用-n代替，代替， 再取共轭得到：再取共轭得到： x*(-n) = xe(n)-xo(n) (2.2.17) 比较两式，比较两式， 得得 ：1( ) ( )()21( ) ( )()2eox nx nxnx nx nxn(2.2.18) (2.2.19) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版在频域，函数在频域，函数X(ej)也有类似的概念和结论：也有类似的概念和结论： X(ej) = Xe(ej)+Xo(ej) (2.2.20) 共轭对称部分共轭对称部分Xe(ej)和共轭反对称部分和共轭反对称部分Xo(ej) 满足：满

10、足： Xe(ej) =X*e(e-j) (2.2.21) Xo(ej) = -X*o(e-j) (2.2.22) 同样有下面公式：同样有下面公式： 1()()()21()()()2jjjejjjoXeX eXeXeX eXe(2.2.23) (2.2.24) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 FT的对称性的对称性 (a) 将序列将序列x(n)分成实部分成实部xr(n)与虚部与虚部xi(n) x(n) = xr(n) + jxi(n) 进行进行FT，得：，得： X(e j) = Xe(e j) + Xo(e j) 式中式中 xr(n)和和xi(

11、n)都是实数序列。都是实数序列。Xe(ej) 具有共轭对称性，其实部是偶函数，虚部是奇函数。具有共轭对称性，其实部是偶函数，虚部是奇函数。 Xo(ej) 具有共轭反对称性质，其实部是奇函数，虚部是偶函数。具有共轭反对称性质，其实部是奇函数，虚部是偶函数。 结论结论： x(n) = xr(n) + jxi(n) X(e j) = Xe(e j) + Xo(e j)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 (b) 将序列分成共轭对称部分将序列分成共轭对称部分xe(n)和共轭反对称部和共轭反对称部分分xo(n)，即：，即： x(n) = xe(n)+xo(

12、n) (2.2.25) 由由(2.2.18)式和式和(2.2.19)式：式： 1( ) ( )()21( ) ( )()2eox nx nxnx nx nxn 将上面两式分别进行将上面两式分别进行FT，得：，得： FTxe(n)=1/2X(ej)+X*(ej)= ReX(ej)= XR(ej) FTxo(n)=1/2X(ej) -X*(ej)= jImX(ej)= jXI(ej) 因此对因此对(2.2.25)式进行式进行FT得到：得到： X(ej) = XR(ej)+jXI(ej) (2.2.26)结论结论： x(n) = xe(n) + xo(n) X(ej) = XR(ej) + jXI(

13、ej)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版x(n) = xe(n) + xo(n) X(ej) = XR(ej) + jXI(ej)x(n) = xr(n) + jxi(n) X(e j) = Xe(e j) + Xo(e j)利用利用FT的对称性的对称性，可得以下四个结论：可得以下四个结论：（1） x(n)为为实序列实序列（xi(n)=0），得），得X(ej) = Xe(ej)为共轭为共轭对称函数，即对称函数，即 X(ej) = X*(e-j)（2） x(n)为为实偶序列实偶序列（ xi(n)=0且且x(n)= x(-n)，x0(n)=0），）

14、，得得X(ej)为实偶函数，即为实偶函数，即 X(ej) = X(e-j)（3） x(n)为为实奇序列实奇序列（xi(n)=0且且x(n)= -x(-n)，xe(n)=0），），得得X(ej)为纯虚奇对称函数，即为纯虚奇对称函数，即 X(ej) = X*(e-j)=-X(e-j)（4） x(n)为为实因果序列：实因果序列：x(n)= xe(n) +xo(n) ，0, 00, )(0),(2)(nnnxnnxnxee或：或：0, 00, )0(0),(2)(nnxnnxnxox e(n)=1/2x(n)+ x(-n) x o(n)=1/2x(n) - x(-n)第第1章章 时域离散信号和时域离散

15、系统时域离散信号和时域离散系统 数字信号处理第三版对对实因果序列实因果序列，只要知道，只要知道XR(ej) ，就可求得，就可求得x(n)，过程如下：过程如下：已知已知：XR(ej)=FTxe(n) xe(n) x(n) X(ej) 已知已知XI(ej)和和 x(0) ：jXI(ej) xo(n) x(n) X(ej)v 对实对实因果因果序列：其傅里叶变换序列：其傅里叶变换X(ej)的实部包含了的实部包含了X(ej)或或x(n)的全部信息，即的全部信息，即X(ej) 中有冗余信息。中有冗余信息。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版例例 2.2.3

16、x(n)=anu(n)， 0a1， 求其偶函数求其偶函数xe(n) 和和奇函数奇函数xo(n)。 解：解： x(n) = xe(n)+xo(n)(0),01( ),021(),02xnx nnxnn( )ex n 1,01,021,02nnnanan第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版5. 时域卷积定理时域卷积定理 设设 y(n)=x(n)*h(n), 则则 Y(e j)=X(e j)H(e j) (2.2.32) 定理说明，定理说明， 两序列卷积的两序列卷积的FT，服从相乘的关系。，服从相乘的关系。 对对LTI系统，其输出的系统，其输出的FT等

17、于输入信号的等于输入信号的FT乘以单位脉乘以单位脉冲响应的冲响应的FT。 因此因此求系统的输出求系统的输出信号，信号， (1)可以在时域用卷积公式可以在时域用卷积公式(1.3.7)； (2)可以在频域按照可以在频域按照 (2.2.32)式，求出输出的式，求出输出的FT，再作逆，再作逆FT求出输出信号。求出输出信号。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版6. 频域卷积定理频域卷积定理 设设y(n) = x(n)h(n)，则：，则：(2.2.33) 定理说明，定理说明，在时域两序列相乘，对应频域为卷积关系。在时域两序列相乘，对应频域为卷积关系。 7.

18、帕斯维尔帕斯维尔(Parseval)定理定理 定理说明定理说明，信号时域的总能量等于频域的总能量信号时域的总能量等于频域的总能量。这里频域总能量是指这里频域总能量是指|X(e j)|2在一个周期中的积分再乘以在一个周期中的积分再乘以1/(2)。(2.2.34) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2.3 周期序列的离散傅里叶级数周期序列的离散傅里叶级数 及傅里叶变换表示式及傅里叶变换表示式 2.3.1 周期序列的离散傅里叶级数周期序列的离散傅里叶级数 设设 是以是以N为周期的周期序列，为周期的周期序列， 由于周期性，由于周期性， 可以展成傅里叶级

19、数：可以展成傅里叶级数：( )x n2( )jknNkkx na e(2.3.1)式中式中ak是傅里叶级数的系数。是傅里叶级数的系数。-k1 X(z)存在的条件是存在的条件是|z-1|1， 由由X(z)表达式表明，极点是表达式表明，极点是z=1，单位圆上的，单位圆上的Z变换不存变换不存在，或者说收敛域不包含单位圆。因此其傅里叶变换不存在，在，或者说收敛域不包含单位圆。因此其傅里叶变换不存在，更不能用更不能用(2.5.4)式求式求FT。该序列的该序列的FT不存在，但如果引进奇异函数不存在，但如果引进奇异函数()，其傅里，其傅里叶变换可以表示出来叶变换可以表示出来(见表见表2.3.2)。该例说明该

20、例说明一个序列的傅里叶变换不存在，在一定收敛域一个序列的傅里叶变换不存在，在一定收敛域内内Z变换是存在的。变换是存在的。 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 2.5.2 序列特性对收敛域的影响序列特性对收敛域的影响 序列的特性决定其序列的特性决定其Z变换收敛域，了解序列特性与收变换收敛域，了解序列特性与收敛域的一些关系，有助于敛域的一些关系，有助于Z变换的使用。变换的使用。 x(n) n1nn2 x(n)= 0 其它其它其其Z变换为：变换为：21( )( )nnn nX zx n z1. 有限长序列有限长序列 如序列如序列x(n)满足：满足：

21、第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版x(n)为有界序列为有界序列，由于是有限项求和，除，由于是有限项求和，除0与与两点是两点是否收敛与否收敛与n1、n2取值有关外，整个取值有关外，整个z平面均收敛。平面均收敛。 若：若：n10，出现，出现z-n项，则收敛域不包括项，则收敛域不包括z=0点；点；如果是因果序列，收敛域包括如果是因果序列，收敛域包括z=点。点。v有限长序列的收敛域表示如下：有限长序列的收敛域表示如下：21( )( )nnn nX zx n zn10， n20， 0|z|:n10， 00， 0|z|:第第1章章 时域离散信号和时域离散系

22、统时域离散信号和时域离散系统 数字信号处理第三版例例 2.5.2 求求x(n)=RN(n)的的Z变换及其收敛域。变换及其收敛域。解解： 1101( )( )1NNnnNnnzX zRn zzz这是一个因果的有限长序列，因此收敛域为这是一个因果的有限长序列，因此收敛域为0z。 从从X(z)的分母看到的分母看到z=1似乎是似乎是X(z)的极点，但同时分子的极点，但同时分子多项式在多项式在z=1时也有一个零点，极零点对消，时也有一个零点，极零点对消，X(z)在单位圆在单位圆上仍存在。上仍存在。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2. 右序列右序列 右序

23、列是在右序列是在nn1时，序列值不全为零，而在时，序列值不全为零，而在nn1 时，时，序列值全为零的序列。序列值全为零的序列。 v第一项为有限长序列，设第一项为有限长序列，设n1-1，其收敛域为，其收敛域为0|z|。 v第二项为因果序列，其收敛域为第二项为因果序列，其收敛域为Rx-|z|，Rx-是第二项是第二项最小的收敛半径。最小的收敛半径。将两收敛域相与，其收敛域为将两收敛域相与，其收敛域为Rx- |z|。如果是因果序列，如果是因果序列， 收敛域定为收敛域定为Rx- |z|。 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版v 如果如果n20, z=0点收

24、敛，点收敛，z=点不收敛，其收敛域为点不收敛，其收敛域为0|z|0，则收敛域为，则收敛域为0|z|n2， 序列值全为零的序列。序列值全为零的序列。 其其Z变换为：变换为： 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版4. 双边序列双边序列 一个双边序列可以看作一个左序列和一个右序列之一个双边序列可以看作一个左序列和一个右序列之和，其和，其Z变换为：变换为：X(z)的收敛域是的收敛域是X1(z)和和X2(z)收敛域的公共收敛区域。收敛域的公共收敛区域。v 如果如果Rx+Rx-，其收敛域为，其收敛域为Rx-|z|Rx+ ，这是一个环状域，这是一个环状域v 如

25、果如果Rx+ Rx- ，两个收敛域没有公共区域，两个收敛域没有公共区域，X(z)没有收敛没有收敛域，域， 因此，因此，X(z)不存在。不存在。 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版例例2.5.5x(n)=a|n|, a为实数，求为实数，求x(n)的的Z变换及其收变换及其收敛域。敛域。解解第一部分收敛域为第一部分收敛域为|az|1，得，得|z|a|1; 第二部分收敛域为第二部分收敛域为|az1|a|。如果。如果|a|1, 两部分的公共收敛域为两部分的公共收敛域为|a|z|a|1, 其其Z变换如下式变换如下式:如果如果|a|1，则无公共收敛域，因此，

26、则无公共收敛域，因此X(z)不存在。当不存在。当0aa，求其逆，求其逆Z变变换换x(n)。解解： 用留数定理求解，用留数定理求解， 要先找出要先找出F(z)的极点，的极点，极点有：（极点有：（1）z=a （2）当）当n0时，时，z=0也是极点也是极点其中极点其中极点z=0与与n的取值有关：的取值有关：n0时，时，n=0不是极点。不是极点。 n0时，时，z=0是一个是一个n阶极点。阶极点。 因此要分成因此要分成n0和和n0两种情况求两种情况求x(n)。 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版（1）n0 时，只有一个极点：时，只有一个极点：( )Re

27、( ), ()nz anx ns F z azzazaa （2）n0时，增加时，增加z=0的的n阶极点，不易求留数，采用留数辅阶极点，不易求留数，采用留数辅助定理求解，检查助定理求解，检查(2.5.10)式是否满足，由于式是否满足，由于n0，只要，只要N-M0，(2.5.10)式就满足。式就满足。 本例满足本例满足(2.5.10)式。式。N-M-n1 (2.5.10)图图 2.5.4 例中例中na，根据前面分析的序，根据前面分析的序列特性对收敛域的影响知道，列特性对收敛域的影响知道，x(n)一定是因果的右序列，这一定是因果的右序列，这样样n0部分一定为零，无需再求。部分一定为零，无需再求。第第

28、1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版例例2.5.7已知已知, 求其求其逆变换逆变换x(n)。解该例题没有给定收敛域，为求出唯一的原序解该例题没有给定收敛域，为求出唯一的原序列列x(n)，必须先确定收敛域。分析，必须先确定收敛域。分析X(z)， 得到其极点得到其极点分布如图分布如图2.5.5所示。图中有两个极点：所示。图中有两个极点：z=a和和z=a1，这样收敛域有三种选法，它们是这样收敛域有三种选法，它们是 （1） |z|a1|，对应的，对应的x(n)是因果序列；是因果序列；（2） |z|a|，对应的，对应的x(n)是左序列；是左序列；（3） |a|

29、z|a1|:这种情况的原序列是因果的右序列，无须求这种情况的原序列是因果的右序列，无须求n0时的时的x(n)。当。当n0时，时，F(z)在在c内有两个极点：内有两个极点：z=a和和z=a1，因此因此2111( )(1)(1)naF zzazaz211()()naza za za第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版最后表示成：最后表示成：x(n)=(anan)u(n)。 1( )Res ( ), Res ( ),x nF z aF z a12211(1)(1)()()()(1)()()nnz az aazazzazazaaza za zannaa第

30、第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版（2） 收敛域为收敛域为|z|a|:这种情况原序列是左序列，无须计算这种情况原序列是左序列，无须计算n0情况。实情况。实际上，当际上，当n0时，围线积分时，围线积分c内没有极点，因此内没有极点，因此x(n)=0。n0时，时，c内只有一个极点内只有一个极点z=0，且是，且是n阶极点，改求阶极点，改求c外外极点留数之和。极点留数之和。n0时时, 1( )Res ( ), Res ( ),x nF z aF z a 122111(1)(1)()()()()()()nnz az aazazzazaa za zaa za

31、za ()nnnnaaaa 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版最后将最后将x(n)表示成封闭式：表示成封闭式：x(n)=(anan)u(n1)（3） 收敛域为收敛域为|a|z|a1|:这种情况对应的这种情况对应的x(n)是双边序列。根据被积函数是双边序列。根据被积函数F(z)，按，按n0和和n0两种情况分别求两种情况分别求x(n)。n0时，时，c内只有内只有1个极点：个极点：z=a, x(n)=ResF(z), a=ann0时，时，c内极点有内极点有2个，其中个，其中z=0是是n阶极点，改求阶极点，改求c外极点留数，外极点留数，c外极点只有外极

32、点只有z=a1， 因此因此x(n)=ResF(z), a1=an最后将最后将x(n)表示为表示为即即 x(n)=a|n|0( )0nnanx nan第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2.5.4 Z 变换的性质和定理变换的性质和定理1. 线性线性 设设 X(z)=ZTx(n), Rx- |z| Rx+ Y(z)=ZTy(n), Ry- |z| Ry+ 则则 ZT a x(n)+b y(n) =aX(z)+bY(z), R m-|z|R m+ (2.5.15)其中：其中： Rm+= min Rx+ ,Ry+ Rm-= max Rx- ,Ry-即即Z

33、变换的收敛域变换的收敛域(Rm-，Rm+)是是X(z)和和Y(z)的公共收敛域，的公共收敛域，若无公共收敛域则若无公共收敛域则Z变换不存在。变换不存在。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2. 序列移位序列移位 设设 X(z)=ZTx(n), R x-|z|R x+ 则则 ZTx(n-n0)= z-n0X(z), R x-|z|R x+ (2.5.16)3. 乘指数序列乘指数序列 设设 X(z)=ZTx(n), R x-|z|R x+ y(n)=anx(n), a为常数为常数 则则 Y(z)=ZTanx(n) =X(a-1 z) |a|R x-

34、|z| |a|R x+ (2.5.17)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 4. 序列乘序列乘n 设设 ( ) ( )( )( )xxxxX zZT x nRzRdX zZT nx nzRzRdz 则则(2.5.18) 5. 复序列的共轭复序列的共轭 设设6. 初值定理初值定理 设设 x(n)是因果序列，是因果序列，X(z)=ZTx(n)(0)lim( )xxX z(2.5.20) 则则第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版7. 终值定理终值定理 若若x(n)是因果序列，其是因果序列，其Z变换的极点

35、，除可以变换的极点，除可以有一个有一个一阶极点在一阶极点在z=1上上，其它极点均在单位圆内，则，其它极点均在单位圆内，则 ：终值定理也可用终值定理也可用X(z)在在z=1点的留数表示：点的留数表示：如果单位圆上如果单位圆上X(z)无极点，则无极点，则x()=0。 (2.5.22) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版9. 复卷积定理复卷积定理如果如果 ZTx(n)= X(z)， R x-|z|R x+ ZTy(n)= Y(z), R y-|z|R y+ w(n) = x(n)y(n)则则W(z)的收敛域的收敛域 (2.5.24) (2.5.25)

36、8. 序列卷积序列卷积 设设 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版10. 帕斯维尔帕斯维尔(Parseval)定理定理 那么那么 v 平面上，平面上，c 所在的收敛域为：所在的收敛域为：11max(,)min(,)xxyyRvRRR设设 (2.5.27)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 2.5.5 利用利用Z变换解差分方程变换解差分方程 用用Z变换求解差分方程，将差分方程变成了代数方程，变换求解差分方程，将差分方程变成了代数方程，使求解过程简单。使求解过程简单。 N阶线性常系数差方程为：阶线性常

37、系数差方程为：利用线性和序列移位性利用线性和序列移位性对于对于N阶差分方程，求其解必须已知阶差分方程，求其解必须已知N个初始条件。个初始条件。 设设x(n)是因果序列（是因果序列（x(n)=0,nmax(|a|,|b|)1111( )2(),0nnny nbabnab式中第一项为零输入解，第二项为零状态解。式中第一项为零输入解，第二项为零状态解。11111bzbazaba|,*2bzwhenbbn第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2.6.2 用系统函数的极点分布分析系统的因果性和稳定性用系统函数的极点分布分析系统的因果性和稳定性 因果因果(可实

38、现可实现)系统系统其单位脉冲响应其单位脉冲响应h(n)一定满足当一定满足当n0时，时，h(n)=0，那么其系统函数，那么其系统函数H(z)的收敛域一定包含的收敛域一定包含点，点，即即点不是极点，极点分布在某个圆的圆内，点不是极点，极点分布在某个圆的圆内，收敛域在某收敛域在某个圆外。个圆外。( )nh n 稳定系统稳定系统要求要求 ，对照，对照Z变换定义，系统稳定变换定义，系统稳定要求要求收敛域包含单位圆收敛域包含单位圆。 如果如果系统因果且稳定系统因果且稳定，收敛域包含，收敛域包含点和单位圆，那点和单位圆，那么收敛域可表示为：么收敛域可表示为： r|z|， 0r1 v 系统因果且稳定，系统因果

39、且稳定，H(z)的极点集中在单位圆的内部。的极点集中在单位圆的内部。v 具体系统的因果性和稳定性可由系统函数的具体系统的因果性和稳定性可由系统函数的极点分布极点分布确定。确定。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版例例2.6.1 已知已知 ，分析其因果性，分析其因果性和稳定性。和稳定性。211( ),01(1)(1)aH zaazaz解：解：H(z)的极点为的极点为z=a，z=a-1，如图所示。，如图所示。 (1)收敛域收敛域a-1|z|，对应的系统是因果系统，但由于收敛，对应的系统是因果系统，但由于收敛域不包含单位圆，因此是不稳定系统。单位脉冲响

40、应：域不包含单位圆，因此是不稳定系统。单位脉冲响应：h(n)=(an-a-n)u(n)(见例题见例题2.5.7)，这是一个因果序列，但不收敛。，这是一个因果序列，但不收敛。 (2)收敛域收敛域0|z|a，对应的系统是非因果且不稳定系统。其，对应的系统是非因果且不稳定系统。其单位脉冲响应单位脉冲响应h(n)=(a-n-an)u(-n-1)(见例题见例题2.5.7)，这是一个非因，这是一个非因果且不收敛的序列。果且不收敛的序列。 (3)收敛域收敛域a|z|a-1，对应的系统是一个非因果系统，但由于，对应的系统是一个非因果系统，但由于收敛域包含单位圆，因此是稳定系统。其单位脉冲响应收敛域包含单位圆，

41、因此是稳定系统。其单位脉冲响应h(n)=a|n|，这是一个收敛的双边序列，如图，这是一个收敛的双边序列，如图2.6.1(a)所示所示。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版下面分析本例这种系统的下面分析本例这种系统的可实现性可实现性：在在H(z)H(z)的三种收敛域中，前二种系统不稳定，不能选用；的三种收敛域中，前二种系统不稳定，不能选用；在最后一种收敛域中，系统稳定但非因果，还是不能具体在最后一种收敛域中，系统稳定但非因果，还是不能具体实现。实现。严格地讲，这种系统是无法具体实现的。严格地讲，这种系统是无法具体实现的。但是我们利用数字系统或计算机

42、的存贮性质，可以近似实但是我们利用数字系统或计算机的存贮性质，可以近似实现第三种情况。现第三种情况。方法方法：将图：将图2.6.1(a)2.6.1(a)的的h(n)h(n)从从-N-N到到N N截取一段，再向右移，截取一段，再向右移，形成如图形成如图2.6.1(b)2.6.1(b)所示的所示的h(n)h(n)序列，将序列，将h(n)h(n)作为具体实现作为具体实现的系统单位脉冲响应。的系统单位脉冲响应。N N愈大，愈大，h(n)h(n)表示的系统愈接近表示的系统愈接近h(n)h(n)系统。系统。具体实现时，预先将具体实现时，预先将h(n)h(n)存贮起来，备运算时应用。存贮起来，备运算时应用。

43、v这种非因果但稳定系统的近似实现性，是数字信号处理技这种非因果但稳定系统的近似实现性，是数字信号处理技术比模拟信息处理技术优越的地方。术比模拟信息处理技术优越的地方。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版2.6.3 利用系统的极零点分布分析系统的频率特性利用系统的极零点分布分析系统的频率特性 将将(2.6.2)式因式分解，得：式因式分解，得： 式中：式中：A=b0/a0，影响传输函数的幅度大小；，影响传输函数的幅度大小； cr是是H(z)的的零点零点，dr是其是其极点极点。 零点零点cr和极点和极点d 的分布影响系统的特性。的分布影响系统的特性。

44、下面用几何方法来研究系统零极点分布对系统频率特下面用几何方法来研究系统零极点分布对系统频率特性的影响。性的影响。将将(2.6.4)式分子分母变为正幂次，得：式分子分母变为正幂次，得：00( )( )( )MiiiNiiibzY zH zX za z(2.6.4) 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版11()11()( )()()()()MrN MrNrrMjrjjN MrNjrrzcH zAzzdecH eAeed设系统稳定，将设系统稳定，将z=e j，得到频率响应函数：，得到频率响应函数： (2.6.5) (2.6.6) 若若N=M，则：，则：

45、(2.6.7)第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版jrrjrrc Becd Bed 和和 分别称为零点矢量和极点矢量，将它们用分别称为零点矢量和极点矢量，将它们用极坐标表：极坐标表：rc B rd B将将 和和 表示式代入表示式代入(2.6.7)式，得：式，得：rc B rd Brc B rd B 在在z平面上，平面上，ej-cr用一条由零点用一条由零点cr指向单位圆上指向单位圆上ej点点B的向量的向量 表示；同样表示；同样ej-dr用由极点指向用由极点指向ej点点B的的向量向量 表示，如图所示。表示，如图所示。第第1章章 时域离散信号和时域离散

46、系统时域离散信号和时域离散系统 数字信号处理第三版()111111()()()()NrjjjrNrrNrjrNrrNNrrrrc BH eAH eed Bc BH eAd B (2.6.8)(2.6.9) 式中：式中：系统的传输特性或信系统的传输特性或信号的频率特性由号的频率特性由(2.6.8)式式和和(2.6.9)式确定。式确定。当频率当频率从零变化到从零变化到2时，这些向量的终点时，这些向量的终点B沿单位圆逆时针旋转一周，沿单位圆逆时针旋转一周，按照上述两式，可以分别按照上述两式，可以分别估算出系统的幅度特性和估算出系统的幅度特性和相位特性。相位特性。所有零点矢量相角之和所有零点矢量相角之

47、和所有极点矢量模之积所有零点矢量模之积)()()(argarg)()()()()(zHzHMNBdBcMNzHzHBdBceHrrrrrrrrj图图2.6.2 频响的几何表示法频响的几何表示法下图表示了有一个零下图表示了有一个零点和二个极点的频率特性。点和二个极点的频率特性。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版 按照按照(2.6.8)式，由零极点分布，可以确定零极点位置对式，由零极点分布，可以确定零极点位置对系统特性的影响。系统特性的影响。当当B点转到点转到极点附近极点附近时，极点矢量长度最短，因而幅度时，极点矢量长度最短，因而幅度特性可能特性可

48、能出现峰值出现峰值，且极点愈靠近单位圆，极点矢量长度愈，且极点愈靠近单位圆，极点矢量长度愈短，峰值愈高愈尖锐。短，峰值愈高愈尖锐。如果极点在单位圆上，则幅度特性为如果极点在单位圆上，则幅度特性为 ，系统不稳定。，系统不稳定。对于零点，情况相反，当对于零点，情况相反，当B点转到点转到零点附近零点附近，零点矢量，零点矢量长度变短，幅度特性将长度变短，幅度特性将出现谷值出现谷值，零点愈靠近单位圆，谷值，零点愈靠近单位圆，谷值愈接近零。当零点处在单位圆上时，谷值为零。愈接近零。当零点处在单位圆上时，谷值为零。结论：结论： 极点位置影响频响的峰值位置及尖锐程度，极点位置影响频响的峰值位置及尖锐程度， 零

49、点位置影响频响的谷点位置及形状。零点位置影响频响的谷点位置及形状。第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 数字信号处理第三版例例2.6.2 已知已知H(z)=z-1，分析其频率特性。，分析其频率特性。 解：解： 由由H(z)=z-1，极点为，极点为z=0，幅度特性，幅度特性|H(e j)|=1，相，相位特性位特性()=-，频响如图，频响如图2.6.3所示。所示。 图图2.6.3 H(z)=z-1的频响的频响 用几何方法确定：从用几何方法确定：从=0转到转到=2时，极点矢量的时，极点矢量的长度始终为长度始终为“1 1”。 该例说明该例说明：处于原点处的零点或极点，由于零点矢：处于原点处的零点或极点，由于零点矢量长度或者是极点矢量长度始终为量长度或者是极点矢量长度始终为“1 1” ，因此，因此原点处的原点处的零极点不影响系统的频率特性零极点不影响系统的频率特性。