teste_som

//Bloco de ciração do filtro
n   = 5; // Ordem do filtro
Fs      = 48000; // Frequência de amostragem
Fc = 1000;   // Frequência de corte
Freq = Fc/Fs;
hz = iir(n,'lp','butt',[Freq 0],[0 0]);//Butterworth - transformação bilinear

[y1,Fs,bits]=wavread("/home/edmau/Downloads/audiocheck.net_sin_1500Hz_-3dBFS_5s.wav");

//H = poly(hz,'x','coeff');
H = horner(hz,1/%z);

numerador = round(2^20.*coeff(numer(H)));
denominador = round(2^20.*coeff(denom(H)));


[audio1,audio2]=filter(numerador,denominador,y1);

sound(audio1,Fs,bits);