实验5+利用FFT实现快速卷积.docVIP

实验5 利用FFT实现快速卷积 专业班级： 信息111 学号： 201112030121姓名： 吕 治 报告日期： 2014.5.17 一、实验目的： 1、加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理。 2、掌握循环卷积和线性卷积两者之间的关系。 二、实验原理： 数字滤波器根据系统的单位脉冲响应是有限长或无限长可分为有限长单位脉冲响应（FTR系统）和无限长单位脉冲响应（IIR系统）。 对于FIR滤波器来说，除了可以使用数字网络实现以外，还可以使用FFT来实现。 一个信号序列通过FIR滤波器时，其输出应该是和的卷积。 当是一个有限长序列，即是FIR滤波器，且时 应用FFT实现数字滤波器实际上就是用FFT实现有限长序列的线性卷积。 这种方法就是先将输入信号通过FFT变换为它的频谱采样值，然后再和FIR滤波器的频响采样值相乘，可实现存放在存储器中，最后再将乘积经过快速傅里叶变换的反变换IFFT还原为时域序列，即得到输出序列。 下面以FFT求有限长序列间的卷积及求有限长序列与较长序列间的卷积为例来讨论FFT快速卷积算法。 （1）序列和的长度差不多，设的长度为，的长度为，要求 用FFT完成这一卷积的步骤如下： ①为使两有限长序列的线性卷积可用循环卷积替代而不发生混乱，必须保证，若采用基2-FFT完成卷积，则需要满足（r为整数） ②用补零的方法使和变成长度为N的序列 ③用FFT变换求出和的N点离散傅里叶变换。 ④完成和的乘积， ⑤利用FFT求出的反变换； （2）当长度很长时，即，通常不允许等全部采集齐之后再进行卷积，否则使输出相对于输入有较长的延迟，另外，若太大，则要补上太多的零，很不经济，且FFT要进行很长的运算，为此采用分段卷积的方法，本实验不作要求。 三、实验内容： 1、数字滤波器的脉冲响应，取8输入序列可选下列几种情况： ①，取8 ②，取8 ③，取8 2、实验前预先计算好的值 3、实验前，预先编制一个应用FFT实现数字滤波器实验的程序 编写一个应用FFT实现数字滤波器的function文件如下所示： function [ny,y]=getConvByFFT(nx,x,nh,h) n1=length(nx); n2=length(nh); N=n1+n2-1; n3=pow2(round(log2(N))); ny=0:n3-1; y1=zeros(1,n3); y2=y1; y1(1:n1)=x; y2(1:n2)=h; nx=ny; nh=ny; x=y1; h=y2; xk=fft(x); hk=fft(h); yk=xk.*hk; y=ifft(yk); stem(ny,y,fill),grid on 四、实验结果 五、实验总结 通过此次实验的练习，加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理，并且掌握了循环卷积和线性卷积两者之间的关系。