用fft作谱分析1.doc

文档介绍：西安电子科技大学课程论文 数字图像处理
快速抽取傅里叶变换
班级：070821
作者：董文凯
学号：07082001

西安电子科技大学课程论文 数字图像处理
快速抽取傅里叶变换
班级：070821
作者：董文凯
学号：07082001
时间：2011-06-30
用FFT作谱分析
实验目的
1、进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一
种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。
2、熟悉FFT算法原理和FFT子程序的应用。
3、学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出
现的分析误差及其原因，以便在实际中正确应用FFT。
实验步骤
1、复习DFT的定义、性质和用DFT作谱分析的有关内容。
2、复习FFT算法原理与编程思想，并对照DIT—FFT运算流图和程序框图，
读懂本实验提供的FFT子程序。
3、编制信号产生子程序，产生以下典型信号供谱分析：

，
，
， 其它n

，
，
， 其它n



应当注意，如果给出的是连续信号，则首先要根据其最高频率确定采样速率以及由频率分辨率选择采样点数N，然后对其进行软件采样(即计算，)，产生对应序列。对信号，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。
4、实验程序如下：
N=8;
n=0:N-1;
X1n=[1 1 1 1 0 0 0 0];
X2n=[1 2 3 4 4 3 2 1];
X3n=[4 3 2 1 1 2 3 4];
X4n=cos(pi*n/4);
X5n=sin(pi*n/8);
X1kn=fft(X1n,N);
X2kn=fft(X2n,N);
X3kn=fft(X3n,N);
X4kn=fft(X4n,N);
X5kn=fft(X5n,N);
M=16;
m=0:M-1;
X1m=[1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0];
X2m=[1 2 3 4 4 3 2 1 0 0 0 0 0 0 0 0];
X3m=[4 3 2 1 1 2 3 4 0 0 0 0 0 0 0 0];
X4m=cos(pi*m/4);
X5m=si