Home > PR > operadorCSF.m

operadorCSF

PURPOSE ^

fs = frecuencia de muestreo

SYNOPSIS ^

function opCSF = operadorCSF(fs,M);

DESCRIPTION ^

 fs = frecuencia de muestreo
 M  = tamaņo del bloque (y por tanto del operador)

 function opCSF = operadorCSF(fs,M);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

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