% Analisis de la probabilidad de error en deteccion escogiendo los valores
% mas optimos de cada uno de los parametros que intervienen en el calculo
% de los fingerprints.

% En esta prueba tomamos con punto de partida el valor optimo del numero de bandas espectrales y vamos variando, primeramente, el numero de muestras por frame y
% luego el tiempo de captura

%Fecha de creacion: 1 de Abril de 2007


% N variable
% time_detect variable

num_band=9;     % valor optimo
over_lap=31/32;
finf=300;
fsup=2000;
alpha=0.18;                  % umbral para considerar deteccion correcta


fi=fopen('p_global_numband_N_timedetect','w');     %abrimos un fichero para almacenar los resultados

% Calculo de los fingerprints de los anuncios empleados en la prueba

p=1;

for N=[11960:-1000:2960]
    for time_detect=[8000 12000 16000 20000 24000]
          
        if (N<time_detect)
              
          f_22_23_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\22_23_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_22_23_02_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\22_23_02_8000_8'),N,over_lap,num_band,finf,fsup);
          f_22_23_03_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\22_23_03_8000_8'),N,over_lap,num_band,finf,fsup);
          f_22_23_04_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\22_23_04_8000_8'),N,over_lap,num_band,finf,fsup);
          f_01_07_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\01_07_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_08_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\08_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_09_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\09_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_10_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\10_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_11_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\11_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_12_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\12_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_13_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\13_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_14_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\14_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_17_19_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\17_19_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_18_23_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\18_23_01_8000_8'),N,over_lap,num_band,finf,fsup);
          f_20_23_01_8000_8=fingerprint_actualizable(wavread('Anuncios\Procesado\20_23_01_8000_8'),N,over_lap,num_band,finf,fsup);

          %simulamos como si se cogiera un intervalo de tiempo de una emision de
          %radio y se analizara si se corresponde con algun anuncio de la base de datos

          anuncio_detectado=wavread('Anuncios\Procesado\22_23_01_8000_8');
          t=1;
          s=1;
          key=1;

          for repeticion=1:1000                                             % repeticion para cada valor de N
                 l=round((length(anuncio_detectado)-time_detect)*rand(1));         %valor escogido al azar para tomar un trozo al azar del anuncio leido

                 t1=cputime;

                 finger_detectado=fingerprint_actualizable(anuncio_detectado(l:l+time_detect-1,1),N,over_lap,num_band,finf,fsup);

                 aux_times(repeticion)=cputime-t1;

                % La primera vez que entramos en el bucle, preparamos una variable,
                % 'probab', para almacenar las probabilidades de cada anuncio de ser
                % considerado como detectado

                 if key==1
                    key=0;
                    [m,n]=size(finger_detectado);
                    num_bits=m*n;                          %numero de bits que se comparan en la busqueda
                    probab_alarm=zeros(1,m*n);
                    probab_no_alarm=zeros(1,m*n);
                 end

                 j=1;
                 distancia(j)=min(busqueda(finger_detectado,f_22_23_01_8000_8));
                 if distancia(j)~=0
                     probab_alarm(1,distancia(j))=probab_alarm(1,distancia(j))+1;
                 end
                 dist_alarm(p,t)=distancia(j);                % almacena distancia en deteccion correcta para calculo de probabilidad de alarma
                 t=t+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_22_23_02_8000_8));
                 if distancia(j)~=0
                     probab_alarm(1,distancia(j))=probab_alarm(1,distancia(j))+1;
                 end
                 dist_alarm(p,t)=distancia(j);
                 t=t+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_22_23_03_8000_8));
                 if distancia(j)~=0
                      probab_alarm(1,distancia(j))=probab_alarm(1,distancia(j))+1;
                 end
                 dist_alarm(p,t)=distancia(j);
                 t=t+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_22_23_04_8000_8));
                 if distancia(j)~=0
                      probab_alarm(1,distancia(j))=probab_alarm(1,distancia(j))+1;
                 end
                 dist_alarm(p,t)=distancia(j);
                 t=t+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_01_07_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);          % almacena distancia en deteccion correcta para calculo de probabilidad de falsa alarma
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_08_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_09_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_10_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_11_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_12_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_13_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_14_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_17_19_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_18_23_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;
                 j=j+1;
                 distancia(j)=min(busqueda(finger_detectado,f_20_23_01_8000_8));
                 if distancia(j)~=0
                     probab_no_alarm(1,distancia(j))=probab_no_alarm(1,distancia(j))+1;
                 end
                 dist_no_alarm(p,s)=distancia(j);
                 s=s+1;

               [h,q]=min(distancia);
          
              if (h/num_bits<alpha)                  %si el porcentaje de bits incorrectos es superior al alpha% se considera que el trozo de audio que se esta 
                                                     %estudiando no se corresponde con ningun anuncio de nuestra base de anuncios
                  fprintf(fi,'Anuncio detectado: ANUNCIO NUMERO %1.0f con una probabilidad de deteccion correcta del %2.2f\n',q,(1-h/num_bits)*100);
              else
                  fprintf(fi,'NO SE HA DETECTADO NINGUN ANUNCIO. La probabilidad de deteccion mas alta ha sido de %2.2f\n',(1-h/num_bits)*100');
              end
          
          end

        times(p)=mean(aux_times);     % Guardamos el tiempo empleado en calcular el finprint para cada valor de N


        figure;
        plot((1:m*n)/m/n,probab_alarm/sum(probab_alarm),'rx')
        hold on;
        plot((1:m*n)/m/n,probab_no_alarm/sum(probab_no_alarm),'yx')
        plot((1:m*n)/m/n,zeros(1,m*n),'kx')
        title(['N = ',num2str(N),'  NumBand = ',num2str(num_band-1),'  DetectionTime = ',num2str(time_detect/8000),' seg.']);
        xlabel('Tasa de error de bit entre fingerprints');
        ylabel('pdf');

        %A�adimos una linea vertical que simula el umbral
        plot([alpha alpha],[0 max(probab_no_alarm/sum(probab_no_alarm))]) 
        hold off;


        % Para evitar problemas con 'wblfit'
        if(min(dist_alarm(p,:))==0)      
            dist_alarm(p,:)=dist_alarm(p,:)+10^(-6);
        end

        % Calculo de los parametros de la distribucion que mas se aproxima a cada
        % uno de los dos casos:
        [a,b]=wblfit(dist_alarm(p,:));
        [c,d]=normfit(dist_no_alarm(p,:));


        % Probabilidad de producirse una deteccion incorrecta

        prob1(p)=1-wblcdf(alpha,a(1,1)/(m*n),a(1,2));
        prob2(p)=normcdf(alpha,c/(m*n),d/(m*n));


        fprintf(fi,'Tama�o de la ventana de Hanning: %1.0f muestras\n',N);
        fprintf(fi,'Numero de bandas: %2.0f muestras\n',num_band-1);
        fprintf(fi,'Overlap: %1.4f muestras\n',over_lap);
        fprintf(fi,'Tiempo de captura: %1.1f segundos\n',time_detect/8000);
        fprintf(fi,'Tiempo en calcular el fingerprint: %2.4f   Numero de anuncios: %1.0f\n',times(p),j);
        fprintf(fi,'Numero de bits por fingerprint detectado: %3.0f x%3.0f\n',m,n);
        fprintf(fi,'Pn: %2.4e   Pf: %2.4e\n\n\n',prob1(p),prob2(p));


        p=p+1;
        end
    end
end



fi=fopen('p_global_numband_N_timedetect','r');
type p_global_numband_N_timedetect<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="es">
	<head>
            <title>e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.</title>

                <meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>
<meta name="description" content="Buscador de proyectos ETSI" />
<meta name="author" content="Daniel Callejo"/>
<meta name="author" content="Federico L&oacute;pez"/>
<meta name="author" content="Juli&aacute;n P&eacute;rez"/>
<meta name="author" content="Jos&eacute; Mar&iacute;a Vidal Vidal"/>
<base href="https://biblus.us.es/bibing/proyectos/" /><link rel="shortcut icon" href="./favicon.ico" />
<link rel="stylesheet" href="./publico/css/estructura.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/bootstrap.css" type="text/css" media="screen"/>
<link rel="stylesheet" href="./publico/css/enlaces-capas.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/pestanas.css" type="text/css"  /> 
<link rel="stylesheet" href="./publico/css/formularios.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/lightbox.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/autocompleter.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/resultados.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/milista.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/explorador.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/ultimos.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/contactar.css" type="text/css"  />
<!-- <link rel="stylesheet" href="./css/print.css" type="text/css" media="print" /> -->
<script src="./publico/javascript/jquery-1.9.1.min.js" type="text/javascript"></script>
<script src="./publico/javascript/bootstrap.js" type="text/javascript"></script>

<script src="./publico/javascript/mootools-1.2.4-core.js" type="text/javascript"></script>
<script src="./publico/javascript/mootools-1-1.2.4.2-more.js" type="text/javascript"></script>
<script src="./publico/javascript/targetbl.js" type="text/javascript"></script>
<script src="./publico/javascript/mensajes.js" type="text/javascript"></script>
<script src="./publico/javascript/form.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Observer.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.Request.js" type="text/javascript"></script>
<script src="./publico/javascript/domready.js" type="text/javascript"></script>

<!-- Google tag (gtag.js) -->
<script async src="https://www.googletagmanager.com/gtag/js?id=G-2HP38G406P"></script>
<script>
  window.dataLayer = window.dataLayer || [];
  function gtag(){dataLayer.push(arguments);}
  gtag('js', new Date());

  gtag('config', 'G-2HP38G406P');
</script>



<script type="text/javascript">
function seleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=1
}
function deseleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=0
} 
</script>
	</head>
    <body style="min-width: 960px; background-color: #e2e2e2; height: 100%">	
		<a name="arriba"></a>
<div>
<div id="header">
    <table id="Tabla_01" height="105" border="0" cellpadding="0" cellspacing="0">
    <tr>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_01.jpg" width="3" height="105" alt="">
        </td>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_02.jpg" width="100" height="105" alt="">
        </td>
        <td width="100%" align="right"  style="min-width: 613px; background-repeat: no-repeat; background-position: top right; font-size: medium; font-weight: bold; color: white;">
            <table width="100%" height="105" border="0">
                <tr>
                    <td style=" vertical-align: top; width: 547px">
                        <!-- <a style="text-decoration: none" href="./"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a> -->
                        <a style="text-decoration: none" href="https://biblus.us.es/bibing/proyectos/"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a>
                     </td>
                     <td  background="./publico/imagenes/cabecera/banner_new2_04.jpg" style="background-repeat: no-repeat; background-color: #c6d4e8;"><br><br><br>
                         &nbsp;
                        </td>
                    </tr>
             </table>
        </td>
    </tr>
    <tr style="border-top: 1px solid #ffffff;">
        <td colspan="3" id="menu_sup" style="font-size: small; text-align: right; background-color: #9c1f2f; vertical-align: bottom">
                                     

                                <a href="https://biblus.us.es/bibing/proyectos/" class="enlace " style=" color: #ffffff">Inicio&nbsp;&nbsp;|</a>
                                <a href="https://biblus.us.es/bibing/proyectos/acerca/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Acerca de&nbsp;&nbsp;|</a>
                                
                                <a href="https://biblus.us.es/bibing/proyectos/contactar/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Contactar&nbsp;&nbsp;</a>
                         <!--       <a href="https://biblus.us.es/bibing/proyectos/ayuda/" style=" color: #ffffff" class="enlace  ">&nbsp;&nbsp;Ayuda&nbsp;&nbsp;|</a> -->
                                <!-- <form id="idioma" class=""  method="post" action="">
                                <div class="inline" style="color: white">
                                <input type="submit" class="tipo_letra " style=" color: #ffffff"name="idioma" value="Castellanossss" />&nbsp;&nbsp;
                                </div>	
                                </form> -->

                                </td>

    </tr>
    </table>
</div>
<br><br>
<div class="colmask rightmenu">
<div class="panel panel-default" style="margin:20px">
            <div class="panel-heading"><h4>OPTIMIZACIÓN DE UN SISTEMA DE AUDIO FINGERPRINTING PARA LA DETECCIÓN DE ANUNCIOS EN TIEMPO REAL</h4><h5 style="color:grey"><span class="glyphicon glyphicon-user""></span> : <b>Álvarez Dujat Des Allimes, Alejandro</b><br><span class="glyphicon glyphicon-book""></span>  :<b> Ingenier&iacute;a Telecomunicación</b></h5></div><div class="panel-body"><div class="contenido_proyecto"><div class="titulo2_explorador">Contenido del proyecto:</div><div class="ruta"><a href="use/abreproy/11429">Directorio ra&iacute;z</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/11429/direccion/C%C3%B3digos+de+Matlab+%252F">Códigos de Matlab </a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/11429/direccion/C%C3%B3digos+de+Matlab+%252Fprueba_global_NumBand.m%252F">prueba_global_NumBand.m</a></div><div class="listing"></div></div></div></div>
</div>
</div>

        
	</body>
 
 
<!-- Inicio Google Analytics -->
<!--
	<script type="text/javascript">
		var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www.");
		document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E"));
	</script>
	
	<script type="text/javascript">
		try {
		var pageTracker = _gat._getTracker("UA-7789620-1");
		pageTracker._trackPageview();
	} 
		catch(err) {}
	</script>
-->
<!-- Fin Google Analytics -->

</html>