《数字信号处理实验六.docVIP

实验六?利用FFT实现快速卷积 一．实验目的 1.加深理解FFT在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT进行数字信号处理。 2.掌握圆周卷积和线性卷积两者之间的关系。 二．实验内容及结果 ⑴给定两个序列x(n)=[2,1,1,2],h(n)=[1,-1,-1,1]。首先直接在时域计算两者的线性卷积；其次用FFT快速计算二者的线性卷积，验证结果。 ⑵数字滤波器的脉冲响应为 ,N可自定，本实验取N=17 输入序列x(n)可选下列几种情况： ①x(n)=R(n),N可自取16 ②x(n)=cos(,N=16 ③x(n)=( R(n), N=16 ⑶实验前，预先编制一个应用FFT实现数字滤波器的通用程序。 ⑷上机独立调试，并打印或记录实验结果。 ⑸将实验结果与预先笔算的结果比较，验证其正确性。 解: (1)程序： x=[2,1,1,2]; h=[1,-1,-1,1]; y=conv(x,h); subplot(2,1,1); stem(y); title(时域上线性卷积); n=[0:1:3]; m=[0:1:3]; N1=length(n); N2=length(m); N=N1+N2-1; X=fft(x,N); H=fft(h,N); Y2=X.*H; y2=ifft(Y2,N); if all(imag(x)==0)(all(imag(h)==0)) y=real(y); end subplot(2,1,2); stem(y2); title(FFT快速卷积); (2) ①x(n)=R (n),N 可自取16 程序： clc clear all n=[0:1:15]; m=[0:1:16]; N1=length(n); N2=length(m); x=ones(1,N1); h=(-1/2).^m; N=N1+N2-1; X=fft(x,N); H=fft(h,N); Y=X.*H; y=ifft(Y,N); if all(imag(x)==0)(all(imag(h)==0)) y=real(y); end x1=0:N-1; stem(x1,y); title(FFT快速卷积); ②x(n)=cos(,N=16 程序： clc clear all n=[0:1:15]; m=[0:1:16]; N1=length(n); N2=length(m); x=cos(2*pi*n/16); h=(-1/2).^m; N=N1+N2-1; X=fft(x,N); H=fft(h,N); Y=X.*H; y=ifft(Y,N); if all(imag(x)==0)(all(imag(h)==0)) y=real(y); end x1=0:N-1; stem(x1,y); title(FFT快速卷积); ③x(n)=( R(n), N=16 程序： clc clear all n=[0:1:15]; m=[0:1:16]; N1=length(n); N2=length(m); x=(1/3).^n; ; h=(-1/2).^m; N=N1+N2-1; X=fft(x,N); H=fft(h,N); Y=X.*H; y=ifft(Y,N); if all(imag(x)==0)(all(imag(h)==0)) y=real(y); end x1=0:N-1; stem(x1,y); title(FFT快速卷积); 三．心得体会 通过本次实验，我加深理解了FFT在实现数字滤波（或快速卷积）中的重要作用、掌握了圆周卷积和线性卷积两者之间的关系，受益匪浅。