% Archivo: robustez.m
%
% Script para calcular la BER entre dos trozos de audio
%

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


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

audio1=readwav('jimmy_3sg_5512.wav');
f1=fingerprint_princ(audio1);
[m1,n1]=size(f1);

% cargamos el bloque modificado con el que vamos a comparar

audio2=readwav('jimmy_3sg_comp_5512.wav');
f2=fingerprint_princ(audio2);
[m2,n2]=size(f2);

% para tener en cuenta el hecho de que por la distorsión
% puede que los dos trozos sean de distinta duración

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