MATLAB基础与应用教程(人民邮电出版社-蔡旭辉)第七章b.pptVIP

7.3 离散傅里叶变换（DFT） 一、DFT的定义 设x(n)是一个长度为M的有限长序列， 则定义x(n)的N点离散傅里叶变换为 X(k)的离散傅里叶逆变换为 式中， ，N称为DFT变换区间长度， 二、DFT和Z变换的关系 设序列x(n)的长度为N， 其Z变换和DFT分别为： 实际上， 任何周期为N的周期序列 都可以看作长度为N的有限长序列x(n)的周期延拓序列， 即 7.3 离散傅里叶变换（DFT） 定义DFT： 用类似于例7.9中的方法，可把（7.3）式写成矩阵乘法运算。 Xk=xn*Wnk 其中，xn为序列行向量，Wnk是一 N×N 阶方阵， 而 称为旋转因子。 Wnk 的MATLAB表示： Wnk = WN.^([0:N-1]’* [0:N-1]) 用矩阵乘法计算N点DFT的程序如下: ? MATLAB程序q73a.m %用矩阵乘法计算N点DFT clear;close all xn=input(请输入序列x= ); N = length(xn); % n=0:N-1; k=n; nk=n*k; % 生成N×N方阵 WN=exp(-j*2*pi/N); % 旋转因子 Wnk=WN.^nk; % 产生旋转因子矩阵 Xk=xn*Wnk; % 计算N点DFT 这种方法计算DFT概念清楚、编程简单，但占用内存大、运行速度低，所以不实用。MATLAB基础部分提供了fft、ifft、fft2、ifft2等等快速计算傅里叶变换的函数，使DFT的运算速度量提高了若干数量级，在后面的例题中均直接调用这些函数。 例7.15 序列的离散傅立叶变换 求复正弦序列 余弦序列 正弦序列 的离散傅立叶变换，分别按N =16和N =8进行计算。绘出幅频特性曲线， 进行比较讨论。 解：直接产生序列x1n、x2n和x3n，调用fft函数求解 在截取16点时，得到的是完整的余弦波形；而截取8点时，得到的是半截的余弦波形，当然有大量的谐波成分。 例7.16 验证N点DFT的物理意义 （1） 绘出幅频曲线和相频曲线。 （2）计算并图示x(n)的8点DFT。 （3）计算并图示x(n)的16点DFT。 解: 序列x(n)的Ｎ点DFT的物理意义是 在[0，2?]上进行Ｎ点等间隔采样。 程序先密集采样，绘制出幅频曲线图。然后再分别做8点和16点DFT来验证这个采样关系。 程序略。 频域采样定理 频域采样是指对有限长序列的傅氏变换在频域离散化得到X(k)的过程。 本节讨论两个基本问题： a：频域采样（DFT）不失真的条件，即由 X (k)不失真地恢复x(n)的条件 b：用X (k)表示X (z)和 的插值公式(内插公式) 图3-3-1 时域恢复示意图 结论：若序列长度为L，频域采样点数（DFT的长度）为N，且L≤N，则频域采样后可不失真地恢复原序列 ；但若LN，则频域采样后不能不失真地恢复原序列 。 例7.17 频域与时域采样对偶性 （1）产生三角波序列 （2）对M = 40，计算x(n)的64点DFT，并图示x(n)和 X(k) = DFT[x(n)]，k = 0, 1, …, 63。 （3）对（2）中所得X(k)在 [0，2?] 上进行32点抽样得 （4）求 的32点IDFT，即 （5）绘出 的波形图，评述它与x(n)的关系。 % 第七章例7.17程序q717 % 时域与频域采样的对偶性验证 clear;close all M=40;N=64;n=0:M; %产生M长三角波序列x(n) xa=0:floor(M/2); xb= ceil(M/2)-1:-1:0; xn=[xa,xb]; Xk=fft(xn,64);%64点FFT[x(n)] X1k=Xk(1:2:N);%隔点抽取Xk得到X1(K) x1n=ifft(X1k,N/2);%32点IFFT[X1(k)]得到x1(n) nc=0:3*N/2; % 取97点为观察区 xc=x1n(mod(nc,N/2)+1);%x1(n)的周期延拓序列 subplot(3,2,1);stem(n,xn,.) title(40点三角波序列x(n));xlabel(n);ylabel(x(n)) k=0:N-1; subplot(3,2,3);stem(k,abs(Xk),.) title(64点DFT[x(n)]);xlabel(k);ylabel(|X(k)|) k