Series de Fourier (LAJS)

1 Introducción

La manera más común de representar la serie de Fourier de una función [math]f(x) [/math] es mediante

[math] f(x) \sim\frac{a_0}{2 \pi} + \sum_{n=1}^{N} \left( a_n \frac{1}{\sqrt{\pi}} \cos(n x) + b_n \frac{1}{\sqrt{\pi}} \sin (n x)\right) [/math]

siendo [math] a_0,a_n,b_n [/math] los coeficientes de Fourier, que se calculan mediante las integrales:

[math] a_0 = \int_{-\pi}^{\pi} \frac{1}{\sqrt{\pi}} f(x) dx \quad a_n = \int_{-\pi}^{\pi} \frac{1}{\sqrt{\pi}} \cos(n x) f(x) dx \quad b_n = \int_{-\pi}^{\pi} \frac{1}{\sqrt{\pi}} \sin(n x) f(x) dx [/math]

Es decir, la función [math]f(x) [/math] se expresa como una combinación lineal de funciones seno y coseno con coeficientes constantes. Sin embargo, en algunos contextos reales, como la modelización de fenómenos naturales, no es razonable suponer que dichos coeficientes son conocidos de forma exacta. En algunos casos, la información que conocemos sobre el sistema puede ser intrínsecamente aleatoria, lo que motiva la introducción de modelos probabilísticos. Una forma natural de incorporar esta aleatoriedad consiste en sustituir los coeficientes de Fourier que conocemos por variables aleatorias. Este planteamiento permite definir funciones aleatorias cuyos valores, aun conservando una estructura armónica, presentan fluctuaciones controladas en función de la distribución estadística de los coeficientes.

2 Modelo matemático

Comenzamos definiendo la ecuación general que utilizaremos:

[math] f_{\sigma}(x) = \frac{a_0}{2} + \sum_{n=1}^{N} \left( a_n \cos(n\pi x) + b_n \sin(n\pi x) \right) [/math]

donde los coeficientes [math] a_n [/math] y [math] b_n [/math] serán variables aleatorias independientes e idénticamente distribuidas que se generarán a partir de distribuciones de probabilidad en lugar de calcularlos mediante integrales. Estudiaremos dos escenarios diferentes:

• Distribución uniforme: los coeficientes seguirán una distribución uniforme en el intervalo [math] [-1,1] [/math], es decir, [math]a_n, b_n \sim \mathcal{U}(-1, 1) [/math]

• Distribución normal: los coeficientes seguirán una distribución normal estándar, [math] a_n, b_n \sim \mathcal{N}(0, 1) [/math]

3 Cálculo de la esperanza, varianza y covarianza

3.1 Esperanza

Para empezar, en el caso que los coeficientes [math] a_n, b_n [/math] siguen una distribución [math] \mathcal{N}(0, 1) [/math], la esperanza se calcula:

[math] E[f_{\sigma}(x)] = E\left[ \frac{a_0}{2} + \sum_{n=1}^{N} \left( a_n \cos(n\pi x) + b_n \sin(n\pi x) \right) \right] = \frac{E[a_0]}{2} + \sum_{n=1}^{\infty} \big( E[a_n] \cos(n \pi x) + E[b_n] \sin(n \pi x) \big) = 0 [/math]

ya que [math] E[a_0]= E[a_n]= E[b_n]=0 [/math] para todo [math] n \in \{1,\dots , N\} [/math].

Por otro lado, en el caso de la distribución uniforme, esto es, [math]a_n, b_n \sim \mathcal{U}(-1, 1) [/math], la esperanza es