clear;
j = sqrt(-1);
H=zeros(1,16);
H(1) = 1;
H(2) = exp(-j*pi*15/16);
H(3) = 0.7*exp(-j*pi*30/16);
% H(4) = exp(-j*pi*45/16);
% H(5) = exp(-j*pi*60/16);
% H(13) = conj(H(5));
% H(14) = conj(H(4));
H(15) = conj(H(3));
H(16) = conj(H(2));

h = real(ifft(H));

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

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

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

% Usando la funcion fir2

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

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