实验二 应用FFT对信号进行频谱分析.docxVIP

实验二 使用FFT对信号进行频谱分析 姓 名: 刘竟伦 班 级: 11电子A 学 号: 1115105037 日 期 2014 年 3月15日 指导老师 戴在平 华侨大学信息科学和工程学院电子工程系 一、实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 (2) 熟悉FFT算法原理和FFT子程序的使用。 （3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法， 了解可能出现的分析误差及其原因， 以便在实际中正确使用FFT。 二、实验原理和方法 (1) 复习DFT的定义、 性质和用DFT作谱分析的有关内容。 (2) 复习FFT算法原理和编程思想， 并对照DIT-FFT运算流图和程序框图， 读懂本实验提供的FFT子程序。 (3) 编制信号产生子程序， 产生以下典型信号供谱分析用： 应当注意，如果给出的是连续信号xa（t），则首先要根据其最高频率确定采样速率fs以及由频率选择采样点数N，然后对其进行软件采样（即计算择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行分析，否则有可能产生较大的分析误差。请实验者根据DFT的隐含周期性思考这个问题。 (4) 编写主程序。 三、上机实验内容及实验结果图 （1）对2中所给出的信号逐个进行谱分析。下面给出针对各信号的FFT变换区间N以及对连续信号x6（t）的采样频率fs。 x1(n),x2(n),x3(n),x4(n),x5(n),:N=8,16 x6(t):fs=64(hz),N=16,32,64 x1(n)=R4(n) 对应x1(n)的matlab源程序代码： n=[0:7]; x=[1 1 1 1 0 0 0 0] f1=fft(x,8) f2=fft(x,16) subplot(2,2,1) stem(n,x); axis([0 8 0 2]) xlabel(n) ylabel(xl(n)) title(xl的波形 ) subplot(2,2,4) k=[0:15] stem(k,abs(f2)); axis([0 16 0 5]) xlabel(k) ylabel(|xl(k)|) title(xl(n)的16点fft) subplot(2,2,3) k=[0:7] stem(k,abs(f1)); axis([0 10 0 5]) xlabel(k) ylabel(|xl(k)|) title(xl(n)的8点fft) 得到如下波形图： 2. 对应x2(n)的matlab源程序代码： n=[0:7]; x=[1 2 3 4 4 3 2 1] f1=fft(x,8) f2=fft(x,16) subplot(2,2,1) stem(n,x); axis([0 8 0 4]) xlabel(n) ylabel(x2(n)) title(x2的波形) subplot(2,2,4) k=[0:15] stem(k,abs(f2)); axis([0 16 0 20]) xlabel(k) ylabel(|x2(k)|) title(x2(n)的16点fft) subplot(2,2,3) k=[0:7] stem(k,abs(f1)); axis([0 10 0 20]) xlabel(k) ylabel(|x2(k)|) title(x2(n)的8点fft) 3. 对应x3(n)的matlab源程序代码： n=[0:7]; x=[4 3 2 1 1 2 3 4] f1=fft(x,8) f2=fft(x,16) subplot(2,2,1) stem(n,x); axis([0 8 0 4]) xlabel(n) ylabel(x3(n)) title(x3的波形) subplot(2,2,4) k=[0:15] stem(k,abs(f2)); axis([0 16 0 20]) xlabel(k) ylabel(|x3(k)|) title(x3的16点fft) subplot(2,2,3) k=[0:7] stem(k,abs(f1)); axis([0 8 0 20]) xlabel(k) ylabel(|x3(k)|) title(x3(n)的8点fft) 4.x4（n）=cos