数字信号处理 实验6 利用FFT实现快速卷积

《数字信号处理 实验6 利用FFT实现快速卷积》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、x=[2,1,1,2];h=[1-1-11];N=7;XK=fft(x,N);HK=fft(h,N);YK=XK.*HK;yn=ifft(YK,N);ifall(imag(x)==0)&(all(imag(h)==0))yn=real(yn);endy=conv(x,h);y1=y-yn;n=0:N-1;subplot(3,1,1);stem(n,y);ylabel('时域计算');subplot(3,1,2)stem(n,yn);ylabel('FFT快速');subplot(3,1,3);stem(n,y1)N=N

2、1+N2-1;XK=fft(xn,N);HK=fft(hn,N);YK=XK.*HK;yn=ifft(YK,N);ifall(imag(xn)==0)&(all(imag(hn)==0))yn=real(yn);endx=0:N-1;stem(x,yn);n=[0:1:15];m=[0:1:16];N1=length(n);N2=length(m);xn=cos(2*pi*n/N1);hn=(-0.5).^m;N=N1+N2-1;XK=fft(xn,N);HK=fft(hn,N);YK=XK.*HK;yn=ifft(Y

3、K,N);ifall(imag(xn)==0)&(all(imag(hn)==0))yn=real(yn);endx=0:N-1;stem(x,yn);n=[0:1:15];m=[0:1:16];N1=length(n);N2=length(m);xn=(1/3).^n;hn=(-0.5).^m;N=N1+N2-1;XK=fft(xn,N);HK=fft(hn,N);YK=XK.*HK;yn=ifft(YK,N);ifall(imag(xn)==0)&(all(imag(hn)==0))yn=real(yn);endx

4、=0:N-1;stem(x,yn);