基于Matla DFT 及FFT 频谱分析_ 刘小群.docVIP

. . 　　 本文档下载自 HYPERLINK 文档之家，如果内容可能不完整，您可以点击以下网址继续阅读或下载：  HYPERLINK /doc/0a0732b9d15abe23482f4dec.html /doc/0a0732b9d15abe23482f4dec.html 基于Matlab的DFT及FFT频谱分析_刘小群 山西电子技术2010年第4期 文章编号:1674-4578(2010)04-0048-02 软件技术 基于Matlab的DFT及FFT频谱分析 刘小群,周云波 1 2 (宝鸡文理学院,陕西宝鸡721007) 摘 要:DFT及FFT是数字信号处理的重要内容。DFT是FFT的基础,FFT是DFT的快速算法,在MATLAB中可以利用函数FFT来计算序列的离散傅里叶变换DFT。基于此首先介绍了Matlab软件;其次给出了基于Matlab软件的DFT和FFT频谱分析的方法,利用Matlab软件方法,使得设计方便、快捷,大大减轻了工作量;最后结合实例给出了仿真结果。 关键词:Matlab;DFT;FFT;频谱分析 中图分类号:TP311.52:TN929.5 文献标识码:A n1=[0:1:19]; xa1=sin(2*p*in1/k);subplot(2,2,1)plot(n1,xa1) xlabel(ct/Tc);ylabel(cx(n)c);xk1=fft(xa1);xk1=abs(xk1);subplot(2,2,2)stem(n1,xk1) xlabel(ckc);ylabel(cX(k)c);n2=[0:1:15]; xa2=sin(2*p*in2/k);subplot(2,2,3)plot(n2,xa2) xlabel(ct/Tc);ylabel(cx(n)c);xk2=fft(xa2);xk2=abs(xk2);subplot(2,2,4)stem(n2,xk2) xlabel(ckc);ylabel(cX(k)c); 计算结果示于图1,上方两个图分别是N=20时的截取信号和DFT结果,由于截取了两个半周期,频谱出现泄漏;下方两个图分别是N=16时的截取信号和DFT结果,由于截取了两个整周期,得到单一谱线的频谱。上述频谱的误差主要是由于时域中对信号的非整周期截断产生的频谱泄漏。 0 引言 随着信息时代和数字世界的到来,数字信号处理己成为当今一门极其重要的学科和技术领域,数字信号处理在通信、语音、图像、自动控制、医疗和家用电器等众多领域得到了广泛的应用。任意一个信号都具有时域与频域特性,信号的频谱完全代表了信号,因而研究信号的频谱就等于研究信号本身。通常从频域角度对信号进行分析与处理,容易对信号的特性获得深入的了解。因此,信号的频谱分析是数字信号处理技术中的一种较为重要的工具[1]。 1 Matlab简介 Matlab语言是当今科学界最具影响力、也是最有活力的软件。它起源于矩阵运算,并已经发展成一种高度集成的计算机语言。它是由美国MathWorks公司于1982年推出的软件产品,取名来源于MatrixLaboratory,简称为/Matlab0。Matlab是一个完整的、可扩展的、高性能数值计算的可视化软件,是一种进行科学工程计算的交互式程序设计语言。经过不断的扩充和完善,它已成为各类科学研究与工程应用中的标准工具,它集数值分析,矩阵运算,信号处理与图形显示于一体,构成一个界面友好的用户环境。它的信号处理工具箱包含了各种经典的和现代数字信号处理技术。是一个非常优秀的算法研究与辅助设计的工具 [2] 。 2 序列的DFT变换 DFT是FFT的基础,FFT是DFT的快速算法,在MAT-LAB中可以利用函数fft来计算序列的离散傅里叶变换DFT。 例如:对连续的单一频率周期信号按采样频率的幅度谱。此时离散序列 k=8; 3 FFT频谱分析 在实际应用中,FFT是最常见的数字信号处理算法,它在各种数字信号处理系统中扮演重要的角色。在信号处理过程中,频域分析往往比时域分析方便和高效,FFT是时域和频域转换的基本运算。例如,通过计算信号序列的FFT可以直接分析它的数字频谱;在计算数字滤波器的输出响应时需要进行卷积运算,它可以通过先进行FFT运算而后将信号频谱与系统频响相乘再进行IFFT变换得到:匹配滤波 fs=8fa 采样,截取长度N,分别选N=20和N=16,观察其DFT结果 x(n)=sin(2Pnfa/fs)= sin(2Pn/8),即k=8。程序如下: 收稿日期:2010-05-28 作者简介:刘小群(1977-),女,陕西兴平人,讲师,硕士研究生,研究方向为通信与信息系统。 /doc/0a0732b9d15abe23482f4dec.html第4期 刘小群,等:基于Matla