% Archivo: robustez_dda.m
%
% Script para calcular la BER entre dos trozos de audio con preprocesado
% DDA. Se diferencia de los anteriores en que se cargan los archivos a
% 44100 Hz
%

clear all;


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

audio11=readwav('vangelis44100_3sg.wav');
preproc1=preprocdda2(audio11);

audio1=decimate(preproc1,2);     %se pasa de 11025 Hz a 5512 Hz

f1=fingerprint_princ(audio1);
[m1,n1]=size(f1);

% cargamos el bloque modificado con el que vamos a comparar

audio21=readwav('vangelis44100_3sg_noise.wav');
preproc2=preprocdda2(audio21);

audio2=decimate(preproc2,2);     %se pasa de 11025 Hz a 5512 Hz

f2=fingerprint_princ(audio2);
[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