试验6FFT算法的应用.docVIP

实验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,xm