clear;
j = sqrt(-1);
H=zeros(1,17);
H(1) = 1;
H(2) = exp(-j*pi*16/17);
H(3) = exp(-j*pi*32/17);
H(4) = exp(-j*pi*48/17);
H(5) = exp(-j*pi*64/17);
H(14) = conj(H(5));
H(15) = conj(H(4));
H(16) = conj(H(3));
H(17) = conj(H(2));

h = real(ifft(H));

[Hf,Wf] = freqz(h,1,500); 

Ff = Wf/(2*pi);
F=[0:8]/17;
plot(Ff,abs(Hf),F,abs(H(1:9)),'o');zoom;grid;
hold;

h = h.*hamming(17)';
Hf1 = freqz(h,1,500); 
plot(Ff,abs(Hf1),'r');

% Usando la funcion fir2

B = fir2(16,[0 1 2 3 4 5 6 7 8 8.5]*2/17,[1 1 1 1 1 0 0 0 0 0],boxcar(17));
Hf2 = freqz(B,1,500);

plot(Ff,abs(Hf2),'m');