Asumo que la frecuencia (0,0) esta en el centro
0001 % Asumo que la frecuencia (0,0) esta en el centro 0002 0003 function operador = formatea_operador(op) 0004 0005 M = size(op,1); 0006 0007 modulo = abs(op); 0008 modulo = modulo-min(min(modulo)); 0009 modulo(M/2+1,M/2+1) = 0; 0010 0011 0012 %[f1,f2] = freqspace(M,'meshgrid'); 0013 %win = fspecial('gaussian',M,9); 0014 %win = win ./ max(win(:)); 0015 %h = fwind2(modulo,win); 0016 0017 h = fsamp2(modulo); 0018 0019 operador = freqz2(h,M); 0020 0021 %angulo = angle(operador); 0022 %operador = abs(operador); 0023 %operador = modulo-min(min(operador)); 0024 %operador(M/2+1,M/2+1) = 0; 0025 %operador = operador*exp(sqrt(-1)*angulo); 0026 0027 0028 %figure,mesh(abs(operador));pause,close; 0029 0030 operador = fftshift(operador); 0031