Grupo 4

TEMA: DISEÑO DEL FILTRO

DIGITAL IIR MEDIANTE EL
EMPLEO DE MATLAB
Diseño del filtro digital IIR médiate el empleo de MATLAB
• La herramienta de procedimiento de señal de
MATLAB incluye una variedad de archivos M para
diseñar filtros digitales tantos IIR como FIR.
• El proceso de diseño de filtros digitales IIR incluye
dos etapas. En la primera, el orden N del filtro y el
factor de escalamiento de frecuencia Wn se
determinan a partir de las especificaciones dadas.
Con base en estos parámetros y los rizos
especificados, se termina los coeficientes de la
función de transferencia en la siguiente etapa. A
continuación se describe la implementación de
MATLAB de estas dos etapas.
 Estimación de orden (primera etapa)

• Para estimar el orden del filtro digital IIR que

se va a diseñar mediante el método de la
transformación bilineal, los archivos M que se
utilizan son: buttord para los filtros
butterworth, cheblord para los chebyshes tipo
1, cheb2ord para los shebyshev tipo 2 y
ellipord para los elípticos. El uso de los
archivos M para estimar el orden se ilustra en
los siguientes ejemplo.
 EJEMPLO 1.1 orden mínima del filtro digital IIR pasaalta
Chebyshev tipo 2. Determinar el orden mínimo de la función de
transferencia de un filtro pasaalta digital Chebyshev tipo2 que
opera a una tasa de muestreo de 4 kHz con la siguiente
especificaciones:

• Borde de banda de paso 1000 Hz

• Borde de banda supresora de 40db
• Borde normalizado de la banda de paso Wp es de
2*1000/4000=0.5
• Borde normalizado de banda supresora Ws es de
2*600/4=0.3
En el diseño que utiliza el método de la transformación
bilineal se recurre a la instrucción.
[N,Wn]=cheb2ord(0.5,0.3,1.40);
El cual produce N=5 y Wn =(0.3224)
Diseño de Filtro (segunda etapa)
• Para diseñar filtros IIR basados en la transformación bilineal, el
Signal Processing Toolbox de MATLAB incluye funciones para cada
una de las cuatro técnicas de aproximación de magnitud. De
manera especifica, se dispone de los siguientes archivos M: butter
para diseñar filtros Butterworth, cheby1 para diseñar filtros
Chebyshev tipo1, cheby2 para diseñar filtros Chebyshev tipo2 y
ellip para el diseño de filtros elípticos. Es posible determinar los
coeficientes del numerador y del denominador de la función de
transferencia a partir de los últimos datos mediante el empleo de la
función zp2tf. De manera alternativa, la función zp2sos puede
utilizarse para determinar los factores de segundo orden del
numerador y del denominador de la función de transferencia.
• Después de que se ha calculado la función de transferencia, la
respuesta en frecuencia se calcula con el uso del archivo M freqz.
• En el ejemplo siguiente se ilustra el diseño de filtros pasabajas
digitales utilizando MATLAB.
EJEMPLO 1.3 Diseño de un filtro pasabajas IIR elíptico.
Determine la función de transferencia y grafique la ganancia de la
respuesta de un filtro pasabajas IIR elíptico con las siguientes
especificaciones:
• Borde de la banda de paso FP=800 Hz
• Borde de la banda supresora FS =1 kHz
• Rizo pasabanda de 0.5 dB, atenuación mínima de la banda
supresora de 40 dB
• Tasa de muestreo FT =4 kHz
• De acuerdo con las ecuaciones de frecuencias angulares
normalizadas y se llega a los bordes normalizados de la banda
dados por:
• ;
• Cuando se ejecuta el programa, pide que se escriba las
especificaciones del filtro. Primero se calcula el orden mínimo N del
filtro y la frecuencia de corte WN deseada necesaria para cumplir
con las especificaciones dadas. En el caso del filtro elíptico, WN=
WP= 0.4
La grafica muestra la ganancia de la respuesta
que se genera también los coeficientes del
numerador y el denominador de la función de
transferencia.
Diseño asistido por computadora por
computadora de filtros digitales IIR
• Los algoritmos de diseño de filtros digitales IIR. Se
basan en el diseño de un prototipo de un filtro
digital IIR. Estos algoritmo se usan en aplicaciones
que requieren filtros con una repuestas de
magnitud selectiva de frecuencia con
característica pasabajas, pasabanda o supresora
de banda.
• En esta sección primero se expone la idea básica
de las técnicas de diseño iterativo basadas en
computadoras y después describe una aplicación
especifica para ecualizar el retazo de grupo de
filtros digitales IIR.
 IDEA BASICA.
• El objetivo es determinar de manera iterativa los
coeficientes de la función de transferencia de manera
que se minimice la diferencia entre respuestas en
frecuencia de la función de transferencia digital H (
) y D( ) para todos los valores de w sobre sus
intervalos cerrado de( ) la cual se especifica como
una función de error ponderado E(w) dado por
• =( )
• =función de ponderación positiva especificada
por el usuario
• =respuesta en frecuencia de la función de
transferencia digital H(Z)
• =respuesta en amplitud deseada
• MEDIDAS DE APROXIMACION
• Criterio Chebyshev o minimax, se emplea para
minimizar el valor absoluto máximo del error de
ponderación E (W)
• donde R es un conjunto de bandas de frecuencias
disjuntas en el intervalo ( ),sobre el cual se define
la respuesta en frecuencia deseada. En aplicaciones de
filtrado, R está compuesta por las bandas de pasos y las
bandas de supresión del filtro que se va a diseñar. Por
ejemplo ,en un diseño de filtro pasa baja R es la unión
disjuntas de los intervalos de frecuencias ([ ] y [ ]),
donde y , son respectivamente el borde de la
banda de paso él y borde de la banda de supresión.
• Criterio del mínimo p, se usa para minimizar la
integral P-ésima potencia de la función de error
ponderada E (w): E=
• Sobre el intervalo de frecuencia especificado R
con P un entero positivo. Menudo y por
simplicidad se utiliza el criterio de mínimos
cuadrados que se obtiene de la anterior ecuación
con P=2.
• El filtro FIR resultante exhibe grande errores de
pico cerca de los bordes de banda debido al
fenómeno Gibbs. Por consiguiente =1, por lo
común no se utiliza. En la practica la medida del
error integral de la ecuación se aproxima por
medio de la suma finitada por:
• Donde Wi,(1 ) es una densa rejilla elegida
adecuadamente de frecuencias angulares digitales. El
criterio mínimo de cuadrados que se obtiene de la
ecuación de la suma fintada con p=2 se usa a menudo
por simplicidad.
• En el caso del diseño del filtro IIR, H ( ) y D ( ) se
sustituyen por sus funciones de magnitud. La función
deseada H(z) se supone que es una función real y
racional de z con ordenes fijos de los polinomios del
numerador y del denominador. Los paralelos del filtro
ajustable sin los coefientes de los polinomios del
numerador y del denominador o los polos y ceros de la
función de transferencia.
ECUALIZACION DEL RETRASO DE GRUPOS DE FILTROS DEL Y
DIGITALES IIR
• Para llegar a un filtro digital IIR selectivo de frecuencia con
un retraso de grupo constante un procedimiento practico
que se sigue a menudo es diseñar, primero un filtro IIR que
cumpla con las especificaciones de las respuestas de
magnitud y diseñar, después una selección pasa todas de
manera que el retraso de grupo total del filtro digital IIR en
cascada con sección pasa todas tenga un retraso de grupo
constante en la banda de paso.
• El ecualizador de retraso pasa todas suele diseñar mediante
el empleo de un método de optimización asistido por
computadora. Sea H(z) la función transferencia de un filtro
digital IIR con un retraso de grupo dado por (W). El
objetivo es diseñar una sección pasa todas estables con una
función de transferencia
|

A (Z)=
• con un retraso de grupo a(w)de manera que el retraso
de grupo total

• Del sistema en cascada sea aproximadamente

constante en la banda de paso del filtro. Además para
garantizar la estabilidad de sección pasatodas,
necesitamos asegurar que los coeficientes de la
función de la transferencia pasa todas satisfagan las
restricciones.
• |<1 , |<1 +
• Para)=el diseño del ecualizador de retraso pasa todas
puede formularse como un problema de optimización
mínimax en la cual se minimiza el valor absoluto
máximo de error ()
• en la banda de paso del filtro los parámetros
ajustables en el procedimiento de
optimización son el retraso deseado , y
los coeficientes de la función de transferencia
pasa todas.
 resumen
• El problema del diseño de filtros digitales se relaciona con
el desarrollo de una función de trasferencia adecuada que
cumpla con las especificaciones de respuesta en frecuencia
,la cuales , en este capítulo, se restringen a las
especificaciones suelen indicarse en términos de las
frecuencias deseadas del borde de las bandas de paso y de
supresión, así como de las desviaciones tolerables respecto
de los niveles de magnitud(ganancia) que se desean de las
bandas de paso y supresión. Este capítulo considero del
diseño de filtros digitales de respuesta al impulso infinita
(IIR), causales y estables.
• El diseño de filtro IIR suele efectuarse al transformar un
función de transferencia analógica prototipo mediante un
mapeo adecuado de la variable de frecuencia compleja en la
variable compleja z .El método de la trasformación bilineal ,
que se utiliza ampliamente y que se estudia en este capítulo,
se basa en dicho procedimiento.
• En este capítulo también se explica algunos de
los algoritmos del diseño de filtros digitales
IIR que están disponibles e el signal processing
toolbox de la MATLAB como funciones. En
particular, incluye el diseño de filtros digitales
IIR con respuesta de magnitud de
butterworth, chebyshev y elíptica. Por último,
se explica la idea básica del diseño de filtros
digitales IIR mediante el empleo el empleo de
técnicas iterativas asistidas por computadora y
describe una aplicación especifica de este
procedimiento en el diseño de ecuaciones
de retraso de grupo.
 EJERCICIOS CON MATLAB

M9.1 diseñe un filtro pasa bajas butterworth digital

que opere una tasa de muestreo de 100kHz con una
frecuencia de corte de 10 kHz a 0.4 dB y una
atenuación mínima de la banda supresora de 50
dB a 30 kHz mediante el empleo del método de la
transformación bilineal. Determine el orden del
prototipo del filtro analógico mediante el uso del
archivo M buttap de MATLAB. Trasforme la función
de trasferencia digital que se desea empleando el
archivo M bilineal. Grafique las curvas de ganancia
de las respuestas y de las de fase utilizando
MATLAB .indique los pasos utilizados en le
diseño.
• M 9.2 modifique el programa 9.3 para
diseñar un filtro pasa bajás butterworth
digital mediante el método de la
transformación bilineal. los datos de entrada
que requiere el programa modificado deben
ser los bordes de la banda de paso y de la
banda supresora que se desean y de la
desviación máxima de la banda de paso, así
como la atenuación mínima de banda de
supresión en dB.con el uso del programa
modificado ,diseñe el filtro pasa bajas
butterworth digital del ejercicio M9.1.