conv_dft.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

% Name:conv_dft
% Linear convolution Via the DFT
x=input ('Type in the first sequence = ');
h=input ('Type in the second sequenc = ');
% Determine the length of the result of convolution
L=length(x)+length(h)-1;
XE=fft(x,L);
HE=fft(h,L);
% Determine the IDFT of the product
y1=ifft(XE.*HE);
% Plot the sequence generated by DFT-based convolution
% and the error from direct linear convolution
k=0:L-1;
subplot(211),stem(k,y1),xlabel('Time index n'),ylabel('Amplitude')
title('Result of DFT-based linear convolution')
y2=conv(x,h);error=y1-y2;
subplot(212),stem(k,abs(error)),xlabel('Time index n'),ylabel('Amplitude')
title('Magnitude of error sequence')