数字信号处理上机实验

(5)、用FFT分别实现xa(n)（p＝8，q＝2）和 积和线性卷积。

xb(n)（a＝0.1，f＝0.0625）的16点圆周卷

当N=16时，圆周卷积不等于线性卷积。要用FFT计算圆周卷积需在序列后面补零使其长度L=16+16-1=31。

M=16; N=16; nx=0:15; nh=0:15;

xn=exp(-(nx-8).*(nx-8)/2);

hn=exp(-0.1*nh).*sin(2*pi*0.0625*nh); L=pow2(nextpow2(M+N-1)); XK=fft(xn,L); HK=fft(hn,L); YK=XK.*HK; yn=ifft(YK,L); ny=0:L-1;

subplot(3,1,1);stem(nx,xn,'.');title('x(n)'); subplot(3,1,2);stem(nh,hn,'.');title('h(n)'); subplot(3,1,3);stem(ny,abs(yn),'.');title('y(n)');

图5 xa(n)、xb(n)及其线性卷积波形