实验6-FFT算法的应用讲课稿.pdfVIP

实 验 6 - FFT 算 法 的 应 用 精品文档 实验 6 FFT 算法的应用 实验目的： 加深对离散信号的 DFT 的理解及其 FFT 算法的运用。 实验原理： N 点序列的 DFT 和 IDFT 变换定义式如下： N 1 kn 1 N 1 kn X [ k ] x [ n ]W N ， x[ n] X [ k ]W N n 0 N k 0 2 j nk kn N 利用旋转因子 WN e 具有周期性，可以得到快速算法（ FFT）。 在 MATLAB 中，可以用函数 X=fft （x ，N ）和 x=ifft （X ，N ）计算 N 点序列 的 DFT 正、反变换。 例 1 对连续的单一频率周期信号 按采样频率 采样，截取长度 N分别 N N 选 =20 和 =16 ，观察其 DFT结果的幅度谱。 解 此时离散序列 ，即 k=8。用 MATLAB 计算并作图，函数 fft 用于计算离散傅里叶变换 DFT，程序如下： k=8; n1=[0:1:19]; xa1=sin(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) xlabel(t/T);ylabel(x(n)); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) xlabel(k);ylabel(X(k)); n2=[0:1:15]; xa2=sin(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) xlabel(t/T);ylabel(x(n)); xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4) stem(n2,xk2) 收集于网络，如有侵权请联系管理员删除 精品文档 xlabel(k);ylabel(X(k)); 计算结果示于图 2.1 ，(a) 和(b) 分别是 N=20 时的截取信号和 DFT结果，由 于截取了两个半周期，频谱出现泄漏； (c) 和 (d) 分别是 N=16 时的截取信号和 DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主 要是由于时域中对信号的非整周期截断产生的频谱泄漏。 实验内容： （1） 2N 点实数序列 2 1 2 cos( 7n)