Series de Fourier

Las series de Fourier son una herramienta matemática que se puede emplear,
entre otras aplicaciones, para representar señales periódicas como una
sumatoria infinita de términos senoidales. Para tal efecto, considerar lo
siguiente:
Sea 𝑓(𝑡) una señal continua y periódica, de periodo 𝑇, integrable en el intervalo
[−𝑇/2, 𝑇/2], y que posiblemente contenga un número finito de discontinuidades
(Condiciones de Dirichlet). Definiendo su frecuencia fundamental como:

2𝜋
𝜔0 = = 2𝜋𝑓 (1)
𝑇
Es posible aproximar dicha señal mediante la siguiente expansión:
∞

𝑓(𝑡) = 𝑎0 + ∑[𝑎𝑛 cos(𝑛𝜔0 𝑡) + 𝑏𝑛 sen(𝑛𝜔0 𝑡)] (2)

𝑛=1

Donde la expresión (2) se le conoce como expansión en serie trigonométrica de

Fourier para 𝑓(𝑡). El significado primordial de dicha ecuación consiste en que
una señal periódica se puede descomponer como una sumatoria infinita de
términos senoidales, ponderados por coeficientes 𝑎0 , 𝑎𝑛 y 𝑏𝑛 (conocidos como
coeficientes de Fourier), los cuales aumentan su frecuencia fundamental a razón
de múltiplos enteros de esta (mediante el término 𝑛). Estos términos de
frecuencia múltiple entera, se les denomina armónicos. Los coeficientes citados,
se calculan de la siguiente forma:
1 𝑇/2
𝑎0 = ∫ 𝑓(𝑡)𝑑𝑡 (3)
𝑇 −𝑇/2
2 𝑇/2
𝑎𝑛 = ∫ 𝑓(𝑡)cos⁡(𝑛𝜔0 𝑡)𝑑𝑡 (4)
𝑇 −𝑇/2
2 𝑇/2
𝑏𝑛 = ∫ 𝑓(𝑡)sen⁡(𝑛𝜔0 𝑡)𝑑𝑡 (5)
𝑇 −𝑇/2
Paréntesis: Algunas expresiones útiles
Integración por partes

∫ 𝑢𝑑𝑣 = 𝑢𝑣 − ∫ 𝑣𝑑𝑢

Donde se eligen las funciones 𝑑𝑣 y 𝑑𝑢, de tal forma, que 𝑑𝑣 sea fácil de
integrar y que 𝑢 desaparezca o involucre una forma conveniente

Derivadas e integrales de funciones trigonométricas

Ejemplo: Calcular la expansión en series de Fourier para f(t)

Se observa que:
𝑇 = 2𝜋
2𝜋
𝜔0 = = 1⁡[𝑟𝑎𝑑/𝑠]
𝑇
Solución:
• Para el término de DC 𝑎0

1 𝑇/2
𝑎0 = ∫ 𝑓(𝑡)𝑑𝑡
𝑇 −𝑇/2
1 𝜋
𝑎0 = ∫ 𝑡𝑑𝑡
2𝜋 −𝜋
+𝜋
1 𝑡2
𝑎0 = [ ]
2𝜋 2 −𝜋
1 𝜋 2 (−𝜋)2
𝑎0 = [ − ]
2𝜋 2 2

𝑎0 = 0

• Para el término 𝑎𝑛
𝑇
2 2
𝑎𝑛 = ∫ 𝑓(𝑡) cos(𝑛𝜔0 𝑡) 𝑑𝑡
𝑇 −𝑇
2
𝜋𝜋
1
𝑎𝑛 = ∫ 𝑡cos⁡(𝑛𝜔0 𝑡)𝑑𝑡
𝜋 −𝜋
1
Eligiendo 𝑑𝑣 = cos(𝑛𝜔0 𝑡) → ⁡⁡𝑣 = 𝑛𝜔 𝑠𝑒𝑛(𝑛𝜔0 𝑡) y 𝑢 = 𝑡⁡ → ⁡⁡𝑑𝑢 = 𝑑𝑡
0

+𝜋
1 𝑡 1
𝑎𝑛 = [ 𝑠𝑒𝑛(𝑛𝜔0 𝑡) − ∫ 𝑠𝑒𝑛(𝑛𝜔0 𝑡) 𝑑𝑡]
𝜋 𝑛𝜔0 𝑛𝜔0 −𝜋
+𝜋
1 𝑡 1
𝑎𝑛 = [ 𝑠𝑒𝑛(𝑛𝜔0 𝑡) + 𝑐𝑜𝑠(𝑛𝜔0 𝑡)]
𝜋 𝑛𝜔0 (𝑛𝜔0 )2 −𝜋

1 𝜋 1 −𝜋 1
𝑎𝑛 = {[ 𝑠𝑒𝑛(𝑛𝜔0 𝜋) + 2
𝑐𝑜𝑠(𝑛𝜔0 𝜋)] − [ 𝑠𝑒𝑛(−𝑛𝜔0 𝜋) + 𝑐𝑜𝑠(−𝑛𝜔0 𝜋)]}
𝜋 𝑛𝜔0 (𝑛𝜔0 ) 𝑛𝜔0 (𝑛𝜔0 )2
1 𝜋 1 𝜋 1
𝑎𝑛 = {[ 𝑠𝑒𝑛(𝑛𝜔0 𝜋) + 2
𝑐𝑜𝑠(𝑛𝜔0 𝜋)] − [ 𝑠𝑒𝑛(𝑛𝜔0 𝜋) + 𝑐𝑜𝑠(𝑛𝜔0 𝜋)]}
𝜋 𝑛𝜔0 (𝑛𝜔0 ) 𝑛𝜔0 (𝑛𝜔0 )2

𝑎𝑛 = 0
• Para el término 𝑏𝑛

2 𝑇/2
𝑏𝑛 = ∫ 𝑓(𝑡)sen⁡(𝑛𝜔0 𝑡)𝑑𝑡
𝑇 −𝑇/2

1 𝜋
𝑏𝑛 = ∫ 𝑡sen⁡(𝑛𝜔0 𝑡)𝑑𝑡
𝜋 −𝜋
 −1
Eligiendo 𝑑𝑣 = sen(𝑛𝜔0 𝑡) → ⁡⁡𝑣 = 𝑛𝜔 𝑐𝑜𝑠(𝑛𝜔0 𝑡) y 𝑢 = 𝑡⁡ → ⁡⁡𝑑𝑢 = 𝑑𝑡
0

+𝜋
1 −𝑡 −1
𝑏𝑛 = [ 𝑐𝑜𝑠(𝑛𝜔0 𝑡) − ∫ 𝑐𝑜𝑠(𝑛𝜔0 𝑡) 𝑑𝑡]
𝜋 𝑛𝜔0 𝑛𝜔0 −𝜋
+𝜋
1 −𝑡 1
𝑏𝑛 = [ 𝑐𝑜𝑠(𝑛𝜔0 𝑡) + 𝑠𝑒𝑛(𝑛𝜔0 𝑡)]
𝜋 𝑛𝜔0 (𝑛𝜔0 )2 −𝜋

1 −𝜋 1 −(−𝜋) 1
𝑏𝑛 = {[ 𝑐𝑜𝑠(𝑛𝜔0 𝜋) + 2
𝑠𝑒𝑛(𝑛𝜔0 𝜋)] − [ 𝑐𝑜𝑠(−𝑛𝜔0 𝜋) + 𝑠𝑒𝑛(−𝑛𝜔0 𝜋)]}
𝜋 𝑛𝜔0 (𝑛𝜔0 ) 𝑛𝜔0 (𝑛𝜔0 )2

1 −𝜋 1 𝜋 1
𝑏𝑛 = {[ 𝑐𝑜𝑠(𝑛𝜔0 𝜋) + 2
𝑠𝑒𝑛(𝑛𝜔0 𝜋)] − [ 𝑐𝑜𝑠(𝑛𝜔0 𝜋) − 𝑠𝑒𝑛(𝑛𝜔0 𝜋)]}
𝜋 𝑛𝜔0 (𝑛𝜔0 ) 𝑛𝜔0 (𝑛𝜔0 )2

1 −2𝜋 2
𝑏𝑛 = [ 𝑐𝑜𝑠(𝑛𝜔0 𝜋) + 𝑠𝑒𝑛(𝑛𝜔0 𝜋)]
𝜋 𝑛𝜔0 (𝑛𝜔0 )2

Considerando que 𝜔0 = 1
1 −2𝜋 2
𝑏𝑛 = [ 𝑐𝑜𝑠(𝑛𝜋) + 2 𝑠𝑒𝑛(𝑛𝜋)]
𝜋 𝑛 𝑛
Donde 𝑠𝑒𝑛(𝑛𝜋) = 0
−2
𝑏𝑛 = 𝑐𝑜𝑠(𝑛𝜋)
𝑛
Se sabe que 𝑐𝑜𝑠(𝑛𝜋) = (−1)𝑛
−2
𝑏𝑛 = (−1)𝑛
𝑛
2
𝑏𝑛 = (−1)𝑛+1
𝑛

Finalmente

∞
𝟐
𝒇(𝒕) = ∑ (−𝟏)𝒏+𝟏 𝐬𝐞𝐧(𝒏𝒕)
𝒏
𝒏=𝟏

Empleando el script “ejemplo_serie_Fourier.m” se generaron las siguientes

aproximaciones empleando la serie de Fourier calculada de forma analítica
Simplificaciones por simetría
Existen ocasiones en las cuales la señal a la cual se le desea calcular su expansión
presenta alguna clase de simetría (par o impar). En caso de presentarse alguna
de las mencionadas condiciones, y dadas las propiedades de la serie de Fourier,
se pueden emplear las siguientes simplificaciones.
• Si la señal 𝑓(𝑡) es par en el intervalo [−𝑇/2, 𝑇/2]
∞

𝑓(𝑡) = 𝑎0 + ∑[𝑎𝑛 cos(𝑛𝜔0 𝑡)] (6)

𝑛=1
Con

2 𝑇/2
𝑎0 = ∫ 𝑓(𝑡)𝑑𝑡 (7)
𝑇 0
4 𝑇/2
𝑎𝑛 = ∫ 𝑓(𝑡)cos⁡(𝑛𝜔0 𝑡)𝑑𝑡 (8)
𝑇 0

𝑏𝑛 = 0 (9)
 • Sila señal 𝑓(𝑡) es impar en el intervalo [−𝑇/2, 𝑇/2]
∞

𝑓(𝑡) = ∑[𝑏𝑛 sen(𝑛𝜔0 𝑡)] (10)

𝑛=1
Con
4 𝑇/2
𝑏𝑛 = ∫ 𝑓(𝑡)sen⁡(𝑛𝜔0 𝑡)𝑑𝑡 (11)
𝑇 0
𝑎0 = 0 (12)
𝑎𝑛 = 0 (13)

Ejemplo: Calcular la expansión en series de Fourier para la siguiente señal

simétrica

Solución:
i) En el intervalo de [−𝜋, 𝜋] la función es igual a 𝑓(𝑡) = |𝑡|, con 𝜔0 = 1 y
T=2𝜋. De esta forma, se puede emplear la siguiente propuesta

𝑡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡0 < 𝑡 < 𝜋

𝑓(𝑡) = {
−𝑡⁡⁡⁡⁡⁡ − 𝜋 < 𝑡 < 0
Ahora bien, se plantean las correspondientes integrales. Para este caso, la
integral completa en el intervalo [−𝜋, 𝜋] se divide en dos, cada una con un
distinto argumento de integración, de acuerdo con sus correspondientes
límites
1 0 𝜋
𝑎0 = [∫ (−𝑡)𝑑𝑡 + ∫ 𝑡𝑑𝑡]
𝑇 −𝜋 0

2 0 𝜋
𝑎𝑛 = [∫ −𝑡cos⁡(𝑛𝜔0 𝑡)𝑑𝑡 + ∫ 𝑡cos⁡(𝑛𝜔0 𝑡)𝑑𝑡]
𝑇 −𝜋 0
 2 0 𝜋
𝑏𝑛 = [∫ −𝑡sen⁡(𝑛𝜔0 𝑡)𝑑𝑡 + ∫ 𝑡sen⁡(𝑛𝜔0 𝑡)𝑑𝑡]
𝑇 −𝜋 0

Con este ejemplo, se muestra el procedimiento a seguir si se desea obtener

la expansión en series de Fourier para una señal definida a trozos, sin
importar cuantas definiciones posea en el intervalo [-T/2, T/2]. A
propósito de tal comentario, se debe aclarar que dicho intervalo no es
estricto, la única condición a cumplir consiste en que la diferencia de
ambos límites debe formar una banda de valor T, pudiéndose emplear
[0,T], [-T/2,T/2], [T/4,5T/4] etcétera; en cualquier caso, es imperativo
elegir los intervalos de tal forma que el número de integrales a calcular sea
el mínimo. En el caso de que tales operaciones sean resueltas empleando
alguna paquetería computacional, el a resolver no es un gran problema.

Paréntesis: Solución de integrales empleando Matlab

Las integrales se resuelven de forma muy simple empleando el comando
“int”, donde sus parámetros se enuncian a continuación

Para el caso de las integrales mencionadas, se empleó el script

“ejemplos_serieF_integrales.m”, entregando el siguiente resultado
Ahora bien, empleando las propiedades de simetría (señal par) se pueden aplicar
las ecuaciones de (6)-(9)
𝑏𝑛 = 0

2 𝑇/2
𝑎0 = ∫ 𝑓(𝑡)𝑑𝑡
𝑇 0

1 𝜋 1 𝑡2 1 𝜋 2 02
𝑎0 = ∫ 𝑡𝑑𝑡 = [ ] = [ − ]
𝜋 0 𝜋 2 𝜋 2 2
𝜋
𝑎0 =
2

4 𝑇/2
𝑎𝑛 = ∫ 𝑓(𝑡)cos⁡(𝑛𝜔0 𝑡)𝑑𝑡
𝑇 0
2 𝜋
𝑎𝑛 = ∫ 𝑡cos⁡(𝑛𝜔0 𝑡)𝑑𝑡
𝜋 0
1
Eligiendo 𝑑𝑣 = cos(𝑛𝜔0 𝑡) → ⁡⁡𝑣 = 𝑛𝜔 𝑠𝑒𝑛(𝑛𝜔0 𝑡) y 𝑢 = 𝑡⁡ → ⁡⁡𝑑𝑢 = 𝑑𝑡
0

𝜋
2 𝑡 1
𝑎𝑛 = [ 𝑠𝑒𝑛(𝑛𝜔0 𝑡) − ∫ 𝑠𝑒𝑛(𝑛𝜔0 𝑡)𝑑𝑡]
𝜋 𝑛𝜔0 𝑛𝜔0 0
𝜋
2 𝑡 1
𝑎𝑛 = [ 𝑠𝑒𝑛(𝑛𝜔0 𝑡) + 𝑐𝑜𝑠(𝑛𝜔0 𝑡)]
𝜋 𝑛𝜔0 (𝑛𝜔0 )2 0

Recordar que 𝜔0 = 1
2 𝜋 1 0 1
𝑎𝑛 = {[ 𝑠𝑒𝑛(𝑛𝜋) + 2
𝑐𝑜𝑠(𝑛𝜋)] − [ 𝑠𝑒𝑛(𝑛0) + 𝑐𝑜𝑠(𝑛0)]}
𝜋 𝑛 (𝑛) 𝑛 (𝑛)2
2 1 1
𝑎𝑛 = [ 2 (−1)𝑛 − 2 ]
𝜋 𝑛 𝑛
2
𝑎𝑛 = [(−1)𝑛 − 1]
𝜋𝑛2
∞
𝜋 2
𝑓 (𝑡) = + ∑ { 2 [(−1)𝑛 − 1] cos(𝑛𝑡)}
2 𝜋𝑛
𝑛=1
Series exponenciales de Fourier
Tomando como punto de partida las series trigonométricas, se puede realizar
una generalización de estas empleando funciones exponenciales complejas. Para
ello se debe considerar la identidad de Euler.

𝑒 ±𝑗𝜃 = cos(𝜃) ± 𝑗𝑠𝑒𝑛(𝜃) (14)

Lo que conlleva a las siguientes definiciones:

1
cos(𝜃) = [𝑒 𝑗𝜃 + 𝑒 −𝑗𝜃 ] (15)
2
1 𝑗𝜃 (16)
sin(𝜃) = [𝑒 − 𝑒 −𝑗𝜃 ]
2𝑗

Aplicando las ecuaciones (15) y (16) en (2)

∞
1 1 1 1
𝑓(𝑡) = 𝑎0 + ∑ [𝑎𝑛 ( 𝑒 𝑗𝑛𝜔0𝑡 + 𝑒 −𝑗𝑛𝜔0𝑡 ) + 𝑏𝑛 ( 𝑒 𝑗𝑛𝜔0𝑡 − 𝑒 −𝑗𝑛𝜔0𝑡 )] (17)
2 2 2𝑗 2𝑗
𝑛=1

Reagrupando
∞
1 1
𝑓(𝑡) = 𝑎0 + ∑ [ 𝑒 𝑗𝑛𝜔0𝑡 (𝑎𝑛 − 𝑗𝑏𝑛 ) + 𝑒 −𝑗𝑛𝜔0𝑡 (𝑎𝑛 + 𝑗𝑏𝑛 )] (18)
2 2
𝑛=1
Ahora, se proponen las siguientes igualdades
𝑐0 = 𝑎0 (19)
1 (20)
𝑐𝑛 = (𝑎𝑛 − 𝑗𝑏𝑛 )
2

La serie descrita por la ecuación (18) se puede plantear como:

𝑓(𝑡) = 𝑐0 + ∑[𝑐𝑛 𝑒 𝑗𝑛𝜔0𝑡 + 𝑐𝑛∗ 𝑒 −𝑗𝑛𝜔0𝑡 ] (21)

𝑛=1

Donde los coeficientes se pueden calcular de la siguiente forma:

1 𝑇/2 (22)
𝑐0 = ∫ 𝑓(𝑡)𝑑𝑡
𝑇 −𝑇/2
1 𝑇/2
𝑐𝑛 = ∫ 𝑓(𝑡)𝑒 −𝑗𝑛𝜔0𝑡 𝑑𝑡 (23)
𝑇 −𝑇/2

Tomando en cuenta al complejo conjugado

1 𝑇/2
𝑐𝑛∗ = ∫ 𝑓(𝑡)𝑒 𝑗𝑛𝜔0𝑡 𝑑𝑡 = 𝑐−𝑛 (24)
𝑇 −𝑇/2

Se plantea la siguiente expresión, conocida como la serie exponencial de Fourier

∞

𝑓(𝑡) = ∑ [𝑐𝑛 𝑒 𝑗𝑛𝜔0𝑡 ] (25)

𝑛=−∞
Con coeficientes generalizados
1 𝑇/2
𝑐𝑛 = ∫ 𝑓(𝑡)𝑒 −𝑗𝑛𝜔0𝑡 𝑑𝑡 (26)
𝑇 −𝑇/2
La serie exponencial de Fourier tendrá coeficientes de valores complejos; así las
cosas, se cumplirán las siguientes propiedades:
1) El coeficiente generalizado 𝑐𝑛 es una función de la variable discreta 𝑛𝜔0
(armónicos), donde n es un entero y 𝜔0 representa la frecuencia
fundamental de la señal periódica f(t)
𝑐𝑛 = 𝑐𝑛 (𝑛𝜔0 ) (27)

2) El coeficiente generalizado cumple la siguiente condición

𝑐𝑛 (−𝑛𝜔0 ) = [𝑐𝑛 (𝑛𝜔0 )]∗ (28)

3) Al ser una función compleja, tendrá un módulo y un argumento

𝑐𝑛 = |𝑐𝑛 |𝑒 𝑗Φ𝑛 (29)

Donde
|𝑐𝑛 |: Espectro de Amplitud (función par)

Φ𝑛 : Espectro de Fase (función impar)

4) La señal periódica f(t) también se puede escribir como:

∞

𝑓(𝑡) = 𝑐0 + ∑ 2|𝑐𝑛 |cos⁡(𝑛𝜔0 + Φ𝑛 ) (30)

𝑛=1

Ejemplo: Calcular la expansión en series de Fourier para f(t)

 Se observa que:
𝑇 = 2𝜋
2𝜋
𝜔0 = = 1⁡[𝑟𝑎𝑑/𝑠]
𝑇

Calculando
1 𝑇/2
𝑐𝑛 = ∫ 𝑓(𝑡)𝑒 −𝑗𝑛𝜔0𝑡 𝑑𝑡
𝑇 −𝑇/2

1 𝜋 −𝑗𝑛𝜔 𝑡
𝑐𝑛 = ∫ 𝑡𝑒 0 𝑑𝑡
2𝜋 −𝜋

−1
Eligiendo 𝑑𝑣 = 𝑒 −𝑗𝑛𝜔0𝑡 → ⁡⁡𝑣 = 𝑗𝑛𝜔 𝑒 −𝑗𝑛𝜔0𝑡 y 𝑢 = 𝑡⁡ → ⁡⁡𝑑𝑢 = 𝑑𝑡
0

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

−1 𝜋 −𝑗𝑛𝜔 𝜋 1 −𝑗𝑛𝜔0 𝜋
−𝜋 𝑗𝑛𝜔 𝜋 1
𝑐𝑛 = {[ 𝑒 0 + 𝑒 ] − [ 𝑒 0 + 𝑒 𝑗𝑛𝜔0𝜋 ]}
2𝜋 𝑗𝑛𝜔0 (𝑛𝜔0 )2 𝑗𝑛𝜔0 (𝑛𝜔0 )2
−1 𝜋 −𝑗𝑛𝜔 𝜋 1 𝜋 𝑗𝑛𝜔 𝜋 1
𝑐𝑛 = { 𝑒 0 +
2
𝑒 −𝑗𝑛𝜔0𝜋 + 𝑒 0 − 𝑒 𝑗𝑛𝜔0𝜋 }
2𝜋 𝑗𝑛𝜔0 (𝑛𝜔0 ) 𝑗𝑛𝜔0 (𝑛𝜔0 )2
−1 𝜋 1
𝑐𝑛 = { [𝑒 𝑗𝑛𝜔0𝜋 + 𝑒 −𝑗𝑛𝜔0𝜋 ] + [𝑒 −𝑗𝑛𝜔0𝜋 − 𝑒 𝑗𝑛𝜔0𝜋 ]}
2𝜋 𝑗𝑛𝜔0 (𝑛𝜔0 )2
−1 𝜋 1
𝑐𝑛 = { [2cos⁡(𝑛𝜔0 𝜋)] + [−2𝑗𝑠𝑒𝑛(𝑛𝜔0 𝜋)]}
2𝜋 𝑗𝑛𝜔0 (𝑛𝜔0 )2

Sustituyendo 𝜔0 = 1
−1 2𝜋 (−1) (−2)
𝑐𝑛 = [cos⁡(𝑛𝜔0 𝜋)] + [𝑗𝑠𝑒𝑛(𝑛𝜔0 𝜋)]
2𝜋 𝑗𝑛𝜔0 2𝜋 (𝑛𝜔0 )2
−1 1
𝑐𝑛 = [cos⁡(𝑛𝜋)] + 2 [𝑗𝑠𝑒𝑛(𝑛𝜋)]
𝑗𝑛 𝜋𝑛
 𝑗
𝑐𝑛 = cos⁡(𝑛𝜋)
𝑛
∞
𝑗
𝑓(𝑡) = ∑ [ cos⁡(𝑛𝜋)𝑒 𝑗𝑛𝜔0 𝑡 ]
𝑛
𝑛=−∞