Matlab code for digital filter example with Handel audio.

digitalFilter.m

load handel.mat
y = y';
Wp = 1000/Fs; %%Frecuencia de corte normalizada
Ws = 2000/Fs; %%Frecuencia de atenuación normalizada
[n,Wn] = buttord(Wp,Ws,3,60); %%orden de butterworth
[b,a] = butter(n,Wn,'low');   %%coeficientes de butterworth 
yh = filter(b, a, y);       %%Filtro de la señal
L=length(y);
NFFT = 2^nextpow2(L); % Siguiente potencia de 2 del vector y
f = Fs/2*linspace(0,1,NFFT/2+1);

figure(1)
yf=fft(y,NFFT)/L;
plot(f,2*abs(yf(1:NFFT/2+1)))
title('Espectro de la señal')
xlabel('Frecuencia [Hz]')
grid on

figure(2)
yf2=fft(yh,NFFT)/L;
plot(f,2*abs(yf2(1:NFFT/2+1)))
title('Espectro de la señal filtrada')
xlabel('Frecuencia [Hz]')
grid on

figure(3)
[H,W] = freqz(b,a,1000);
plot(Fs.*W./pi,abs(H))
title('Respuesta en frecuencia del filtro')
xlabel('Frecuencia [Hz]')
grid on
sound(yh)

sound (y)