DFT和FFT频谱分析.docx

上海第二工业大学 实验报告 DFT和FFT频谱分析 课程名称 数字信号处理 姓 名 刘万军 班 级 13通信A1 学 号 20134821333 2015年12月25日 DFT和FFT频谱分析 实验目的 掌握DFT频谱分析的原理与编程方法。 理解FFT算法的编程思想。 熟练掌握利用FFT对信号做频谱分析，包括正确地进行参数选择，画频谱及读频谱图。 利用FFT频谱分析进行快速卷积。 实验环境 Windows xp2操作系统 安装MATLAB2007a软件 实验原理 1.离散傅里叶变换（DFT） 离散傅里叶变换（Discrete?Fourier?Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式，将信号的时域采样变换为其DTFT的频域采样。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作其周期延拓的变换。在实际应用中通常采用快速傅里叶变换计算DFT。 DFT的定义是针对任意的离散序列中的有限个离散抽样的，它并不要求该序列具有周期性。 由DFT求出的离散谱是离散的周期函数，周期为、离散间隔为。离散谱关于变元k的周期为N。 在时域和频域范围内的N点分别是各自的主值区间或主值周期。 X(n)的的N点离散傅里叶变换： X(K)的傅里叶反变换： 2.快速傅里叶变换（FFT) FFT是一种DFT的高效算法，称为快速傅立叶变换（fast Fourier transform）。 有限长序列通过离散傅里叶变换(DFT)将其频域离散化成有限序列.但其计算量太大(与 N的平方成正比, 4N2次实数相乘和 2N(2N-1)次实数相加), 很难实时地处理问题, 因此引出了快速傅里叶变换(FFT)。 FFT并不是一种新的变换形式,它只是 DFT的一种快速算法.并且根据对序列分解与选取方法的不同而产生了 FFT 的多种算法。 8点DNT-FFT蝶形运算图： 程序设计 （1）根据公式设计DFT原理程序如下： clear all xn=input(x(n)=); M=length(xn); N=input(变换区间N=); xn=[xn zeros(1,N-M)]; n=0:N-1;k=0:N-1; nk=n*k; wn=exp(-j*2*pi/N); wnK=wn.^nk; xk=xn*wnK subplot(211);stem(k,abs(xk),.);grid on; subplot(212);plot(k,abs(xk));grid on; （2）FFT谱分析程序如下： clear clc fs=1000 t=0:1/fs:0.25; N=512; f1=50;f2=600; s=sin(2*pi*f1*t)+sin(2*pi*f2*t); x=s+randn(size(t)); Y=fft(x,N); PSD=Y.*conj(Y)/N; f=fs/N*(0:N/2-1); k=0:N/2-1; subplot(311);plot(x); subplot(312);plot(k,abs(Y(1:N/2))); subplot(313);plot(f,PSD(1:N/2)); 实验小结 通过数字信号处理课程设计，我自己学习很多知识。MATLAB类的编程书籍，各类信号处理类的书籍等，不仅丰富了自己的知识范围，而且对与自己所学的知识有了更深的了解和认识。 ? 在整个课程设计的过程中，我不仅学习很多知识，也帮助其他小组分析设计题目，不仅锻炼自己的学习能力，也加强了我与其他同学合作的能力。查找资料的过程中我也增强自己学习的能力。提高了我遇到问题分析问题解决问题的能力。这些都是我在以后的学习、生活和工作中的宝贵财富。