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