试验FFT算法的应用.DOC

实验4 FFT算法的应用 一 、实验目的： 加深对离散信号的DFT的理解及其FFT算法的运用。 二 、实验原理： N点序列的DFT和IDFT变换定义式如下： ， 利用旋转因子具有周期性，可以得到快速算法（FFT）。 在MATLAB中，可以用函数X=fft（x，N）和x=ifft（X，N）计算N点序列的DFT正、反变换。 三、实验内容 1. 对连续的单一频率周期信号 按采样频率采样，截取长度N分别选N =20和N =16，观察其DFT结果的幅度谱。 解 此时离散序列： 即N=8。 用MATLAB计算并作图，函数fft用于计算离散傅里叶变换DFT，程序如下： q=8; k1=[0:1:19]; xa1=sin(2*pi*k1/q); subplot(2,2,1) plot(k1,xa1) xlabel(k);ylabel(x1(k)); xm1=fft(xa1);xm1=abs(xm1); subplot(2,2,2) stem(k1,xm1) xlabel(m);ylabel(X1(m)); k2=[0:1:15]; xa2=sin(2*pi*k2/q); subplot(2,2,3) plot(k2,xa2) xlabel(k);ylabel(x2(k)); xm2=fft(xa2);xm2=abs(xm2); subplot(2,2,4) stem(k2,xm2) xlabel(m);ylabel(X2(m)); ?????? ??图2.1 计算结果示于图2.1，(c) 和(d) 分别是N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱；而(a) 和(b)分别是N=20时的截取信号和DFT结果，由于截取了两个半周期，不再是单一谱线的频谱，而是含有丰富的谐波分量，与N=16时的差别很大；上述频谱的误差主要是由于有限长序列的离散傅里叶变换是把有限长序列作为周期序列的一个周期来表示的，隐含有周期性意义。 2. 分别计算16点序列 的16点和32点DFT,绘出幅度谱图形。 (a) (b) (c) (d)