《DSP第十一讲》PPT课件

第十一讲,3.2.5DFT的共轭对称性3.3频域抽样理论-抽样Z变换3.4.1用DFT计算线性卷积,要点,为什么要定义圆周对称？DFT对称性的特点？频域抽样提出的背景？,3.2.5DFT的共轭对称性,与DTFT对称性的区别DTFT以(-,+)为变换空间，所以在讨论对称性质中，以原点为对称中心，序列的移位范围无任何限制，因为无论如何不会移出变换区间；DFT以(0,N-1)为变换空间，所以在讨论对称性质中，序列的移位会移出变换区间，所以要在区间(0,N-1)上定义有限长序列的共轭对称序列和反对称序列；DFT以(0,N-1)为变换空间，所以在讨论对称性质中，将会得出其对称中心为n=N/2。,1.有限长序列的共轭对称分量与共轭反对称分量有限长序列的共轭对称分量与共轭反对称分量分别定义为：,由于,所以,这表明长为N的有限长序列可分解为两个长度相同的两个分量。,1.有限长共轭对称序列与共轭反对称序列,上式已给出有限长共轭序列对称共轭反对称序列的对称中心为n=N/2，任意有限长序列其圆周共轭对称分量与圆周共轭反对称分量可简写为：,共轭对称与共轭反对称序列示意图,有限长序列x(n)的实、虚分解及其DFT对称性,则有：,证明：,2.DFT的共轭对称性,有限长序列x(n)的对称分量分解及其DFT对称性,参见式（3.2.12）,实、纯虚序列的对称特性,当x(n)为实序列时，则X(k)=Xep(k)又据Xep(k)的对称性：,当x(n)为纯虚序列时，则X(k)=Xop(k)又据Xop(k)的对称性：,(1)X(k)共轭对称，即X(k)=X*(N-k)k=0,1,N-1(2)如果x(n)是实偶对称序列，即x(n)=x(Nn)，则X(k)实偶对称，即X(k)=X(Nk)(3)如果是奇对称序列，即x(n)=x(Nn)，则X(k)纯虚奇对称，即X(k)=X(Nk),实序列的对称特性小结,序列DFT,有限长序列的共轭对称性总结1：复数序列的共轭对称性,序列DFT,有限长序列共轭对称性总结2:实数序列的共轭对称性,有限长序列共轭对称性总结3:纯虚序列的共轭对称性,序列DFT,例3.2.2假设x1(n)和x2(n)都是N点的实数序列，设想用一次N点DFT运算来计算它们各自的DFT:,共轭对称性的应用,解：,由(3.2.17)、(3.2.18)和(3.2.19)式得到：所以，由X(k)可以求得两个实序列x1(n)和x2(n)的N点DFT：,3.3频率域采样频域采样定理讨论：时域采样:对一个频带有限的信号,根据抽样定理对其进行抽样,所得抽样信号的频谱是原带限信号频谱的周期延拓,因此,完全可以由抽样信号恢复原信号。频域采样:对一有限序列(时间有限序列)进行N点DFT所得x(k)就是序列傅氏变换的采样.所以DFT是频域N点抽样的结果。,问题:,能否由频域抽样X(k)恢复序列x(n)能否由频域抽样X(k)恢复x(z)或若能恢复其条件是什么？如何推导频域内插恢复公式?,回忆时域内插恢复公式!,X(n)为M点的有限长序列。,IDFTX(k)=XN(n),FT,DTFT,DFS,讨论之前先明确一些概念,x(n),一.由频域抽样恢复原序列,一.由频域抽样恢复原序列（续）,频域抽样时域以N点为周期进行延拓的主值区间,原序列,X(z)在单位圆上的N点等间隔采样所得到的周期序列X(k)的IDFS是原序列x(n)以N为周期进行延拓的周期序列。,x(n)为无限长序列混叠失真x(n)为有限长序列，长度为M,由于时域抽样造成频域周期延拓，同样，频域抽样造成时域周期延拓。,分两种情况讨论周期延拓是否造成混叠失真：,频域采样定理,若序列长度为M，则只有当频域采样点数:时，才有即可由频域采样不失真地恢复原信号，否则产生时域混叠现象。,1.由X(k)恢复X(Z),则：,内插公式与内插函数,用频域采样恢复的内插公式,2.,记住此公式，第七章数字滤波器的设计中，我们将会看到，该公式提供了一种有用的滤波器结构和滤波器设计途径。,【例3.3.1】长度为26的三角形序列x(n)如图3.3.1(a)所示。编写MATLAB程序验证频域采样理论。解解题思想：先计算x(n)的32点DFT，得到其频谱函数X(ej)在频率区间0，2上等间隔32点采样X32(k)，再对X32(k)隔点抽取，得到X(ej)在频率区间0，2上等间隔16点采样X16(k)。最后分别对X16(k)和X32(k)求IDFT,得到：绘制x16(n)和x32(n)波形图验证频域采样理论。,MATLAB求解程序ep331.m如下：%数字信号处理(第三版）第3章例3.3.1程序ep331.%频域采样理论验证M=26;N=32;n=0:M;xa=0:M/2;xb=ceil(M/2)-1:-1:0;xn=xa,xb;%产生M长三角波序列x(n)Xk=fft(xn,512);%512点FFTx(n)X32k=fft(xn,32);%32点FFTx(n)x32n=ifft(X32k);%32点IFFTX32(k)得到x32(n)X16k=X32k(1:2:N);%隔点抽取X32k得到X16(k)x16n=ifft(X16k,N/2);%16点IFFTX16(k)得到x16(n)以下绘图部分省略。,图3.3.1频域采样定理验证,3.4DFT的应用举例,3.4.1用DFT计算线性卷积3.4.2用DFT进行信号的谱分析,3.4.1用DFT计算线性卷积,0kL-1,则由时域循环卷积定理有Y(k)=DFTy(n)=H(k)X(k),0kL-1,如果,1.用DFT计算循环卷积,由此可见，循环卷积既可在时域直接计算，也可在频域计算。由于DFT有快速算法FFT，当N很大时，在频域计算的速度快得多，因而常用DFT(FFT)计算循环卷积。,图3.4.1用DFT计算循环卷积的原理框图,背景及意义：在实际应用中，为了分析LSI系统或者对序列进行滤波处理时，需要计算两个序列的线性卷积。为了提高运算速度，也希望用DFT(FFT)计算线性卷积。而与DFT对应的是循环卷积，为此需导出线性卷积和循环卷积之间的关系以及循环卷积与线性卷积相等的条件。假设h(n)和x(n)都是有限长序列，长度分别是N和M。它们的线性卷积和循环卷积分别表示如下：,2.循环卷积与线性卷积,长度为N+M-1,长度为L,其中，LmaxN,M,可以看出，上式中,线性卷积与循环卷积的关系,图3.4.2线性卷积与循环卷积,线性卷积与循环卷积图示,x(n)=1，1，1，1h(n)=1，1，1，1，1线性卷积y(n)=1，2，3，4，4，3，2，16点圆周卷积X(n)=1，1，1，1，0，0h(0-m)=1，0，1，1，1，1y(0)=3h(1-m)=1，1，0，1，1，1y(1)=3h(2-m)=1，1，1，0，1，1y(2)=3h(3-m)=1，1，1，1，0，1y(3)=4h(4-m)=1，1，1，1，1，0y(4)=4h(5-m)=0，1，1，1，1，1y(5)=3,-6-5-4-3-2-1012345671，2，3，4，4，3，2，11，2，3，4，4，3，2，13，3，3，4，4，3,直接做圆周卷积,利用圆周卷积与线性卷积的关系,图3.4.3用DFT计算线性卷积框图,MN。若仍选取LNM1，以L为循环卷积区间，并用上述快速卷积法计算线性卷积，则要求对短序列补很多零点，而且长序列必须全部输入后才能进行快速计算。因此要求存储容量大，运算时间长，并使处理延时很大，不能实现实时处理。通常采用分段卷积。若将x(n)均匀分段，每段长度取M，则,于是，h(n)与x(n)的线性卷积可表示为,3.长序列的分段卷积,图3.4.4重叠相加法卷积示意图,用DFT计算分段卷积yk(n)的方法：,（1）i=0；L=NM1；计算并保存H(k)=DFTh(n)L;（2）读入xi(n)=x(n)RM(nkM)，构造变换区间0，L-1上的序列，实际中就是将xi(n)的M个值存放在长度为M的数组中,并计算（3）；（4），n=0,1,2,L1；（5）计算：（6）i=i1，返回(2)。应当说明，一般x(n)是因果序列，假设初始条件y-1(n)=0。,MATLAB信号处理工具箱中提供了一个函数fftfilt，该函数用重叠相加法实现线性卷积的计算。调用格式为：y=fftfilt(h,x，M)。式中,h是系统单位脉冲响应向量；x是输入序列向量；y是系统的输出序列向量（h与x的卷积结果）；M是由用户选择的输入序列x的分段长度，缺省M时，默认输入序列x的分段长度M=512。【例3.4.1】假设h(n)=R5(n)，x(n)=cos(n/10)+cos(2n/5)u(n)，用重叠相加法计算y(n)=h(n)*x(n)，并画出h(n)、x(n)和y(n)的波形。,解h(n)的长度为N=5,对x(n)进行分段，每段长度为M=10。计算h(n)和x(n)的线性卷积的MATLAB程序如下：%例3.4.1重叠相加法的MATLAB实现程序：ep341.mLx=41;N=5;M=10;%Lx为信号序列x(n)长度hn=ones(1,N);hn1=hnzeros(1,Lx-N);%产生h(n)，其后补零是为了绘图好看n=0:L-1;xn=cos(pi*n/10)+cos(2*pi*n/5);%产生x(n)的Lx个样值yn=fftfilt(hn,xn,M);%调用fftfilt