Home > PR > maximo_respuesta.m

maximo_respuesta

PURPOSE ^

function mr = maximo_respuesta(observado,M,H,k1,k2,Nbloques);

SYNOPSIS ^

function mr = maximo_respuesta(observado,M,N,H,k1,k2,Nbloques);

DESCRIPTION ^

 function mr = maximo_respuesta(observado,M,H,k1,k2,Nbloques);

  Esta función se utiliza en costeinvrespfour para obtener el 
  maximo de respuesta en una serie de bloques. 
  Este máximo se utiliza (divido por un Juan's number)
  para establecer el corte al hacer la inversa de la respuesta.

   obsevado: imagen observada
   M: tamaño del bloque
   H: funciones de interaccion (se obtiene en constantes_sin_continua)
   k1: alpha (se obtiene en constantes_sin_continua)
   k2: beta (se obtiene en constantes_sin_continua)
   Nbloques: numero de bloques en los que obtener la respuesta

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % function mr = maximo_respuesta(observado,M,H,k1,k2,Nbloques);
0002 %
0003 %  Esta función se utiliza en costeinvrespfour para obtener el
0004 %  maximo de respuesta en una serie de bloques.
0005 %  Este máximo se utiliza (divido por un Juan's number)
0006 %  para establecer el corte al hacer la inversa de la respuesta.
0007 %
0008 %   obsevado: imagen observada
0009 %   M: tamaño del bloque
0010 %   H: funciones de interaccion (se obtiene en constantes_sin_continua)
0011 %   k1: alpha (se obtiene en constantes_sin_continua)
0012 %   k2: beta (se obtiene en constantes_sin_continua)
0013 %   Nbloques: numero de bloques en los que obtener la respuesta
0014 %
0015 %
0016 %
0017 %
0018 
0019 function mr = maximo_respuesta(observado,M,N,H,k1,k2,Nbloques);
0020 
0021 
0022 [fil,col]=size(observado);
0023 
0024 
0025 r=zeros(M);
0026 
0027 
0028 filas = round(M/2 + rand(Nbloques,1)*(fil-M));
0029 columnas = round(M/2 + rand(Nbloques,1)*(col-M));
0030 
0031 
0032 for i=1:Nbloques
0033    
0034    bl=observado(filas(i)-M/2+1:filas(i)+M/2,columnas(i)-M/2+1:columnas(i)+M/2);
0035 
0036    r = r+ respuestafourier(bl,M,H,k1,k2);
0037 
0038 end
0039 
0040 
0041 respuestamedia = r/Nbloques;
0042 
0043 mr = max(max(respuestamedia));

Generated on Wed 29-Nov-2006 16:19:19 by m2html © 2003