ms13r41.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 26 行

%方法3,用FFT求时域循环卷积

xn1=[5,4,-3,-2];      %建立x1(n)序列
xn2=[1,2,3,0];        %建立x2(n)序列
N=length(xn1);      
n=0:N-1;k=0:N-1; 
Xk1=fft(xn1,N);  %由x1(n)的DFT求X1(k) 
Xk2=fft(xn2,N);  %由x2(n)的DFT求X2(k)
Yk=Xk1.*Xk2;                       %Y(k)=X1(k)X2(k) 
yn=ifft(Yk,N);   %由Y(k)的IDFT求y(n)
yn=abs(yn)
figure(2)
subplot(2,3,1),stem(n,xn1);
title('x1(n)');
subplot(2,3,2),stem(n,xn2); 
title('x2(n)');
subplot(2,3,3),stem(n,yn); 
title('y(n)');
subplot(2,3,4),stem(k,abs(Xk1)); 
title('|X1(k)|');
subplot(2,3,5),stem(k,abs(Xk2)); 
title('|X2(k)|');
subplot(2,3,6),stem(k,abs(Yk)); 
title('|Y(k)|');