T

to

X(e ) k=
x(t)y (t)dt =
x
[n mod N ] X(f )Y (f )df
X [k]
n x[n]
j d(z)
(x y)[n]
X(z)Y
contiene Rx Ry
d
j
para
[n]
1 cualquier
z
x[n]
X(e to . )
Im[x(t)] = 0 X(f ) = X (f )
Transformada Z
1
an u[n] ! t +T
|z|
>
|a|

1
[n]
1
1az
"
x(t) = x (t) Im[X(f )] = 0
1
+
1

+
n
x(t)y
(t)dt
=
X[k]Y
[k]
!
1] 1az1
|z| < |a|
1
[n no ] a u[n
ejno
T t
n (f )
Re[x(t)]
==
0 X(f
) = X
no $ k=
X(z)
x[n]z
x(n)
=
X(z)z n1 dz
[n no ]
z
z excepto 0 o
+

2j
1
2(
+
2k)
1
C
x(t)
=
x
(t)

Re[X(f
)]
=
0
n=
1z k=
cos o
cos(o n)u[n]
1 X[k]
Im[x(t)] =12z
0
=2
X [k]|z| > 1
cos1o +z
an u[n] (|a| < 1)
1
Funci
on
Transformada
z1ae
sinj
o
donde
C
es
una
curva
antihoraria
en la region de convergencia y
$ |z| > 1
sin(o n)u[n]
Formulario
TDS
1 cos

du[n]
1 o +z 2= 0
x(t) = x12z
(t)
Im[X[k]]
aX(f ) + bY (f )
+ k= ( + 2k) ax(t) + by(t)
1ej
que
envuelve
al origen.
Formulario
(x y)(t)
X(f )Y (f )
1 etrica
Se
Procesamiento
Digital
de
nales
Geom
(n + 1)an u[n]
(|a|
< 1)= 0 Serie
Re[x(t)]

X[k]
=
X
[k]
j
2
(1ae Digital
)
x(t)y(t) Secuencia
(X Y )(f )Transformada Z
Procesamiento
de Se
nales
ROC
N2
nte.
N1
+1
r n sin o (n+1)
!
1

N2Re[X[k]]
u[n]
(|r|
1)
k <
x
(t)
X (f )
x(t)
=
x
(t)
=
0
j
2 ej2
sin o
12r
cos

e
+r
=
sio N2 >
N1 y += 1
%
&
ax[n]
+
by[n]
aX(z)
+
bY
(z)
contiene
R x Ry
1 $

X(t)
x(f )
sin o n
k=N1
+2k
n
o
j2f
t
k
n
2
x[n

n
]
z
X(z)
R
,
quiz
a
0o
o
x(t td )
o X(f )e
x
Funcion
Transformada
'
M
j t n
M
+1 + bYM[k]
M
x(t)e
X(f

f
)
1 0 X
n +M
ax(t)
by(t)X
aX[k]
+1
cX(z/zo )
zo x[n]
|zo |Rx
1 sin((M+1)/2)
1 ejM/2
n
x[n] =
X(f
f )e ,+X(f +f )e
(1)
sin(/2)
(t
) =
ej2k /T X[k]
n =
(1)
x(t) cos(kc t + )
d k
0 otrofj2mt/T
n
2
)
*
n
x[n]
z
X(z)
R
1m]
x , quiza 0 o
$X[k
f
dz
1
e
x(t) n=0 1
n=0
x(t)
#
|| X

ejo n
2(

+
2k)
o
1 T
k=
X (z )
Rx
1
X[k]Y [k]
d x(t) x [n]
X
T10 x(t )y()d
X
(j2f )n X(f )
1 n =X[k] 1Y [k], |a| < 1
# dt
x(t)y(t)n
t
x[n]
X(1/z)
1/R
X(0)

(2)
1
x
%
&
= de Poisson
, 1 |a| < 1
(2)
$ F
ormula
tnT
j2f X(f ) + 2 (f )
x()d
1]z
1

T sinc(k/T )
d X(f ) (z)
n=0
n
nX(z)Y
(x

y)[n]
contiene
R x Ry
t x(t)
(j2)
n=0
df
(
)
+
+
j2f t
[no ]
!
!
(t tCuadro
)
e
1
z
d[n]
2
2
Identidades
trigonom
etricas
1: Transformada
z de algunas funciones comunes
jkt
tk
n
ej( t+)
ej (f fc ) 1
ej =e cos +=j sin

a
u[n]
|z| > |a|

1
1az
1 (f
B)k=
= sen(A) cos(B) cos(A) sen(B)e(at) n
(3)/a) 1
k=
cos sen(A
= e +e
a e$
%
&
$
2
Se
n
al
x(n)
Transformada
z,

a u[n
1] 1az
|z|X(z)
< |a| ROC
1 (3)
sen(A sin
B)
1k
e=
esen(A) cos(B) cos(A) sen(B)
k= f T
cos(A
B) = cos(A) cos(B) sen(A) sen(B) k= (t kT ) T1 (4)
=
2j
n
o
Desigualdad
de
Schwarz
] jf(4) z
z exceptoPlano
0 o z
1
= sin(
+ 90 )
cos(A cos
B)
= cos(A)
cos(B)
sen(A) sen(B) sign t [n no(n)
1
1
1
2
1z 1 cos o
u(t)
+
(f
)
cos
(A)
=
(1
+
cos(2A))
(5)
j2f
2
sin = cos(
|z| > 1
*2 90 ) IA EL ECTRICA "
% & cos(o n)u[n]
c **"
ver 2.8 !
2003
INSTITUTO DE INGENIER
12z 1 cos o +z 2
1
b
b
t
sinc f
2 2 + cos*2 1= 1" b 2
sin
z 1 sin o
% t & sin( n)u[n]
*
*
cos
=
(1
+
cos(2A))
(5)
u(n)
|z| > 1
(A)
2
2
2
|z| > 1
1
o

sinc f 12z 1 cos o1+z 2
v()w
|w()| d
*
* 22
cos2()d
= 1/2(1
cos 2)
z 1
sen
=|v()|
(1 d
cos(2A))
(6)
* a
* +(A)
a
a
3
cos = 1/4(3
1cos + cos23)
1
Cola
Gaussiana Serie Geometrica
2
1
sensin2 (A)
= (1cos 2)cos(2A))
an u(n)(6)
|z| > |a|
= 1/2(1
sen(A)
sen(B)
=
(cos(A
B)
cos(A
+
B))
(7)
1
az 1
2
3
N
2
La igualdadsinseda
si v()
K 3)
= 1/4(3
sin =
sin
N1
N2 +1
!
2 w() con K constante.

k
sin( ) = 1
sin cos 1cos sin
azsi 1N2 > N1 y += 1
=
n (7)
sen(A) sen(B)
=
(cos(A
B)
cos(A
+
B))
na
u(n)
|z| > |a|
1
cos(A)
cos(B)
=
(cos(A
B)
+
cos(A
+
B))
(8)
Transformada
cos(
) =2cos cos Zunilateral
sin sin
k=N
(1 az 1 )2
1
tan tan 2
tan( ) =11tan
tan 1
1
+
!
sinsen(A)
sin =
=cos(B)
1/2(cos(A
cos(
1/2
cos( +
cos(A) cos(B)
B)
+ B))+ B))
= )(sen(A
B))+ sen(A
(9)
(an )u( n(8)
1)
|z| < |a|
n + cos(A
2
x[n]z
21/2=cos(
1
az 1
cos cosX(z)
=

)
+
1/2
cos(
+
)
r
n=0 ) +
sin cos = 1
1/2Asin(
1 1/2 sin( + )
az 1
sen(A) cos(B)sen
= (sen(A
B) +cos(A))
sen(A + B))
=
(1
(10)1)
n(an )u( n(9)
|z| < |a|
2
2
2
(1 az 1 )2
r

X
(z)
r
x[n]
u

Algunas funciones utiles

A n1 A1
1
1 z 1 cos !0
no(1 + cos(A))
o=e
cos
=
(11)

p(x)
=
Distribuci
on
de
Gauss
cos(!
n)u(n)
|z| > 1
sen
(1
cos(A))
(10)
0
x[n + n
]

z
X
(z)

x[0]z

x[n

1]z
o
u
o
2
1 2z 1 cos !0 + z 2
!2 22
2 2
1

/2
n
=z n
o re j! d

Q(k)
Cola
Gaussiana
o +1
cos(!)
+
j sen(!)
(12)
z 1 sen !0
x[n no ]
x[1]z
+ + x[no ]
A sin2tXku (z)e1+ =
sen(!
n)u(n)
|z| > 1
0
j!
cos
=
(1 + j!
cos(A))
(11)
1 2z 1 cos !0 + z 2
sinc t =
Sinc
2 't cos(!)
2 = e +e
(13)
1
tF>ormula
0
1Signo
1 az 1 cos !0
2
sign t = j!
an cos(!0 n)u(n)
|z| > |a|
j!
e' 1
= cos(!)
+ej!
j sen(!)
(12)
t<0
e
1 2az 1 cos !0 + a2 z 2
j
sen(!)
=
(14)
1 t >e0j! + 1 =j!0
u(t) =
Escal
2jon
az 1 sen !0
0 t <e0 + e
' =
cos(!)
(13)
an sen(!0 n)u(n)
|z| > |a|
%t&
1 |t| Hoja
< 2 2de formulas MPD/ SISCOM ver 2.8 !
1IA EL
2az 1 cos ! + a2 z 2
c k > 3 laINSTITUTO
DE INGENIER
=

Rectangulo
Cuando 2003
siguiente es una
buena aproximaci
on ECTRICA 0

DTFT

Formulario

X(ej ) =

+
!

( 0

|t| j!
> 2
|t|
|t|

% &
1=
sen(!)
t =

<
|t|2j
>

X1 (ej )X2 (ej )d

Im

x[

Re[

x[

Funcion
a1 x[n] + a2 y[n]
(x y)[n]
x[n]y[n]
x [n]
x[n no ]
x[n]ejno

Transformada
a1 X(ej ) + a2 Y (ej )
X(ej )Y (ej )
#
1
j
j(
))d
2 X(e )Y (e
X (ej )
X(ej )ejno
X(ej(o ) )

n x[n]
x[n]
[n]
[n no ]
1
an u[n] (|a| < 1)
u[n]
(n + 1)an u[n] (|a| < 1)

)
j d X(e
d
X(ej )
1
ejno
$+

r n sin o (n+1)
u[n]
sin o
sin o n
n

x[n] =
e

jo n

'

F
ax[n] + by[n
X[n]
x[(n m) m
WNln x[n]
$N 1
m=0 x[m]
x[n]y[n]
x [n]
x [n mod

k= 2( + 2k)
1
1aej $
1
+
k= (
1ej
1
(1aej )2
1
12r cos
% o ej +r
& 2 ej2
$
+2k
k
2o

(|r| < 1)

1 0nM
0 otro n

X(z) =

+ 2k)

+
!

c 2003 INSTITUTO DE INGENIER IA EL ECTRICA

Hoja de formulas MPD/SISCOM ver 2.8 !

ejkt =

k=

2( o + 2k)

(
)
+
2 !
2
tk

k=

Desigualdad de Schwarz
*"
*

*2
*

"

"

1
b

donde C es u

Secuenc
ax[n] + by
x[n no ]
zon x[n]

Formula de Poisson

2
1
Q(k) (14)ek /2
2k

+
!

n=

sin((M+1)/2) jM/2
e
sin(/2)
k=

j!

Triangulo

x[n] = x [n] Re[X(ej )] = 0

(xm)2
22

X[k] =

X(ej )ejn d

Re[x[n]] = 0 X(ej ) = X (ej )

x[n] = x [n] Im[X(e )] = 0

"

1
=
2

"

1
2

Im[x[n]] = 0 X(ej ) = X (ej )

x1 [n]x2 [n]

n=

x[n] =

x[n]ejn

n=

+
!

nk x[n]
x [n]
x[n]
(x y)[n]
[n]
an u[n]
an u[n
[n no ]
cos(o n)u

3.4.4.Comparacin entre los distintos mtodos de diseo.

Histricamente,
el primero
en aparecer,
losdiseo.
otros dos
desarrollaron
en la de
dcada
3.4.4.Comparacinfue
entre
los distintos
mtodos
de
3.4.4.Comparacin
entreselos
distintos mtodos
diseo.

Filtros FIR. Resumen

Cuadro 2: Propiedades de la transformada z

MTODO
DE70.
LAS
VENTANAS
de Dominio
los
ROC
Dominio n
z

Propiedad

Cuadro 3: Propiedades de la DFT

Propiedad
Dominio temporal
MTODOX(z)
DE LAS VENTANASROC: r2 <
MTODO
DElos
LAS
VENTANAS
|z|
< r1
Histricamente,
fue adecuada
el primerolas
en
aparecer,
otros
dos
desarrollaron
en dependen
la dcada del
No fija de manera
frecuencias
crticas
! se
y!
ya que stas

Conjugacion

Dominio frecuencial
p
s
Notacion
x(n), y(n)
X(k), Y (k)
ROC1
X1 (z)
Histricamente, fue el primero
en aparecer,
los otros dos
se desarrollaron
enaparecer,
la dcada
Periodicidad
x(n)otros
= x(ndos
+ N )se desarrollaron
X(k) = X(k + Nen
) la dcada

Histricamente,
fue
el
primero
en
los
ROC
2
x2 (n)
X2de
(z)
de
los
70.
Linealidad
a1 x1 (n) + a2 x2 (n) a1 X1 (k) + a2 X2 (k)
tipo
ventana y de la longitud seleccionada.
de los
70.
Reflexion temporal
x(N n)
X(N k)
por lo menos ROCde
los 270.
1 \ROC
a1 x1 (n) + a2 x2 (n)
1 X1 (z) + a2 X2 (z)
No afija
de manera adecuada las frecuencias crticas Desplazamiento
!p y !s ya que
stas
dependen
dell))N
temporal
circular
x((n
X(k) e j2kl/N
como
de X(z)
z =de
0 si
No fija
de manera adecuada
laslafrecuencias
!
!s adecuada
ya que stas
Desplazamiento
frecuencial
circular del
x(n)ej2ln/N
X((k stas
l))N dependen del
p y
excepto
Nocrticas
fija
manera
las dependen
frecuencias
crticas
!p y !s ya que
x(n k)
z k X(z)
k > 0 y z = 1 si k < 0
Conjugaci
o
n
compleja
x(n)
X(N k)
tipo
de
ventana
y
de
la
longitud
seleccionada.
MUESTREO
EN FRECUENCIA seleccionada.
Convolucion circular
x1 (n) N x2 (n)
X1 (k)X2 (k)
tipo de 1ventana y de la longitud
|a|r2 < |z| < |a|r1 tipo de ventana y de la longitud seleccionada.
n
a x(n)
X(a z)
Correlacion circular
x(n) N y( n)
X(k)Y (k)
1
Controlamos perfectamente
la anchura
de la zona Multiplicaci
de transicin,
que es igualx1a(n)x
2"/N.
on de dosya
secuencias
X (k) N X2 (k)
2 (n)
1
1
N 1
< |z| <
N
N 1
x( n)
X(z 1 )
X1
1 X
r2
MUESTREO
EN
FRECUENCIAr1
Teorema
de
Parseval
x(n)y(n)
MUESTREO
EN
FRECUENCIA
MUESTREO
ENde
FRECUENCIA
Cuadro 3:
Propiedades
de la DFT
Hay
procedimientos
rpidos
el clculo
los coeficientes, bien basndose
en la N k=0 X(k)Y (k)
ROC para
n=0
x(n)
X(z)

Parte realPropiedad

Re{x(n)}

Notacion

x(n)
x1 (n)

Ventanas:

Linealidad
Desplazamiento en n
Escalado en z

Muestreo en Frecuencia:

Reflexion en n
3

Controlamos
perfectamente
anchura
de
la zona
deperfectamente
transicin,
ya
que
es aigual
a 2"/N.

Controlamos
perfectamente
lalaanchura
la zona
de transicin,
ya que
igual
2"/N.
1
deControlamos
laes
anchura
de
la zona
de transicin,
ya que es igual a 2"/N.
Incluye
ROC
Dominio
temporal
Dominio
frecuencial
FFT,
bien
con
las
ecuaciones
propuestas
anteriormente.
Especialmente
interesante
si la
X(z)
+ X(z)
2
Cuadro
4: Simetr
asVentanas.
en filtros FIR de fase lineal
Diseo
por

Notacion
y(n)
X(k),
Ypara
(k)
Incluye
ROC
Hay
procedimientos
rpidos
para
de
coeficientes,
bien
basndose
Hay1 x(n),
procedimientos
rpidos
el el
clculo
de
los los
coeficientes,
bien
basndose
en ladeen
clculo
Hay
procedimientos
rpidos
para
el
clculo
loslacoeficientes, bien basndose en la
Parte imaginaria
Im{x(n)}
X(z)
X(z)
mayor
parte
de
los
puntos
el
mdulo
de
la
ganancia
son
cero
o
uno.
2 = x(n + N )
Periodicidad
x(n)
X(k) = X(k + N )
Simetra Simetrica h(n) = h(M 1 n)
Antisimetrica h(n) = h(M 1 n)
dX(z)
|z| <
r1
a1 x1z(n)
+ acon
a1 X1 (k) r+2 <
apropuestas
FFT,
bien
las
ecuaciones
anteriormente.
Especialmente
interesante
si la si la Especialmente interesante si la
DerivacionLinealidad
en z
nx(n)
FFT,
bien
con
las
ecuaciones
propuestas
anteriormente.
Especialmente
interesante
FFT,
bien
con
las
ecuaciones
propuestas
anteriormente.
2 x2 (n)
2X
2 (k)
1
X
dz
Como
sePor
tiene
un pobre
control Mdeparla respuesta
fuera
de esos puntos
y el
Reflexion temporal
x(Ninconveniente,
n)
X(N
k)menos
Hr (0) = 2
h(k)
Hr (0) = 0
lo
ROC1 \ROC2
Convoluci
o
n
x
(n)

x
(n)
X
(z)X
(z)
j2kl/N
1
2
1
2
mayor
parte
de
los
puntos
el
mdulo
de
la
ganancia
son
cero
o
uno.
mayor
parte
de
los
puntos
el
mdulo
de
la
ganancia
son
cero o uno.
Desplazamiento temporal circular mayor
x((nparte
l))Nde los puntos
X(k) e el mdulo de la ganancia son cero o uno.
k=0
j2ln/N
no
apto
como
filtro paso bajos
procedimiento
puede
convertirse
en
un
esquema
de
prueba
y
error.
Por
lo
menos
la
intersecci
o
n
de
Desplazamiento frecuencial circular
x(n)e
X((k l))N
1
Como
inconveniente,
tiene
un
pobre
control
de
la
respuesta
fuera
de
esos
puntos
y
el
Correlacion
rx x (n) = x1 (n) x2
( Como
n)
Rx x inconveniente,
= X1 (z)X2 (z 1 ) se

Como
inconveniente,
se
tiene
un
pobre
control
de
la
respuesta
fuera de esos puntos y el
ROC
y
la
de
X
(z
)
1
2
se
tiene
un
pobre
control
de
la
respuesta
fuera
de
esos
puntos
y
el
Conjugacion compleja
x(n)
X(N k)

X
M
2

1 2

1 2

M 3
2

M 1
xx(0)
x2X(z)
(n)puede convertirse
X1 (k)X2 (k) en un esquema
M impar
Hconvertirse
+2
h(k) de
Hr (0)
= Hr ()
=0
1 (n)= Nl
r (0) = h
procedimiento
de
prueba
y
error.
m
procedimiento
puede
prueba
y error.
procedimiento
convertirse
en un esquema de prueba y error. 2en un esquema
k=0
x(n)INz!1
y( n) puede
X(k)Y
(k)
no apto como filtro paso bajos o altos
1x (n)x (n) z 1 1X (k)Por
N lo
X2menos
(k) r1l r2l < |z| < r1u r2u
1
X21 (v)X2
vN dv1
2j
v
N
X1IC
X1
1 N
1
X
1 x(n)y(n)DE1 TCHEBYSHEV
APROXIMACIN
Teorema
Relacion de
Parsevalde Parsevalx1 (n)x2 (n) =
X1 (v)X2
vN1 dv X(k)Y (k)
Cuadro 5: Funciones utilizadas como ventanas.
2j
v
n=0 C
k=0

Convoluci
on circular
Teorema del
valor inicial
Si x(n) es causal
Correlacion circular
Multiplicaci
on de
dos2 (n)
secuencias
Multiplicaci
on
x1 (n)x

Aproximacin de TChebyshev:
n= 1

APROXIMACIN
TCHEBYSHEV
TCHEBYSHEV
Ventana
h(n),a0 frecuencias,
n M 1 Ancho
Permite unDE
control
total de APROXIMACIN
las caractersticasDEdel
filtro en cuanto

Pico lobulo

APROXIMACIN DE TCHEBYSHEV
lobular lateral [dB]
en
Permite
unydecontrol
total de las caractersticas
filtro entotal
cuanto
a frecuencias,
Permite del
un control
de las
caractersticas del filtro en cuanto a frecuencias,
Cuadro 4: Simetras
filtros FIR
fase lineal
ganancias
longitud.
Rectangular
1
4/M
13
M 1
Permite
un
control
total
de
las
caractersticas
del
filtro
en
cuanto
a
frecuencias,
2
n
2
y longitud.
Simetrica h(n) = h(M 1 ganancias
n)
Antisim
etrica h(n) = h(M 1 n) ganancias y longitud.
(triangular)
1
8/M
27
M del
1 filtro,
No existe una forma fcil de optimizar el diseo Bartlett
respecto
a la longitud
2n
ganancias
y
longitud.
1
Hamming
0,54Aproximacin
0,46
cos
8/M a la 32
No existe una forma fcil de optimizar
el existe
diseo una
respecto
a fcil
la longitud
del
filtro,
X
No
forma
de Orden
optimizar
el diseo
respecto
longitud del filtro,
Estimacin
del
M 1 Tchebysheb

Hr (0) = 2
h(k)
Hrexisten
(0) = 0 aproximaciones como la de Kaiser:
aunque
1
2n
Filtrosla Equiripple.
No
existe
una aproximaciones
forma fcil decomo
optimizar
el diseo
respecto
Hanning a la longitud
1 del
cos filtro,
8/M
43
aunque
existen
la de
Kaiser:
k=0
aunque
existen aproximaciones
como
2 de Kaiser:
M 1

Simetra

M
2

M par

no apto como filtro paso bajos

M impar

Hr (0) = h

1
2

aunque existen aproximaciones como la de Kaiser:

Estimacin del Orden.

Aprox. Kaiser

M 3
2

+2

h(k)

Hr (0) = Hr () = 0

k=0

no apto como filtro paso bajos o altos

Marcelino Martnez Sober

Tratamiento Digital de Seales GIET

Cuadro 5: Funciones utilizadas como ventanas.

lunes 5 de marzo de 12

Ventana

h(n), 0 n M

Ancho
lobular

Pico lobulo
lateral [dB]

Filtros Analgicos. Revisin

Filtros de Butterworth

Cuadro 7: Caractersticas de Filtros Paso Bajo Analogicos

log #!1 / + " % 1$

1
Filtros Analgicos.
Revisin
Filtros
Analgicos. Revisin
H (( ) &
N &
2

Filtro

Funcion de Transferencia

Caractersticas

Butterworth
(todo-polos)

|H()|2 =

1
1 + (/C )2N
C : frecuencia de corte

|H()| monotona en las bandas de paso

y rechazo.

|H()|2 =

1 + 2

1
2
TN
2
TN

|H()|2 =

1
1+

2 UN

'e

j ( 2 k )1)

2N

Filtros de Chebyshev I

k & 0 ,1,..., N % 1

0 1 1 7 2 2 1 1 7 2 3 (1 2 6 2 ) 2
2
) 1 &
log10 .pasante
+
$$
ch 11 ''
Respuesta plana en banda
y 6atenuada.
3 7 22
.
+
1
I
NGENIERA
E
LECTRNICA
2
/
,
( 7 26 %
Estimacin
Orden
Solo
Polos.
0
1
6
H
(
)
N
#
4
5
5
4
*
Rizado constante en la banda de paso
y
Tchebyshevmontona decreciente.
1 2 6 2 3 T 2 !# / # c " Respuesta
ch 11 !# s / # c "
0#
2
caracterstica monotona en la banda Nde
log10 . s 22 !# s / # c " 11
1+
Polos distribuidos en una circunferencia
# c) 4 de radio unidad
1 1 71 4 1 1 2 6 2
H /(#
rechazo.
4.5.- Transformaciones
en frecuencia. ,

1 2 6 2 3 T 2 !# / # c "

0 5 6 51

"

"

x 51

150

x =1

200 2
10

UN : funcion elptica Jacobiana

de orden N

N=2
N=3
N=4
N=5

8Presentan una cada monotnica en la banda

corte
8Menor orden
que Butterworth

0.9
0.8

1
2
0.7 T (x), 0el polinomio
que se define como:
10
10
10 de Tchebyshev
10
Rizado constante
tanto en la banda de
N Frecuencia
TN+1(x)=23x3TN(x)-TN-1(x) con T0(x)=1 y T1(x)=x.
0.6Filtro de Butterworth c=1
paso como en la de rechazo.
<9cos!N 3 cos 11 x "
x 51
1
P. Baja ! P.0.5Alta
" nueva frecuencia de
>TN(x)>51
? @x@<1
TN ( x) 4 ;
11
0.8
TN(1)=1
? N,
x =1
9:cosh !N 3 cosh x "
0.4
corte
Todas las races de TN(x) estn en xAB-1,+1C
1

1
0.6
Todas las curvas pasan por
0.3
BN (s)
N=1
0.4
N=2
Respuesta de fase lineal en la banda de
(x)=23x3TN(x)-TN-1(x) con T0(x)=1 y T1(x)=x.
BN (s): funciones de Bessel
0.2 TN+1
N=1
N=3
N=2
0.21
1
paso (aunque se destruye en la
2 converN=4
N=3
H (# 4 # c ) 4
4
B0 (s) = 1
0.1
N=5
sion a digital).
1 2 6 2 3 TN2 !1" 1 20 622 P. Baja !1 P. Banda N=4
?10@x@<1
N=5 0 >TN(x)>51
1
2
B1 (s) = s + 1
10
10
10
10
0
Frecuencia
10
10
10
10
10
T
N(1)=1
?
N,
2
"
frecuencia
inferior
de
la
Frecuencia
BN (s) = (2N 1)BN 1 (s) + s BN 2 (s)
l

Si aplicamos la transformacin
(todoBilineal a las Transformaciones
en frecuencia analgicas
obtenemos las T. en
Tratamiento
Digital
Seales GIET
Tratamiento
Digital de
Sealesde
GIET
Frecuencia en
el dominio
Z.
Cuadro 8: Transformaciones de frecuencia para filtros analogicos.
|H()|

5 dedemarzo
de 12
lunes 5lunes
de marzo
12

Tipo de filtro
deseado

Transformacion

Paso bajo
Paso alto
Paso banda
Supresor de banda
P :
l :
u :

P
s ! 60 s
P
P 0P
s!
s 2
s + l u
s ! P
s(u l )
s(u l )
s ! P 2
s + l u

lunes 5 de marzo de 12

Constantes

H (# 4 # c )

1
1 2 6 2 3 TN2 !1"

u , l

lunes 5 de marzo de 12

u , l

"u frecuencia superior de la

1
12 62

1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2

N=1
N=2
N=3
N=4
N=5

0.1
0 2
10

P. Baja ! Elimina Banda

Tratamiento Digital de Seales GIET
"l frecuencia inferior de la

0P

"u frecuencia superior de la

filtro pasa-baja

0P

Todas las races de TN(x) estn en xAB-1,+1C

banda Marcelino
TodasSober
las curvas pasan por
Marcelino Martnez
Sober Martnez

" frecuencia de corte del

10

Marcelino Martnez Sober

banda

banda

frecuencia de corte del prototipo

frecuencia de corte inferior
frecuencia de corte superior

Transformacion

banda

Nuevas
frecuencias
de corte

" frecuencia de corte del

Digital
de Seales
Cuadro 9: Transformaciones deTratamiento
frecuencia para
filtros
digitales.GIET
Tipo de filtro

0 11 7 2 2
2
log10 .
./
N4
0#
log10 . s 2
/ #c

N
8 Son filtros solo polos
=1
8 Presentan rizado constanteFiltro
endelaButterworth
bandapasante
2
c
1 7 1 4no
1 pasante.
12 6
0
8Presentan
una cada monotnica en la1banda
8Menor orden que Butterworth

<9cos N 3 cos 11 x
TN ( x) 4 ;
11
9:cosh N 3 cosh x

H(s) =

Bessel
polos)

Rizado constante en la banda de recha50

Las transformaciones
digitales
son:
Transf.
en frecuencia
digitales 8 Son filtros solo polos
Filtro de Chebyshev I =1
zo y caracterstica
mon
otona en de
la banTN(x),
el polinomio
Tchebyshev que se define como: 1
8 Presentan rizado constante en la banda pas
100
P. Baja ! P. Baja
" nueva frecuencia de
N=1
da de paso.

TN : polinomio de Chebyshev

Elptico
o Cauer (polos
y ceros)

Polos : s k & ( c ' e

2 ' log 10 !( s / ( c "

Transformaciones en frecuencia en el dominio Z

TN : polinomio de Chebyshev
T0 (x) = 1
T1 (x) = x
TN +1 (x) = 2xTN (x) TN 1 (x)

Chebyshev, Tipo II (polos y

ceros)

Filtros de *Chebyshev
* I

|H()|

Chebyshev,
Tipo I (todopolos)

1 + 2 TN2

1 ) !( / ( c "

|H()|

|H()| =

10

2N

|H()|(dB)

Estimacin Orden Butterworth

Nuevas

Marcelino Martnez Sober

32

k=0
N

lineal de las salidas del sistema todo ceros (ver Proakis).

(k) z k

1 + aN (k) z k
k =1

Desc. Celosa

C (z)
= M
AN (z)

Cm 1 (z) = Cm (z) vm Bm (z) siendo vm = cm (m)

CELOSIA IIR ceros y polos

= =

Recursin Descendente
Obtenemos los coeficientes de reflexin

Estructuras en celosa. (lattice-ladder)

= +

Recursin Ascendente

Calculamos Am(z) y finalmente AM(z)=H(z)

= =

CELOSA IIR con ceros y polos

La salida total del sistema se puede obtener como una combinacin

Estructuras en celosa.
(lattice-ladder)
lineal de
las salidas
del sistema
todo ceros (ver Proakis).
Estructuras
en
celosa.
(lattice)
c (k) z

H (z) =

k) z k
=
(k) z

C M (z)
AN (z)

CELOSA FIR

k=0
N

Tratamiento Digital de Seales GIET

1 consideramos que es una celosa todo-polos y calculamos los

C M (z)

Estructura multietapa
lunes 5 de marzo de 12

k
Dado H (z) = h(k)z
k=0

Definimos el conjunto de filtros:

do polos

=1

(k)z k

Marcelino Martnez Sober

Tratamiento Digital de Seales GIET

Marcelino Martnez Sober

Estructuras en celosa. (lattice)

lunes 5 de marzo de 12

=
coeficientes de reflexin.
robustez
Presentan una gran
AN (z) numrica.

C 1 (z) = Cm (z)
vm Bm (z) siendo
vm = cC
(m)
k
mM
La salida total
del
sistema
se
puede
obtener
como una combinacin
1
+
a
(k)

z
2 Calculamos losmcoeficientes
de la escalera,
vm utilizando
(z)

Modularidad
para
la
implementacin.
N

ksistema
=1
lineal
de
las
salidas
del
todo
ceros
(ver
Proakis).

Constituyen
un
mtodo
para
determinar
la
estabilidad
de
un
sistema
IIR
(Schur-Cohn
Test)

de estructuras no DIRECTAS

Se trata

C (z) = Cm (z) vm Bm (z) siendo vm = cm (m)

Sistemas MA(FIR), AR(IIR solo polos), ARMA(IIR con ceros y polos)
Distinguimosm 13 casos:

con ceros y polos

H (z) =

Etapa genrica

1
=
AN (z)

si intercambiamos y(n) y x(n)

tenemos un sistema FIR y
podemos calcular Km.

y(n) = x(n) + N (k)x(n k)

k =1

Am (z) = m (k)z k
k=0

H (z) = AM (z)

m 1
m (0) = 1
m (m) = K m

Km= COEFICIENTES DE REFLEXIN

(k)y(n
k) + x(n) que
1
consideramos
que es una celosa
todo-polos y calculamos los

1 consideramos
es una celosa
todo-polos
y calculamos
los

coeficientes de reflexin.
coeficientes de reflexin.
CELOSIAdeIIR
Todo
Polosvmde
2 Calculamos los
coeficientes
laEcuaciones
escalera,
utilizando
CM(z)
2recurrencia
Calculamos
los coeficientes de la escalera, vm utilizando CM(z)

f0 (n) = g0 (n) = x(n)

fm (n) = fm 1 (n) + km gm 1 (n 1)
Marcelino Martnez Sober
gm (n) =Digital
km fde
(n) + gm 1 (n 1)
m 1
Tratamiento
Seales GIET
lunes 5 de fmarzo
de =
12 y(n)
M (n)

recurrencia

tal de Seales GIET

n) km gm 1 (n 1)
1 (n) + gm 1 (n 1)
g0 (n)

Cmo obtener los

coeficientes de reflexin a
partir de H(z) ?

Marcelino Martnez Sober

Tratamiento Digital de Seales GIET

lunes 5 de marzo de 12

Marcelino Martnez Sober

<

El sistema es