function fir_freq_demo

N = 16;
H = zeros(1,N);
H(1+ 3) = 1;
H(1+ N-3) = 1;

% h(n) = IDFT(H)
h = ifft(H);
stem(0:N-1,h)

% segnale x(n) casuale
Nx = 100;
%x = randn(1,100);
x = cos(2*pi*(0:Nx-1)*3/16);
y = fir_freq(H,x);
%y = conv(h,x);

subplot(2,1,1)
stem(0:Nx-1,x)
set(gca,'YLim',[-3,3])
subplot(2,1,2)
stem(0:Nx+N-2,imag(y))
set(gca,'YLim',[-3,3])
