离散傅里叶变换课件.ppt

离散傅里叶变换;;　　因此，时域离散系统的研究与应用在许多方面取代了传统的连续时间系统。所以说，DFT不仅在理论上有重要意义，而且在各种信号的处理中亦起着核心作用。 　　本章主要讨论DFT的定义、物理意义、基本性质以及频域采样和DFT的应用举例等内容。;;;;; 　　设变换区间N=8，则  　由此例可见，x(n)的离散傅里叶变换结果与变换区间长度N的取值有关。对DFT与Z变换和傅里叶变换的关系及DFT的物理意义进行讨论后，上述问题就会得到解释。 ;; 　　(3.1.3)式表明序列x(n)的N点DFT是x(n)的Z变换在单位圆上的N点等间隔采样。(3.1.4)式则说明X(k)为x(n)的傅里叶变换X(ejω)在区间[0, 2π]上的N点等间隔采样。这就是DFT的物理意义。由此可见，DFT的变换区间长度N不同，表示对X(ejω)在区间[0, 2π]上的采样间隔和采样点数不同，所以DFT的变换结果不同。上例中， x(n)=R4(n)，DFT变换区间长度N分别取8、16时，X(ejω)和X(k)的幅频特性曲线图如图3.1.1所示。由此容易得到x(n)=R4(n)的4点DFT为X(k)=DFT[x(n)]4=4δ(k)，这一特殊的结果在下面将得到进一步解释。 　　;图3.1.1 R4(n)的FT和DFT的幅度特性关系 ;;;; 图3.1.2 x(n)及其周期延拓序列;;;　　现在解释DFT[R4(n)]4=4δ(k)。根据DFT第二种物理解释可知，DFT[R4(n)]4表示R4(n)以4为周期的周期延拓序列R4((n))4的频谱特性，因为R4((n))4是一个直流序列，只有直流成分（即零频率成分）。 ;3.1.4 用MATLAB计算序列的DFT 　　MATLAB提供了用快速傅里叶变换算法FFT(算法见第4章介绍)计算DFT的函数fft，其调用格式如下: Xk = fft (xn, N)； 调用参数xn为被变换的时域序列向量，N是DFT变换区间长度，当N大于xn的长度时，fft函数自动在xn后面补零。函数返回xn的N点DFT变换结果向量Xk。当N小于xn的长度时，fft函数计算xn的前面N个元素构成的N长序列的N点DFT，忽略xn后面的元素。  　　Ifft函数计算IDFT，其调用格式与fft函数相同，可参考help文件。;　　【例3.1.2】 设x(n)=R4(n)，X(ejω)=FT[x(n)]。分别计算X(ejω)在频率区间［0，2π］上的16点和32点等间隔采样，并绘制X(ejω)采样的幅频特性图和相频特性图。 　　解 由DFT与傅里叶变换的关系知道，X(ejω)在频率区间［0??2π］上的16点和32点等间隔采样，分别是x(n)的16点和32点DFT。调用fft函数求解本例的程序ep312.m如下:　 ;　% 例3.1.2程序ep312.m 　% DFT的MATLB计算 　xn=[1 1 1 1]; %输入时域序列向量xn=R4(n) 　Xk16=fft(xn, 16); %计算xn的16点DFT 　Xk32=fft(xn, 32); %计算xn的32点DFT 　%以下为绘图部分（省略，程序集中有） 程序运行结果如图3.1.3所示。;;3.2 离散傅里叶变换的基本性质 3.2.1 线性性质 　　如果x1(n)和x2(n)是两个有限长序列，长度分别为N1和N2，且 y(n)=ax1(n)+bx2(n) 式中，a、b为常数，取N=max[N1, N2], 则y(n)的N点DFT为 Y(k)=DFT[y(n)]N=aX1(k)+bX2(k) 0≤k≤N－1 　　　　　　　　　　　　　　　　 (3.2.1) 其中X1(k)和X2(k)分别为x1(n)和x2(n)的N点DFT。 ;;　　显然，y(n)是长度为N的有限长序列。观察图3.2.1可见，循环移位的实质是将x(n)左移m位，而移出主值区(0≤n≤N-1)的序列值又依次从右侧进入主值区。“循环移位”就是由此得名的。 　　由循环移位的定义可知，对同一序列x(n)和相同的位移m，当延拓周期N不同时，y(n)=x((n+m))NRn(n)则不同。请读者画出N = M=6，m=2时，x(n)的循环移位序列y(n)波形图。;图3.2.1 x(n)及其循环移位过程;;;;;式中，L称为循环卷积区间长度，L≥max[N，M]。上式显然与第1章介绍的线性卷积不同，为了区别线性卷积，用 　表示循环卷积，用　表示L点循环卷积，即yc(n)=h(n)　 x(n)。观察(3.2.5)式，x((n－m))L是以L为周期的周期信号，n和m的变化区间均是[0,