hablemos sobre filtros digitales....

[cerebro]

Hola gente!
Empece a realizar una breve explicación sobre filtros digitales y queria compartirla con ustedes a ver que opinan, he consultado varios libros, apuntes de clases y también ideas que me han dado vuelta por la cabeza y realizado una especie de minilibro...
Esta es la primera parte y primera version,  en este caso me dedique unicamente a desarrollar contenido teorico sin relacion con los microcontroladores todo en general pero seguro que continuaré trabajando en los momentos libres para abarcar mayor cantidad de temas y ver que tanto se puede meter mano a los filtros en pic's.
Les dejo los temas que he explicado hasta el momento (de manera muy breve, yo creo que cada tema puede llegar a ser un solo libro  pero como guia....).

Sistemas digitales                                                   2
La Transformada Z                                                   3
Propiedades de la Transformada Z                           3
Función transferencia de un sistema                           4
Transformada inversa                                            4
Transformada de Fourier                                           5
Diseño de filtros                                                           6
Diseños de Filtros IIR (respuesta infinita al impulso)   7
técnica de Euler                                                           8
técnica Invariante al impulso                                   9
Transformación Bilineal                                          10
   
Los sistemas digitales es muy pero muy pobre, porque me parecio que no tenia nada de nuevo como para perder tiempo explicando....  dejo el documento en formato DOC para que si alguien esta interesado en agregar cosas lo pueda hacer y compartir con todos. Tiene mucha matemática en juego asi que no duden en preguntar.

Saludos, espero no los aburra... y se copen 

[migsantiago]

Hola

Deberías ponerle tu nombre al documento 

Actualmente yo estoy diseñando un filtro IIR para la ponderación A usando la bilineal Z. Según yo ya lo terminé pero necesito probar su respuesta en frecuencia. Ya tengo los coeficientes de las 'z'. ¿Sabes cómo probar su respuesta de frecuencia usando Matlab?

La ecuación a la que llegué es:



El denominador lo escribí en dos renglones pero es uno mismo  Quedaron las 'z' ordenadas de z^0 a z^-6, pero ya no sé cómo ponerla a prueba en Matlab.

[gera]

Cerebro compadre, esto lo tendrias q haber posteado hace 3 meses, antes de q rinda analisis de señales jaja. Te dejo el sitio web donde mi profe subio todos los apuntes de clase. Quizas puedas sacar algo de material de ahi, esta muy completo, aunq el tipo tiende a irse por las ramas jeje.
http://www.geocities.com/analisis_de_senales/
saludos y felicidades por la iniciativa, me parece muy interesante

[MGLSOFT]

Donde puedo encontrar informacion "entendible para un no universitario" sobre filtros Butterwoth / Bessel que se entienda y que tenga si es posible ejemplos en lenguaje C ???

[cerebro]

hola mig
pero cual es el problema? no podes represantar la funcion transferencia? el comando que suelo utilizar es freqz() para obtener la repuesta en frecuencia de un filtro digital.
Calculaste los coeficientes con la funcion filter() ? edit: conviene usar [NUMd,DENd] = BILINEAR(NUM,DEN,Fs) es mas sencilla.
No le quise poner nombre porque tenia unas imagenes robadas de libros, cuando cambie las imagenes y mejore un poquito el texto se le pongo.
  Saludos

Donde puedo encontrar informacion "entendible para un no universitario" sobre filtros Butterwoth / Bessel que se entienda y que tenga si es posible ejemplos en lenguaje C ???

filtros butterworth y bessel son distintas aproximaciones para a un filtro pasabajos normalizado. butterworth aproxima magnitud y bessel aproxima fase (pero trabaja con retardo de grupo que es la derivada del retardo de fase por ser mas trabajable). y se utilizan estas aproximaciones para disminuir las exigencias sobre el filtro (inecesarias), lo que quiero decir es que podes construir filtros de esas caracteristicas con filtros pasivos con componentes (R C,L C) o con algun filtro activo (sallem key, v. de estado) o bien con un filtro digital para el cual teniendo la funcion transferencia del filtro analogico podes aplicar los metodos que explique y despues en C es facil implementar las ecuaciones en diferencias  .

Te dejo una hoja de datos sobre filtros analogicos de national.

http://www.national.com/an/AN/AN-779.pdf

espero que te sirva 
saludos!

Gracias gera seguro le pego una leida mi idea es llegar hasta filtros adaptativos, si es que puedo...

[migsantiago]

hola mig
pero cual es el problema? no podes represantar la funcion transferencia? el comando que suelo utilizar es freqz() para obtener la repuesta en frecuencia de un filtro digital.
Calculaste los coeficientes con la funcion filter() ? edit: conviene usar [NUMd,DENd] = BILINEAR(NUM,DEN,Fs) es mas sencilla.
No le quise poner nombre porque tenia unas imagenes robadas de libros, cuando cambie las imagenes y mejore un poquito el texto se le pongo.
  Saludos

Ah, ok, entonces lo voy a probar con freqz y la otra que no conocía y te aviso. 

[cerebro]

Hola gente les entrego una segunda parte de mi pequeña aventura en el mundo filtrado, esta vez decidi hacer un ejemplo sobre la primera parte que habia realizado solo teoria. La idea es basicamente analizar un poco la señal en MatLab, obtener los coeficientes por bilinear de un filtro butterworth analizar que sucede con la señal filtrada y luego implementarlo en un micro he optado por un 16f877A y visualizar la salida con un conversor mcp4921. Hice una explicación de todo en un pdf (le agregue mi nombre ya que no tiene imagenes robadas  ) y les paso todos los archivos (pdf,c,txt,m,sch). 

Saludos! 

Cosas para descargar

Documento pdf
Todos los Archivos

[migsantiago]

¿Quién necesita a los dspic si tenemos al poderoso PIC16F877? 

Excelente ejemplo Cerebro.

[beatAurex1]

Muy buen ejemplo. Una pregunta si diseño un filtro pasobajos 2do orden butterworth con fc=3.4Khz  y f muestreo de 8kHz podria ser suficiente con un pic18f4550 a 48 Mhz o ya seria encesario un dspic?.

[CAT_12]

Buenas tardes, primero que nada felicidades por ese trabajo, solicito su ayuda a quien haya consultado este trabajo o al mismo autor para ver si me podria explicar como sacar lo de la funcion de transferencia del archivo en matlab.
Me refiero a lo siguiente que esta en el pdf

%Calculo de la funcion transferencia del filtro

a=[1];                                         %num
b=[1/(251.3^3), 2/(251.3^2), 2/251.3, 1];     %den

H=tf(a,b)                     %obtenemos la funcion transferencia del filtro

Pues en el programa que revise de matlab aparce asi y no se que hizo para poner esos valores en los numeradores y denominadores.
O como los obtengo, pues supongo que son especificos para cada señal.

gracias.