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.