Diseño de filtros digitales
10
1
Asignatura: Sistemas y señales discretos. Tema: Diseño de filtros digitales. Lugar de Ejecución: Instrumentación y control (Edificio de electrónica) Facilitador: Denis Altuve.Objetivos específicos
Realizar un filtro de Wiener-Kolmogoroff. Realizar un filtro de estadístico de mediana.
Introducción teórica
Filtro de Wiener-Kolmogoroff Es un tipo de filtro norecursivo que puede ser usado para eliminar ruido, ecualización, reducción de eco y predicción. Para hacer el diseño de los coeficientes del filtro se requiere obtener primero: Auto-correlación R yy yla Correlación cruzada R yx . Lo cual significa que se requiere conocer la señal de entrada y de salida que se desea filtrar. En Matlab, la función de correlación filtrada se obtiene con la funciónXCORR.
XCORR Estima la función de correlación-cruzada. XCORR(A,B), donde A y B son vectores de longitud M, regresa la secuencia de correlacióncruzada de longitud 2*M-l en un vector columna. XCORR(A),donde A es un vector, es la secuencia de auto-correlación. XCORR(A), donde A es una matriz M-por-N, es un matriz grande con 2*M-l filas cuyas N^2 columnas contienen las secuencias decorrelación-cruzada para todas las combinaciones de las columnas de A. Y para la generación de la matriz hermitiana se usa la función TOEPLITZ. TOEPLITZ(C) es una matriz Toeplitz simétrica (o Hermitiana) teniendo a Ccomo su primera columna y fila.
Materiales y equipos
1 Computadora con sistema operativo Windows 95 o superior 1 Programa MATLAB 5.x o superior
Procedimiento
PARTE I. Filtro deWiener-Kolmogoroff. 1. Escriba el siguiente guión que implementa el filtro: [x,fs,b]=wavread('pds'); n=1:100; %100 muestras t=(n-1)/fs;
Guía
y=x(n+5)+x(n); %5 muestras de atraso (eco) + la original x=x(n);figure(1) subplot(3,1,1) plot(t,x) subplot(3,1,2) plot(t,y) l=length(x); ryy=xcorr(y); ryx=xcorr(y,x); clear x ryy=toeplitz(ryy(l:2*l-1)); ryx=ryx(l:2*l-1); h=inv(ryy)*ryx; num=h'; den=zeros(1,l);...
Regístrate para leer el documento completo.