数字信号处理教学课件PPT信号的取样.ppt

华中科技大学电信系,1,2.5 信号的取样,2.5.1 连续信号的取样（重点）表示模拟信号，或 表示取样信号，T为取样周期。可把取样器看作每T秒闭合一次的电子开关S，开关每接通一次，便得到一个输出取样值。,模拟信号？数字信号？,华中科技大学电信系,2,2.29,华中科技大学电信系,3,1.理想取样：取样器开关S的闭合时间 的极限情况下，取样脉冲序列变成冲激函数序列p(t)：模拟信号经过取样器，其输出的取样信号为：由于d(t-nT)只在t=nT时取非0值，上式即为,华中科技大学电信系,4,2.频谱延拓研究取样信号 与模拟信号 的频谱之间的关系,将p(t)展成傅里叶级数：,其中：,所以：,华中科技大学电信系,5,理想取样信号 的频谱为：,可见,p(t)的傅里叶变换为：,华中科技大学电信系,6,取样信号的频谱包括原信号频谱和无限个经过平移的原信号频谱，这些频谱系数都要乘以系数1/T,2.30,归一化,华中科技大学电信系,7,设原信号是最高频率为 的带限信号，从上图可见，当 或 时平移后的频谱必然相互重叠，重叠部分的频率成分的幅值与原信号不同，这就是混叠现象。如果原信号不是带限信号，则必然存在混叠现象。当原信号是带限信号时，为使取样后的信号频谱不产生混叠，必须有取样频率大于等于两倍信号的最高频率，即：称为奈奎斯特频率，称为折叠频率。,华中科技大学电信系,8,3.频率归一化研究离散时间信号 的频谱 和取样信号 的频谱 之间的关系。假设离散时间信号是模拟信号通过周期性取样得到的，即（2.64）取样信号的频谱为：,（2.65）,华中科技大学电信系,9,离散时间信号 的FT为：（2.66）与（2.65）比较：将频谱延拓的（2.63）代入上式：（2.67）或（2.68）由（2.67）可知，在 的条件下，离散时间信号的频谱 与取样信号的频谱 相等。,，则有,图,华中科技大学电信系,10,由于（为取样频率）是 对 归一化的结果，因此可认为离散时间信号的频谱是取样信号的频谱经频率归一化后的结果。4.信号重建如果取样信号（或离散时间信号）的频谱不存在混叠，则取样信号（或离散时间信号）通过一个理想低通滤波器可以完全恢复出原信号：,(2.72),或,华中科技大学电信系,11,假设取样信号的频谱不存在混叠，所以取样信号经过一个理想低通滤波器(LPF)可以完全重建，而不损失任何信息。由图2.30(b)可知,将取样信号(或离散时间信号)经过一个截止频率为Ws/2的理想低通滤波器,就可将取样信号(或离散时间信号)频谱中的基带频谱取出来，恢复原来的模拟信号。,理想低通滤波器的频率特性为：,华中科技大学电信系,12,其输出信号的频谱为：,华中科技大学电信系,13,从时域理解：模拟信号:,又因为,所以,华中科技大学电信系,14,由模拟信号与取样信号的关系，即由（2.65）和（2.67），得：,华中科技大学电信系,15,所以输出信号为：其中：称为内插函数。重建信号值=原信号抽样点的值与内插函数乘积之和。内插函数的特性：,Sa函数在抽样点nT上值为1;其余抽样点上值为0。,华中科技大学电信系,16,被恢复的信号 在取样点的值恰好等于原来连续信号 在取样时刻 t=nT的值，而取样点之间的部分由各内插函数的波形叠加而成。如图所示：,华中科技大学电信系,17,本节主要是讨论模拟信号，取样信号，离散时间信号，重建的模拟信号，这四种信号的频谱之间的关系。如果保证不存在混叠，取样信号和离散时间信号都可以通过一个LPF重建模拟信号。重建信号值=原信号抽样点的值与内插函数（序列）乘积之和,华中科技大学电信系,18,离散,思考：如果模拟信号不是低通，而是带通信号，则取样信号无混叠对采样频率有什么要求？,2.30,华中科技大学电信系,19,例,对 进行理想取样，取样间隔T=0.25s，得到,，再让,通过理想低通滤波器G(jW)，,G(jW)用下式表示：,设,，要求,（1）写出 的表达式（2）求出理想低通滤波器的输出信号,解：（1）,华中科技大学电信系,20,（2）,华中科技大学电信系,21,如下图所示，,为,的频谱，,为,的频谱。,由图可知，通过 后其输出信号为：,华中科技大学电信系,22,2.5.2 离散时间信号的取样（1）时域表示：可看作是信号调制,华中科技大学电信系,23,（2）频域表示：,其中:,(ws=2p/N),则:,华中科技大学电信系,24,为了不发生混叠,取样频率应满足条件:,2.34,华中科技大学电信系,25,（3）x(n)序列的恢复：,华中科技大学电信系,26,在xp(n)序列的频谱没有混叠失真的情况下，用一个增益为N，截止频率大于wM而小于(ws-wM)的低通滤波器，对xp(n)进行滤波，可恢复出原信号x(n)。恢复的序列为：其中 即为内插序列。,所恢复的序列xr(n)可以由离散时间取样序列的取样值x(kN)与内插序列hr(n-kN)相乘并求和得到。,华中科技大学电信系,27,2.5.3 离散时间信号的抽取和内插(1)离散时间信号的抽取:(Decimation)这种抽取N的整数倍点上的样本的过程称为抽取。也称为减取样(Down Sampling),将xp(n)的非零取样值抽取出来组成一个新的序列xd(n)。,华中科技大学电信系,28,由于xp(n)在N的整数倍点外的取样值均为0，则,华中科技大学电信系,29,取样序列xp(n)和抽取序列xd(n)的频谱只是频率尺度不同。,2.37,华中科技大学电信系,30,(2)离散时间信号的内插:(Interpolating)也叫增取样（Upsampling）是抽取的逆过程，分2步进行：先在xd(n)的每相邻两个序列值之间插入N-1个零值，得到序列xp(n)；然后用一个理想低通滤波器从xp(n)得到内插后的序列x(n),下图中:,华中科技大学电信系,31,2.39,N=2,华中科技大学电信系,32,2.6 Z变换,2.6.1 Z变换的定义1.Z变换的定义 双边Z变换：极坐标形式：因此，序列在单位圆(r=1)上的Z变换等于序列的傅里叶变换。单边Z变换：X(z)一般情况下是一个有理分式,华中科技大学电信系,33,ZT的收敛域：使 收敛的z值，一般为某个环域：；ZT的零点：使 的z值，在Z平面用“o”表示；ZT的极点：使 的z值，在Z平面用“”表示。ZT存在的条件：,序列的ZT的收敛域和 极零点分布图,华中科技大学电信系,34,例2.16 求序列 的Z变换。解：当，即 时，上列级数收敛，且有。本序列的ZT的收敛域如图所示：,华中科技大学电信系,35,2.6.2 几种特殊序列的Z变换及收敛域（1）有限长序列 其ZT为：,华中科技大学电信系,36,由 可知，在 内，X(z)是有限项级数和，只要级数的每一项有界，有限项的和也就有界。一般情况下，其收敛域为：；因为X(z)为有限项，当X(z)展开时，既有z的正幂项，又有z的负幂项，故：，当 时，收敛域为：（此时X(z)无z的正幂项）当 时，收敛域为：,（此时X(z)无z的负幂项）,华中科技大学电信系,37,（2）右边序列 其ZT为：其收敛域是以 为半径 的圆的外部，即,华中科技大学电信系,38,右边序列的特例因果序列因果序列：的右边序列。其收敛域为：因果序列的ZT无正幂项右边序列的ZT如果无正幂项，则该系统是因果的。对于右边序列，如 存在，则该序列为因果序列。,华中科技大学电信系,39,的右边序列称为非因果序列。,华中科技大学电信系,40,（3）左边序列其ZT为：其收敛域是以 为半径 的圆的内部，即,华中科技大学电信系,41,左边序列的特例逆因果序列逆因果序列：的左边序列。其收敛域为：逆因果序列的ZT无负幂项左边序列的ZT如果无负幂项，则该系统是逆因果的。对于左边序列，如 存在，则该序列为逆因果序列。,华中科技大学电信系,42,的右边序列称为非因果序列。,华中科技大学电信系,43,（4）双边序列 双边序列指 时,序列都有非零值的序列，可看作左边序列和右边序列之和。其ZT为：收敛域为X1(z)和X2(z)的收敛域的公共部分，一般情况下为：如果，双边序列 的ZT无收敛域。,华中科技大学电信系,44,ZT的结论：序列ZT为有理分式的收敛域以极点为边界（包括0，）收敛域内不能包括任何极点，可以包含零点；相同的极零点分布可能分别对应不同的收敛域，即：由同一幅极零点分布图可以得到不同的序列，而不同的序列可能有相同的ZT；见p56的图2.48(a)-(d)收敛域口诀：右外、左内、双环。,华中科技大学电信系,45,例2.17 求 的ZT及其收敛域。解：该序列为双边序列，其ZT为：这是一个有理分式，极点为z1=a，z2=b；零点为z1=0，z2=(a+b)/2，收敛域为一个环域：,华中科技大学电信系,46,2.6.3 逆Z变换已知X(z)及其收敛域,反过来求序列x(n)的变换称作逆Z变换。记作:x(n)=Z-1X(z)（1）幂级数法（长除法）如果序列的ZT能表示成幂级数的形式，则序列x(n)是幂级数 中的z-n的系数。这种方法只对某些特殊的ZT有效。判断右边序列还是左边序列:右边序列：展开成负幂级数左边序列：展开成正幂级数,华中科技大学电信系,47,例2.18 求 的逆变换。解：因为 所以有 又由 可知原序列为右边序列，因此,华中科技大学电信系,48,一些常用的幂级数公式：,华中科技大学电信系,49,若Z变换为有理函数，可用长除法将X(z)展开成幂级数。注意：在使用长除法之前，应先根据给定的收敛域确定对应的原序列是右边序列还是左边序列。若为右边序列（或因果序列），可将X(z)展成负幂级数，即商为z的负幂次。若为左边序列（或逆因果序列），可将X(z)展成正幂级数，即商为z的正幂次。例1：求 的逆变换。,解：因为|z|3，说明对应的序列是右边序列，又因为收敛域包含，所以它是一个因果序列。注意：此时除式与被除式都是降幂排列。,华中科技大学电信系,50,华中科技大学电信系,51,b可正可负,华中科技大学电信系,52,（2）部分分式展开法(重点)如果X(z)是一个有理分式（两个多项式之比），分母的阶次大于分子的阶次，且只有单阶极点，则X(z)可表示成部分分式，即：,其中dk（k=1,N）是X(z)的极点。,分解系数Ak的求解有2种方法：方程组法和留数法。,留数法：,华中科技大学电信系,53,X(z)的收敛域为：以最大极点的模为半径的圆的外部，即：,例2.21 用部分分式法求下列Z变换的逆变换：,解：因为|z|2，说明对应的序列是右边序列，又因为收敛域包含，所以它是一个因果序列。X(z)的2个极点z1=2，z2=0.5。将X(z)展成部分分式：,华中科技大学电信系,54,留数法求系数：,华中科技大学电信系,55,例2.22,求下列Z变换的逆变换：,解：由环形收敛域可知这是一个双边序列。分解得到,华中科技大学电信系,56,（3）留数定理法可以由柯西积分公式导出求逆Z变换的公式。柯西积分公式：其中c为反时针方向环绕原点的围线。ZT的定义式：,由n=k时，利用柯西积分公式得：,逆Z变换的公式,对正的n和负的n都成立。,华中科技大学电信系,57,对于有理z变换，这个围线积分可以用留数定理计算。设在有限Z平面上，ak(k=1,2,N)是在围线c内部的极点集，bk(k=1,2,M)是在围线c外部的极点集。围线c内的极点一般对应于一个因果序列，因此，当 n0时：围线c外的极点一般对应于一个逆因果序列，因此，当 n0时,因为当X(z)zn-1在z=处有二阶或二阶以上的零点，即X(z)zn-1的分母多项式的阶数比分子多项式的阶数高二阶或二阶以上时，无穷远处的留数为零，故有：,或：,华中科技大学电信系,58,如果X(z)zn-1是z的有理函数，且z=z0处有s阶极点，即 其中(z)在z=z0处无极点，那么,特别地，s=1时：,推荐使用：幂级数法（长除法）和部分分式展开法，其中部分分式展开法要求重点掌握。,例题见p61的例2.23和例2.24。,华中科技大学电信系,59,2.6.4 Z变换的性质和定理（1）线性：设,，则：,收敛域：,一般情况下，收敛域变小。但在组合ZT可能出现新的零极点抵消的情况时，收敛域可能增大。,如：，收敛域：,但，收敛域：整个Z平面。,华中科技大学电信系,60,例1 已知,求其z变换。,解：,华中科技大学电信系,61,（2）序列移位：,，则：,设：,注意：收敛域在 处有例外，如 收敛域为整个Z平面。但 在 处不收敛，在 处不收敛。,例2 求 的逆变换。,方法1：应用时移性质，,华中科技大学电信系,62,方法2：,华中科技大学电信系,63,（3）乘以指数序列an：,，则：,设：,a为正实数：z平面尺度的缩小或扩大，极零点沿径向移动。,a为模为1的复数：z平面的旋转，极零点沿以原点为中心 的圆周移动。,例3 求 x(n)=rncos(w0n)u(n)的z变换。,解：,然后利用指数相乘性质，得到：,华中科技大学电信系,64,再由线性性质有：,华中科技大学电信系,65,（4）序列的折叠：,，则：,设：,（5）复序列的共轭：,，则：,设：,证明：,对实序列x(n)：,华中科技大学电信系,66,（6）X(z)的微分：,，则：,设：,证明：,对其两端求导得:,华中科技大学电信系,67,例4 求 x(n)=nanu(n)的z变换。,解：,x(n)=nanu(n)=nanu(n),由微分性质和幂指数乘积性质:,华中科技大学电信系,68,（7）初值定理：,证明：,对于因果序列x(n),则,对于逆因果序列x(n),则,显然，,华中科技大学电信系,69,（8）终值定理：,若x(n)是因果序列，且X(z)除在z=1处有一阶极点外，其它全部极点都在单位圆内，则：,证明：,由于(z-1)X(z)抵消了函数X(z)在z=1处的可能极点，故(z-1)X(z)的收敛域将包括单位圆，所以对上式两端求极限得：,华中科技大学电信系,70,（9）序列的卷积：,设,，则：,，,华中科技大学电信系,71,证明：,令m=n-k,则:,W(z)的收敛域是X(z)和Y(z)收敛域的重合部分.,华中科技大学电信系,72,例5.,解：,华中科技大学电信系,73,（10）复卷积定理：,设,，,则:,式中c是v平面收敛域中任一条环绕原点的逆时针方向的闭合曲线，v平面的收敛域为：,z平面的收敛域为：,vk是围线c包含的所有极点。,华中科技大学电信系,74,例2.25 设x(n)=anu(n)，y(n)=bnu(n)，相应的z变换为：,求w(n)=x(n)y(n)的z变换。,解：由复卷积定理,得：,极点：v1=b，v2=z/a,华中科技大学电信系,75,由复卷积定理，v平面的收敛域应为,即,由右图可见，只有极点v1=b在围线c1内，由留数法计算W(z)：,华中科技大学电信系,76,12.帕塞瓦尔定理(parseval),设,，则:,，且：,，,式中c是 收敛域中环绕逆时针方向的围线，v平面的收敛域由下式确定：,华中科技大学电信系,77,几点说明:,Parseval公式的物理意义：在时域中对序列求能量与在频域中对频谱求能量是一致的。,华中科技大学电信系,78,2.6.5 Z变换与拉氏变换的关系,函数的LT为:,(1)模拟信号（连续时间信号）和取样信号的LT的关系：,上式两边取LT得：,其中：,华中科技大学电信系,79,即：连续时间信号经理想取样后，取样信号的LT是原信号的LT在S平面上沿虚轴的周期延拓。,(2)取样信号的LT和离散时间信号的ZT的关系：,则，取样信号的LT为：,离散时间信号的ZT为：,华中科技大学电信系,80,X(n)=Xa(nT),即：当z=esT时，离散时间信号的ZT等于取样信号的LT。,S平面用直角坐标表示为：s=s+jW,Z平面用极坐标表示为：z=rejw,z=esT时,s=0,即S平面的虚轴 r=1,即Z平面单位圆；,s0,即S的左半平面 r1,即Z的单位圆内；,s0,即S的右半平面 r1,即Z的单位圆外；,r=esT,华中科技大学电信系,81,当 时，,当 时，,当 时，,即幅角旋转一周，将整个Z平面映射一次。,当,时,w=WT,华中科技大学电信系,82,小结：,(1)模拟信号（连续时间信号）和取样信号的LT的关系：,(2)取样信号的LT和离散时间信号的ZT的关系：,S平面的每2p/T的条带，映射整个Z平面一次。,华中科技大学电信系,83,2.7 系统函数,描述线性非移变系统的方式：,线性常系数差分方程、单位取样响应、频率响应描述、系统函数。,设x(n)、y(n)和h(n)分别是线性非移变系统的输入、输出和单位取样响应，X(z)、Y(z)和H(z)分别表示相应的Z变换。系统函数定义为,它是单位取样响应h(n)的Z变换。,1、系统函数的定义,华中科技大学电信系,84,2、系统函数与系统差分方程的关系,线性非移变系统可以用线性常系数差分方程描述：,对上式两边求Z变换，利用线性性质和时不变性质，得,因此,可见系统函数分子分母的系数也正是其差分方程的系数。,华中科技大学电信系,85,系统函数还可以进一步分解成：,式中，dk)和cr分别表示H(z)在z平面上的极点和零点。这样，系统函数可以用z平面上的极点、零点和常数A来确定。,例1 根据系统函数:,求该系统的差分方程。,华中科技大学电信系,86,为了求满足该系统输入输出的差分方程，可以将H(z)的分子和分母各因式乘开，而得到如下的形式：,于是，,其差分方程就是,华中科技大学电信系,87,3、系统函数的收敛域与系统的稳定性,为了使h(n)的z变换存在,则：,当|z|=1时，上式变成,这就是系统稳定的充要条件。,因此，若系统函数在单位圆上收敛，则系统是稳定的。H(z)的收敛域包括单位圆 系统是稳定的。,请思考，一个稳定的因果系统的系统函数的收敛域应该是怎么样的？,华中科技大学电信系,88,一个稳定的因果系统的系统函数的收敛域应该是,例2.26（重点）设一个线性非移变系统的系统函数为,试画出极-零点分布图，并确定H(z)的收敛域和稳定性。,华中科技大学电信系,89,解对H(z)的分母进行因式分解得,极点为z1=-1/4,z2=-1/2；零点为z1=0,z2=1/2。,(1)若收敛域是极点z2=-1/2所在的圆的外部区域，且,那么系统是因果的。,华中科技大学电信系,90,(2)若收敛域选的是极点z1=-1/4所在的圆的内部区域，且,那么系统是逆因果的，系统的收敛域为,因为收敛域没有包含单位圆，所以系统是不稳定的。,(3)若收敛域是极点z1=-1/4与z2=-1/2所在的两个圆之间的环域，即,又因为单位圆没有包含在收敛域中，所以系统是不稳定的。,系统的收敛域为,因为该收敛域包含了单位圆，所以系统是稳定的。,那么系统是非因果的。,华中科技大学电信系,91,从上例可知，因果性和稳定性并不存在必然联系。,要使线性时不变系统既因果又稳定，相应系统函数的收敛域必须是：位于最外围极点的外面，又包括单位圆。,也就是说，要求该系统函数的全部极点都在单位圆内。,华中科技大学电信系,92,系统的稳定性判别方法总结：,直接由定义判别：若|x(n)|M，则|y(n)|对于线性非移变系统:由其单位取样响应是否绝对可和由其系统函数的收敛域是否包含单位圆,华中科技大学电信系,93,三个方法的等价性：,若|x(n)|M，则|y(n)|,所以,华中科技大学电信系,94,若H(z)收敛域包含单位圆，,华中科技大学电信系,95,系统的因果性判别方法总结：,直接由定义判别：y(n)是否取决于将来的输入x(n+1)等。对于线性非移变系统:由其单位取样响应 h(n)=0，n0由其系统函数的收敛域是否包含,华中科技大学电信系,96,三个方法的等价性：,y(n)是否取决于将来的输入x(n+1)等。,例,满足,h(n)=0，n0,所以，收敛域包含,华中科技大学电信系,97,4.系统的频率响应,又因为稳定系统的系统函数的收敛域包括单位圆，则对于稳定系统:,华中科技大学电信系,98,例2.27 设一因果IIR系统如图所示，试确定描述系统的差分方程，极零点分布图和频率响应。,解:由图可得系统的差分方程和系统函数分别为:,华中科技大学电信系,99,系统的极点为z1=0.9j，z2=-0.9j,零点为z1=0，z2=1。极零点分布图如右图示。,因为该系统是因果系统,所以系统的收敛域为|0.9j|z|,可见收敛域包含单位圆,故系统是稳定的,可以有:,华中科技大学电信系,100,z=1处是零点，则|H(ejw)|=|H(e0)|=0,又z=0.9j处是极点，则|H(ejw)|在p/2附近升至峰值。,幅度响应和相位响应图为：,华中科技大学电信系,101,例:下图所示表示一个因果的线性非移变系统，求：（1）该系统的系统函数H(z)，并画出极零点图；如果该系统是稳定的，求b的取值范围；（2）求系统对以下输入的稳态响应:,华中科技大学电信系,102,（1）,b的取值范围为：0|b|1。,华中科技大学电信系,103,（2）系统的频率响应为：,该系统对,的响应为：,华中科技大学电信系,104,5.系统极零点对系统频率响应的影响,影响系统特性的是零点cr和极点dk的分布。,结论：极点位置主要影响频响的峰值位置及尖锐程度，零点位置主要影响频响的谷点位置及形状。,华中科技大学电信系,105,2.8 全通系统与最小相位系统,2.8.1 全通系统用于滤波器设计，信道相位均衡等。,华中科技大学电信系,106,华中科技大学电信系,107,极零点互为倒数的共轭对,a(a*)-1,两种方法判断是否全通系统：1.极零点互为倒数的共轭对2.,华中科技大学电信系,108,2.8.2 最小相位系统,要求系统对输入信号的处理具有最小相位延迟。,所有极零点都在单位圆内的系统，即最小相位系统（最小相位延迟系统）。,一个因果系统不一定是最小相位的，但所有最小相位系统都是因果的。,华中科技大学电信系,109,任何系统都可以表示为一个最小相位系统和一个全通系统的级联：,证明：假设H(z)有一个零点z=1/c*在单位圆外,|c|1,其余零点都在单位圆内，则,单位圆外的极点可作类似处理。,华中科技大学电信系,110,2.8.3 相位延迟和群延迟,若信号x(n)经过系统延迟了时间m后为y(n)=x(n-m)，则,可认为y(n)是x(n)经过滤波器 后的输出信号。若滤波器H(ejw)的相位响应为f(w),则,反映了信号的延迟时间,华中科技大学电信系,111,相位延迟,群延迟,表示输入频率为w的单一正弦波的延迟时间,反映了某一频率w邻域内（窄带信号）的延迟性质，或者说反映了对某一频率的包络的延迟时间。,定义,华中科技大学电信系,112,w1=0.85p,w2=0.25p,w3=0.5p.,群延迟的例子,华中科技大学电信系,113,滤波器的群延迟和幅度响应,w1=0.85p,w2=0.25p,w3=0.5p.,华中科技大学电信系,114,t(0.25p)200,t(0.5p)50,输入,输出,