资源描述：
《利用FFT计算线性卷积.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字信号处理实验三实验三利用FFT计算线性卷积实验记录1.比较当序列长度分别为8,16,32,64,256,512,1024时两种方法计算的时间一、直接线性卷积程序代码forn1=3:10;L=2^n1;x=ones(1,L);h=cos(0.2*pi.*x);ticy=conv(x,h);tocend程序结果：Elapsedtimeis0.000104seconds.Elapsedtimeis0.000022seconds.Elapsedtimeis0.000018seconds.Elapsedtimeis0.0000
2、21seconds.Elapsedtimeis0.000031seconds.Elapsedtimeis0.000059seconds.Elapsedtimeis0.000152seconds.Elapsedtimeis0.000508seconds.二、快速卷积程序代码forn0=3:10;L=2^n0;n=0:L;x=heaviside(n)-heaviside(n-L);h=cos(0.2*pi.*n);X=fft(x);H=fft(h);Y=X.*H;ticy=ifft(Y);tocend数字信号处理实验三程序结
3、果Elapsedtimeis0.000041seconds.Elapsedtimeis0.000441seconds.Elapsedtimeis0.000075seconds.Elapsedtimeis0.000070seconds.Elapsedtimeis0.001507seconds.Elapsedtimeis0.000544seconds.Elapsedtimeis0.005577seconds.Elapsedtimeis0.002435seconds.数字信号处理实验三2.考察当L=2048和4096，M=25
4、6时两种方法计算时间1）N=2048时一、直接线性卷积程序代码L=2048;M=256;n0=0:L;x=heaviside(n0)-heaviside(n0-L);n1=0:M;h=cos(0.2*pi.*n1);ticy=conv(x,h);toc程序结果Elapsedtimeis0.011091seconds二、快速卷积程序代码L=2048;M=256;n0=0:L;x=heaviside(n0)-heaviside(n0-L);n1=0:M;h=cos(0.2*pi.*n1);ticX=fft(x,2048);H
5、=fft(h,2048);Y=X.*H;y=ifft(Y);toc程序结果Elapsedtimeis0.005779seconds.数字信号处理实验三2）N=4096时一、直接线性卷积程序代码L=4096;M=256;n0=0:L;x=heaviside(n0)-heaviside(n0-L);n1=0:M;h=cos(0.2*pi.*n1);ticy=conv(x,h);toc程序结果Elapsedtimeis0.002153seconds.二、快速卷积程序代码L=4096;M=256;n0=0:L;x=heavisi
6、de(n0)-heaviside(n0-L);n1=0:M;h=cos(0.2*pi.*n1);ticX=fft(x,2048);H=fft(h,2048);Y=X.*H;y=ifft(Y);toc程序结果Elapsedtimeis0.005149seconds.数字信号处理实验三3.利用重叠相加法计算线性卷积，考察L=2048和4096时的计算时间1）L=2048程序代码L=2048;M=256;x=ones(1,L);n=0:M-1;h=cos(0.2*pi*n);N=L+M-1;y=zeros(1,N);subLe
7、n=M;subNum=L/M;subLenY=2*M-1;ticH=fft(h,2*M-1);forcount=1:subNumx_temp=x(((count-1)*subLen+1):(count*subLen));X=fft(x_temp,2*M-1);Y=X.*H;y_temp=ifft(Y,2*M-1);y(((count-1)*subLen+1):(count*subLen+M-1))=y(((count-1)*subLen+1):(count*subLen+M-1))+y_temp;endtoc程序结果El
8、apsedtimeis0.003825seconds.数字信号处理实验三2）L=4096程序代码L=4096;M=256;x=ones(1,L);n=0:M-1;h=cos(0.2*pi*n);N=L+M-1;y=zeros(1,N);subLen=M;subNum=L/M;subLenY=2*M-1;ticH=ff
