%Método de las series de Fourier
%Filtro pasobajo
%Especificaciones
fp=1500;
fs=3000;
fm=8000;
Ap=1;
As=45;

%Se calculan las frecuencias digitales
Fp=fp/fm;
Fs=fs/fm;

%Cálculo del orden del filtro

Ws=1.91;								%Se obtiene de la tabla 3 de los apuntes
									%En este caso vamos a utilizar la ventana de Hamming
N=ceil(Ws/(Fs-Fp));

%Frecuencia de corte
Fc=1.1*Fp;						

%Filtro pasobajo

n=[-(N-1)/2:(N-1)/2];
hLP=2*Fc*sinc(2*n*Fc);

%Se aplica la ventana de Hamming
h=hLP.*hamming(N)';
[H,F]=freqz(h,1,1000,fm);

plot(F,20*log10(abs(H)));hold on;
plot([fp fs],[-Ap -As],'*r');