Home > PR > costeinvrespfour.m

costeinvrespfour

PURPOSE ^

costerespfour calcula los coeficientes en en dominio de fourier

SYNOPSIS ^

function F=costeinvrespfour(observado,lcuan,H,k1,k2,umbral);

DESCRIPTION ^

 costerespfour calcula los coeficientes en en dominio de fourier
 del operador de regularización.
 
 Devuelve la inversa de la respuesta a los coeficientes de la transformada
 del bloque.

  umbral se utiliza para determinar fijar el valor minimo de la respuesta
      valormin = max(max(respuesta))/umbral;
      Valores menores que este valor en la respuesta se ponen a valormin
 
 USO: C=costerespfour(bloque,ladobloque,H,k1,k2,umbral);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function F=costeinvrespfour(observado,lcuan,H,k1,k2,umbral);
0002   
0003 % costerespfour calcula los coeficientes en en dominio de fourier
0004 % del operador de regularización.
0005 %
0006 % Devuelve la inversa de la respuesta a los coeficientes de la transformada
0007 % del bloque.
0008 %
0009 %  umbral se utiliza para determinar fijar el valor minimo de la respuesta
0010 %      valormin = max(max(respuesta))/umbral;
0011 %      Valores menores que este valor en la respuesta se ponen a valormin
0012 %
0013 % USO: C=costerespfour(bloque,ladobloque,H,k1,k2,umbral);
0014 
0015 
0016 % Modelo no-lineal con interacción con los vecinos (H no es diagonal)
0017 
0018 C=respuestafourier(observado,lcuan,H,k1,k2);
0019 
0020 %figure
0021 %subplot(1,3,1),colormap(gray);imagesc(C.^0.3);axis image
0022 
0023 %C=conv2(C,fspecial('gaussian',5,.3),'same');
0024 
0025 %subplot(1,3,2),colormap(gray);imagesc(C.^0.3);axis image
0026 %subplot(1,3,3),colormap(gray);imagesc(abs(fftshift(cub).^0.3));axis image
0027 %pause,close;
0028 
0029 C(lcuan/2+1,lcuan/2+1) = 1;
0030 
0031 
0032 %umbral = zeros(32);
0033 %for i=24:32
0034 %   v = zeros(1,32);
0035 %   v(i) = 1;
0036 %   f = simrevol(v);
0037 %   f = double(f>0);
0038 %   valores = C.*f;
0039 %   media = sum(sum(valores))/sum(sum(f));
0040 %   umbral(f>0) = media;
0041 %
0042 %   %figure,mesh(umbral);pause;close;
0043 %
0044 %end
0045 
0046 
0047 F = C.*( C > umbral ) + umbral.*( C <= umbral );
0048 
0049 %
0050 %figure(200);
0051 %subplot(1,4,1);
0052 %prod(size(umbral))
0053 %if (prod(size(umbral))>1)
0054 %   mesh(umbral);
0055 %end
0056 %subplot(1,4,2);
0057 %mesh(C);
0058 %subplot(1,4,3);
0059 %mesh(F);
0060 
0061 
0062 F=1./F;
0063 
0064 F = F/max(max(F));
0065 
0066 F(lcuan/2+1,lcuan/2+1) = 1;
0067 
0068 
0069 %subplot(1,4,4);
0070 %mesh(F);
0071 %pause(0.2)
0072 
0073 
0074 
0075 
0076 
0077 %F(lcuan/2+1,lcuan/2+1) = 0;
0078 
0079 
0080 
0081 %figure,mesh(fftshift(F));pause;close;

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