实验3 FFT算法的应用.docVIP

实验3 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结果的幅度谱。 解 此时离散序列：即k=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，(a) 和(b)分别是N=20时的截取信号和DFT结果，由于截取了两个半周期，频谱出现泄漏；(c) 和(d) 分别是N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主要是由于时域中对信号的非整周期截断产生的频谱泄漏。 2.教材P130 2-12 利用MATLAB实现由DFT计算有限序列线性卷积，并与线性卷积直接计算的结果进行比较，其中， 程序如下： x=[1,2,0,1];h=[2,2,1,1]; l=length(x)+length(h)-1; XE=fft(x,l);HE=fft(h,l); y1=ifft(XE.*HE); k=0:l-1; subplot(2,1,1); stem(k,real(y1)); axis([0 6 0 7]); title(线性卷积结果); xlabel(k);ylabel(A); y2=conv(x,h); error=y1-y2; subplot(2,1,2); stem(k,abs(error)); xlabel(k);ylabel(A); title(计算误差); (d) (c) (b) (a) 1