Diseño de Filtros Analogos - Del Libro PDS

Diseno de filtros en tiempo continuo
4
Quiza una de las operaciones mas comunes para el ingeniero en comuni-
caciones y electronica es el filtrado de senales. Ya sea una aplicacion en sis-
temas de comunicacion, control, o instrumentacion, el filtrado aparecera en
alguna fase del sistema. As como el procesamiento de la senal es la operacion
que nos permite extraer o transformar la informacion util al observador, el
filtrado nos permite alterar o seleccionar algun segmento de frecuencias de la
senal a ser filtrada. Esta operacion de filtrado permite atenuar algunas fre-
cuencias mientras otras pueden permanecer sin alterar o con una alteracion
mnima (si no se desea introducir alguna ganancia en amplitud).
Dado que el diseno de una cierta clase de filtros digitales esta basado en
la transformacion de un filtro analogico, en este captulo se estudia el diseno
de filtros en tiempo continuo, o analogicos, que despues seran empleados en
el diseno de filtros digitales con respuesta infinita al impulso (RII).
Por diseno entendemos la obtencion de la funcion de transferencia H(s),
a partir de las especificaciones dadas sobre la respuesta en frecuencia. Son
varios metodos de aproximacion de filtros analogicos, funciones de aproxi-
macion entre ellos: filtros de Butterworth, filtros de Chebyshev, y filtros
Elpticos. Se dice aproximacion porque un filtro ideal dado a partir de sus
propiedades de su respuesta en frecuencia, es decir
(
1, si c
H(j) = (4.1)
0, si > c ,
donde c es la frecuencia de corte, es fsicamente no realizable, o no causal.

Esto se puede demostrar al obtener la respuesta al impulso h(t) de la funcion
de transferencia dada por H(j).
79
80 4.1. ESPECIFICACIONES DE UN FILTRO
4.1. Especificaciones de un filtro

Generalmente la especificacion de un filtro incluye tres regiones: una
banda pasante, una banda de transicion, y una banda de rechazo. En la
figura 4.1 se muestran estas regiones para el caso de un filtro pasa-bajas,
donde k1 y k2 son las ganancias a las frecuencias de corte c y de rechazo
r , respectivamente.
k1
banda banda banda de

pasante de transicin rechazo
k2
c
r
Figura 4.1: Especificaciones de un filtro.
4.2. La aproximacion de Butterworth

En la aproximacion de Butterworth, tambien llamada aproximacion con
amplitud maximamente plana, la amplitud al cuadrado del filtro esta dada
por
1
|H(j)|2 = 2n , (4.2)
1 + c
donde n es el orden de la funcion de transferencia correspondiente y c

esta definida como la frecuencia de corte de
-3 dB. Puede verse que, cuando
= c , la respuesta en magnitud sera 1/ 2 veces la ganancia de corriente
directa, es decir, esto corresponde a una atenuacion de 3.01 dB, si se tiene
la ganancia en dB. Mas aun, podemos hacer el siguiente analisis.
Sea la ganancia Gn (), entonces
Gn () = 20 log |H(j)| = 10 log |H(j)|2n

1 2n

= 10 log 2n = 10 log 1 + , (4.3)
[1 + c

] c
donde log es el logaritmo base 10.
G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

4.2. LA APROXIMACION DE BUTTERWORTH 81
Para << c , tenemos que

2n
0, y entonces
c
Gn () = 10 log |1 + 0| = 0.
Para >> c ,
2n

1 << y entonces
c
2n

Gn () 10 log[ ] = 20 n log ,
c c
es decir, presenta una atenuacion de 20 n dB/decada.
En la figura 4.2 se muestra la respuesta en frecuencia para el filtro de
Butterworth de orden 1 a 3, donde claramente se ve la atenuacion de 20n
db/decada.
0
10
n=1
20
Magnitud en dB
n=2
30
40 n=3
50
60
70 1 0 1
10 10 10
Frecuencia en (rad/s)
Figura 4.2: Respuesta en frecuencia para filtros Butterworth, orden 1 a 3.
Para encontrar la funcion de transferencia decimos que = s/j, ya que

s = j, y con c = 1 (frecuencia de corte normalizada) tenemos
1
|Hn (j)|2 = Hn (j)H(j) = .
1 + (s/j)2n
Para encontrar los polos de la funcion de transferencia necesitamos encontrar
las races del denominador. Estas races estan dadas por la ecuacion (6.8):
(

si n es par, Sk = 16 2n + k
n , k=0,1,. . . (2n-1)
H(j) = k (4.4)
si n es impar, Sk = 16 n , k=0,1,. . . (2n-1),
donde Sk es la k-esima raiz o polo del denominador, con magnitud unitaria

y angulo 6 .
Con el uso de (6.8) se pueden obtener los polinomios de Butterworth del
orden n deseado. En la tabla 4.1 se muestran algunos de estos polinomios.

82 4.2. LA APROXIMACION DE BUTTERWORTH
Tabla 4.1: Algunos polinomios de Butterworth.

orden n Bn (s)
1 (s + 1)
2 (s2 + 2s + 1)
3 (s2 + s + 1)(s + 1)
4 (s2 + 0.765s + 1)(s2 + 1.84s + 1)
.. ..
. .
Ejemplo 4.2.1 Encontrar la funcion de transferencia H 1 (s) para un filtro

de Butterworth normalizado de orden 1.
Solucion:
Como n = 1 tenemos dos races ya que k = 0 y 1, as
0
S0 = 1 6
1
1
S1 = 1 6 .
1
Estas dos races se muestran en la figura 4.3. Tomando los polos del lado
izquierdo, es decir, los polos que daran una funcion de transferencia estable,
nos queda
1 1
H1 (s) = = .
s (1) s+1
Planos
1.5
0.5
Eje Imaginario
0.5
1.5
1.5 1 0.5 0 0.5 1 1.5
Eje Real
Figura 4.3: Races de la funcion de Butterworth para n=1.
Ejemplo 4.2.2 Encontrar la funcion de transferencia H 2 (s) para el filtro

de Butterworth normalizado de orden 2.
Solucion:
Como n = 2 tenemos cuatro races, con k = 0, 1, 2, y 3.
k
Sk = 1 6 + , k = 0, 1, 2, 3.
4 2

Desarrollando para n=2 y cada valor de k se obtiene

S0 = 1 6 ,
4
S1 = 16 + ,
4 2
S1 = 16 + ,
4
3
S1 = 16 + ,
4 2
Estas cuatro races se muestran en la figura 4.4. Tomando los polos del lado
izquierdo, es decir, los polos que daran una funcion de transferencia estable,
nos queda
1
H2 (s) =
[s (0.7071 + j0.7071)][s (0.7071 j0.7071)]
1
H2 (s) = .
s2 + 2s + 1
Planos
1.5
0.5
Eje Imaginario
0.5
1.5
1.5 1 0.5 0 0.5 1 1.5
Eje Real
En las figuras 4.5 y 4.6 se muestran los polos (o races) en el plano-s

para los polinomios de Butterworth de orden n=6 y n=7. Puede observarse
que los polos de Butterworth estan distribuidos en un crculo de radio 1.

Planos
Eje
Imaginario
j
Eje
Real
Planos
Eje
Imaginario
j
Eje
Real
1 1

Tabla 4.2: Polinomios de Butterworth de orden n, Bn (s) = bn sn + bn1 sn1 +

. . . + b1 s + b 0 .
orden Bn (s) orden Bn (s)
n=1 b0 =1.0000 n=2 b0 = 1.0000
b1 =1.0000 b1 = 1.4142
b2 = 1.0000
n=3 b0 = 1.0000 n=4 b0 = 1.0000
b1 = 2.0000 b1 = 2.6131
b2 = 2.0000 b2 = 3.4142
b3 = 1.0000 b3 = 2.6131
b4 = 1.0000
n=5 b0 = 1.0000 n=6 b0 = 1.0000
b1 = 3.2361 b1 = 3.8637
b2 = 5.2361 b2 = 7.4641
b3 = 5.2361 b3 = 9.1416
b4 = 3.2361 b4 = 7.4641
b5 = 1.0000 b5 = 3.8637
b6 = 1.0000
n=7 b0 = 1.0000 n=8 b0 = 1.0000
b1 = 4.4940 b1 = 5.1258
b2 = 10.0978 b2 = 13.1371
b3 = 14.5918 b3 = 21.8462
b4 = 14.5918 b4 = 25.6884
b5 = 10.0978 b5 = 21.8462
b6 = 4.4940 b6 = 13.1371
b7 = 1.0000 b7 = 5.1258
b8 = 1.0000
Para el diseno de filtros de Butterworth haremos uso de tablas que con-

tienen los polinomios de Butterworth, una vez que se haya determinado el
orden del filtro deseado. En la tabla 4.2 se muestran los valores para cons-
truir las funciones de transferencia de filtros Butterworth normalizados de
orden 1 a 8. Conociendo el orden n del filtro, decimos que H(s) = 1/B n (s)
donde Bn (S) son los polinomios de Butterworth.
4.2.1. Transformaciones analogico-analogico

En la seccion anterior se vio como encontrar la funcion de transferencia
de un filtro de Butterworth normalizado de orden n. Si los polinomios de
Butterworth se anotan en forma de tabla es muy facil obtener la funcion de
transferencia del filtro de Butterworth normalizado correspondiente. Como
su nombre lo indica, estos son filtros normalizados, es decir, tienen una
frecuencia de corte de 1 rad/s; sin embargo, en una aplicacion practica, rara

vez necesitaremos un filtro precisamente con esa frecuencia de corte. En esta

seccion, se empleara una transformacion que permita al filtro de orden n, que
inicialmente tiene una frecuencia de corte de 1 rad/s, tener una frecuencia
de corte segun se desee.
Las transformaciones llamadas transformaciones de frecuencia, o trans-
formaciones analogico-analogico permiten no solo cambiar la frecuencia de
corte sino, tambien cambiar el tipo de respuesta en frecuencia del filtro, es
decir, es posible transformar un filtro normalizado pasa-bajas en un filtro
pasa-altas con frecuencia de corte d , o bien a un filtro pasa-banda, o a un
filtro rechazo de banda.
Si en nuestra funcion de transferencia H(s), reemplazamos s por s/ d ,
obtenemos una nueva funcion de transferencia H 0 (s):

H 0 (s) = H(s)

= H(s/d ),
ss/d
donde H 0 (s) es la nueva funcion de transferencia, obtenida a partir de H(s)

con una frecuencia de corte deseada d .
Para obtener la respuesta en frecuencia evaluamos H 0 (s) en s = j, as:
|H 0 (j)| = |H(j/d )|,
y cuando = d ,
|H 0 (j)| = |H(j1)|,
es decir, |H 0 (j)| es igual al valor de la funcion de transferencia del filtro
normalizado en = 1. De esta manera hemos desplazado la frecuencia
de corte de 1 rad/s a d rad/s. Esto tambien se llama escalamiento en
frecuencia, y puede emplearse para transformar un filtro Pb (pasa-bajas) en
Pa (pasa-altas), o de un filtro Pb a PB (pasa-banda), entre otros.
Por medio de la transformacion s = d /s, se obtiene una respuesta en
frecuencia pasa-altas con frecuencia de corte d .
Transformacion pasa-banda
Como se vio en la seccion anterior, la transformacion s/ u transfor-
ma una respuesta en frecuencia normalizada pasa-bajas en una respuesta
en frecuencia no normalizada pasa-bajas. Para obtener una transformacion
pasa-bajaspasa-banda (Pb-PB) se requiere una transformacion doble. Si
tenemos un FPb caracterizado por alguna distribucion de ceros y polos so-
bre el plano s, estos tendran su eje de simetra alrededor del eje real (). Si
desplazamos el eje de simetra a alguna frecuencia 0 , entonces 0 reempla-
zara a . Esto requiere de una transformacion 0 ; sin embargo, esta
transformacion no es util, ya que para tener un filtro realizable se requieren
polos complejos conjugados, as que la transformacion que se utiliza es
s 0
s + . (4.5)
0 s

Si se desea, esto podra considerarse como una doble transformacion, una

pasa-altas ( s0 ) mas una pasa-bajas ( s0 ). Ademas,
s 0 s2 + 2
+ = , (4.6)
0 s 0 s
donde el termino (s2 + 2 ) nos da los polos complejos conjugados, mientras

que el termino 1/s dara los ceros en el origen y en el infinito lo que asegura
una respuesta pasa-banda. El termino 1/ 0 afecta al ancho de banda.
En las figuras 4.7 a 4.9 se muestra la descripcion grafica de la transfor-
macion Pb-PB vista como un desplazamiento del eje real a 0 . En la figura
4.7 se muestra la posicion de un polo simple para un FPb Butterworth de
primer orden. En la figura 4.8 se muestra la nueva posicion del polo cuando
el eje real se ha desplazado a 0 dando un par de polos complejos conjuga-
dos. Es importante notar la inclusion del cero en el origen para garantizar la
respuesta pasa-banda. Finalmente, en la figura 4.9 se muestra la respuesta
en frecuencia correspondiente a la distribucion de polos y ceros de la figura
4.8.
Planos
1.5
j
1
0.5 polo de un
Eje Imaginario
filtro pasabajas
0
0.5
1.5
1.5 1 0.5 0 0.5 1 1.5
Eje Real
Figura 4.7: Posicion del polo pasa-bajas Butterworth de primer orden.
Planos
1.5
j
1
desplazamiento
0.5 del eje real
Eje Imaginario
polos a +0
complejos
0
conjugados
0.5
1.5
1.5 1 0.5 0 0.5 1 1.5
Eje Real
Figura 4.8: Transformacion del eje real a 0 .

Magnitud
3
6
0
10
Frecuencia [rad/s]
Figura 4.9: Respuesta pasa-banda por desplazamiento del eje real.
De lo anterior se puede deducir que el FPB resultante tendra un orden

que es el doble del equivalente pasa-bajas. Cabe hacer notar que la transfor-
macion por lo general no es aritmeticamente simetrica. Si L y U son los
extremos de la banda de paso, entonces el centro de la banda esta dado por
p
0 = L U , (4.7)
y no como 0 = (L + U )/2. Para filtros con un ancho de banda estrecho

la diferencia entre la media geometrica y la media aritmetica es pequena.
En terminos del ancho de banda W = U L y de la frecuencia 0 =
U L , la transformacion requerida queda como
0 s 0 s2 + U L
s= ( + )= . (4.8)
W 0 s s(U L )
Para obtener una transformacion PbPB la transformacion requerida es
la inversa de (4.8) es decir:
s(U L )
s= . (4.9)
s2 + U L
4.2.2. Diseno de filtros Butterworth pasa-bajas

Se ha visto que es posible desplazar la frecuencia de corte de un filtro
normalizado de 1 rad/s a cualquier frecuencia deseada, pero falta determinar
cual es el orden requerido del filtro Butterworth para cumplir con los reque-
rimientos de ganancia en la banda pasante, as como la atenuacion mnima
en la banda de rechazo.
Normalmente, los requerimientos del filtro se especificaran con cuatro
parametros: k1 , k2 , 1 , 2 , donde 1 = c , y 2 = r , corte y rechazo
respectivamente. Empleamos 1 y 2 en lugar de c y r , y reservamos c
para el filtro normalizado. Ademas, como se vera mas adelante, la n depende
de (1 /2 ) en filtros Butterworth y 2 /1 en filtros Chebyshev y parece
mas facil de recordar de este modo.

De la figura 4.1, se tienen las siguientes condiciones sobre la respuesta

en frecuencia deseada,
0 20 log |H(j)| k1 , para toda 1 y

20 log |H(j)| k2 , para toda 2 . (4.10)
Si evaluamos (4.10) con H(j) = 1/Bn (j), se tiene que cuando = 1
se debe satisfacer

1
10 log 2n = k1 y (4.11)
1
1+ c

1
10 log 2n = k2 , (4.12)
2
1+ c
de donde 2n
1
= 10k1/10 1 y (4.13)
c
2n
2
= 10k2/10 1. (4.14)
c
Dividiendo (4.13) entre (4.14)
2n
1 10k1 /10 1
= , (4.15)
2 10k2 /10 1
para que, finalmente, resolviendo para n,
l log[(10k1 /10 1)/(10k2/10 1)] m
n= , (4.16)
2 log(1 /2 )
l m
donde indica redondeo hacia arriba, as que, si n es fraccionario, em-
plearemos el numero entero siguiente superior.
Una vez determinado el valor numerico de n, tomamos la funcion pro-
totipo de la tabla 4.2 de polinomios de Butterwoth. Ahora, para que esta
funcion de transferencia, construida con los polinomios de la tabla 4.2, tenga
una frecuencia de corte 1 , si y solo si k1 = 3dB, debemos desplazar la
frecuencia de corte normalizada a 1 con la transformacion s s/1 . Si k1
es de valor arbitrario debemos hacer un ajuste al escalar la frecuencia de la
funcion prototipo.
Para realizar este ajuste tenemos dos opciones. Si nos interesa mas cum-
plir con el requerimiento de ganancia en 1 , entonces calculamos la frecuen-
cia de escalamiento con
1
esc = k1/10
, (4.17)
(10 1)1/2n

y entonces aplicaremos la transformacion s = s/ esc .

Si, por el contrario, interesa mas exceder el requerimiento en 2 entonces
calculamos la frecuencia de escalamiento con
2
esc = (4.18)
(10k2/10
1)1/2n
para aplicar la transformacion.
Ejemplo 4.2.3 Disenar un filtro con respuesta monotonica con las siguien-
tes especificaciones, frecuencia de corte, f c , de 1000 Hz con -2 dB y una
atenuacion mnima de -20 dB a una frecuencia de rechazo, f r , de 3000 Hz.
La frecuencia en Hz se denotara por f .
Solucion:
Tenemos pues, k1 = 2, 1 = 2fc = 21000 k2 = 20, 2 = 2fr =
23000, sustituyendo
l log[(100.2 1)/(102 1)] m
n= =3
2 log(1000/3000)
entonces, tomamos la funcion prototipo

1
HB3 (s) = ,
(s2 + s + 1)(s + 1)
ahora aplicamos una transformacion PbPb, con
6.2831853 103
esc = = 6.8706913 103
(100.2 1)1/6
y
1
H(s) = | 3
(s2 + s + 1)(s + 1) s=s/6.870691310
para obtener
3.243406044 1011
H(s) =
s3 + 13741.383s2 + 94412799.811s + 3.243406044 1011
En la figura 4.10 se muestra la respuesta en frecuencia de la funcion de
transferencia encontrada. Ademas, en las figuras 4.11 y 4.12 se puede obser-
var, despues de hacer una ampliacion, que se cumple con el requerimiento
de -2 dB en 1000 Hz y se excede, con poco mas de 5 dB, el requerimiento
en 3000 Hz.

10
10
Ganancia en (dB)
20
30
40
50
60
70 1 2 3 4 5
10 10 10 10 10
Frecuencia (Hz)
Figura 4.10: Respuesta en frecuencia del filtro del ejemplo 4.2.3.
0
Ganancia en (dB)
7
3
10
Frecuencia (Hz)
Figura 4.11: Ampliacion en k1 = 2 dB.
10
15
20
Ganancia en (dB)
25
30
35
40
45
3 4
10 10
Frecuencia (Hz)
Figura 4.12: Ampliacion en k2 < 20 dB.

4.2.3. Diseno de filtros Butteworth pasa-banda

El diseno de filtros Butterworth pasa-banda esta basado en la aplicacion
de una transformacion sobre un filtro prototipo pasa-bajas. Los requeri-
mientos tpicos son los que se muestran en la figura 4.13. Donde se cumple
con
20 log |H(j)| k2 , para 1
0 20 log |H(j)| k1 , para L U
20 log |H(j)| k2 , para 2
0
k
1
k2
1 U 2
L
Figura 4.13: Especificaciones para un filtro pasa-banda.
La transformacion necesaria es

HP B (s) = HP b (s) s2 +L U , (4.19)
s= s(
U L )
donde PB significa pasa-banda, y Pb significa pasa-bajas. Es decir, se apli-

cara una transformacion analogico-analogico sobre un filtro pasa-bajas para
obtener la funcion de transferencia del filtro pasa-banda.
Para satisfacer el requerimiento de ganancia k 2 en 1 debemos mantener
la igualdad dentro de la transformacion, esto es,
(j21 + L U )
jr =
j1 (U L )
Resolviendo para r tenemos
21 + L U
r =
1 (U L )
y similarmente para satisfacer k2 en 2 , tenemos
22 + L U
r = .
2 (U L )

Ahora bien, dependiendo del valor numerico de 1 y 2 (donde 1 <

2 ) con respecto al producto L U , la r podra ser positiva o negativa.
Tambien en la mayora de los casos las dos r anteriores no seran iguales
y deberemos seleccionar el caso mas restrictivo. As pues, el valor de r
sera aquel que cumpla con
r = mn{|A|, |B|}, (4.20)
donde
21 + L U
A=
1 (U L )
22 L U
B= . (4.21)
2 (U L )
Notese el cambio de signo en el calculo de B.

Esta r sera empleada en el calculo del orden del filtro pasa-bajos pro-
totipo, empleando (4.16) con 1 = 1, y 2 = r .
El procedimiento es el siguiente:
1. Dadas las especificaciones del FPB deseado encontrar A y B.
2. Determinar r .
3. Calcular el orden n del FPb prototipo.
4. Encontrar (o leer de las tablas) la funcion de transferencia del FPb

normalizado.
5. Aplicar la transformacion dada en (4.19).
Ejemplo 4.2.4 Disenar un FPB con las siguientes especificaciones: respues-

ta en frecuencia maximamente plana, con -3 dB a 50 Hz y 2kHz, y con al
menos 20 dB de atenuacion en 20Hz y 4.5kHz
Solucion:
tenemos los siguientes datos
1 = 2(20) = 125.6637 rad/s,

2 = 2(4500) = 2.8274 104 rad/s,
L = 2(50) = 314.1593 rad/s,
U = 2(2000) = 1.2566 104 rad/s.
Calculando A y B tenemos que A = 2.5640 y B = 2.2963; as que r =

2.2963. Entonces,
l log[(100.3 1)/(102 1)] m l m
n= = 2.7667 = 3,
2 log(1/2.2963)

94 4.3. LA APROXIMACION DE CHEBYSHEV
de modo que nuestra funcion prototipo sera

1
HP b (s) = ,
s3 + 2s2 + 2s + 1
aplicando la transformacion (4.19) tenemos
1
HP B (s) = 2 +3.9478106 2 +3.9478106 2 6 .
{[ ss(1.2534910 4) ]
3 + 2[ ss(1.2534910 4) ]
2 + 2[ ss(1.2534910
+3.947810
4 ) ] + 1}
Reduciendo terminos, la solucion es
1.83926132639871 1012 s3
H(s) = ,
s6 + b 5 s5 + b 4 s4 + b 3 s3 + b 2 s2 + b 1 s + b 0
donde
b5 = 2.450442269800038 104
b4 = 3.120768911624457 108
b3 = 2.032740492883779 1012
b2 = 1.232030183398072 1015
b1 = 3.819125666120342 1017
b0 = 6.152890838882032 1019
En la figura 4.14 se muestra la respuesta en frecuencia para la funcion

de transferencia encontrada, puede observarse que se cumple ampliamente
con las especificaciones del filtro.
20
Ganancia en (dB)
40
60
80
100
120 1 2 3 4
10 10 10 10
Frecuencia (Hz)
Figura 4.14: Respuesta en frecuencia del FPB del ejemplo 4.2.4.
4.3. La aproximacion de Chebyshev

Otra funcion de aproximacion esta basada en los polinomios de Chebys-
hev tambien conocidos como Tschebysheff. Al igual que en los filtros de But-
terworth se emplearon los polos estables de un polinomio de Butterworth,
para los filtros Chebyshev se emplearan los polos estables de un polinomio

4.3. LA APROXIMACION DE CHEBYSHEV 95
de Chebyshev. Estos filtros exhiben una cada mas acentuada que los de
Butterworth, a cambio de un rizo que aparece en la banda pasante o en la
banda de rechazo segun se trate de Chebyshev tipo I o Chebyshev tipo II,
respectivamente.
4.3.1. Los polinomios de Chebyshev

Antes de definir la funcion de transferencia de un filtro de Chebyshev,
conviene revisar algunas caractersticas importantes de estos polinomios.
Los polinomios de Chebyshev estan definidos por
(
cos(n cos1 x) si |x| 1
Tn (x) = (4.22)
cosh(n cosh1 x) si |x| > 1
donde
T0 (x) = 1, es decir, n = 0
T1 (x) = x, es decir, n = 1
Tn = Tn (x). (4.23)
Los polinomios de Chebyshev, al igual que las funciones sinusoidales,

son ortogonales, pero con una funcion de ponderacion sobre el intervalo
1 x 1.
Para el caso de las funciones coseno, tenemos

Z 0,
si m 6= n
cos(m) cos(n)d = /2, si m = n 6= 0 (4.24)
0
, si, n = m = 0,
haciendo = cos1 x, se tiene que
1
d = d(cos1 x) = ,
1 x2
cambiando lmites
si 0, entonces x 1
si , entonces x 1,
as pues

Z 1 Tm (x)Tn (x) 0,
si m 6= n
p dx = /2, si m = n 6= 0 (4.25)
1 (1 x2 )
, si n = m = 0
por lo que la funcion de ponderacion es
1
p ,
(1 x2 )

que, como puede verse, crece a medida que x se acerca a 1.

Los polinomios de Chebyshev oscilan entre 1 dentro del intervalo 1 <
x < 1, mientras que fuera de este intervalo crecen rapidamente hacia + o
. Esto puede observarse en la figura 4.15 con T 5 (x), donde se observa,
ademas, que el numero de oscilaciones es igual a n 1.
T (x)
0
n
1
5
1.5 1 0.5 0 0.5 1 1.5
x
Figura 4.15: Polinomio de Chebyshev de orden 5, T5 (x).
Usando la identidad trigonometrica
cos[(n + 1)] + cos[(n 1)] = 2 cos cos n, (4.26)
y cambiando por cos1 x tenemos
cos[(n + 1) cos 1 x] + cos[(n 1) cos 1 x] = 2x cos(n cos1 x)

Tn+1 (x) + Tn1 (x) = 2xTn (x), (4.27)
esto es, tenemos una formula de recursion de tres terminos que permite
encontrar Tn+1 , conociendo los dos anteriores Tn (x) y Tn1 (x). Se obtiene
T0 (x) = 1
T1 (x) = x
T2 (x) = 2x2 1
T3 (x) = 4x3 3x
T4 (x) = 8x4 8x2 + 1 (4.28)
Notese que para n > 1 el primer coeficiente es 2 n1 , y que los polinomios

son alternativamente funciones pares e impares de x. Tambien notese que
T0 (1) = T1 (1) = 1, y
Tn (1) = 1, para toda n.

4.3.2. Los filtros de Chebyshev

La funcion de transferencia de Chebyshev esta dada por
1
|H(j)|2 = . (4.29)
1+ 2 [T n (/c )]
2
Para || c
1
|H(j)|2 = |H()H()| 1.
1 + 2
En exactamente = c ,
1
|H()|2 = ; ya que Tn (1) = 1.
1 + 2
Cuando |x| > 1, esperamos un incremento rapido en T n (x), por lo tanto
el efecto en la funcion de transferencia es hacer que la banda de transicion
sea muy estrecha.
Para || = r = 2 , la frecuencia de rechazo, deseamos una ganancia
1/A2 , es decir,
1 1
= 2, (4.30)
1 + 2 Tn2 (2 /1 ) A
donde 2 = r , y 1 = c . Entonces

A2 1
|Tn (2 /c )| = , (4.31)

y resolviendo para n,
A2 1
cosh1
n= . (4.32)
cosh1 2
1
Notese que en el denominador aparece el termino ( 2 /1 ), mientras que en

el calculo de n para los filtros de Butterworth fue ( 1 /2 ).
Para encontrar la funcion de transferencia H(s) necesitamos encontrar
las races del polinomio
L2 (s) = 1 + 2 Tn2 () = 0.
Recuerdese que Tn () = cosh(n cosh1 ), y haciendo = s/j, entonces
1 + 2 [cosh(n cosh1 s/j)]2 = 0, y despejando tenemos

q
cosh(n cosh1 s/j) = 1/2 = j/ (4.33)
si s = + j es un cero de L2 (s), podemos escribir
cosh[n cosh1 (j + )] = j/, (4.34)

haciendo
u + jv = cosh1 (j + ) (4.35)
j + = cosh(u + jv), (4.36)
de (4.33)
cosh[n(u + jv)] = j/, (4.37)
y de la bien conocida identidad trigonometrica
cos( + ) = cos cos sensen;
podemos escribir (4.37) como
j/ = cos(nv) cosh(nu) + jsenh(nu)sen(nv) = j + . (4.38)
Comparando ambos lados de esta ecuacion,
= sen(nv)senh(nu) (4.39)
= cos(nv) cosh(nu), (4.40)

como j/ tiene parte real cero y parte imaginaria 1/, entonces
cos(nv) cosh(nu) = 0
sen(nv)senh(nu) = 1/.
Como para una nu real, cosh(nu) es diferente de cero, entonces cos(nv) debe
ser cero, y entonces

v= + k /n, k = 1, 2, . . . (4.41)
2

/2+k
(ya que cos(/2) = 0, cos = 0) para sen(nv)senh(nu) = 1/, y
n

como sen(nv) = sen 2 + k /n = 1, entonces
1
u= senh1 (1/) . (4.42)
n
De (4.39) (4.40) (4.41) y (4.42), se obtiene
1 1 2k + 1
k = senh senh1 sen ,
n 2n
1 1 2k + 1
k = cosh senh1 cos , k = 1, 2, . . . , (n 1),
n 2n
y podemos ver que
k2 2k
+ = 1,
senh2 (u) cosh2 (u)

es decir, los ceros de L2 (s) estan localizados sobre una elipse. Los polos de
L2 (s) estan dados por
1 1 2k + 1 1 1 2k + 1
pk = c senh senh1 sen + jc cosh senh1 cos .
n 2n n 2n
Finalmente, la funcion de transferencia normalizada se puede formar
como
K K
Hn (s) = Qn = , (4.43)
i=1 (s p i ) V n (s)
donde
Vn (s) = sn + bn1 sn1 + . . . + b1 s + b0 , (4.44)
y ademas se puede comprobar que
(
b0 , si n es impar
K= b0 ,
(4.45)
1+2
si n es par.
En las tablas 4.3, 4.4, y 4.5 se muestran los polinomios para filtros
Chebyshev tipo I para diferentes ordenes y niveles de rizo aceptable en la
banda pasante.

Tabla 4.3: Polinomios de Chebyshev de orden n, para un rizo de 0.5 dB ( =

0.34931140018895).
orden Kn y Vn (s) orden Kn y Vn (s)
n=1 Kn = 2.86277516124319 n=2 Kn = 1.43138758062160
b0 = 2.86277516124319 b0 = 1.51620262694593
b1 = 1.42562451364020
n=3 Kn =0.71569379031080 n=4 Kn = 0.35784689515540
b0 = 0.71569379031080 b0 = 0.37905065673648
b1 = 1.53489545855561 b1 = 1.02545527713700
b2 = 1.25291297268055 b2 = 1.71686620544884
b3 = 1.19738565671119
n=5 Kn = 0.17892344757770 n=6 Kn = 0.08946172378885
b0 =0.17892344757770 b0 = 0.09476266418412
b1 = 0.75251811034289 b1 = 0.43236692045929
b2 =1.30957474478165 b2 = 1.17186133283420
b3 = 1.93736749474810 b3 = 1.58976350135366
b4 =1.17249093365202 b4 = 2.17184462271973
b5 = 1.15917610630976
n=7 Kn = 0.04473086189442 n=8 Kn = 0.02236543094721
b0 = 0.04473086189442 b0 = 0.02369066604603
b1 = 0.28207222652099 b1 = 0.15254443887580
b2 = 0.75565110403040 b2 = 0.57356040094292
b3 = 1.64790292616448 b3 = 1.14858936912993
b4 = 1.86940790812050 b4 = 2.18401538172037
b5 = 2.41265095657871 b5 = 2.14921726034592
b6 = 1.15121757854778 b6 = 2.65674980658908
b7 = 1.14608010766183

Tabla 4.4: Polinomios de Chebyshev de orden n, para un rizo de 1 dB ( =

0.50884713990959).
n=1 Kn =1.96522672836027 n=2 Kn = 0.98261336418014
b0 = 1.96522672836027 b0 = 1.10251032805385
b1 = 1.09773432856393
n=3 Kn = 0.49130668209007 n=4 Kn = 0.24565334104503
b0 = 0.49130668209007 b0 = 0.27562758201346
b1 = 1.23840917357824 b1 = 0.74261937310676
b2 = 0.98834120988476 b2 = 1.45392476228002
b3 = 0.95281137931914
n=5 Kn =0.12282667052252 n=6 Kn = 0.06141333526126
b0 = 0.12282667052252 b0 = 0.06890689550337
b1 = 0.58053415132206 b1 = 0.30708063841820
b2 = 0.97439607307168 b2 = 0.93934552954141
b3 = 1.68881597917823 b3 = 1.20214038896508
b4 = 0.93682013127199 b4 = 1.93082492260129
b5 = 0.92825096024867
n=7 Kn = 0.03070666763063 n=8 Kn = 0.01535333381531
b0 = 0.03070666763063 b0 = 0.01722672387584
b1 = 0.21367139021183 b1 = 0.10734472587535
b2 = 0.54861981077436 b2 = 0.44782572369701
b3 = 1.35754480295084 b3 = 0.84682432066278
b4 =1.42879430819559 b4 = 1.83690238444740
b5 = 2.17607847362704 b5 = 1.65515567030966
b6 = 0.92312347346067 b6 = 2.42302641920325
b7 = 0.91981130587012

Tabla 4.5: Polinomios de Chebyshev de orden n, para un rizo de 2 dB ( =

0.76478310157921).
n=1 Kn = 1.30756027157908 n=2 Kn = 0.65378013578954
b0 = 1.30756027157908 b0 = 0.82306042667169
b1 = 0.80381643012779
n=3 Kn =0.32689006789477 n=4 Kn = 0.16344503394738
b0 = 0.32689006789477 b0 = 0.20576510666792
b1 = 1.02219033985978 b1 = 0.51679810179633
b2 = 0.73782157715775 b2 = 1.25648193319480
b3 = 0.71621495822804
n=5 Kn =0.08172251697369 n=6 Kn = 0.04086125848685
b0 = 0.08172251697369 b0 = 0.05144127666698
b1 = 0.45934912106488 b1 = 0.21027055620056
b2 = 0.69347695849584 b2 = 0.77146177104735
b3 = 1.49954326711260 b3 = 0.86701492190888
b4 = 0.70646056806109 b4 = 1.74585874591746
b5 = 0.70122570676988
n=7 Kn =0.02043062924342 n=8 Kn = 0.01021531462171
b0 = 0.02043062924342 b0 = 0.01286031916675
b1 = 0.16612634957256 b1 = 0.07293732194127
b2 = 0.38263807763868 b2 = 0.35870427523376
b3 = 1.14459656515260 b3 = 0.59822138558167
b4 =1.03954580250709 b4 = 1.57958072373214
b5 = 1.99366531716017 b5 = 1.21171207766297
b6 = 0.69809070637012 b6 = 2.24225292847038
b7 = 0.69606454940671

Ejemplo 4.3.1 Disenar un filtro pasa-bajas Chebyshev normalizado con

rizo de 2dB y una atenuacion de 20 dB o mayor a una frecuencia de 1.5
rad/s en adelante.
Solucion:
De las figuras 4.16 y 4.17 observamos que en = c = 1 se cumple la
2
1/(1+ )
2
1/(A )
1 r
Figura 4.16: Respuesta en frecuencia pasa-bajas Chebyshev n impar.
1/(1+2)
1/(A2)
1 r
Figura 4.17: Respuesta en frecuencia pasa-bajas Chebyshev n par.
siguiente relacion:
1/2
1
20 log |H(j1)| = 20 log
1 + 2

1
= 10 log = 2,
1 + 2
mientras que en = r
1 1/2

20 log |H(j1.5)| = 20 log 2

A
10 log(1/A2 ) = 20.

Resolviendo para en (4.46) obtenemos un valor para de 0.76478, y para

A en (4.46) obtenemos A = 10. Estos dos valores son necesarios en el calculo
del orden n del filtro, ya que, aplicando (4.32) tenemos

l cosh1 A2 1 m

n =
cosh1 2
1

l cosh1 0.76478
99 m
=
cosh1 1.3
1
l m
= 4.3 = 5
Entonces requerimos de un filtro prototipo de Chebyshev de quinto orden

cuya funcion de transferencia esta dada por:
0.081722
H5 (s) = .
s5 + 0.70646s4 + 1.49954s3 + 0.69347s2 + 0.45934s + 0.08172
La respuesta en frecuencia de H5 (s) se muestra en la figura 4.18. Puede
observarse que se cumple el requerimiento de -2 dB a la frecuencia de corte
de 1 rad/s.
10
10
Ganancia en (dB)
20
30
40
50
60
0
10
Frecuencia (Hz)
Figura 4.18: Respuesta en frecuencia Chebyshev del ejemplo 4.3.1.

4.4. CONSTRUCCION DE TABLAS DE DISENO DE FILTROS
NORMALIZADOS 105
4.4. Construccion de tablas de diseno de filtros
normalizados
Antes de continuar con algunos ejemplos de diseno de filtros Chebyshev,
conviene introducir una forma de construir tablas de diseno que pueden ser
muy utiles cuando no se tiene una tabla completa. Para esto haremos uso
de el ambiente de programacion Matlab. Con esto se podran obtener los
valores de filtros normalizados o no, tipo Butterworth, Chebyshev, Elpticos,
entre otros.
4.4.1. Tabla de filtros de Butterworth usando Matlab

Los polinomios de Butterworth de la tabla 4.1 se pueden obtener usando
la siguiente instruccion:
[num,den]=butter(n,1,s)
donde n es el orden del filtro, 1 = c y s indica analogico. As pues, para

diferentes valores de n tenemos los resultados que se muestran en la tabla
4.6.
4.4.2. Tabla de filtros de Chebyshev usando Matlab

Los polinomios de Chebyshev son obtenidos usando la siguiente instruc-
cion:
[num,den]=cheby1(n,e,1,s),
donde n es el orden del filtro, e es el rizo permisible en la banda de paso,

1 = c y s indica analogico. As pues, para diferentes valores de n tenemos
las funciones de transferencia que se muestran en las tablas 4.7 a 4.9.

106 NORMALIZADOS
Tabla 4.6: Construccion de HBn (s)

orden n
Resultado de Matlab
n=1

1
num = 0 1
s+1
den = 1 1
n=2
)
num = 0 0 1 1

den = 1 1.4142 1 s2 + 2s + 1
n=3
)
num = 0 0 0 1 1
den = 1 2.0 2.0 1 s3 + 2s2 + 2s + 1
n=4
)
num = 0 0 0 0 1
HB4 (s)
den = 1 2.6131 3.4142 2.6131 1
1
=
s4 + 2.6131s3 + 3.4142s2 + 2.6131s + 1
..
n=5 .

NORMALIZADOS 107
Tabla 4.7: Construccion de HChn (s) con rizo = 0.5 dB.

orden n
Resultados de Matlab
n=1
)
num = 0 2.8628 2.8628
den = 1 2.8628 s + 2.8628
n=2
)
num = 0 0 1.4314 1.4314
den = 1 1.4256 1.5162 s2 + 1.4256s + 1.5162
n=3
)
num = 0 0 0 0.71569379
den = 1 1.2529129 1.5348954 0.71569379
0.7156
=
s3 + 1.25291s2 + 1.5348s + 0.7156
n=4
)
num = 0 0 0 0 0.35784689
den = 1 1.197385 1.7168662 1.0254552 0.379050
0.3578
=
s4 + 1.1973s3 + 1.7168s2 + 1.0254s + 0.379
..
n=5 .

108 NORMALIZADOS
Tabla 4.8: Construccion de HChn (s) con rizo = 1 dB.

orden n
n=1
)
num = 0 1.96522 1.96522
den = 1 1.96522 s + 1.96522
n=2
)
num = 0 0 0.98261 0.98261
den = 1 1.09773 1.102510 s2 + 1.09773s + 1.10251
n=3
)
num = 0 0 0 0.4913066
den = 1 0.9883412 1.238409 0.4913066
0.4913066
=
s3 + 0.9883412s2 + 1.238409s + 0.4913
n=4
)
num = 0 0 0 0 0.2456533
den = 1 0.9528113 1.453924 0.7426193 0.275627
0.2456533
=
s4 + 0.9528113s3 + 1.453924s2 + 0.742619s + 0.27562
..
n=5 .

NORMALIZADOS 109
Tabla 4.9: Construccion de HChn (s) con rizo = 2 dB.

orden n
n=1
)
num = 0 1.307560 1.30756
den = 1 1.307560 s + 1.30756
n=2
)
num = 0 0 0.65378 0.65378
den = 1 0.803816 0.823060 s2 + 0.803816s + 0.82306
n=3
)
num = 0 0 0 0.326890
den = 1 0.737821 1.0221903 0.32689
0.326890
=
s3 + 0.737821s2 + 1.0221903s + 0.32689
n=4
)
num = 0 0 0 0 0.1634450
den = 1 0.7162149 1.2564819 0.516798 0.2057651
0.1634450
=
s4 + 0.7162149s3 + 1.256481s2 + 0.5167s + 0.2057651
..
n=5 .

110 NORMALIZADOS
Ejemplo 4.4.1 Disenar un filtro pasa-bajas Chebyshev que tenga un rizo
de -2 dB en la banda pasante y una frecuencia de corte de 60 rad/s, con una
atenuacion de 20 dB al menos en 80 rad/s.
Solucion:
El procedimiento general es cambiar las especificaciones a un filtro (pro-
totipo) pasa-bajas normalizado. Para esto calculamos la frecuencia crtica
2
cr = 1
= 80/60 = 1.333. Luego disenamos un filtro Chebyshev norma-
lizado con 2 dB de rizo en 1 rad/s, y 20 dB de atenuacion en 1.33 rad/s.
Entonces obtenemos n = 5, as que de H C5 (s) del ejemplo 4.3.1 solo necesi-
tamos aplicar una transformacion pasa-bajas pasa-bajas con s = s/60.
En la figura 4.19 se muestra la respuesta en frecuencia para este caso.
10
5
Ganancia en (dB)
10
15
20
25
30
35
40 0 1 2
10 10 10
Frecuencia (rad/s)
Figura 4.19: Respuesta en frecuencia Chebyshev del ejemplo 4.4.1.
Ejemplo 4.4.2 Disenar un filtro pasa banda (FPB) Chebyshev con las si-
guientes especificaciones: -2dB de rizo en la banda de frecuencias de 200 Hz
a 3.5 kHz, y -30 dB o mayor a frecuencias menores de 50 Hz y mayores de
10kHz.
Solucion:
Comenzamos encontrando la frecuencia crtica para determinar el orden del
filtro prototipo.
1 = 2(50),
L = 2(200),
U = 2(3500),
2 = 2(10000),
21 + L U
A =
1 (U L )

4.5. FILTROS ELIPTICOS 111
= 4.2424
22 L U
B =
2 (U L )
= 3.009.
As
pues, cr = mn{|A|, |B|} = 3.009 y = 0.76478 mientras que A r =
10k2/20 = 31.6277, y

l cosh1 ( 31.622 1/0.76478) m
n =
cosh1 (3.009)
l m
= 2.4998 = 3.
La funcion prototipo sera

0.3268
Hn (s) = .
s3 + 0.7378s2 + 1.022s + 0.3268
Finalmente, aplicamos la transformacion PbPB con

H(s) = Hn (s) s2 +L U
s s(
U L )
2.913957075 1012 s3
= ,
s6 + b 5 s5 + b 4 s4 + b 3 s3 + b 2 s2 + b 1 s + b 0
donde
b0 = 2.11044155 1022 ,
b1 = 1.1683170 1019 ,
b2 = 1.44354926 1016 ,
b3 = 3.75949469 1012 ,
b4 = 5.223647 108 ,
b5 = 1.5298369 104 .
En la figura 4.20 se muestra la respuesta en frecuencia para el filtro

pasa-banda del ejemplo 4.4.2.
4.5. Filtros elpticos

Una clase de filtros que presentan una cada mas pronunciada en la
banda de transicion son los filtros elpticos. Estos presentan rizo tanto en la
banda pasante como en la banda de rechazo, y un pequeno aumento en la
complejidad del calculo de la funcion de transferencia.

112 4.5. FILTROS ELIPTICOS
10
10
Ganancia en (dB)
20
30
40
50
60
2 3
10 10
Frecuencia (Hz)
Figura 4.20: Respuesta en frecuencia Chebyshev pasa-banda del ejemplo 4.4.2.
Los filtros elpticos se basan en las propiedades de la funcion elptica,

la cual fue analizada primero por Jacobi en 1829 y que ha sido estudiada
extensivamente en algunos tratados de matematicas. Para nuestros proposi-
tos es importante reconocer que la funcion elptica jacobiana, denotada por
sn(u), es una funcion doblemente periodica de la variable compleja u, es
decir, periodica en Re{u}, y en Im{u}, y analtica en el plano-u, excepto en
los polos simples de la funcion. El patron basico de dos ceros y dos polos se
repite infinitamente a lo largo de los dos ejes. El desarrollo matematico de
los filtros elpticos es algo mas complicado que los polinomios de Chebyshev.
Un filtro elptico se basa en una transformacion de una funcion
1
T 2 (w) = , (4.46)
1+ 2 sn2 (w, k)
desde el plano complejo-w al plano complejo-s, donde
Z z=sn(w,k) dt
w = sn1 (z, k) = p (4.47)
0 (1 t2 )(1 k 2 t2 )
es una integral elptica de primera clase, z = sn(w, k) es una funcion elptica,

k es el modulo y 0 < k < 1. Puede verse, a partir del patron de ceros y polos
que para valores reales de w, T 2 (w) tiene un maximo en w = 2mK (m es
un entero) y un mnimo en w = (2m + 2)K (este en el comportamiento que
se desea en la banda de paso) donde K es una constante dada por:
Z 1 dt
K = K(k) = p .
0 (1 t2 )(1 k 2 t2 )
Es importante notar que a lo largo de Im(w) = K 0 , T 2 (w) = 0 en w =

2mK = jK 0 y T 2 (w) tiene un maximo en w = (2m + 1)K = jK 0 ; y ese es
el comportamiento en amplitud deseado en la banda de rechazo. Para mas

detalles sobre las propiedades de las funciones elpticas se puede consultar

la bibliografa indicada al final del captulo.
El filtro elptico normalizado tiene una respuesta en frecuencia dada por
1
|Hn (j)|2 = , (4.48)
1 + 2 Rn2 ()
donde Rn () es la funcion elptica. Como en los filtros Chebyshev, tene-

mos dos casos, uno para n par y el otro para n impar. De igual manera es
conveniente emplear filtros normalizados como una base para obtener filtros
pasa-bajas, pasa-altas, pasa-banda, y rechazo de banda.
Para el filtro elptico, = 1 es la media geometrica de 1 y 2 , es decir,
(1 2 )1/2 = 1.
El parametro que representa lo pronunciado de la banda de transicion esta de-

finido por la razon
r = 2 /1 .
La funcion de transferencia para el filtro elptico pasa-bajas normalizado

esta dada como sigue:
(n1)/2
H0 Y s2 + A0i
Hn (s) = , n impar, (4.49)
(s + s0 ) i=1
s2 + B1i + B0i
n/2
Y s2 + A0i
Hn (s) = H0 , n par. (4.50)
i=1
s2 + B1i + B0i
El filtro queda especificado determinando el orden n, H 0 , el polo simple
s0 , y los coeficientes A0i , B1i B0i . Estos parametros son determinados a
partir de las especificaciones de diseno, , A, y r , o por su equivalente Rp ,
Rr y r , donde

1
Rp = 20 log = 20 log |Hn (j1 )|,
(1 + 2 )1/2
Rs = 20 log(1/A2 ) = 20 log |Hn (j2 )|,
Rp = rizo en la banda de paso en dB,
Rr = rizo en la banda de rechazo en dB.
4.5.1. Diseno de filtros elpticos analogicos en Matlab

Para el diseno de filtros elpticos se hara uso del ambiente de computo
Matlab donde la orden de entrada tiene la siguiente forma:
[num,den] = ellip(N,Rp,Rs,Wn,s),

114 4.5. FILTROS ELIPTICOS
donde,
num = es el numerador, un polinomio en potencias decrecientes de s,
den = es el denominador, un polinomio en potencias decrecientes de s,
N = orden del filtro,
Rp = rizo en la banda de paso,
Rs = riso en la banda de rechazo,
Wn = frecuencia natural del filtro.
Si Wn es un vector de dos elementos Wn=[W1 W2], entonces se obtiene un
filtro pasa-banda de orden 2N, con una banda de paso W1 < W < W2.
Para obtener un diseno pasa-altas se emplea:
[num,den] = ellip(N,Rp,Rs,Wn,high,s).
Si, nuevamente, Wn es un vector de dos elementos, entonces se obtiene un

filtro rechazo-de banda con
[num,den] = ellip(N,Rp,Rs,Wn,stop).
El orden N del filtro y la frecuencia natural Wn, se obtienen con:
[N, Wn] = ellipord(Wp, Ws, Rp, Rs, s),
donde
Wp = la frecuencia (final) de la banda de paso y
Ws = la frecuencia (inicial) de la banda de rechazo.
Tambien se puede obtener el resultado en forma de ceros, polos y ganan-
cia, empleando
[Z,P,K] = ellip(...),
donde los ceros estaran en un vector Z de N columnas, y los polos en un

vector P de N columnas. La ganancia K es un escalar. Para otras formas en
el resultado ver la gua del usuario de Matlab.
Ejemplo 4.5.1 Disenar un filtro pasa-bajas elptico con las siguientes es-
pecificaciones: frecuencia de corte de 1000 Hz con rizo de 3 dB, frecuencia
de rechazo a 2000 Hz con -60 dB de atenuacion.
Solucion:
Bajo el ambiente de Matlab definimos los siguientes valores y obtenemos
el orden del filtro introduciendo las siguientes lneas:
Wp = 2*pi*1000;
Ws = 2*pi*2000;
Rp = 3;
Rs = 60;
[N,Wn] = ellipord(Wp,Ws,Rp,Rs,s);
obteniendo el resultado siguiente:

N = 5
Wn = 6.2832e+003.
Con este resultado introducimos:
[num,den] = ellip(N,Rp,Rs,Wn,s)
Ahora obtenemos la respuesta en frecuencia con:
[H,w] = freqs(num,den);\\
Mag = abs(H);
En lugar de emplear la funcion plot emplearemos la funcion semilogx para

obtener una mejor presentacion en la grafica de la respuesta en frecuencia.
Ademas, calculamos la ganancia en dB haciendo 20*log10( ). As pues
introducimos la siguiente orden
semilogx(w/(2*pi), 20*log10(Mag))
Si deseamos que aparezcan las unidades en los ejes, entonces introducimos:
xlabel(Frecuencia (Hz))
ylabel(Ganancia en (dB))
Finalmente, colocamos una rejilla para una mejor lectura de las ganancias
y frecuencias en los puntos importantes, con:
grid
La figura 4.21 muestra la grafica de la respuesta en frecuencia, y puede

observarse que se cumple con las especificaciones de diseno.
10
10
20
Ganancia en (dB)
30
40
50
60
70
80
90
100 1 2 3 4 5
10 10 10 10 10
Frecuencia (Hz)
Figura 4.21: Respuesta en frecuencia del ejemplo 4.5.1.

116 BIBLIOGRAFIA
Bibliografa
[1] Gabel, R. A., Roberts, R. A., Senales y Sistemas Lineales, LIMU-
SA,1975.
[2] Haykin, S., Van Veen, B., Senales y Sistema, LIMUSA Wiley, 2001.
[3] Mitra, S. K., Digital Signal Processing: a computer based approach,

McGraw-Hill, 1998.
[4] Neff H. P. Jr., Continuous and Discrete Linear Systems, Harper & Row
Publishers, 1984.
[5] Oppenheim, A. V., Young, I. T., Senales y Sistemas, Prentice Hall

Hispanoamericana, 1994.
[6] Papoulis A., Signal Analysis, McGraw-Hill ISE, 1977.
[7] El-Sharkawy M., Digital Signal Processing Applications with the Moto-
rolas DSP56002 Processor, Prentice Hall PTR, 1996.
Problemas
Problema 4.1 La funcion del sistema H n (s) representa un filtro de But-
terworth normalizado a 1 rad/s de orden-n. Para n = 5,
a) Escriba H5 (s) en forma polinomial y en forma factorizada.
b) Cual es la ganancia |H5 (s)| en = 1?, cual es la ganancia en dB?
c) Repita los incisos anteriores 4.1a) y 4.1b) para un filtro Chebyshev

tipo I, con = 0.7647831.

Problema 4.2 Dado que G(s) = 1/(s2 + 2s + 1) representa un filtro
Butterworth normalizado de segundo orden,
a) Dibuje 20 log |G(j)| para desde 0 a 100. A que frecuencia, en

rad/s, tiene una magnitud de -3 dB?, a cual tiene -20 dB?
b) Aplique una transformacion PbPA, s 10/s, a la G(s) para

obtener un nuevo filtro H1 (s), y trace la grafica de 20 log |H1 (j)|
para de 0 a 100. Tiene el nuevo filtro el comportamiento que se
esperaba? a que frecuencia en rad/s tiene la magnitud 3dB abajo?
cual esta 20 dB abajo?
c) Aplique la transformacion s 5s/(s 2 + 50) a G(s), para obtener un

nuevo filtro H2 (s). Que tipo de filtro resulta y cuales son las frecuen-
cias de interes? Dibuje 20 log |H2 (j)| para verificar sus conclusiones.

BIBLIOGRAFIA 117
Problema 4.3 Compruebe que Tn (1) = (1)n .
Problema 4.4 Calcule T5 (x) y T6 (x).
Problema 4.5 Demuestre que de

A2 1
|Tn (2 /c )| = ,

resolviendo para n, se obtiene

A2 1
cosh1
n= .
cosh1 2
1
Problema 4.6 Disene a) un filtro Butterworth, y b) un filtro Chebyshev

pasa-bajas analogicos que tengan una atenuacion de 3 dB a una frecuencia
de corte de 100 rad/s y una atenuacion de 25 dB o mayor para frecuencias
arriba de 250 rad/s. Dibuje 20 log |H(j)| para los dos filtros y muestre
que se satisfacen los requerimientos a las frecuencias crticas.
Problema 4.7 Cual es el orden n de un filtro Chebyshev pasa-bajas que

tiene una banda de paso desde 0 a 200 Hz con un rizo aceptable de 1dB y una
banda de rechazo monotonica con -40 dB en y mas alla de 250 Hz? Repita
el diseno para un filtro Butterworth y compare las ns. A que conclusiones
puede llegar?
Problema 4.8 Se desea disenar un filtro pasa-bajas analogico que tenga

-0.5 dB a una frecuencia de corte de 75 Hz y que tenga mas de 20 dB de
atenuacion a frecuencias mayores a 150 Hz. Encuentre H(s) que satisfaga
esos requerimientos y dibuje 20 log |H(j)| y arg H(j) para
a) Filtro Butterworth (respuesta maximamente plana).
b) Filtro Chebyshev tipo I (equirizo en la banda de paso).
c) Filtro elptico.
Compare el orden de cada filtro y comente sus resultados.
Problema 4.9 Disene un filtro analogico pasa-banda que satisfaga las si-
guientes especificaciones:
1. Frecuencias de paso inferior y superior de 100 Hz y 3.8 kHz, respecti-

vamente, con -3 dB de atenuacion.
2. Atenuacion en la banda de rechazo de 20 dB a 20Hz y 8 kHz.
3. Sin rizo en la banda de paso ni en la banda de rechazo.

118 BIBLIOGRAFIA
Compruebe sus resultados dibujando 20 log |H(j)| y arg H(j).
Problema 4.10 Repita el Problema 4.9 para -3 dB a 50 Hz y 20 kHz, y

una atenuacion de 20 dB en 10 Hz y 60 kHz.
Problema 4.11 Disene a) un filtro Butterworth, y b) un filtro Chebyshev

pasa-altas que permitan el paso de senales con frecuencias mayores a 100Hz
con una atenuacion no mayor a 2 dB y tenga una atenuacion mayor de 20
dB a frecuencias menores de 20 Hz.
Problema 4.12 Se encontro que un filtro tiene la siguiente funcion de

transferencia:
s
H(s) =
s+1
a) Cual es el orden del filtro?
b) Trace la grafica de la respuesta en frecuencia en magnitud. Tome solo
tres puntos de frecuencia: en = 0.1 rad/s, 1 rad/s, y 10 rad/s.
c) Que tipo de filtro es este?: FPb (filtro pasa-bajas), FPA (filtro pasa-
altas), o FPB (filtro pasa-banda). Podra decirse que es un filtro normali-
zado?
Problema 4.13 Disene un filtro pasa bajas Chebyshev tipo 1 con las si-
guientes especificaciones:
a) Frecuencia de corte de 100 Hz, con -2 dB de rizo.
b) Frecuencia de rechazo de 1500 Hz con -30 dB de atenuacion.
Problema 4.14 Se tiene un filtro analogico Butterworth normalizado de

segundo orden, es decir,
1
Ha (s) = .
s2 + 2s + 1
Realizar una transformacion analogico-analogico para obtener un filtro analogi-
co Butterwoth con respuesta en frecuencia Pasa Altas, con una frecuencia
de corte de 500 Hz.
Problema 4.15 Se tiene un filtro cuya funcion de transferencia esta dada

por
s
H(s) =
s + 10
a)Encuentre la ganancia del filtro, |H(j)| (en magnitud), y la fase theta(j),
si = 10 rad/s.
b) Cual es la frecuencia de corte de -3dB.

Diseño de Filtros Analogos - Del Libro PDS

Cargado por

Copyright:

Formatos disponibles

Diseño de Filtros Analogos - Del Libro PDS

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Diseño de Filtros Analogos - Del Libro PDS

Cargado por

Copyright:

Formatos disponibles

Diseno de filtros en tiempo continuo

donde c es la frecuencia de corte, es fsicamente no realizable, o no causal.

4.1. Especificaciones de un filtro

banda banda banda de

Figura 4.1: Especificaciones de un filtro.

4.2. La aproximacion de Butterworth

donde n es el orden de la funcion de transferencia correspondiente y c

Gn () = 20 log |H(j)| = 10 log |H(j)|2n

donde log es el logaritmo base 10.

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Para << c , tenemos que

Figura 4.2: Respuesta en frecuencia para filtros Butterworth, orden 1 a 3.

Para encontrar la funcion de transferencia decimos que = s/j, ya que

donde Sk es la k-esima raiz o polo del denominador, con magnitud unitaria

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Tabla 4.1: Algunos polinomios de Butterworth.

Ejemplo 4.2.1 Encontrar la funcion de transferencia H 1 (s) para un filtro

Figura 4.3: Races de la funcion de Butterworth para n=1.

Ejemplo 4.2.2 Encontrar la funcion de transferencia H 2 (s) para el filtro

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Desarrollando para n=2 y cada valor de k se obtiene

Figura 4.4: Races de la funcion de Butterworth para n=2.

En las figuras 4.5 y 4.6 se muestran los polos (o races) en el plano-s

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Figura 4.5: Races de la funcion de Butterworth para n=6.

Figura 4.6: Races de la funcion de Butterworth para n=7.

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Tabla 4.2: Polinomios de Butterworth de orden n, Bn (s) = bn sn + bn1 sn1 +

Para el diseno de filtros de Butterworth haremos uso de tablas que con-

4.2.1. Transformaciones analogico-analogico

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

vez necesitaremos un filtro precisamente con esa frecuencia de corte. En esta

donde H 0 (s) es la nueva funcion de transferencia, obtenida a partir de H(s)

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Si se desea, esto podra considerarse como una doble transformacion, una

donde el termino (s2 + 2 ) nos da los polos complejos conjugados, mientras

Figura 4.7: Posicion del polo pasa-bajas Butterworth de primer orden.

Figura 4.8: Transformacion del eje real a 0 .

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Figura 4.9: Respuesta pasa-banda por desplazamiento del eje real.

De lo anterior se puede deducir que el FPB resultante tendra un orden

y no como 0 = (L + U )/2. Para filtros con un ancho de banda estrecho

4.2.2. Diseno de filtros Butterworth pasa-bajas

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

De la figura 4.1, se tienen las siguientes condiciones sobre la respuesta

0 20 log |H(j)| k1 , para toda 1 y

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

y entonces aplicaremos la transformacion s = s/ esc .

para aplicar la transformacion.

entonces, tomamos la funcion prototipo

ahora aplicamos una transformacion PbPb, con

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

Figura 4.10: Respuesta en frecuencia del filtro del ejemplo 4.2.3.

Figura 4.11: Ampliacion en k1 = 2 dB.

Figura 4.12: Ampliacion en k2 < 20 dB.

G. Miramontes, ISBN 968-5923-15-9 PDS: Introduccion con teora y practica

4.2.3. Diseno de filtros Butteworth pasa-banda

Figura 4.13: Especificaciones para un filtro pasa-banda.

1 + 2 [cosh(n cosh1 s/j)]2 = 0, y despejando tenemos

cosh[n cosh1 (j + )] = j/, (4.34)

j/ = cos(nv) cosh(nu) + jsenh(nu)sen(nv) = j + . (4.38)

Tabla 4.3: Polinomios de Chebyshev de orden n, para un rizo de 0.5 dB ( =

Tabla 4.4: Polinomios de Chebyshev de orden n, para un rizo de 1 dB ( =

Tabla 4.5: Polinomios de Chebyshev de orden n, para un rizo de 2 dB ( =

Resolviendo para en (4.46) obtenemos un valor para de 0.76478, y para