% Archivo: falsopositivo_dda.m
%
% Script para el análisis de falso positivo. Se diferencia de
% los anteriores en que se llama a la función preprocdda, que es la que
% hace el preprocesado en sí mismo y en que se trabaja a 44100 Hz

clear all;

%Proceso completo del análisis de falso positivo

Fs=44100;        %frecuencia de muestreo
alpha=0.35;        %umbral

% cargamos el bloque de 3 segundos y obtenemos su huella

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

audio1=decimate(preproc1,2); 
f1=fingerprint_princ(audio1);
[m,n]=size(f1);

% cargamos el archivo con el que vamos a comparar, entero

audio21=readwav('1123075521_5512_20min.wav');

preproc2=preprocdda2(audio21.');

audio2=decimate(preproc2,2);
f2=fingerprint_princ(audio2);

%se comparan ambas huellas

d=busqueda(f1,f2);
d=d./(m*n);             %BER de cada comparación

% sacamos la media y la desviacion estandar

media=mean(d)
desv=std(d)



