Factorización Espectral de Procesos

 

    La densidad espectral de potencia de un proceso puede expresarse de forma general como

donde

es la parte continua del espectro, y

corresponde a la parte discreta.



    Comenzamos con la parte continua del espectro

    En primer lugar, evaluamos si se cumple la condición de Paley-Wiener:

    Si se cumple esta condición necesaria y suficiente, podemos factorizar la parte continua del espectro de la forma siguiente:

donde

- K0 es una constante positiva.
- Hmin(z) es un filtro causal y estable de fase mínima (sus polos y ceros se encuentran dentro de la circunferencia unidad). Tiene inversa causal y estable.
- Hmin*(1/z*)  responde a la parte anticausal. Es un filtro de fase máxima.

    Podemos enunciar las siguientes propiedades de un proceso regular:




    Consideramos ahora la parte discreta de la densidad espectral

    Según el teorema de Descomposición de Wold, un proceso aleatorio general puede ser descrito como suma de dos procesos,

x(n) = xp(n) + xr(n)

donde xr(n) es un proceso aleatorio regular y xp(n) un proceso predecible, siendo xr(n) ortogonal a xp(n),

E{xr(m)xp*(n)} = 0

    Supongamos que queremos estimar el valor de una muestra de un proceso en el instante n, utilizando las muestras anteriores de dicho proceso

obtenemos así un error, diferencia entre el valor real y el estimado:

    Decimos que un proceso es predecible si la potencia del error

es igual a cero   se2 = 0.

    Un proceso estocástico se puede predecir si y sólo si su densidad espectral es la forma:

    Esta densidad de potencia está relacionada con procesos periódicos

donde ai = |ai| ejfi;   |ai|  y fi son variables aleatorias.

    Para que el proceso sea predecible, se ha de cumplir

    Consideramos un filtro de la forma

por lo que la densidad espectral de potencia del error es


    Los ceros del filtro coinciden con las frecuencias de las deltas del espectro predecible.

    Podemos obtener la factorización siguiente

con K0 = 1. Como existe esta factorización espectral, se ha de cumplir la condición de Paley-Wiener para G(ejw)

lo que significa que G(ejw) no puede ser cero, excepto en unos puntos, para que la integral no sea infinita. Por el contrario, para obtener potencia de error cero, Sxp(ejw) si puede ser cero excepto en unos puntos concretos, pues no tiene porqué cumplir la condición anterior.


Volver al temario