%Método de las series de Fourier
%Filtro pasoalto%Especificaciones
fp=5000;
fs=4000;
fm=14000;
Ap=1;
As=45;

%Se calculan las frecuencias digitales
FpHP=fp/fm;
FsHP=fs/fm;

%Se calcula las especificaciones del filtro pasabajo
Fp=FsHP;
Fs=FpHP;

%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));			%N es el número de puntos (el orden es N-1)

%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
hLP=hLP.*hamming(N)';

%Se calcula el filtro pasaalto
delta=zeros(1,N);
delta((N)/2)=1;
hHP=delta-hLP;

%Calculamos el filtro pasoalto
[H,F]=freqz(hHP,1,1000,fm);

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