PROBLEMAS

7.1-1 Dado el proceso aleatorio

𝑋(𝑡) = 𝐴0 𝑐𝑜𝑠(𝜔0 𝑡 + Θ)
Donde 𝐴0 y 𝜔0 son constantes y Θ es una variable aleatoria uniformemente distribuida en el intervalo (0, 𝜋).

(a) ¿Es 𝑋(𝑡) estacionario en sentido amplio?

(b) Hallar la potencia de 𝑋(𝑡) usando la fórmula de la potencia media

(c) Hallar el espectro de potencia de 𝑋(𝑡) usando el espectro de densidad de potencia y calcular la potencia a
partir de la fórmula de potencia ¿Coinciden ambas potencias?

Solución:
∞ 𝜋 1 𝐴0
(a) 𝐸[𝑋(𝑡)] = ∫−∞ 𝑥𝑓Θ (𝜃)𝑑𝜃 = ∫0 𝐴0 𝑐𝑜𝑠(𝜔0 𝑡 + 𝜃) 𝜋 𝑑𝜃 = [𝑠𝑒𝑛(𝜔0 𝑡 + 𝜋) − 𝑠𝑒𝑛(𝜔0 𝑡)]
𝜋

Como se puede observar que la media depende del tiempo, entonces 𝑋(𝑡) no es estacionario en el sentido
amplio.

(b) 𝐸[𝑋 2 (𝑡)] = 𝐸[𝐴0 2 𝑐𝑜𝑠 2 (𝜔0 𝑡 + Θ)]

Por identidad trigonométrica se sabe:

1
𝑐𝑜𝑠 2 (𝛼) = (1 + 𝑐𝑜𝑠(2𝛼))
2
luego

𝐴0 2 𝐴0 2 𝐴0 2 𝜋 1 𝐴0 2
= 𝐸[1 + 𝑐𝑜𝑠(2𝜔0 𝑡 + 2Θ)] = + ∫ 𝑐𝑜𝑠(2𝜔0 𝑡 + 2𝜃) 𝑑𝜃 = +0
2 2 2 0 𝜋 2

De este modo, la potencia media de 𝑋(𝑡)

1 𝑇 1 𝑇 𝐴0 2 𝐴0 2
𝑃𝑋𝑋 = lim ∫ 𝐸[𝑋 2 (𝑡)] 𝑑𝑡 = lim ∫ 𝑑𝑡 =
𝑇→∞ 2𝑇 −𝑇 𝑇→∞ 2𝑇 −𝑇 2 2

(c) En primer lugar se debe hallar la Función de densidad espectral de 𝑋(𝑡) para ello aplico la transformada
de Fourier que se denota así:
𝑇 𝑇
𝑋𝑇 (𝜔) = ∫ 𝑥𝑇 (𝑡) 𝑒 −𝑗𝜔𝑡 𝑑𝑡 = ∫ 𝑥(𝑡) 𝑒 −𝑗𝜔𝑡 𝑑𝑡
−𝑇 −𝑇

Reemplazando valores
𝑇
𝑋𝑇 (𝜔) = ∫ 𝐴0 𝑐𝑜𝑠(𝜔0 𝑡 + Θ) 𝑒 −𝑗𝜔𝑡 𝑑𝑡
−𝑇

Por la expresión de Euler, funciones trigonométricas en variable compleja Z: b

𝑒 𝑗𝑍 + 𝑒 −𝑗𝑍
𝑐𝑜𝑠(Z) =
2
Luego
𝑇
𝑒 𝑗(𝜔0 𝑡+Θ) + 𝑒 −𝑗(𝜔0 𝑡+Θ) −𝑗𝜔𝑡
𝑋𝑇 (𝜔) = 𝐴0 ∫ [ ]𝑒 𝑑𝑡
−𝑇 2

𝐴0 𝑗Θ 𝑇 𝑗(𝜔 −𝜔)𝑡 𝐴0 𝑇
𝑋𝑇 (𝜔) = 𝑒 ∫ 𝑒 0 𝑑𝑡 + 𝑒 −𝑗Θ ∫ 𝑒 −𝑗(𝜔0 +𝜔)𝑡 𝑑𝑡
2 −𝑇 2 −𝑇
Desarrollando por partes:

𝑇 𝑇
𝑗(𝜔0 −𝜔)𝑡
𝑒 𝑗(𝜔0 −𝜔)𝑡 1
∫ 𝑒 𝑑𝑡 = [ ] = [𝑒 𝑗(𝜔0 −𝜔)𝑇 − 𝑒 −𝑗(𝜔0 −𝜔)𝑇 ]
−𝑇 𝑗(𝜔0 − 𝜔) −𝑇 𝑗(𝜔0 − 𝜔)
2 𝑒 𝑗(𝜔0 −𝜔)𝑇 − 𝑒 −𝑗(𝜔0 −𝜔)𝑇
= [ ]
(𝜔0 − 𝜔) 2𝑗

𝑇 𝑇
−𝑗(𝜔0 +𝜔)𝑡
𝑒 −𝑗(𝜔0 +𝜔)𝑡 1
∫ 𝑒 𝑑𝑡 = [ ] = [𝑒 −𝑗(𝜔0 +𝜔)𝑇 − 𝑒 𝑗(𝜔0 +𝜔)𝑇 ]
−𝑇 −𝑗(𝜔0 + 𝜔) −𝑇 −𝑗(𝜔0 + 𝜔)
2 𝑒 𝑗(𝜔0 +𝜔)𝑇 − 𝑒 −𝑗(𝜔0 +𝜔)𝑇
= [ ]
(𝜔0 + 𝜔) 2𝑗

Se utiliza la funcione trigonométrica en variable compleja Z:

𝑒 𝑗𝑍 − 𝑒 −𝑗𝑍
𝑠𝑒𝑛(Z) =
2𝑗
Luego:
𝑇
2𝑠𝑒𝑛[(𝜔0 − 𝜔)𝑇]
∫ 𝑒 𝑗(𝜔0 −𝜔)𝑡 𝑑𝑡 =
−𝑇 (𝜔0 − 𝜔)
𝑇
2𝑠𝑒𝑛[(𝜔0 + 𝜔)𝑇]
∫ 𝑒 −𝑗(𝜔0 +𝜔)𝑡 𝑑𝑡 =
−𝑇 (𝜔0 + 𝜔)

Reemplazando en:

𝐴0 𝑗Θ 2𝑠𝑒𝑛[(𝜔0 − 𝜔)𝑇] 𝐴0 2𝑠𝑒𝑛[(𝜔0 + 𝜔)𝑇]

𝑋𝑇 (𝜔) = ( 𝑒 )( ) + ( 𝑒 −𝑗Θ ) ( )
2 (𝜔 0 − 𝜔) 2 (𝜔0 + 𝜔)

𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]

𝑋𝑇 (𝜔) = 𝐴0 𝑇𝑒 𝑗Θ + 𝐴0 𝑇𝑒 −𝑗Θ
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
A continuación, se determina

|𝑋𝑇 (𝜔)|2 = 𝑋𝑇 (𝜔)𝑋𝑇∗ (𝜔)

Se aplica la propiedad de un numero complejo y su conjugado

|𝑍|2 = 𝑍𝑍̅ = (𝑎 + 𝑗𝑏)(𝑎 − 𝑗𝑏) = 𝑎2 + 𝑏 2

Además:

𝑒 𝑗Θ = 𝑐𝑜𝑥(Θ) + 𝑗𝑠𝑒𝑛(Θ) 𝑒 −𝑗Θ = 𝑐𝑜𝑥(Θ) − 𝑗𝑠𝑒𝑛(Θ)

Luego reemplazando las fórmulas de Euler 𝑒 𝑗Θ y 𝑒 −𝑗Θ para separar en parte real e imaginaria
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
𝑋𝑇 (𝜔) = 𝐴0 𝑇 [𝑐𝑜𝑠(Θ) + 𝑗𝑠𝑒𝑛(Θ)] + 𝐴0 𝑇 [𝑐𝑜𝑠(Θ) − 𝑗𝑠𝑒𝑛(Θ)]
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
Se ordena la parte real y la compleja

𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]

= 𝐴0 𝑇𝑐𝑜𝑠(Θ) [ + ] + 𝑗𝐴0 𝑇𝑠𝑒𝑛(Θ) [ − ]
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇 (𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇

Usando la propiedad de un numero complejo y su conjugado se reemplaza en:

 2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
|𝑋𝑇 (𝜔)|2 = (𝐴0 𝑇𝑐𝑜𝑠(Θ) [ + ])
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
+ (𝐴0 𝑇𝑠𝑒𝑛(Θ) [ − ])
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇

Sea:
𝑠𝑒𝑛[(𝜔−𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔+𝜔0 )𝑇]
𝑃 = 𝐴0 𝑇 𝑄= (𝜔−𝜔0 )𝑇
𝑅= (𝜔+𝜔0 )𝑇

= 𝑃2 𝑐𝑜𝑠(Θ)2 [𝑄 + 𝑅]2 + 𝑃2 𝑠𝑒𝑛(Θ)2 [𝑄 − 𝑅]2

= 𝑃2 𝑐𝑜𝑠(Θ)2 [𝑄 2 + 2𝑄𝑅 + 𝑅 2 ] + 𝑃2 𝑠𝑒𝑛(Θ)2 [𝑄 2 − 2𝑄𝑅 + 𝑅 2 ]

= 𝑃2 [𝑄 2 + 𝑅 2 ][𝑐𝑜𝑠(Θ)2 + 𝑠𝑒𝑛(Θ)2 ] + 2𝑄𝑅𝑃2 [𝑐𝑜𝑠(Θ)2 − 𝑠𝑒𝑛(Θ)2 ]

= 𝑃2 [𝑄 2 + 𝑅 2 ] + 2𝑄𝑅𝑃2 [𝑐𝑜𝑠(2Θ)] = 𝑃2 [𝑄 2 + 𝑅 2 + 2𝑄𝑅𝑐𝑜𝑠(2Θ)]

Entonces:
2 2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
|𝑋𝑇 (𝜔)|2 2
= (𝐴0 𝑇) [( ) +( )
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
+ 2( )( ) 𝑐𝑜𝑠(2Θ)]
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇

Calculando el valor esperado

𝜋 2 2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
𝐸[|𝑋𝑇 (𝜔)|2 ] = ∫ (𝐴0 𝑇)2 [( ) +( )
0 (𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇] 1
+ 2( )( ) 𝑐𝑜𝑠(2Θ)] 𝑑Θ
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇 𝜋
2 2
2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
= (𝐴0 𝑇) [( ) +( ) ]
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇] 1 𝜋
+ (𝐴0 𝑇)2 [2 ( )( )] [ ] ∫ 𝑐𝑜𝑠(2Θ) 𝑑Θ
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇 𝜋 ⏟0
0

Luego queda:
2 2
𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
𝐸[|𝑋𝑇 (𝜔)|2 ] = (𝐴0 𝑇) [(2
) +( ) ]
(𝜔 − 𝜔0 )𝑇 (𝜔 + 𝜔0 )𝑇
2 2
𝐸[|𝑋𝑇 (𝜔)|2 ] 𝐴20 𝜋 𝑇 𝑠𝑒𝑛[(𝜔 − 𝜔0 )𝑇] 𝑇 𝑠𝑒𝑛[(𝜔 + 𝜔0 )𝑇]
= { ( ) + ( ) }
2𝑇 2 𝜋 (𝜔 − 𝜔0 )𝑇 𝜋 (𝜔 + 𝜔0 )𝑇

Por (Lathi, 1968, pag. 24)

2
𝑇 𝑠𝑒𝑛(𝛼𝑇)
lim [ ] = 𝛿(𝛼)
𝑇→∞ 𝜋 𝛼𝑇

El espectro de densidad de potencia es

𝐸[|𝑋𝑇 (𝜔)|2 ] 𝐴20 𝜋

𝒢𝑋𝑋 (𝜔) = lim = [𝛿(𝜔 − 𝜔0 ) + 𝛿(𝜔 + 𝜔0 )]
𝑇→∞ 2𝑇 2
Por último, se utiliza este resultado para obtener la potencia media
Utilizando el concepto:

La función impulso unitario 𝛿(𝑥) puede definirse mediante su propiedad de integral:

∞
𝜙(𝑥0 ) = ∫ 𝜙(𝑥) 𝛿(𝑥 − 𝑥0 )𝑑𝑥
−∞

Donde 𝜙(𝑥) es cualquier función continua en el punto 𝑥 = 𝑥0 ; 𝛿(𝑥) puede interpretarse como una función con
amplitud infinita, área unitaria y duración cero.

1 ∞ 1 ∞ 𝐴20 𝜋
𝑃𝑋𝑋 = ∫ lim 𝒢𝑋𝑋 (𝜔) 𝑑𝜔 = ∫ [𝛿(𝜔 − 𝜔0 ) + 𝛿(𝜔 + 𝜔0 )] 𝑑𝜔
2𝜋 −∞ 𝑇→∞ 2𝜋 −∞ 2

𝐴20 ∞ ∞
𝐴20 𝐴2
= [∫ [1]𝛿(𝜔 − 𝜔0 ) 𝑑𝜔 + ∫ [1]𝛿(𝜔 + 𝜔0 ) 𝑑𝜔] = [1 + 1] = 0
4 −∞ −∞ 4 2

Por lo tanto, se comprueba que el valor obtenido concuerda entre la potencia media y el espectro de potencia.

7.1.12 Determinar cuál de las funciones siguientes puede o no puede ser válida como espectro de densidad de
potencia. Para aquellas que no sean válidas, explicar por qué.
𝜔2
(a) 𝜔6 +3𝜔2 +3
(b) 𝑒𝑥𝑝[−(𝜔 − 1)2 ]

𝜔2 𝜔4
(c) 𝜔4 +1
− 𝛿(𝜔) (d) 1+𝜔2 +𝑗𝜔6

Solución:

Utilizando las Propiedades del espectro de densidad de potencia

(1) 𝒢𝑋𝑋 (𝜔) ≥ 0

(2) 𝒢𝑋𝑋 (−𝜔) = 𝒢𝑋𝑋 (𝜔) 𝑋(𝑡) real

(3) 𝒢𝑋𝑋 (𝜔) es real

(a) Es válida, porque cumple con la propiedad 1, 2 y 3

(b) No es válida, no cumple con la propiedad 2; no es una función par

(c) No es válida, por definición la función impulso unitario o función Delta de Dirac
∞ 𝜔=𝑎
𝛿(𝜔) = {
0 𝜔≠𝑎
En la función puede tomar valor negativo en alguna parte

(d) No es válido, es un numero imaginario y no real (propiedad 3)

7.1.13 Repetir el problema anterior para las siguientes funciones:

𝑐𝑜𝑠(3𝜔) 1
(a) 1+𝜔2
(b) (1+𝜔2 )2

|𝜔| 1
(c) 1+2𝜔+𝜔2
(d)
√1+3𝜔2

Solución:

(a) No es válida, toma valores negativos porque el coseno esta acotado entre -1 y 1 (propiedad 1)

(b) Es válido, porque cumple con la propiedad 1, 2 y 3.

(c) No es válida, no cumple con la propiedad 2; no es una función par.

(d) No es válida, no cumple con la propiedad 1; no es real

7.1.16 Si 𝑋(𝑡) es un proceso estacionario, hallar el espectro de potencia de 𝑌(𝑡) = 𝐴0 + 𝐵0 𝑋(𝑡) en términos del
espectro de potencia de 𝑋(𝑡) si 𝐴0 y 𝐵0 son constantes.

Solución:

𝑅𝑌𝑌 (𝑡, 𝑡 + 𝜏) = 𝐸[𝑌(𝑡)𝑌(𝑡 + 𝜏)] = 𝐸[{𝐴0 + 𝐵0 𝑋(𝑡)}{𝐴0 + 𝐵0 𝑋(𝑡 + 𝜏)}]

= 𝐸[𝐴0 2 + 𝐴0 𝐵0 𝑋(𝑡 + 𝜏) + 𝐴0 𝐵0 𝑋(𝑡) + 𝐵0 2 𝑋(𝑡)𝑋(𝑡 + 𝜏)]

= 𝐸[𝐴0 2 ] + 𝐸[𝐴0 𝐵0 𝑋(𝑡 + 𝜏)] + 𝐸[𝐴0 𝐵0 𝑋(𝑡)] + 𝐸[𝐵0 2 𝑋(𝑡)𝑋(𝑡 + 𝜏)]

= 𝐴0 2 + 𝐴0 𝐵0 𝐸[𝑋(𝑡 + 𝜏)] + 𝐴0 𝐵0 𝐸[𝑋(𝑡)] + 𝐵0 2 𝐸[𝑋(𝑡)𝑋(𝑡 + 𝜏)]

= 𝐴0 2 + 2𝐴0 𝐵0 𝑋̅ + 𝐵0 2 𝑅𝑋𝑋 (𝑡, 𝑡 + 𝜏)

Luego el espectro de potencia se define como:
∞
𝒢𝑌𝑌 (𝜔) = ∫ 𝑅𝑌𝑌 (𝑡, 𝑡 + 𝜏)𝑒 −𝑗𝜔𝜏 𝑑𝜔
−∞

Reemplazando valores:
∞
𝒢𝑌𝑌 (𝜔) = ∫ [𝐴0 2 + 2𝐴0 𝐵0 𝑋̅ + 𝐵0 2 𝑅𝑋𝑋 (𝑡, 𝑡 + 𝜏)]𝑒 −𝑗𝜔𝜏 𝑑𝜔
−∞
∞ ∞
2
𝒢𝑌𝑌 (𝜔) = ∫ [𝐴0
⏟ + 2𝐴0 𝐵0 𝑋̅] 𝑒 −𝑗𝜔𝜏 𝑑𝜔 + 𝐵0 2 ∫ 𝑅𝑋𝑋 (𝑡, 𝑡 + 𝜏) 𝑒 −𝑗𝜔𝜏 𝑑𝜔
−∞ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 −∞

La transformada de Fourier para una constante (Tabla)

𝛼
Si: 𝑥(𝑡) = 𝑋(𝜔) = 𝛼𝛿(𝜔)
2𝜋

Reemplazando se tiene:

𝒢𝑌𝑌 (𝜔) = [𝐴0 2 + 2𝐴0 𝐵0 𝑋̅](2𝜋𝛿(𝜔)) + 𝐵0 2 𝒢𝑋𝑋 (𝜔)