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

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

1、•(N4uq)yj=NH：(N'ux)出二齐x-1NN+1N二N(

2、A4U)UJ9)S：(£T£)】o【dqns•C®¥!ldd.)PQBIx：(uAu)iuo)s(乙T£)2[dqns乂蘆44海闭)弊1引人t(X6u)uiojs：(lT£))o[dqns•I-N-O=ut(q'x)AUOO=Xpusi(uX)[B9J=uX((o二二(q)弘UI!)Z)茅(o二二(x)託UI!)[[EJ!：(曲人)出!=吠•NH*NX=NA•(N4U)W=NH•(N^)yj=NXAN•[II-I-I]=M•irrrrMYK二XK.*HK;yn=ifft(Y

3、K,N);ifall(imag(xn)==O)&(all(imag(hn)==O))yn=real(yn);endx=O:N-l;stem(x,yn);n二[0:1:15];m=[0:l:16];Nl=length(n);N2=length(m);xn=cos(2*pi*n/N1);hn=(-0.5).Am;N=N1+N2-1;XK=ffi(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-

4、l;stem(x,yn);1(?0.80.60.40.2-0.2-0.4-0.6竿egoG-e-oee-e-e-e-e101520253035n=[0:l:15];m=[0:l:16];Nl=length(n);N2=length(m);xn=(l/3).An;hn=(-0.5).Am;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=0:N-l;stem(x,yn

5、);1.20.80.6-0.4-0.2oooooooooooe-e-e-e-e-e-e-e-e-e-e-°-205101520253035