% Archivo: falsopositivo_filt.m
%
% Script que realiza el análisis de falso positivo, se diferencia del
% anterior en la inclusión de la orden filter, que filtra el audio
% antes de calcular la huella

clear all;

%Proceso completo del análisis de falso positivo

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

% cargamos el bloque de 3 segundos y obtenemos su huella

audio1=readwav('pictures5512_3sg.wav');
audio1filt=filter([0.99 -0.99],[1 -0.98],audio1);
f1=fingerprint_princ(audio1filt);
[m,n]=size(f1);

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

audio2=readwav('1123075521_20min.wav');
audio2filt=filter([0.99 -0.99],[1 -0.98],audio2);
f2=fingerprint_princ(audio2filt)

%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)

dfitttool