TC01001 – Funções especiais em

telecomunicações

Série de Fourier

Ronaldo de Freitas Zampolo

FCT-ITEC-UFPA

2022-4
 Agenda

Introdução
Série de Fourier de funções com periodicidade 2π
Convergência da série de Fourier
Série de Fourier de funções de perı́odo arbitrário
Série de Fourier de funções não-periódicas em faixa limitada
Série de Fourier complexa
Propriedades da série de Fourier

TC01001 FET 2 / 34 Série de Fourier

Introdução
 Introdução

• Jean Baptiste Joseph Fourier

@ 21/03/1768, Auxerre
+ 16/05/1830, Paris
• Matemático e fı́sico
• 1822: Théorie Analytique de la
Chaleur

TC01001 FET 4 / 34 Série de Fourier

 Introdução

• Ideia básica: qualquer função f (x) periódica em 2π pode ser representada por
uma série trigonométrica infinita na forma
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

TC01001 FET 5 / 34 Série de Fourier

 Introdução

• Ideia básica: qualquer função f (x) periódica em 2π pode ser representada por
uma série trigonométrica infinita na forma
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Objetivo inicial: resolução de equações diferenciais

TC01001 FET 5 / 34 Série de Fourier

 Introdução

• Ideia básica: qualquer função f (x) periódica em 2π pode ser representada por
uma série trigonométrica infinita na forma
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Objetivo inicial: resolução de equações diferenciais

• Atualmente é ferramenta essencial para análise de formas de onda

TC01001 FET 5 / 34 Série de Fourier

 Introdução

• Ideia básica: qualquer função f (x) periódica em 2π pode ser representada por
uma série trigonométrica infinita na forma
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Objetivo inicial: resolução de equações diferenciais

• Atualmente é ferramenta essencial para análise de formas de onda
• Algumas áreas em que a análise de Fourier tem grande aplicação:
• Processamento de sinais
• Análise de circuitos
• Sistemas de controle
• Telecomunicações

TC01001 FET 5 / 34 Série de Fourier

Série de Fourier de funções com
periodicidade 2π
 Ortogonalidade de funções trigonométricas

• Condição de ortogonalidade de funções (compare com o produto escalar)

Z b
ψn(x)ψm(x)dx = 0, n 6= m
a

As funções ψn e ψm são ditas ortogonais no intervalo entre a e b.

TC01001 FET 7 / 34 Série de Fourier

 Ortogonalidade de funções trigonométricas

• Condição de ortogonalidade de funções (compare com o produto escalar)

Z b
ψn(x)ψm(x)dx = 0, n 6= m
a

As funções ψn e ψm são ditas ortogonais no intervalo entre a e b.

• Considerando que m ∈ Z e n ∈ Z, (convém deduzir, estudante)
Z π 
cos mxdx = 0 m 6= 0 Z π


 0 m 6= n
Z−ππ cos mx cos nxdx = π m = n 6= 0
−π
sin mxdx = 0

2π m=n=0

Z π −π (
Z π
cos mx sin nxdx = 0 0 m 6= n
sin mx sin nxdx =
−π −π π m = n 6= 0

TC01001 FET 7 / 34 Série de Fourier

 Ortogonalidade de funções trigonométricas

• Condição de ortogonalidade de funções (compare com o produto escalar)

Z b
ψn(x)ψm(x)dx = 0, n 6= m
a

As funções ψn e ψm são ditas ortogonais no intervalo entre a e b.

• Considerando que m ∈ Z e n ∈ Z, (convém deduzir, estudante)
Z π 
cos mxdx = 0 m 6= 0 Z π


 0 m 6= n
Z−ππ cos mx cos nxdx = π m = n 6= 0
−π
sin mxdx = 0

2π m=n=0

Z π −π (
Z π
cos mx sin nxdx = 0 0 m 6= n
sin mx sin nxdx =
−π −π π m = n 6= 0

• Logo, as funções 1, cos x, sin x, cos 2x, sin 2x, . . . formam um conjunto ortogonal no
intervalo entre −π e π.
TC01001 FET 7 / 34 Série de Fourier
 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

• Procedimento para determinação dos coeficientes an :

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

• Procedimento para determinação dos coeficientes an :
• Multiplica f (x) por cos mx

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

• Procedimento para determinação dos coeficientes an :
• Multiplica f (x) por cos mx
∞
a0 X
f (x) cos mx = cos mx + [an cos nx cos mx + bn sin nx cos mx]
2 n=1

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

• Procedimento para determinação dos coeficientes an :
• Multiplica f (x) por cos mx
∞
a0 X
f (x) cos mx = cos mx + [an cos nx cos mx + bn sin nx cos mx]
2 n=1

• Integra no intervalo de −π a π

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Série de Fourier:
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Função periódica: f (x + 2π) = f (x)

• Procedimento para determinação dos coeficientes an :
• Multiplica f (x) por cos mx
∞
a0 X
f (x) cos mx = cos mx + [an cos nx cos mx + bn sin nx cos mx]
2 n=1

• Integra no intervalo de −π a π
Z π Z π ∞  π π 
a0 X Z Z
f (x) cos mx dx = cos mx dx+ an cos nx cos mx dx + bn sin nx cos mx dx
−π 2 −π n=1 −π −π

TC01001 FET 8 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Procedimento para determinação dos coeficientes an (cont.):

• Avalia as integrais resultantes
Z π Z π
a0
f (x) cos mx dx = cos mx dx+
−π 2 −π
X∞  Z π Z π 
+ an cos nx cos mx dx + bn sin nx cos mx dx
n=1 −π −π
Z π ∞ π
a0 X Z
= cos mx dx + an cos nx cos mx dx+
2 −π n=1 −π
∞
X Z π
+ bn sin nx cos mx dx
n=1 −π

TC01001 FET 9 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Procedimento para determinação dos coeficientes an (cont.):

• Avalia as integrais resultantes (cont.)
Z π Z π ∞ π
a0 X Z
f (x) cos mx dx = cos mx dx + an cos nx cos mx dx+
−π 2 −π n=1 −π
X∞ Z π
+ bn sin nx cos mx dx
n=1 −π

TC01001 FET 10 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Procedimento para determinação dos coeficientes an (cont.):

• Avalia as integrais resultantes (cont.)
Z π Z π ∞ π
a0 X Z
f (x) cos mx dx = cos mx dx + an cos nx cos mx dx+
−π 2 −π n=1 −π
X∞ Z π
+ bn sin nx cos mx dx
n=1 −π

( Rπ
π a0
−π dx = a0 π, m=0
Z
2
f (x) cos mx dx = Rπ
−π am −π cos mx cos mx dx = am π, m>0

TC01001 FET 10 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Procedimento para determinação dos coeficientes an (cont.):

• Avalia as integrais resultantes (cont.)
Z π Z π ∞ π
a0 X Z
f (x) cos mx dx = cos mx dx + an cos nx cos mx dx+
−π 2 −π n=1 −π
X∞ Z π
+ bn sin nx cos mx dx
n=1 −π

( Rπ
π a0
−π dx = a0 π, m=0
Z
2
f (x) cos mx dx = Rπ
−π am −π cos mx cos mx dx = am π, m>0

• Determina expressão para an

Z π
1
an = f (x) cos nx dx
π −π
TC01001 FET 10 / 34 Série de Fourier
 Determinação dos coeficientes de Fourier

• Seguindo a procedimento semelhante,

Z π
1
bn = f (x) sin nx dx
π −π

TC01001 FET 11 / 34 Série de Fourier

 Determinação dos coeficientes de Fourier

• Seguindo a procedimento semelhante,

Z π
1
bn = f (x) sin nx dx
π −π

• Resumo:
• Seja f (x) uma função periódica com perı́odo 2π, i.e., f (x + 2π) = f (x)
• Pode-se representar f (x) por uma série trigonométrica do tipo
∞
a0 X
f (x) = + [an cos nx + bn sin nx]
2 n=1

• Onde os coeficientes an e bn são determinados pelas expressões

Z π Z π
1 1
an = f (x) cos nx dx bn = f (x) sin nx dx
π −π π −π

TC01001 FET 11 / 34 Série de Fourier

 Expansão de funções em série de Fourier

• Função quadrada periódica

(
−k −π < x < 0
f (x) = f (x + 2π) = f (x)
k 0<x<π

TC01001 FET 12 / 34 Série de Fourier

 Expansão de funções em série de Fourier

• Determinação dos coeficientes de Fourier

Z π π
1 1
Z
a0 = f (x) dx bn = f (x) sin nx dx
π −π π −π
Z 0 Z π Z 0
1 1 1
= (−k) dx + k dx = (−k) sin nx dx+
π −π π 0 π −π
Z π
=0 1
+ k sin nx dx
π 0
(
4k
1 π nπ , n é ı́mpar
Z
an = f (x) cos nx dx =
π −π 0, n é par
Z 0 Z π
1 1
= (−k) cos nx dx + k cos nx dx
π −π π 0
=0

TC01001 FET 13 / 34 Série de Fourier

 Expansão de funções em série de Fourier

• Representação usando série de Fourier

∞ ∞
a0 X X 4k
f (x) = + [an cos nx + bn sin nx] = sin(2n − 1)x
2 n=1 n=1
π(2n − 1)

• Serie truncada
N
X 4k
fN (x) = sin(2n − 1)x
n=1
π(2n − 1)

TC01001 FET 14 / 34 Série de Fourier

 Expansão de funções em série de Fourier

Contribuição de cada componente espectral

N =1 N =6

N =2 N = 26

N =3 N = 50

TC01001 FET 15 / 34 Série de Fourier

Convergência da série de Fourier
 Condições de Dirichlet

Teorema 1 Se uma função periódica f (x) de perı́odo 2π é limitada e contı́nua por

partes, e tem um número finito de máximos e mı́nimos em cada perı́odo, então a
série trigonométrica
∞
a0 X
+ [an cos nx + bn sin nx]
2 n=1

com
1 π
Z
an = f (x) cos nx dx n = 0, 1, 2, . . .
π −π
Z π
1
bn = f (x) sin nx dx n = 1, 2, 3, . . .
π −π
converge para f (x) onde f (x) é contı́nua, e converge para a média dos limites à
esquerda e à direita de f (x) nos pontos de descontinuidade.

TC01001 FET 17 / 34 Série de Fourier

 Condições de Dirichlet

• O teorema anterior estabelece condições suficientes para convergência

TC01001 FET 18 / 34 Série de Fourier

 Condições de Dirichlet

• O teorema anterior estabelece condições suficientes para convergência

• As condições mı́nimas necessárias para convergência não são conhecidas

TC01001 FET 18 / 34 Série de Fourier

 Condições de Dirichlet

• O teorema anterior estabelece condições suficientes para convergência

• As condições mı́nimas necessárias para convergência não são conhecidas
• Pode-se assumir que todas as funções periódicas de interesse prático podem ser
representadas pela série de Fourier

TC01001 FET 18 / 34 Série de Fourier

Série de Fourier de funções de perı́odo
arbitrário
 Mudança de intervalo

• A representação em série de
Fourier pode ser generalizada para
qualquer função periódica com
perı́odo 2L:

f (t) = f (t + 2L)

• Para tanto, deve-se fazer a seguinte

mudança de variável:
L
t= x
π 
L
f (t) =f x ≡ F (x)
π

TC01001 FET 20 / 34 Série de Fourier

 Mudança de intervalo

• A representação em série de • Dessa forma,

Fourier pode ser generalizada para
qualquer função periódica com f (t) =f (t + 2L)
   
perı́odo 2L: L L
f x =f x + 2L
π π
f (t) = f (t + 2L)
 
Lx + 2Lπ
=f
π
 
L
• Para tanto, deve-se fazer a seguinte =f (x + 2π)
π
mudança de variável:
F (x) =F (x + 2π)
L
t= x
π 
L
f (t) =f x ≡ F (x)
π

TC01001 FET 20 / 34 Série de Fourier

 Mudança de intervalo

• A representação em série de • Dessa forma,

Fourier pode ser generalizada para
qualquer função periódica com f (t) =f (t + 2L)
   
perı́odo 2L: L L
f x =f x + 2L
π π
f (t) = f (t + 2L)
 
Lx + 2Lπ
=f
π
 
L
• Para tanto, deve-se fazer a seguinte =f (x + 2π)
π
mudança de variável:
F (x) =F (x + 2π)
L
t= x
π 
L • Essa transformação de variável
f (t) =f x ≡ F (x)
π deve ser aplicada à série
trigonométrica e às expressões para
determinação de an e bn

TC01001 FET 20 / 34 Série de Fourier

 Mudança de intervalo

Série trigonométrica:
∞ h
a0 X π π i
+ an cos n t + bn sin n t
2 n=1
L L

Coeficientes de Fourier:
Z L
1  π 
an = f (t) cos n t dt n = 0, 1, 2, . . .
L −L L
Z L  π 
1
bn = f (t) sin n t dt n = 1, 2, 3, . . .
L −L L
Observações:
• π/L é a frequência fundamental de f (t)
• nπ/L é denominado harmônico de f (t) (múltiplo da frequência fundamental)

TC01001 FET 21 / 34 Série de Fourier

 Exercı́cios

1. Sobre os exemplos 1.3.1, 1.3.2, e 1.3.3 no livro de referência:

• Inicialmente, tentar resolver sozinho e, depois, comparar sua resolução com
a apresentada no livro
• Aproximar as funções usando a série de Fourier truncada (aumente
gradualmente o número de coeficientes usados) no Jupyter/Google Colab
• Observar a qualidade das reconstruções em função do número de
coeficentes usados e das caracterı́sticas das funções
2. Repita exercı́cio anterior para a função f (t) periódica a seguir:

−t

 −1 ≤ t < 1
f (t) = t2 − 4 1 ≤ t < 3

−1 3≤t<4


f (t) =f (t + 5)

TC01001 FET 22 / 34 Série de Fourier

 Série de Fourier de funções pares e ı́mpares

• Função par • Sendo g(t) uma função par,

Z ǫ Z ǫ
f (t) = f (−t) g(t)dt = 2 g(t)dt
−ǫ 0

• Função ı́mpar • Se f (t) for uma função par,

f (t) = −f (−t) 2
Z L  π 
an = f (t) cos n t dt bn = 0
L 0 L
• Produto entre funções:
• Se f (t) for uma função ı́mpar,
• par × par, função par
• ı́mpar × ı́mpar, função par 2
Z L  π 
• par × ı́mpar, função ı́mpar an = 0 bn = f (t) sin n t dt
L 0 L
• Sendo g(t) uma função ı́mpar,
Z ǫ
g(t)dt = 0
−ǫ
TC01001 FET 23 / 34 Série de Fourier
Série de Fourier de funções não-periódicas
em faixa limitada
 Representação de funções em faixa limitada

• Ideia principal: usar série de Fourier para representar funções não periódicas em
um determinado intervalo
• Estratégia: considerar que fora do intervalo de interesse a função se repete, ou
seja, é periódica
• Estratégia melhorada: a partir da função original, formar funções periódicas
pares (half-range cosine expansion) ou ı́mpares (half-range sine expansion)

TC01001 FET 25 / 34 Série de Fourier

 Estratégia melhorada - exemplo

• Função original • Função extendida par

• Função extendida ı́mpar

TC01001 FET 26 / 34 Série de Fourier

Série de Fourier complexa
 Série de Fourier complexa

Seja f (t) expressa pela série Onde

trigonométrica L
an bn 1
Z
−i(n Lπt
) dt
∞ h cn = + = f (t)e
a0 X  π   π i 2 2i 2L −L
f (t) = + an cos n t + bn sin n t an bn 1
Z L
2 n=1 L L c−n = − = f (t)e i(n L
πt
) dt
2 2i 2L −L
Usando a fórmula de Euler a0
c0 =
 π  ei(n Lπ t) + e−i(n Lπ t) 2
cos n t = Pode-se expressar a série como
L 2
 π  ei(n Lπ t) − e−i(n Lπ t) ∞
sin n t = i(n L
πt
)
X
L 2i f (t) = cn e
n=−∞
A série pode, então, ser expressa como
Onde
∞ h
a0 X
i(n L
πt
) −i(n L
πt
)
i
L
1
Z
f (t) = + cn e + c−n e −i(n L
πt
)
2 n=1 cn = f (t)e dt
2L −L

TC01001 FET 28 / 34 Série de Fourier

Propriedades da série de Fourier
 Propriedades da série de Fourier

• Teorema de Parceval ou teorema da energia

Z L ∞
1 2 1 2 1 X
[f (t)] dt = a0 + (a2n + b2n )
2L −L 4 2 n=1

• Soma de potências de inteiros recı́procas

 
f (x) =f (x + 2π) π 1 1 1
= 1 − + − + ...
( 4 3 5 7
−k − π < x < 0 ∞
f (x) = X (−1)n+1
k 0<x<π = (Leibniz, 1673)
  n=1
2n − 1
4k 1 1
= sin x + sin 3x + sin 5x + . . .
π 3 5
Para x = π/2,
π   
4k 1 1
f =k = 1 − + + ...
2 π 3 5
TC01001 FET 30 / 34 Série de Fourier
 Propriedades da série de Fourier

Mais alguns exemplos:

• Exemplo 01
2 ∞
π 1 1 X 1
= 1 + 2 + 2 + ··· =
8 3 5 n=1
(2n − 1)2

Dica: usar o teorema de Parseval no exemplo anterior.

• Exemplo 02

f (x) = x2 para − 1 < x < 1, e f (x + 2) = f (x)

Mostre que
∞ ∞
X (−1)n+1 π2 X 1 π2
= =
n=1
n2 12 n=1
n 2 6
∞ n+1 3 ∞
X (−1) π X 1 π4
= =
n=1
n3 32 n=1
n 4 90
TC01001 FET 31 / 34 Série de Fourier
 Propriedades da série de Fourier

• Integração da série de Fourier

• Premissa: integrais de senos e cossenos são simples
• Ideia: integrar os termos da série trigonométrica ao invés da função f (t)
∞ h
a0 X  π   π i
f (t) = + an cos n t + bn sin n t
2 n=1
L L
Z (X ∞ h
)
a0  π   π i
Z Z
f (t) dt = dt + an cos n t + bn sin n t dt
2 n=1
L L
∞  
a0 π π
Z X Z   Z  
= dt + an cos n t dt + bn sin n t dt
2 n=1
L L
∞
" #
sin n Lπ t cos n Lπ t



a0 X
= t+ an π − bn π +C
2 n=1
n L n L
∞  
a0 X an L  π  bn L  π 
= t+ sin n t − cos n t + C
2 n=1
nπ L nπ L
TC01001 FET 32 / 34 Série de Fourier
 Propriedades da série de Fourier

• Diferenciação da série de Fourier

• Ideia: mesmo procedimento usado para integração
• Cuidado: é necessário atenção com a convergência da série
• Análise: considere que f (x) é contı́nua no intervalo −L < t < L e que sua
derivada pode ser expandida em outra série de Fourier
∞ h
′ a′0X
′
 π 
′
 π i
f (t) = + an cos n t + bn sin n t
2 n=1
L L

onde
L
1 π
Z  
′ ′
an = f (t) cos n t dt
L −L L
1 nπ
= [f (L) − f (−L)] cos nπ + bn usar integração por partes !!
L L
′ 1 nπ
bn = [f (L) − f (−L)] sin nπ − an obtido pelo mesmo procedimento
L L
TC01001 FET 33 / 34 Série de Fourier
 Propriedades da série de Fourier

• Diferenciação da série de Fourier

• Análise (cont.): diferenciando diretamente a série trigonométrica para f (t)
( ∞
)
′ d a 0
X h  π   π i
f (t) = + an cos n t + bn sin n t
dt 2 n=1
L L
∞ h
X nπ  π  nπ  π i
= −an sin n t + bn cos n t
n=1
L L L L

onde
′ nπ
an = bn
L
′ nπ
bn = − an
L
A equivalência do resultado acima com o do slide anterior ocorre quando

f (L) = f (−L) condição head equals tail (a ser satisfeita!)

TC01001 FET 34 / 34 Série de Fourier