实验二应用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,16x6(t):fs=64(hz),N=16,32,64x1(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(pi/4*n)对应x4(n)的matlab源程序代码：n=[0:7];x=cos(0.25*pi*n)f1=fft(x,8)subplot(2,2,1)stem(n,x);axis([0 8 -4 4])xlabel(n)ylabel(x4(n))title(x4的波形)n=[0:15]x=cos(0.25*pi*n)f2=fft(x,16)subplot(2,2,2)stem(n,x);axis([0 16 -4 4])ylabel(x4(n))title(x4的波形)subplot(2,2,4)k=[0:15]stem(k,abs(f