SoloCodigo

hola  
La respuesta en frecuencia no la puedes ver con abs, por que solo te da la magnitud de un valor complejo, la respuesta en frecuencia la puedes visualizar con freqz si tu filtro es digital y tienes los coeficientes en Z (aunque 1024 me parecen muchos coeficientes para un filtro )

o con freqs si tu filtro es analogico y esta definido por funciones en S.

Hola

Mira a ver si tienes instalados los toolbox necesarios, creo que el de freqz es "signal proccessing toolbox" , yo tengo la version 6.5 y si me viene la funcion freqz

Respecto a la respuesta en frecuencia, yo he usado freqz para obtener la respuesta de un filtro de 1D , en imagenes (2D) no se por que motivo, he usado abs(fft2(filtroM)), quizas no haya un freqz para el caso de filtros matriciales.
Tambien puedes usar la herramienta fvtool para estudiar la respuesta en frecuencia del filtro

http://www.mathworks.com/access/helpdesk/h...nal/fvtool.html

salu2  

Hola

He estado mirando y concretamente el archivo .m de la funcion freqz esta en el directorio filterdesign, y este es el toolbox

http://www.mathworks.com/products/filterdesign/tryit.html

Hay una verison trial que no se si te servira.Si no es asi dejame tu correo y te mandare todo el directorio del toolbox , aunque no se si funcionara metiendolo a pelo.

Otra solucion es bajarse una "copia de seguridad" del matlab mediante p2p  

salu2  

Saludos ...

bueno yo busque la funcion freqz y esta en la carpeta signal (quiza porque tengo matlab 7.0) pero si consigues la carpeta con los archivos, creo que solo bastaria con copiarla al directorio "toolbox" de matlab (a menos que internamente llame a otro toolbox)

bueno,  lo del filtrado...

si h es tu señal de tiempo discreto resultado del proceso de filtrado entonces si  
puedes usar

H=abs(fft(h));
plot(H)

y asi podras ver los componentes de frecuencia de h, (aunque las frecuencias estarian en el rango de 0-1 porque estan normalizadas y solo te sirve la mitad del espectro porque esta desde -w a +w)....

pero si h contiene lo coeficientes del propio filtro(por eso la duda de un filtro tan grande)
la respuesta en frecuencia la puedes ver con freqz

Bueno, si no recuerdo mal, la transformada de fourier devuelve un espectro de frecuencias positivas y negativas.
Esto es, la primera mitad del array de puntos devueltos por la fft es el correspondiente a las frecuencias positivas, el valor del medio es el correspondiente a una frecuencia de 0Hz y la segunda mitad del array es la parte correspondiente a frecuencias negativas....  Por tanto sólo interesa la primera mitad del array.
Además hay que realizar una serie de cálculos para poder visualizar correctamente la amplitud de los picos que corresponden a la frecuencia de la señal.  Prueba con lo siguiente:
Vamos a muestrear una señal con varias componentes (una de 50Hz y otra de 20Hz) a una frecuencia de 1000Hz.

t=[0:0.001:1];
h=2*sin(2*pi*50*t)+5*sin(2*pi*20*t);

H=fft(h);
ptos=H.*conj(H)/1024;
frec=1000*(0:512)/1024;
plot(frec,ptos(1:513));


Y teoricamente deberían aparecer dos picos bastante pronunciados a las frecuencias de 50 y 20Hz respectivamente.
Espero que te sirva de ayuda...
Y espero que esté bien el ejemplo....

Perdon, se me olvidó incluir lo siguiente en la llamada a la función fft:
fft (h,1024);

1024 son los puntos devueltos por la fft...
Lo puse así para que se pareciera a lo que querías hacer, variará según los puntos devueltos por la fft.

Lo de la frecuencia es simplemente un vector de puntos que coincidan con las muestras obtenidas de la señal:
Si hemos muestreado a 1000Hz, pues por Nyquist sabemos que la frecuencia de la señal muestreada variará entre 0 y 500Hz.  

Si lo has muestreado a 500Hz, pues tu vector de salida de frecuencia irá de 0 a 250Hz.

y para equiespaciar este vector de frecuencias basta con multiplicar el vector de frecuencia que va de 0 a 250 (si muestreas a 500Hz) por la frecuencia de muestreo y dividirlo por el número de muestras.

Un saludo.