% Archivo: robustez_filt.m
%
% Script para calcular la BER entre dos trozos de audio, tambiéns e
% diferencia del anterior en la inclusión de la orden filter
%

clear all;
Fs=44100/8;        %frecuencia de muestreo


% cargamos el bloque de 3 segundos original y obtenemos su huella

audio1=readwav('vangelis_3sg_5512.wav');
audio1filt=filter([0.99; -0.99],[1; -0.98],audio1);
f1=fingerprint_princ(audio1filt);
[m1,n1]=size(f1);

% cargamos el bloque modificado con el que vamos a comparar

audio2=readwav('vangelis_3sg_lsc-1_5512.wav');
audio2filt=filter([0.99; -0.99],[1; -0.98],audio2);
f2=fingerprint_princ(audio2filt);
[m2,n2]=size(f2);

if m1<m2
    d=busqueda(f1,f2);
    d=d./(m1*n1);
    BER=min(d)

else if m1==m2
    BER=sum(sum(abs(f2-f1)))/(m1*n1)
    
    else
        d=busqueda(f2,f1);
        d=d./(m2*n2);
        BER=min(d)
    end
end
