close all;
%% Calculo de los coeficientes del filtro FIR ideal por las series
%de Fourier

fc = 5000;
fm = 15000;

Fc = fc/fm;

% Ponemos un numero muy elevado de puntos
M = 201;
m = -(M-1)/2:(M-1)/2;

hi = 2*Fc*sinc(2*Fc*m);

figure, fig1=gcf, plot(m, hi,'-o'), hold;

% Truncamos hi con N puntos
N = 27;


% Valido para M y N par
% n = (M-N)/2+1:(M+N)/2;
% h = hi(n);   

% Valido para N impar
n = (M+1)/2-(N-1)/2:(M+1)/2+(N-1)/2;
h = hi(n); 


plot(n-(M+1)/2, h,'-r*');

[Hi,f] = freqz(hi,1,1000,fm);
% Hidb = 20*log10(abs(Hi));

[H,f] = freqz(h,1,1000,fm);
% Hdb = 20*log10(abs(H));

figure,fig2=gcf, plot(f,abs(Hi),'r',f,abs(H),'b');

F=[-500:500]/2000;

Wmag = sin(pi*F*N)./sin(pi*F);

figure, fig3=gcf, plot(F,Wmag);  %% Respuesta frecuancial de la ventana
                       %rectangular
		    
w = hamming(N);

h = h.*w';

figure(fig1), plot(n-(M+1)/2, h, '^m');