ESCUELA POLITCNICA DEL EJRCITO
EXTENSIN LATACUNGA
Escudero Vsconez Miguel Dario
Ruiz Osorio Cristian Mauricio
Ingeniera Electrnica e Instrumentacin, Quinto, Escuela Politcnica del Ejrcito Extensin
Latacunga, Mrquez de Maenza S/N Latacunga, Ecuador.
email :[Link]@[Link]
ruiz.cristian92@[Link]
Fecha de presentacin: 12 de junio del 2014
ESTRUCTURA EN CELOSA
RESUMEN
En el mbito del procesamiento digital de seales es comn utilizar un tipo de estructura de filtro
en el cual se pueda analizar y modelar la seal en base a sus caractersticas, estas prestaciones
nos brinda la representacin celosa, este filtro analiza los coeficientes de reflexin y el anlisis
hacia adelante y hacia atrs para analizar la seal travs de l. Estas estructuras son muy
complejas entenderlas debido a su extenso uso y propiedades que la conforman es por ellos que
en el presente documento presentamos otra estructura de filtro FIR denominada realizacin en
celosa. El uso de los filtros en celosa est muy extendido en las aplicaciones de tratamiento de
voz y en la implementacin de filtros adaptativos.
ABSTRACT
In the field of digital signal processing is common to use a type of filter structure in which to analyze
and model the signal based on its features, these features gives us the lattice representation, this
filter analyzes the reflection coefficients and analysis forward and backward to analyze the signal
there through. These structures are very complex to understand because of its widespread use and
properties that make it is for them that in this paper we present another FIR filter structure called
lattice realization. The uses of lattice filters are widely used in voice processing applications and the
implementation of adaptive filters.
PALABRAS CLAVE
Celosa
FIR
Filtro
DESARROLLO
1. Generalidad de Celosa
polinomio.
1.1.
matemtica, definimos
Ventajas
Nmero reducido de coeficientes
permite que grandes bloques de
datos puedan ser modelados en
tiempo real.
1.2.
al
Am ( z ) y { y (n)}
filtro
es
la
secuencia de salida, tenemos:
Procesamiento
digital
de
seales de voz.
Filtros adaptativos.
Tratamiento
de
seales
geofsicas
y (n)=x (n)+ m (k) x (nk )
k=1
Suponga ahora que tenemos un
filtro de orden m= 1. La salida de tal
filtro es
2. Desarrollo
Considerando una secuencia de filtros FIR
cuyas funciones de sistema son:
y (n)=x (n)+ai(l)x (n1)
Esta salida tambin puede obtenerse
a partir del filtro en celosa de primer
orden o de una sola etapa, mostrado
en la Figura 1, excitando ambas
entradas con x(n) y seleccionando la
salida de la rama superior.
As,
la
salida
es
exactamente
y (n)=x (n)+ai(l)x (n1)
Hm ( z ) =Am ( z ) m=0,1,2. . , M 1
donde, por definicin, Am (z) es el polinomio
m
si
seleccionamos
parmetro
Am ( z ) =1+ m (k ) z m 1
k
K1
k = 1 (1) .
El
de la celosa se
denomina coeficiente de reflexin.
k=1
Ao ( z)=1.
m ( 0 )=1 .
Usos comunes
Adems
conveniencia
Si {x(n)} es la secuencia de entrada
Desventajas
El filtro de estructura celosa puede
ser FIR o IIR por tanto tambin
posee sus desventajas.
1.3.
Por
La respuesta al
impulso unitario del filtro m es
hm(0)=1 y hm(k)=am(k ), k=1, 2,... , m
A continuacin, considere un filtro
FIR para el que m = 2. En este caso,
la salida de una estructura en forma
directa es
y (n)=x (n)+a 2(l) x (n1)+ a 2(2) x (n2)
.
El
subndice
Am (z)
indica
del
el
polinomio
grado
del
Conectando en cascada dos etapas
en celosa como se muestra en la
Figura 2, es posible obtener la
misma salida. Por tanto, la salida de
la primera etapa es:
fi (n)=x (n)+ K 1 x (nl)
gl(n)=K 1 x (n)+ x (nl)
La salida de la segunda etapa es
f 2 ( n ) =f 1 ( n ) + K 2 g 1 ( nl )
g2 ( n )=K 2 f 1 ( n )+ g 1 ( nl )
Si centramos nuestra atencin en f2 (n) y
sustituimos f/ (n) y g/(n 1) obtenemos
K x (nl)+ x (n2)}
f 2(n)=x( n)+ K x (nl)+ K 2
Figura 2. Filtro en celosa de dos
solas etapas
2.1. ESTRUCTURA EN CELOSA Y
CELOSA ESCALONADA PARA SISTEMAS
IIR
Un sistema FIR tiene una funcin de
H ( z )= A N ( z )
transferencia
mientras
que
un
H ( z )=1 / A N (z) ,
sistema
por
lo
IIR
que
intercambiando la entrada por la
f 2(n)=x( n)+ K 1(1+ K 2) x (n1)+ K 2 x( n2)salida se puede obtener un sistema
a partir del otro. Por ejemplo un
sistema IIR todo polos de primer
Ahora la Ecuacin es idntica a la
orden se puede representar a travs
salida de la realizacin de la forma
de:
directa del filtro FIR dada por , si
f N ( n ) =x ( n )
igualamos los coeficientes, es decir,
2 ( 2 )=K 2, 2 ( 1 )=K 1(1+ K 2)
o, lo que es equivalente:
f m1 ( n ) =f m ( n )k m gm 1 ( n1 ) M =N , N1, , 1
gm ( n )=k m f ( n ) + g m1 ( n1 ) M =N , N 1, , 1
m1
y ( n )=f 0 ( n )=g0 (n)
Y un sistema FIR de primer orden:
g1 ( n )=k m y ( n ) + y (n1)
g0(n) = x(n)
Figura 1. Filtro en celosa de una sola
etapa
Los polos son el resultado de la
realimentacin introducida por la
solucin
de
{ f m (n)}
en
orden
descendente. Mientras que los
coeficientes de un sistema FIR son
iguales a los de un sistema IIR pero
en orden inverso; cabe recalcar que
la estructura en celosa todo polo
tiene una trayectoria todo cero con
g0 (n)
entrada
y salida
gn ( n ) ,
trayectoria idntica a la todo ceros
de una estructura en celosa.
Figura 4. Estructura en celosa
escalonada de un sistema polosceros
Figura 3 Estructura en celosa para un
sistema IIR todo polos
Los parmetros K1, K2,, KN son
los mismos para estructuras todo
polos y todo ceros, se diferencian
solo por la interconexin de sus
grafos; si
|K m|<1
para todo m las
estructuras en celosa todo polos
son estables. En la practica a sido
utilizado en las modelaciones del
tranco
vocal
humano
y
la
estratificacin de la tierra.
Un sistema todo ceros tiene como
salida una combinacin lineal de
salidas retardadas de un sistema
todo polos; una estructura en
celosa todo polos con parmetros y
aadiendo una parte escalonada
tomando
como
salida
una
combinacin lineal
gm (n)
y ( n )= v m g m (n)
m=0
representa los parmetros
que nos sirven para determinar los
ceros del sistema, la funcin de
transferencia es :
M
H ( z )=
g (z)
Y (z)
= V m m X ( z ) y F 0 ( z )=G0 (z )
X ( z) m=0
X ( z)
M
G m ( z ) F0 ( z ) M
B (z)
H ( z )= V m
= V m m
A N ( z)
G0 ( z ) F N (z) m=0
m=0
M
V m Bm ( z )
H ( z )=
m=0
A N (z )
Que da como resultado:
es una
combinacin
lineal
de
salidas
presentes y pasadas; la salida del
sistema es:
M
{ vm}
C m ( z )= V m Bm ( z )
m =0
Los
coeficientes
del
polinomio
C M ( z) , sirven para determinar los
coeficientes de ponderacin de la
escalera
{ vm}
del polinomio
y los coeficientes
A N ( z)
determinan
los parmetros de la celosa
{k m }
Los parmetros de la escalera estn
dados por:
en cascada, donde el filtro describe
el conjunto de ecuaciones siguiente:
m 1
Cm ( z )= V k Bk ( z ) +V m Bm ( z )=C m1 ( z ) +V m B m ( F
z )0 ( z ) G0 ( z ) X ( z )
k=0
Estos
parmetros
se
calculan
recursivamente a partir de los
polinomios
m=1,2,..,M,
inversos
como
B m (z)
m (m)
para
todo m, los parmetros se pueden
determinar
mirando
que
V m=C m (z )
m=0,1,,M;
dando
Fm ( z ) Fm 1 ( z ) K m z 1 Gm 1 ( z ),
m=1,2,,M 1
Gm ( z ) K m Fm 1 ( z ) z 1G m 1( z )
m=1,2,,M 1
como resultado :
Cm 1 ( z )=C m1 ( z ) +v m Bm ( z )
Estos filtros en celosa escalonada
requieren un mnimo de memoria
aunque no pocas multiplicaciones,
otra de sus ventajas es que son
filtros muy estables y robustos ante
los efectos de palabras de longitud
finita; por lo que son muy utilizadas
en aplicaciones prcticas tales como
procesamiento de voz, filtrado
adaptivo
y
procesamiento
de
seales geofsicas.
-
Ejercicio.
Obtenga los coeficientes de la
celosa correspondiente al filtro
FIR
con
funcin
de
transferencia.
3
1
1
A( z ) 1 z 1 z 2 z 3
4
2
4
Resolucin
La estructura en celosa se usa
ampliamente en procesado digital
de la voz y en la realizacin de
filtros adaptativos. Un sistema en
celosa presenta una serie de etapas
Donde Km es el parmetro de celosa
de la etapa
m-sima,
tambin
denominados
coeficientes
de
reflexin por ser idnticos a los
coeficientes
de
reflexin
introducidos
en
el
test
de
estabilidad de Schr-Cohn.
Las Ecuaciones m=1,2,,M 1 y se
describen el comportamiento de la
etapa m-sima, donde las entradas
son
Fm_1(z)
y
Gm_1(z),
proporcionndolas salidas Fm(z) y
Gm(z). En conjunto, las Ecuaciones
F0 ( z ) G0 ( z ) X ( z )
a m=1,2,,M
1 son un conjunto de ecuaciones
recursivas que describen el filtro en
celosa.
Como vemos en la figura inferior,
primera etapa, la entrada x(n) est
conectada a f0(n) y g0(n), y la salida
f(n) de la ltima etapa se considera
la salida del filtro
Dado que el sistema tiene dos
salidas, FM(Z) y GM(Z), y una nica
entrada, X(z), podemos diferenciar
dos funciones de transferencia:
AM ( z )
FM ( z ) FM ( z )
,
X ( z)
F0 ( z )
A0 ( z ) B0 ( z ) 1
Am ( z ) Am 1 ( z ) K m z 1 Bm 1 ( z ),
m=1,2,,M 1
Bm ( z ) K m Am 1 ( z ) z 1 Bm 1 ( z ),
B( z )
GM ( z) GM ( z)
,
X ( z)
G0 ( z )
m=1,2,,M 1
Como partimos de los coeficientes
del filtro FIR para la realizacin en
forma directa, tenemos el polinomio
A(z) que es:
3 1 1 2 1 3
z z z
4
2
4
1
3 (0) 3 (1) z 3 (2) z 2 3 (3) z 3
A3 ( z ) 1
Adems,
sabemos
que
los
coeficientes del filtro de salida B(z)
son inversos a los de A(z) por lo
que:
Figura 5. Filtro en celosa de M-1
etapas
3 ( z ) 3 (0) 3 (1) z 1 3 (2) z 2 3 (3) z 3
1 1 1 3 2
z z z 3
4 2
4
y por tanto
3 (1) 3 (2), 3 (0) 3 (3),
3 (3) 3 (0), 3 (2) 3 (1),
Figura 6. Estructura de cada etapa.
por lo que dividiendo las Ecuaciones
F0 ( z ) G0 ( z ) X ( z )
a m=1,2,,M
1 por X(z), tenemos:
Deseamos
determinar
correspondientes parmetros
filtro de celosa {Ki}. Para
K i i (i )
los
del
ello
sabemos que
. Dado que
el grado del polinomio A(z) es tres,
tendremos una celosa de tres
etapas, de la cual podremos obtener
inmediatamente el parmetro
K 3 3 (3) 1 / 4
Al repetir la recursin descendente,
obtenemos:
A1 ( z )
.
Para obtener el parmetro K2
necesitaremos el polinomio A2(z). La
relacin
recursiva
general
se
determina fcilmente a partir de las
Ecuaciones
1
Am ( z ) Am 1 ( z ) K m z Bm 1 ( z ),
y
por
A2 ( z ) K 2 B2 ( z )
2
1 K 2
lo
que
1 1
z
2
finalmente
K 1 1 (1) 1 / 2
con lo que los
coeficientes de la estructura celosa
resultan
Bm ( z ) K m Am 1 ( z ) z Bm 1 ( z ),
donde:
Am ( z ) Am1 ( z ) K m z 1 Bm1 ( z )
Am1 ( z ) K m Bm ( z ) K m A m1 ( z )
K1 = 1/2,
K2 = 1/3,
1/4
K3 =
La estructura en celosa del sistema
FIR propuesto es la representada en
la Figura
K m , Bm
Donde
si conocemos
y A(z)
Am1 ( z ) :
podemos resolver
Am1 ( z )
Am ( z ) K m Bm ( z )
2
1 K m
La cual es precisamente la
recursin descendente usada en el
test de estabilidad de Schr-Cohn.
Mediante
la
recursin
descendiente, con m = 3, se
obtiene:
A2 ( z )
A3 ( z ) K 3 B3 ( z )
2
1 K3
2 1 1 2
z z
3
3
Por lo que:
K 2 2 (2) 1 / 3 yB2 ( z )
1 2 1
z z 2
3 3
Ejercicio 2
Obtenga los coeficientes de
reflexin correspondientes al
filtro FIR con funcin de
transferencia
H ( z) 2
7 1
z z 2
2
Resolucin
Para
aplicar
la
recursin
descendente mediante la que se
obtendrn
los
coeficientes
en
celosa,
tambin
denominados
coeficientes
de
reflexin,
el
coeficiente am (0) debe definirse
como
1
por
conveniencia
matemtica, luego tomaremos:
H ( z ) 2.H ' ( z )
H ' (z) 1
recursiva
que:
7 1 1 2
z z
4
2
As, obtendremos los coeficientes en
celosa de H'(z) y aplicaremos un
factor de ganancia 2 a la salida de
la estructura resultante.
Otra peculiaridad que debe tenerse
en cuenta es que, en el caso de que
K2(0)
hubiera
sido
1,
nos
hubiramos encontrado con K2= - 1
= 2(2).
Ha de tenerse presente que siempre
que un parmetro de celosa es
Km 1
es una indicacin de que el
polinomio Am-1 (z) tiene una raz en
la circunferencia de radio unidad.
As, siempre que se obtiene un
Km 1
parmetro de celosa
se
rompe la ecuacin recursiva y no se
podr
seguir
la
recursividad
descendente.
En estos casos, dicha raz puede ser
factorizada y extrada de A m-1 (z),
continuando el proceso iterativo
para el sistema de orden reducido.
Siguiendo con el caso que nos
ocupa, dado H'(z) tenemos que los
polinomios A2(z) y B2(z) se definen
como
A2 ( z ) 1
7 1 1 2
z z
4
2
1 7
B2 ( z ) z 1 z 2
2 4
Por tanto, K2 =2(2) = - 1/2.
Siguiendo
con
la
ecuacin
descendente
A1 ( z )
tenemos
A2 ( z ) K 2 B2 ( z )
1 K 22
7 1 1 2 1 1 7 1 2
z z z z
4
2
2 2 4
1
1
4
7 1
z
2
y por lo tanto K1=-7/2. La
representacin final del diagrama de
bloques
de
H(z)
segn
una
estructura de celosa se representa
en la Figura, cabe remarcar en esta
realizacin el factor de ganancia dos
de la salida del mismo.
CONCLUSIONES
Las principales aplicaciones de Celosa es en Procesado digital de la voz,
Implementacin de filtros adaptativos y Tratamiento de seales
geofsicas.
Las estructuras en celosa tanto FIR como IIR se caracterizan por los
mismos coeficientes de reflexin, ki, diferencindose nicamente en su
interconexin.
Una etapa de la estructura en celosa de filtro FIR est compuesta por
dos lneas, una Directa y la otra con retardo donde un sumador en cada
lnea adiciona la seal proveniente de la otra, ponderada por una
ganancia. Cada lnea tiene una entrada y una salida propia.
El sistema slo polos ser estable si sus polos se encuentran en el
interior de la circunferencia de radio unidad lo cual implica que |km|< 1
para todo m.
Los algoritmos de conversin de parmetros entre el sistema en forma
directa bm(k) de un sistema FIR y los parmetros de la estructura en
celosa, ki, se aplican tambin a la estructura slo polos.
BIBLIOGRAFIA Y/O ENLACES
[Link]
[Link]
[Link]
[Link]