Home > PR > formatea_operador.m

formatea_operador

PURPOSE ^

Asumo que la frecuencia (0,0) esta en el centro

SYNOPSIS ^

function operador = formatea_operador(op)

DESCRIPTION ^

 Asumo que la frecuencia (0,0) esta en el centro

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

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