数字信号处理实验离散傅里叶变换和其特性验证.docVIP

数字信号处理实验报告 实验名称： 离散傅里叶变换及其特性验证 学号： 姓名： 评语： 成绩： 一、实验目的 1、掌握离散时间傅立叶变换（DTFT）的计算方法和编程技术。 2、掌握离散傅立叶变换（DFT）的计算方法和编程技术。 3、理解离散傅立叶变换（DFT）的性质并用MATLAB进行验证。 二、实验原理与计算方法 1、离散时间傅立叶变换 如果序列x(n)满足绝对可和的条件，即，则其离散时间傅立叶变换定义为： （1） 如果x(n)是无限长的，则不能直接用MATLAB由x(n)计算X(ej?)，但可以用它来估计X(ej?)表达式在[0,π]频率区间的值并绘制它的幅频和相频（或实部和虚部）曲线。 如果x(n)是有限长的，则可以用MATLAB对任意频率?处的X(ej?)进行数值计算。如果要在[0,π]间按等间隔频点估计X(ej?)，则（1）式可以用矩阵－向量相乘的运算来实现。假设序列x(n)在(即不一定在[0, N-1])有N个样本，要估计下列各点上的X(ej?)： 它们是[0,π]之间的(M+1)个等间隔频点，则（1）式可写成： （2） 将{x(nl)}和{X(ej?k)}分别排列成向量x和X，则有： X=Wx （3） 其中W是一个(M+1)×N维矩阵： 将{k}和{n}排成列向量，则 在MATLAB中，把序列和下标排成行向量，对（3）式取转置得： 其中nTk是一个N×(M+1)维矩阵。用MATLAB实现如下： k=[0:M]; n=[n1:n2]; X=x*(exp(-j*pi/M)).^(n’*k); 2、离散傅立叶变换 一个有限长序列的离散傅立叶变换对定义为： （4） （5） 以列向量x和X形式排列x(n)和X(k)，则式（4）、（5）可写成： X=WNx 其中矩阵WN由下式给出： 可由下面的MATLAB函数dft和idft实现离散傅立叶变换运算。 function [Xk] = dft(xn,N) % Computes Discrete Fourier Transform % ----------------------------------- % [Xk] = dft(xn,N) % Xk = DFT coeff. array over 0 = k = N-1 % xn = N-point finite-duration sequence % N = Length of DFT % n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN .^ nk; % DFT matrix Xk = xn * WNnk; % row vector for DFT coefficients function [xn] = idft(Xk,N) % Computes Inverse Discrete Transform % ----------------------------------- % [xn] = idft(Xk,N) % xn = N-point sequence over 0 = n = N-1 % Xk = DFT coeff. array over 0 = k = N-1 % N = length of DFT % n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = ex