fs = frecuencia de muestreo M = tamaņo del bloque (y por tanto del operador) function opCSF = operadorCSF(fs,M);
0001 % fs = frecuencia de muestreo 0002 % M = tamaņo del bloque (y por tanto del operador) 0003 % 0004 % function opCSF = operadorCSF(fs,M); 0005 % 0006 function opCSF = operadorCSF(fs,M); 0007 0008 [fx,fy]=freqspace(M,'meshgrid'); 0009 0010 Fx=fx*fs/2; 0011 0012 Fy=fy*fs/2; 0013 0014 F=sqrt(Fx.^2+Fy.^2); 0015 0016 % Reordenacion 1D del dominio 0017 0018 f=zigzag(F); 0019 0020 [IAF,CSF]=iafet(f,1,0,0,[0 0 0],0); 0021 0022 CSF = dezigzag(CSF); 0023 0024 operador = 1./CSF; 0025 0026 opCSF = ifftshift(operador); 0027 0028 opCSF(1,1)=0; 0029 0030 opCSF = opCSF/max(max(opCSF)); 0031