function c = conv_by_fft( a, b )

% dimensioni delle sequenze
N = numel(a);
M = numel(b);
L = N+M-1;

N_FFT = 2^ceil(log2(L));

% zero padding
a(N_FFT)=0; 
b(N_FFT)=0;

% convoluzione con la dft
temp = ifft(fft(a).*fft(b));

% e rimozione degli zeri in eccesso
c = temp(1:L);
