Y a aqui estan sorry.
:oops:
Bueno Primero algo de teoria muy sencilla.
Los filtros están hechos para pasar algunas frecuencias de ciertas señales o detener otras.
Ambos filtros tanto analógicos como filtros digitales pueden ser considerados como una caja negra, con una señal de entrada y otra de salida
Tipos de filtros digitales.
Este tipo de filtros son usados para remplazar un filtro analógico, procesan la señal en el dominio del tiempo,
donde las señales analógicas primero deben de ser muestreadas y digitalizadas en intervalos discretos usando un
convertidor analógico digital.
Ahí que tener precaución con este tipo de filtros debido a que la señal analógica es muestreada
ocasionando el efecto de alaising, este ocurre cuando la señal analógica tiene un espectro de energía
a frecuencias cercanas a media frecuencia de muestreo, ocasionando que sea imposible recobrar
la señal original cuando se quiere convertir de vuelta a analógica.
Para
prevenir este efecto se debe filtrar la señal con una frecuenta mayor a la de filtrado.
Ahí dos tipos de filtros digitales : Filtros de respuesta al impulso finito (FIR) que son filtros
no recursivos, y filtros de respuesta al impulso infinito( IIR) los cuales son recursivos debido a que
parte de la señal de salida se usa para realimentar la entrada.
Se dice que los filtros IIR son sensibles a los coeficientes del filtro usado para multiplicar factores,
debido a que al truncar los coeficientes puede resultar en una realimentación positiva causando
oscilación, aparte de tener una respuesta de fase no lineal.
Método de diseño.
La mayoría de los filtros IIR son diseñados usando el modelo de un filtro analógico, entre los cuales tenemos a los tipo Butterworth, Chebyshev, Cauer(eliptico) y Chebyshev Inverso
Los ejemplos solo seran del tipo Butterworth
Una forma es usando la transformacion Bilineal usada para convertir la respuesta en frecuencia analógica a una respuesta en el dominio digital.
Donde para pasar del plano s --> a Z tenemos
\s = \frac{1-z^{-1}}{1+z^{-1}}
Los filtros Basicos IIR estan basados en una estructura bicuadratica con elementos de retraso denotados como \z^{-1}
el cual se encuentra en la función de transferencia de forma general denotada como:
\H(z)=\frac{A0+A1z^{-1}+A2z^{-2}...}{1-B1z^{-1}-B2z^{-2}...}
...1
Y su esctructura Bicadratica a manera de bloques es
(http://img402.imageshack.us/img402/138/reograma.jpg)
Para el diseño de filtros de mas alto orden se diseñan mediante etapas bi-cuadráticas en cascada, donde cada una da una respuesta de segundo orden.
Por ejemplo si se requiere diseñar un filtro de orden 4 se usarían dos etapas bi-cuadráticas en serie.
Si se quisiera diseñar un filtro pasabajos de orden impar este requeriría de una o mas etapas de segundo orden , seguidas por la etapa de orden 1. Donde la etapa de primer orden es simplemente un retraso con un coeficiente de retroalimentación como el que se muestra a continuación.
(http://img704.imageshack.us/img704/2499/reograma1.jpg)
Donde la sección de primer orden es lo mismo que las secciones bi-cuadráticas con coeficientes A1,A2 y B2
Filtros Butterworth
Predeformación
Desafortunadamente la transformación bilineal es una aproximación que no produce exactamente la respuesta en frecuencia requerida y se produce una distorsión por lo cual se modifica la transformación bilineal antes de aplicar la transformación a H(s).
Esto lo podemos resumir en la siguiente tabla
(http://img402.imageshack.us/img402/5982/tablado.jpg)
Diseño de un filtro Pasa-Bajos
Diseñar un filtro IIR pasa-bajos de segundo orden con una banda de paso de 3.4 Khz y una frecuencia de muestreo de 8kHz.
Por lo que para la función de transferencia analógica de un filtro butterworth de 2do orden esta dado por
(http://img709.imageshack.us/img709/2226/2do.jpg)
La frecuencia de corte normalizada debido a la predeformación esta dada por.
(http://img168.imageshack.us/img168/7725/98869134.jpg)
Y la nueva transformación bilineal con la predeformación viene dada por la siguiente ecuación.
(http://img709.imageshack.us/img709/4756/wd1.jpg)
Usando las formulas de la tabla anteior
tenemos que
(http://img8.imageshack.us/img8/8991/18454731.jpg)
Y sustituyendo en la funcion de transferencia de segundo orden y reacomodando como en 1
(http://img193.imageshack.us/img193/1706/78341624.jpg)
A0=0.7157374
A1=1.4314748
A2=0.7157374
B1= 1.3489677
B2= 0.05139818
Donde /z^{-1}
representa n- retrasos de tiempo
Por lo que la ecuación en diferencias es:
(http://img402.imageshack.us/img402/3801/63500665.jpg)
Metodo Dos
Diseño Mediante Matlab
Matlab nos brinda la facilidad de poder realizar todo el calculo algebraico y proporcionarnos únicamente los coeficientes ya sea de un filtro analógico o digital mediante los siguientes comandos.
[A,B] = BUTTER(N,Wn,'tipo de filtro')
A, B son los coeficientes del filtro digital
N= Orden del filtro
Wn= frecuencia de corte normalizada
Usando el comando visto y los datos de diseño del ejemplo anterior tenemos.
N=2
(http://img19.imageshack.us/img19/5125/66431269.jpg)
y
(http://img8.imageshack.us/img8/9060/coef.jpg)
Que son los ceoficientes que obtuvimos anteriormente.
Obtención de la respuesta en frecuencia del filtro
Para la respuesta en frecuencia del filtro digital empleamos el comando freqz
FREQZ(A,B,F,Fs)
F (in Hz), frecuencia a probar(Hz)
Fs = frecuencia de muestreo (en Hz).
Freqz(A,B,5000,8000)
(http://img34.imageshack.us/img34/9799/respc.jpg)
Implementación de Filtros Digitales.
El Diagrama de bloques básico de un filtro digital es el siguiente
(http://img340.imageshack.us/img340/8527/bloques.jpg)
Acondicionamiento de la señal
Para acoplar la señal de entrada se ha hecho uso de amplificadores operacionales para agregarle un ofsset de DC y reducirla a tensiones en un rango de 0 [V] para –Vcc y 5 [V] para +Vcc.
(http://img25.imageshack.us/img25/6724/acond.jpg)
Donde el potenciómetro RV1 ajusta el offset de DC y RV2 ajusta la amplitud de la señal.
Y para la reconstrucción de la Señal de Salida se hará uso de un DAC de 8 bits (DAC0800). Cuyo circuito es el siguiente.
(http://img402.imageshack.us/img402/998/dacg.jpg)
Tomando en cuenta que la señal de entrada es de 20 Vpp tenemos que nuestra tensión de entrada en cualquier momento es
(http://img511.imageshack.us/img511/1218/formulas.jpg)
Programación del Filtro Digital hasta aqui le dejo pendiente y al rato subo lo demas