%FIRGIBBS	Filtro passa-baixas FIR e efeito (fenômeno) de Gibbs
%
%Miguel Arjona Ramírez 17/3/2014

%  L = 10;
%  L = 100;
L = 1000;
n = (0:L-1)';
M = (L-1)/2;
w0 = pi/2;
Nf = 4096;
Nm = Nf/2+1;
f = (0:Nm-1)'/Nf;
h2 = sin((n-M)*w0)./(pi*(n-M));
%  h4 = (1-cos((n-M)*w0))./(pi*(n-M)); %não é passa-baixas
H2 = fft(h2,Nf);
M2 = abs(H2(1:Nm));
%  H4 = fft(h4,Nf);
%  M4 = abs(H4(1:Nm));
Wm = sin(pi*f*L)./sin(pi*f);	%magnitude da TFTD da janela retangular
%
figure(1)
plot(f,M2)
grid
%  figure(2)
%  plot(f,M4)
figure(2)
plot(f,Wm)
grid
%
figure(3)
stem(n,h2)
grid
%
figure(4)
plot(f,20*log10(M2))
grid
%
%  wh = hamming(L);
%  h2w = h2.*wh;
%  H2w = fft(h2w,Nf);
%  M2w = abs(H2w(1:Nm));
%  Wh = fft(wh,Nf);
%  MWh = abs(Wh(1:Nm));
%  %
%  figure(5)
%  plot(f,20*log10(M2w))
%  grid
%  %
%  figure(6)
%  plot(f,20*log10(MWh))
%  grid