实验三离散傅立叶变换及其特性验证.doc

实验名称：实验三 离散傅立叶变换及其特性验证 一、实验目的 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）式取转置得： 2、离散傅立叶变换 一个有限长序列的离散傅立叶变换对定义为： （4） （5） 以列向量x和X形式排列x(n)和X(k)，则式（4）、（5）可写成： X=WNx 其中矩阵WN由下式给出： （1）线性性质： 注意：若x1(n)和x2(n)分别是N1点和N2点的序列，则选择N3= max (N1, N2)，将它们作N3点DFT处理。 （2） 周期性：离散傅立叶变换(DFT)是周期序列DFS取主值区间形成的，因此序列及其DFT 具有特性和。通常将结果间的量值表示在k的负值区间。 （3）对称性：实序列的离散傅立叶变换可以表示为，其中实部为偶对称，虚部为奇对称，幅值为偶对称，相位为奇对称。 如果序列是实偶对称序列，则也是实偶对称，即；如果序列是实奇对称序列，则是虚奇对称，即；如果序列是虚偶对称序列，则也是虚偶对称，即；如果序列是虚奇对称序列，则是实奇对称，即。 根据上述关系，对于实序列，则有；对于纯虚序列，则有。 三、实验内容 （1）将实指数函数抽样，作DFT，并作出实部、虚部和幅频、相频特性曲线。 用到的M文件： function [Xk] = dft(xn,N) 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; WNnk = WN .^ nk; % DFT matrix Xk = xn * WNnk; 程序和运行的结果： >> n=[0:127]; Ts=1./32; t=n.*Ts; Xn=exp(-t); Xk=dft(Xn,128); subplot(4,1,1); stem(n,real(Xk),'r');title('实部图形') subplot(4,1,2); stem(n,imag(Xk),'b');title('虚部图形') subplot(4,1,3); stem(n,abs(Xk),'og');title('幅度图形') subplot(4,1,4); stem(n,angle(Xk),'m');title('相位图形') 抽样程序与结果 n=[0:0.5:63]; Ts=1./64; t=n.*Ts; Xn=exp(-t); Xk=dft(Xn,127); stem(n,Xn) （2）将图3-2中的两个连续函数经抽样，作DFT，验证前述的四种奇偶特性，并作出幅频和相频特性曲线。 第一个图的程序及结果：