Volver a la página principal


SUPERFICIE DE ERROR





    La función de coste que utilizamos en la minimización es el error cuadrático medio definido como

z = E{e(n) e*(n)}
que en notación vectorial adquiere la forma
z(w) = rd2(0) - wHp - pH w +wHR w
donde w es el vector de coeficientes del filtro.

    Utilizando las igualdades    p = R wopt   y    RH = R,  con  wopt el vector óptimo de coeficientes que minimizan el error cuadrático medio, obtenemos que el error cuadrático medio dependiente de w es

z(w) = rd2(0) - wHR wopt- woptHR w +wHR w - woptHR w +  woptHR w

    Como el error mínimo es

z min =  rd2(0) -  woptHR w

    Sustituyendo obtenemos

z(w) = zmin  - wHR wopt- woptHR w +wHR w  + woptHR w


z(w) = zmin  +  (w - wopt)HR (w - wopt

 

    La función de coste presenta una función cuadrática; es un paraboloide con un mínimo en el valor w = wopt, como es de esperar, y admite una única solución en el filtrado de Wiener FIR,  propiedad que no siempre es compartida por el filtrado IIR.
    Existen contornos de igual error cuadrático medio que toman la forma de elipses.
 

    Después de la traslación ocasionada por (w - wopt), vamos a efectuar una rotación para cambiar los ejes de la representación. Para ello, expresamos la matriz de autocorrelación en términos de sus autovalores y autovectores

R = Q L QH

donde L es la matriz diagonal de autovalores l1, l2, ..., lM , y Q es la matriz de autovectores asociados a los autovalores.

    Si definimos

Z = QH (w - wopt)

podemos expresar el error cuadrático medio en forma canónica como

z(w) = zmin  +  ZHLZ


    Los nuevos ejes de la representación los determinan los autovectores de la matriz de autocorrelación R, y la excentricidad de las elipses depende de los autovalores. Así, cuando la distribución de autovalores es uniforme, y éstos son similares, los contornos de igual error se asemejan a circunferencias, y cuanto mayor sea la diferencia entre estos autovalores, más excéntricas serán las elipses.



    De este estudio, podemos extraer las siguientes conclusiones: