# Série de Fourier

Série de Fourier é uma forma de série trigonométrica usada para representar funções infinitas e periódicas complexas dos processos físicos, na forma de funções trigonométricas simples de senos e cossenos. Isto é, simplificando a visualização e manipulação de funções complexas. Foi criada em 1807 por Jean Baptiste Joseph Fourier (1768-1830).
A forma geral da série é:
$T(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[a_{n}\cdot \cos \left({\frac {n\pi t}{L}}\right)+b_{n}\cdot \operatorname {sen} \left({\frac {n\pi t}{L}}\right)\right]$
em que os coeficiente $a_{0}$, $a_{n}$ e $b_{n}$ são números que variam de acordo com a função que será representada, de período fundamental $2L$. Esses coeficientes são as amplitudes de cada onda em série, que são calculadas com as seguintes fórmulas:
$a_{0}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\,dt$
$a_{n}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\cos \left({\frac {n\pi t}{L}}\right)\,dt$
$b_{n}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\,\operatorname {sen} \left({\frac {n\pi t}{L}}\right)\,dt$
A Série de Fourier é importante na técnica de compactação digital, como por exemplo: para reproduzir músicas digitais por streaming, para ver imagens online de rápido carregamento, e no cancelamento de ruído nos fones de ouvido.

## História

A série surgiu na tentativa de Fourier solucionar um problema físico, que gerou novas fronteiras na matemática. Durante o estudo da propagação de calor em corpos sólidos. Admitindo que essa propagação deveria ser por ondas de calor, levando em consideração que a forma mais simples de uma onda é uma função senoidal. Assim Fourier demonstra através da transformada que qualquer função complexa, pode ser decomposta em uma combinação infinita de senoides, dividida como uma soma de senos e cossenos.
A ideia de decompor funções arbitrárias em termos de funções trigonométricas simples movimentou grandes nomes da matemática começando por volta de 1750 com L. Euler (1707-1783) e D. Bernoulli (1700-1782), seguindo com J. d'Alembert (1717-1783) e J. L. Lagrange (1736-1813).
Mais tarde Jean-Baptiste Joseph Fourier (1768-1830) estudou sistematicamente tais séries infinitas, na tentativa de resolver a equação do calor. Em 1811, em sua Théorie mathématique de la chaleur (Teoria matemática de condução do calor), Fourier explicitou os coeficientes de tais séries (que ficaram conhecidos como coeficientes de Fourier, embora Euler já conhecesse o formato dos mesmos) e escreveu as séries de senos e cossenos de várias funções. De um ponto de vista moderno, os resultados de Fourier, embora muito importantes a forma da série que recebeu o seu nome, são informais, em boa parte devido à falta de uma definição concisa de funções e integrais até o início do século XIX.
P. G. Dirichlet (1805-1859) foi um dos primeiros a reconhecer que nem toda função poderia ser representada por uma série de Fourier (fato que Fourier acreditava), obtendo uma condição suficiente para a validade da representação a partir da série estudada. Em um trabalho de 1829, Dirichlet dá a primeira demonstração rigorosa de que a série de Fourier de uma função f converge, em cada ponto x, para a média aritmética dos limites laterais de f nesse ponto. Nesse trabalho Dirichlet dá origem ao conceito de função como hoje se é conhecido.
Aparentemente por influência de Dirichlet, G. B. Riemann (1826-1866) se interessou pelo estudo das séries trigonométricas, sendo levado a estudar a integral que leva hoje o seu nome e publicando em 1854 um trabalho intitulado "Sobre a representação de funções por meio de séries trigonométricas".
Em 1876 du Bois-Reymond (1818-1896) construiu função cuja série de Fourier divergia em um dado ponto, e mais tarde ele mesmo construiu uma função cuja série divergia em um conjunto denso. Exemplos mais simples foram dados por L. Fejér (1880 -1959) em 1909. Vale citar também que em 1861 K. Weierstrass (1815-1897) deu o primeiro exemplo de função contínua sem derivada em ponto algum, sendo tal função definida por uma série trigonométrica que converge uniformemente (portanto uma série de Fourier).
Não esqueçamos de citar G. Cantor (1845-1918), o qual teve grande influência pelo trabalho de Dirichlet e investigou o problema da unicidade da representação de funções por séries trigonométricas. Tais influências foram decisivas para a definição de números reais como sequência de números racionais e para a criação da Teoria dos Conjuntos, o que mostra o quão importante para o desenvolvimento da fundamentação teórica da matemática foi a teoria das séries de Fourier, podendo esta, então, ser considerada uma das teorias mais importantes da Análise.
Embora o objetivo inicial de Joseph Fourier fosse resolver a equação do calor, depois que o método foi estudado e encontrado ele foi sendo usado para resolver muitos problemas matemáticos e físicos e, em especial, os que continham equações diferenciais. Esse estudo, conhecido hoje por Série de Fourier, tem aplicação direta nas áreas da engenharia elétrica, análise de vibrações, acústica, óptica, processamento de sinais, processamento de imagens, econometria. De encontro com esse aspecto está que, durante a elaboração da Série de Fourier, Joseph Fourier não estava focado em entender o calor enquanto entidade física (explicitando que não era o foco dele, mas que na pressão de dar uma explicação ele diria que a condução do calor se daria como radiação) mas para ele o importante era descrever matematicamente o comportamento do calor.

## Funções ortogonais

Um conjunto ortogonal de funções é uma generalização de um conjunto ortogonal de vetores, seja um conjunto ortogonal de vetores um conjunto de vetores perpendiculares entre si. As propriedades de um conjunto ortogonal de funções são derivadas por analogia dos conjuntos ortogonais de vetores.
Seja V ou um certo V(k) um vetor com n coordenadas definimos a fórmula que relaciona a norma do vetor com cada coordenada como:
$\left\|V\right\|^{2}=\sum _{k=1}^{n}V(k)^{2}$ .
Definição 1: A condição para que dois vetores sejam ortogonais é que o produto interno entre ambos seja igual a zero, e o produto interno é definido como:
$(V_{1},V_{2})=\sum _{k=1}^{n}V_{1}(k)V_{2}(k)$
Precisamos definir o que é um conjunto ortonormal de vetores, para isto usaremos primeiramente um conjunto de 3 vetores ortogonais. Seja $V_{n}$ um conjunto de 3 vetores, ou seja, teremos $V_{1},V_{2},V_{3}$ onde cada vetor possui 3 coordenadas. Calculamos as componentes de um vetor que faz parte do conjunto ortonormal de vetores segundo a fórmula:
$\phi _{n}(k)=\left\|V_{n}\right\|^{-1}V_{n}(k)$
em que n representa o vetor e k representa a sua respectiva coordenada. Note que o vetor é unitário e sua direção não mudou, portanto, neste exemplo, o conjunto ortonormal de vetores são os 3 $\phi _{n}$ que representam cada vetor $V_{1},V_{2},V_{3}$.
Agora veja, podemos definir qualquer vetor, por exemplo $F$, do espaço tridimensional como uma combinação linear do conjunto ortonormal, ou seja, o conjunto ortonormal forma uma base do espaço. Portanto:
$F=c_{1}\phi _{1}(k)+c_{2}\phi _{2}(k)+c_{3}\phi _{3}(k)$.
Podemos encontrar o coeficiente $c_{n}$, basta usar as propriedades do produto interno:
seja $(\phi _{m},\phi _{n})=\delta _{mn}$, onde $\delta _{mn}$ é o delta de Kronecker. Portanto:
$C_{n}=(F,\phi _{n})=\sum _{k=1}^{3}F(k)\phi _{n}(k)$.
Supondo agora que temos uma função contínua por partes onde $f(k)$ representa a ordenada associada à abscissa k, se a função for definida para x dentro de um intervalo, agora não conseguimos mais obter a norma do vetor simplesmente usando a soma trivial, é natural pensar portanto que a soma do quadrado das componentes será feito por integração:
$\left\|f\right\|={\begin{Bmatrix}\int _{a}^{b}f(x)^{2}dx\end{Bmatrix}}^{\frac {1}{2}}$
usando o mesmo princípio definimos o produto interno de duas funções como:
$(f_{1},f_{2})=\int _{a}^{b}f_{1}(x)f_{2}(x)dx$
também podemos definir o conjunto ${\begin{Bmatrix}\phi _{n}(x)\end{Bmatrix}}$. O conjunto ortonormal sobre o intervalo pode ser calculado da mesma forma que em vetores, como segue: $\phi _{n}(x)=\left\|f_{n}\right\|^{-1}f_{n}(x)$
e o produto interno entre eles também, lembrando é claro que agora temos que cobrir todo o intervalo e não usamos mais soma e sim a integral dentro do intervalo:
$\int _{a}^{b}\phi _{m}(x)\phi _{n}(x)dx$.
A ideia agora é escrever a função f(x) como uma combinação linear do conjunto ortonormal de funções, ou seja, $f(x)=c_{1}\phi _{1}(x)+c_{2}\phi _{2}(x)+...+c_{n}\phi _{n}(x)$
lembrando que x precisa estar dentro do intervalo. Quando estávamos lidando com vetores podíamos descobrir os coeficientes da série, e com funções não será diferente, podemos usar a mesma ideia proposta quando estávamos trabalhando com vetores: obter o produto interno de ambos os lados envolvendo $\phi _{n}(x)$, obtendo desta forma o coeficiente $c_{n}$ através da definição de produto interno. Portanto:
$c_{n}=\int _{a}^{b}f(x)\phi _{n}(x)dx$ (produto interno entre as funções)
Agora já estamos aptos a definir a série de Fourier generalizada através das funções ortogonais e do produto interno entre a função e o conjunto de funções ortogonais que define cada coeficiente da série:
$f(x)=\sum _{n=1}^{\infty }\phi _{n}(x)\int _{a}^{b}f(\xi )\phi _{n}(\xi )d\xi$.

## Definição

Uma série trigonométrica é uma série da forma
$T(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[a_{n}\cdot \cos \left({\frac {n\pi t}{L}}\right)+b_{n}\cdot \operatorname {sen} \left({\frac {n\pi t}{L}}\right)\right]$
Seja $f$ uma função periódica de período $2L$, ou seja $f(t+2L)=f(t)$ para todo $t$, a qual satisfaz às seguintes condições, conhecidas como as condições de Dirichlet:
* A função é unívoca/injetora (de um para um), e contínua exceto em um número finito de descontinuidade ordinárias dentro do período 2L;
* A função tem um número finito de máximos e mínimos dentro do período 2L;
* A função é absolutamente integrável, ou seja, a integral $\int _{0}^{2L}|f(t)|\,dt$ converge;
Então define-se a Série de Fourier da função $f$ como a série trigonométrica dada pelos coeficientes:
$a_{0}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\,dt$
$a_{n}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\cos \left({\frac {n\pi t}{L}}\right)\,dt$, $n\geq 1$
$b_{n}={\frac {1}{L}}\int _{c}^{c+2L}f(t)\,\operatorname {sen} \left({\frac {n\pi t}{L}}\right)\,dt$ , $n\geq 1$
para n inteiro. Observamos aqui que, como $f$ periódica de período $2L$, o intervalo de integração pode ser qualquer intervalo de comprimento $2L$, sendo que geralmente são utilizados $[0,2L]$ ou $[-L,L]$.
Os coeficientes $a_{n}$, $n\geq 0$ e $b_{n}$, $n\geq 1$ são conhecidos como coeficientes de Fourier.

## Forma Harmônica

Também podemos expressar a série de Fourier da função $f(t)$ com período $T$ como
    $f(t)=A_{0}+\sum _{n=1}^{\infty }A_{n}\cdot \cos(w_{n}t-\theta _{n})$
em que
    $A_{0}={\frac {a_{0}}{2}}$, $A_{n}={\sqrt {a_{n}^{2}+b_{n}^{2}}}$ para $n\geq 1$
    $\theta _{n}=\arctan({\frac {b_{n}}{a_{n}}})$
    $w_{n}={\frac {2n\pi }{T}}$
A igualdade entre a forma harmônica e a forma trigonométrica se dá por meio da identidade trigonométrica do cosseno da diferença:
    $\cos(a-b)=\cos(a)\cos(b)+\operatorname {sen}(a)\operatorname {sen}(b)\,$
Aplicando a identidade na forma harmônica, tem-se que:
    $f(t)=A_{0}+\sum _{n=1}^{\infty }[A_{n}\cos(\theta _{n})\cos(w_{n}t)+A_{n}\sin(\theta _{n})\sin(w_{n}t)]$
Comparando os termos com os da representação trigonométrica, tem-se que:
    $A_{0}={\frac {a_{0}}{2}},a_{n}=A_{n}\cos(\theta _{n})$ e $b_{n}=A_{n}\sin(\theta _{n})$

## Forma complexa

Usando-se as expressões para as funções trigonométricas provenientes da fórmula de Euler: $e^{ix}=cos(x)+isen(x)$
$\operatorname {sen} z={\frac {e^{iz}-e^{-iz}}{2i}}=-i{\frac {e^{iz}-e^{-iz}}{2}}$
$\cos z={\frac {e^{iz}+e^{-iz}}{2}}$
Substituindo-as na forma harmônica da série de Fourier:
$f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[a_{n}{\frac {e^{iw_{n}t}+e^{-iw_{n}t}}{2}}-ib_{n}{\frac {e^{iw_{n}t}-e^{-iw_{n}t}}{2}}\right]$
em que
$w_{n}={\frac {n\pi }{L}}$
$f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[{\frac {a_{n}-ib_{n}}{2}}e^{iw_{n}t}+{\frac {a_{n}+ib_{n}}{2}}e^{-iw_{n}t}\right]$
Considerando a paridade das funções cosseno (par) e seno (ímpar) é possível notar que
$a_{-n}=a_{n}$
$b_{-n}=-b_{n}$ e $b_{0}=0$
em que foi usado
$cos(-x)=cos(x)$ e $sen(-x)=-sen(x)$
$w_{-n}={\frac {(-n)\pi }{L}}=-w_{n}$
Como notação é tomado: $c_{n}={\frac {a_{n}-ib_{n}}{2}}$ , tal que $c_{0}={\frac {a_{0}-ib_{0}}{2}}={\frac {a_{0}}{2}}$ dado por
$c_{0}={\frac {1}{2L}}\int _{-L}^{L}f(t)\,dt$
Para calcular os coeficientes com índices diferentes de zero utilizamos a fórmula de Euler:
${\frac {1}{2L}}\int _{-L}^{L}f(t)\,(cos(w_{n}t)-isen(w_{n}t))dt$
Logo,
$c_{n}={\frac {1}{2L}}\int _{-L}^{L}f(t)\,e^{-iw_{n}t}dt$
Para a expressão da função obtemos
$f(t)=c_{0}+\sum _{n=1}^{\infty }c_{n}\,e^{iw_{n}t}+\sum _{n=1}^{\infty }c_{-n}\,e^{-iw_{n}t}$
que pode ser simplificada em
$f(t)=\sum _{n=-\infty }^{\infty }c_{n}\,e^{iw_{n}t}$

## Exemplos

| Domínio do tempo $f(x)$ | Domínio da Frequência ${\begin{aligned}&a_{0}\\&a_{n}\quad {\text{for }}n\geq 1\\&b_{n}\quad {\text{for }}n\geq 1\end{aligned}}$ | Observações | Referência |
| $f(x)=A\left|\sin \left({\frac {2\pi }{T}}x\right)\right|\quad {\text{for }}0\leq x<T$ | ${\begin{aligned}a_{0}=&{\frac {4A}{\pi }}\\a_{n}=&{\begin{cases}{\frac {-4A}{\pi }}{\frac {1}{1-n^{2}}}&\quad n{\text{ par}}\\0&\quad n{\text{ ímpar}}\end{cases}}\\b_{n}=&0\\\end{aligned}}$ | Onda completa retificada (seno) | $^{:p.193}$ |
| $f(x)={\begin{cases}A\sin \left({\frac {2\pi }{T}}x\right)&\quad {\text{for }}0\leq x<T/2\\0&\quad {\text{for }}T/2\leq x<T\\\end{cases}}$ | ${\begin{aligned}a_{0}=&{\frac {2A}{\pi }}\\a_{n}=&{\begin{cases}{\frac {-2A}{\pi }}{\frac {1}{1-n^{2}}}&\quad n{\text{ par}}\\0&\quad n{\text{ ímpar}}\end{cases}}\\b_{n}=&{\begin{cases}{\frac {A}{2}}&\quad n=1\\0&\quad n>1\end{cases}}\\\end{aligned}}$ | Meia onda retificada (seno) | $^{:p.193}$ |
| $f(x)={\begin{cases}A&\quad {\text{for }}0\leq x<D\cdot T\\0&\quad {\text{for }}D\cdot T\leq x<T\\\end{cases}}$ | ${\begin{aligned}a_{0}=&2AD\\a_{n}=&{\frac {A}{n\pi }}\sin \left(2\pi nD\right)\\b_{n}=&{\frac {2A}{n\pi }}\left(\sin \left(\pi nD\right)\right)^{2}\\\end{aligned}}$ | $0\leq D\leq 1$ | |
| $f(x)={\frac {Ax}{T}}\quad {\text{for }}0\leq x<T$ | ${\begin{aligned}a_{0}=&A\\a_{n}=&0\\b_{n}=&{\frac {-A}{n\pi }}\\\end{aligned}}$ | | $^{:p.192}$ |
| $f(x)=A-{\dfrac {Ax}{T}}\quad {\text{for }}0\leq x<T$ | ${\begin{aligned}a_{0}=&A\\a_{n}=&0\\b_{n}=&{\frac {A}{n\pi }}\\\end{aligned}}$ | | $^{:p.192}$ |
| $f(x)={\frac {4A}{T^{2}}}\left(x-{\frac {T}{2}}\right)^{2}\quad {\text{for }}0\leq x<T$ | ${\begin{aligned}a_{0}=&{\frac {2A}{3}}\\a_{n}=&{\frac {4A}{\pi ^{2}n^{2}}}\\b_{n}=&0\\\end{aligned}}$ | | $^{:p.193}$ |
Para a função $g(t)$ $2\pi$-periódica dada por
$g(t)={\begin{cases}-1,\,-\pi <t<0\\+1,\,0\leq t<\pi \end{cases}}$
encontramos os seguintes coeficientes de Fourier
$a_{0}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(t)\,dt=0$ ,
$a_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(t)\cos \left({\frac {n\pi t}{\pi }}\right)\,dt=0$
$b_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(t)\,\operatorname {sen} \left({\frac {n\pi t}{\pi }}\right)\,dt={\frac {2}{\pi }}\cdot {\frac {1-(-1)^{n}}{n}}$
Logo, a série de Fourier de $g(t)$ é dada por
$\sum _{n=1}^{\infty }{\frac {2}{\pi }}\cdot {\frac {1-(-1)^{n}}{n}}\cdot \operatorname {sen}(nt)$
A função onda quadrada escrita na forma trigonométrica da Série de Fourier é dada por:
${\displaystyle h(t)={\frac {4}{\pi }}(sen(\pi \,t)+{\tfrac {1}{3}}sen(3\pi \,t)+{\tfrac {1}{5}}sen(5\pi \,t)+...)}$ com ${\displaystyle \omega n=n\pi }$
Encontrando os seguintes coeficientes: ${\displaystyle C_{0}={\frac {a_{0}}{2}}=0}$ e ${\displaystyle C_{n}={\frac {0-i2{\tfrac {1-(-1)^{n}}{\pi n}}}{2}}=i{\frac {(-1)^{n}-1}{\pi n}},n\neq 0}$
A função onda quadrada pode ser escrita na forma exponencial como
${\displaystyle h(t)=...+{\frac {2}{3\pi }}ie^{-i3\pi t}+{\frac {2}{\pi }}ie^{-i\pi t}-{\frac {2}{\pi }}ie^{i\pi t}-{\frac {2}{3\pi }}ie^{i3\pi t}-...}$

## Diagramas de espectro

Os coeficientes $c_{n}$ da expansão de uma função em Série de Fourier Complexa são números complexos. É comum representar os $c_{n}$ na forma de módulo e fase, isto é:
    $c_{n}=|c_{n}|e^{i\phi _{n}}$
Pode-se representar $c_{n}$ através de dois diagramas: um para as suas amplitudes (conhecido como espectro de amplitudes) e outro para as suas fases (conhecido como espectro de fases):
* Espectro de amplitudes: é um diagrama onde grafa-se os valores das amplitudes $|C_{n}|$ dos coeficientes de Fourier versus $n\omega _{0}$, isto é, um gráfico da forma: $|C_{n}|$ × $n\omega _{0}$.
* Espectro de fase: é um diagrama onde grafa-se os valores das fases $\phi _{n}$ dos coeficientes de Fourier versus $n\omega _{0}$, isto é, um gráfico da forma: $\phi _{n}$ × $n\omega _{0}$.

## Convergência

Nas aplicações da engenharia, normalmente se presume que a Série de Fourier converge em todos os pontos, exceto nas descontinuidades. Isso se deve ao fato de que essas funções são mais comportadas do que as que os matemáticos fornecem em contrapartida. Em particular, se $s$é contínua e derivada de $f(x)$, que pode apresentar descontinuidade, é integrável ao quadrado, então a Série de Fourier $f$converge para $f(x)$. Se uma função é integrável ao quadrado no intervalo $[xo,xo+T]$, então a série de Fourier converge para a função em praticamente todos os pontos. A convergência das séries de Fourier também depende do número finito de máximos e mínimos em uma função que é popularmente conhecida como uma das condições de Dirichlet para as séries de Fourier. É possível definir coeficientes de Fourier para funções ou distribuições mais gerais, nesses casos a convergência em norma ou convergência fraca é usualmente de interesse.

## Propriedades
### Teorema de Parseval

Considerando uma função $f(t)$, com período $T$ e representável por uma série de Fourier, a seguinte identidade é válida:
    ${\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}\,dt=\sum _{n=-\infty }^{\infty }|C_{n}|^{2}$
Sendo $c_{n}$ os coeficientes complexos da série analisada. Esta integral se relaciona com o conceito de potência média de um sinal.

#### Demonstração do Teorema de Parseval

Considerando que
${\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}\,dt={\frac {1}{T}}\int _{0}^{T}f(t){\overline {f(t)}}\,dt$
Como $f(t)=\sum _{n=-\infty }^{\infty }C_{n}e^{iw_{n}t}$, temos
${\overline {f(t)}}={\overline {\sum _{n=-\infty }^{\infty }C_{n}e^{iw_{n}t}}}={\sum _{n=-\infty }^{\infty }{\overline {C_{n}e^{iw_{n}t}}}}={\sum _{n=-\infty }^{\infty }{\overline {C_{n}}}e^{-iw_{n}t}}$
Substituindo esta expressão para ${\overline {f(t)}}$na definição de Potência Média dada por
${\overline {P_{f}}}={\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}dt$
temos,
${\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}dt={\frac {1}{T}}\int _{0}^{T}f(t){\overline {f(t)}}\,dt={\frac {1}{T}}\int _{0}^{T}f(t){\Bigl{\sum _{n=-\infty }^{\infty }{\overline {C_{n}}}e^{-iw_{n}t}}{\Bigr ]}dt$
${\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}dt={\frac {1}{T}}{\sum _{n=-\infty }^{\infty }{\Bigl{\overline {C_{n}}}\int _{0}^{T}f(t)e^{-iw_{n}t}}dt{\Bigr ]}$
Como $C_{n}={\frac {1}{T}}\int _{0}^{T}f(t)e^{-iw_{n}t}dt$,temos:
${\frac {1}{T}}\int _{0}^{T}|f(t)|^{2}dt=\sum _{n=-\infty }^{\infty }{\overline {C_{n}}}C_{n}=\sum _{n=-\infty }^{\infty }|C_{n}|^{2}$

### Funções pares e ímpares

As séries de Fourier transparecem em seus coeficientes a paridade das funções:
* Se f(t) for uma função ímpar, $a_{0}=0$ e $a_{n}=0$ para todo $n\geq 1$ e, além disso $b_{n}={\frac {2}{L}}\int _{0}^{L}f(t)\,\operatorname {sen} \left({\frac {n\pi t}{L}}\right)\,dt$ para $n\geq 1$.
* Se f(t) for uma função par, $b_{n}=0$ para todo $n\geq 1$ e, além disso, $a_{n}={\frac {2}{L}}\int _{0}^{L}f(t)\,\operatorname {cos} \left({\frac {n\pi t}{L}}\right)\,dt$ para $n\geq 0$
Algumas considerações sobre essas funções são importantes:
* A soma e a diferença de duas funções pares é par;
* A soma e a diferença de duas funções ímpares é ímpar;
* A soma e a diferença de uma função par e uma função ímpar não é nem par nem ímpar;
* O produto e o quociente de duas funções pares é par;
* O produto e o quociente de duas funções ímpares é par;
* O produto e o quociente de uma função par e uma função ímpar é ímpar.
Essas afirmações são úteis para o estudo da Transformada de Fourier, pois auxiliam no processo de integração.

### Fenômeno de Gibbs

A convergência das somas parciais de séries de Fourier de funções suaves por partes em torno de um salto, encontra-se oscilações cujas amplitudes não convergem para zero. Tal fenômeno envolve o fato que mesmo com a adição de mais termos a soma parcial, a ultrapassagem nos saltos das descontinuidades não desaparece. Temos a convergência ponto a ponto, no entanto, ao avaliarmos a diferença entre a função e a soma parcial, sempre encontramos um ponto onde esse valor é aproximadamente 8,9% da amplitude do salto.

## Aplicações

As séries de Fourier são capazes de representar uma família de funções periódicas envolvendo tanto funções contínuas como não contínuas, diferentemente das séries de Taylor, que só representam funções contínuas e deriváveis. Algumas das funções representáveis podem ter significado físico, como os sinais musicais ou elétricos. Essas séries também são amplamente utilizadas na resolução de equações diferenciais parciais, como a equação da onda, do calor e problemas envolvendo a Equação de Laplace.
Para funções não periódicas a série de Fourier não está definida. Faz-se, então, o uso da transformada de Fourier, que possui uma possibilidade de aplicação mais ampla.
O uso da série nos sistemas de compressão como MP3, JPG, AVI, mantendo os sons e imagens com boa qualidade, ocupando apenas uma fração do espaço de armazenamento ou meio de transmissão que seus originais sem perdas necessitam. Em outras palavras, a série tornou imagens e músicas digitais em coisas práticas e, com rápido compartilhamento.
Usando as ferramentas e técnicas da espectroscopia por exemplo, astrônomos deduzem a composição química de uma estrela pela análise das componentes frequenciais (espectro) a partir das luzes emitidas pela estrela; engenheiros otimizam o projeto de um sistema de telecomunicações usando informações sobre as componentes espectrais de um sinal de dados que um sistema transportará.
Na medicina e biologia, em estudos de morfologia dentofacial, a aplicação de funções de Fourier para comparar a forma da arcada dentária, os perfis faciais e a forma mandibular entre diferentes populações humanas e gêmeos, forneceu novas informações sobre influências genéticas e ambientais na variação destas estruturas. Diferenças no tamanho e forma dos arcos dentários entre diferentes populações humanas forneceram evidências indiretas das diferenças genéticas subjacentes entre os grupos. Enquanto os aborígines australianos tendiam a ter maior arcada dentária maxilar e mandibular do que japoneses ou caucasianos, as diferenças mais significativas eram aspectos da forma do arco. os caucasianos e os japoneses exibiram arcos maxilares e mandibulares mais afunilados do que os australianos. Além disso, comparações de amplitudes de Fourier descrevendo arcos dentários em gêmeos forneceram evidências da influência de fatores genéticos no tamanho do arco e na relação entre a profundidade / largura do arco.

### Sinais Discretos

Vamos considerar que o espectro da função $f(t)$ é composto apenas por frequências inferiores a $w_{c}$, onde $w_{c}$é chamado de frequência de corte. Mostraremos que se conhecermos apenas os valores de $f(t)$ para $t=kT$, k ∈ Z, onde T é o período de amostragem e $w_{a}:=2w_{c}$ é a frequência de amostragem, então podemos reconstruir exatamente $f(t)$em todos instantes de tempo. Considere $f(t)$ uma função real, deﬁniremos $f_{T}(t)$ uma versão discretizada deste sinal da seguinte forma:
$f_{T}(t)=\sum _{k=-\infty }^{\infty }f(kT)\delta (t-kT)$
Assim $f_{T}(t)$ é um trem de Dirac’s cujas amplitudes coincidem com o valor da função $f(t)$ nos pontos de amostragem $kT$.
A ﬁm de calcularmos a transforma de Fourier de $f_{T}(t)$, observamos que:
$f_{T}(t)=\sum _{k=-\infty }^{\infty }f(kT)\delta (t-kT)=\sum _{k=-\infty }^{\infty }f(t)\delta (t-kT)=f(t)\sum _{k=-\infty }^{\infty }\delta (T-kT)=f(t)\delta _{T}(t)$
onde $\delta _{T}(t)=\sum _{k=-\infty }^{\infty }\delta (t-kT)$é uma função periódica cuja série de Fourier é dada por:
$\delta _{T}(t)=\sum _{k=-\infty }^{\infty }\delta (t-kT)=\sum _{k-\infty }^{\infty }C_{n}e^{iw_{n}t}$ e $C_{n}={\frac {1}{T}}\int _{-{\frac {T}{2}}}^{\frac {T}{2}}\delta _{T}(t)e^{-iw_{n}t}dt={\frac {1}{T}}$ assim,
$\delta _{T}(t)={\frac {1}{T}}\sum _{n=-\infty }^{\infty }e^{iw_{n}t}$ e, portanto:
$f_{T}(t)=f(t)\delta _{T}(t)=f(t){\frac {1}{T}}\sum _{n=-\infty }^{\infty }e^{iw_{n}t}={\frac {1}{T}}\sum _{n=-\infty }^{\infty }f(t)e^{iw_{n}t}$ e finalmente:
$F_{T}(w)={\mathcal {F}}\{f_{T}(T)\}={\frac {1}{T}}{\mathcal {F}}\{\sum _{n=-\infty }^{\infty }f(t)e^{iw_{n}t}\}={\frac {1}{T}}\sum _{n=-\infty }^{\infty }F(w-w_{n})$ onde se usou a propriedade do deslocamento no eixo $w$.
Observamos que se a frequência de amostragem $w_{a}$for superior a $2w_{c}$, então $F_{t}(w)={\frac {1}{T}}F(w)$ no intervalo $[-w_{c},w_{c}]$ e, portanto, toda a informação de $f(t)$ é preservada. De fato, neste caso, podemos escrever:
$f(t)={\frac {1}{2\pi }}\int _{-\infty }^{\infty }F(x)e^{iwt}dw={\frac {1}{2\pi }}\int _{-w_{c}}^{w_{c}}F(x)e^{iwt}dw={\frac {1}{2\pi }}\int _{-w_{c}}^{w_{c}}TF_{T}(w)e^{iwt}dw$
Como $F_{T}(w)$ pode ser calculada apenas com base nos pontos de amostragem, $f(t)$ pode ser reconstruída. Se $w_{a}<2w_{c}$, então existe superposição espectral, o que impede a reconstrução da $f(t)$. Este resultado é conhecido como teorema da amostragem de Nyquist-Shannon ou teorema cardinal da interpolação.

### Passagem do discreto para o contínuo

Suponhamos uma função não periódica. Definimos:
${\begin{cases}f_{T}(t)=f(t),&{\text{-T/2 }}\leq t<{\text{ T/2}}\\f_{T}(t)=f(t+T),&\forall t\in \mathbb {R} \end{cases}}$
Onde:
* $f(t)$ é uma função não periódica;
* $f_{T}(t)$ é uma função T-periódica que coincide com $f(t)$ no intervalo citado.
Dessa forma, é possível representar a função $f(t)$ em Série de Fourier, pois a função $f_{T}(t)$ carrega informação sobre a $f(t)$.
Exemplificando, temos a função não periódica $f(t)=e^{-|t|}$:
* Queremos aproximar o máximo $f_{T}(t)$ de $f(t)$, ou seja, obter o limite quando $T\longrightarrow \infty$;
* $T$ representa o período da função $f_{T}(t)$, dessa forma quando $T$ cresce a frequência fundamental $w_{n}t$ decresce.
Assim, a função $f_{T}(t)$ possui Série de Fourier da forma:
$f_{T}(t)=\textstyle \sum _{n=-\infty }^{\infty }\displaystyle C_{n}e^{iw_{n}t}$
Cujo $C_{n}$ é possível calcular sabendo que $C_{n}={\frac {1}{T}}\int \limits _{-T/2}^{T/2}e^{-|t|}e^{-iw_{n}t}dt={\frac {1}{T}}\int \limits _{-T/2}^{T/2}e^{-|t|}(\cos(w_{nt})-i\operatorname {sen}(w_{n}t))dt={\frac {2}{T}}\int \limits _{0}^{T/2}e^{-|t|}\cos(w_{nt})dt$
$C_{n}={2 \over T}{\frac {1-(-1)^{n}e^{-T \over 2}}{1+w_{n}^{2}}}$
Como quando $T$ cresce a frequência fundamental $w_{n}t$ decresce num gráfico $|TC_{n}|vsW_{n}$ as raias ficam cada vez mais próximas num diagrama de espectro, ou seja, é possível observar um aumento da densidade das linhas espectrais, e podemos escrever $C_{n}$ como:
$C_{n}={\frac {1}{T}}\int \limits _{-T/2}^{T/2}f_{T}(\tau )e^{-iw_{n}\tau }d\tau$
Agora, podemos definir:
$F_{T}(w)=\int \limits _{-T/2}^{T/2}f({\tau })e^{-iw\tau }d\tau$
E também $f_{T}(t)$ em termos de $F_{T}(w)$:
$f_{T}(t)={\frac {1}{2\pi }}\textstyle \sum _{n=-\infty }^{\infty }\displaystyle F_{T}(w_{n})e^{iw_{n}t}\Delta w$
Que fazendo o limite de $T\longrightarrow \infty$ a soma de Riemann converge para a integral $f(t)={\frac {1}{2\pi }}\int \limits _{-T/2}^{T/2}F(w)e^{iwt}dw$ , onde $F(w)=\int \limits _{-\infty }^{\infty }f(t)e^{-iwt}dt$

## Aplicações na Matemática

A Série de Fourier é utilizada na área de Matemática Aplicada, sendo amplamente utilizada na transformada de Fourier. Essa transformada é de extrema importância, pois é capaz de diminuir, por exemplo, a quantidade de informações que precisamos guardar para recriar uma onda. Desse modo, conseguimos recriar a onda apenas guardando as frequências que a compõem, em vez de armazenar cada ponto da amplitude dessa onda. Como essa transformada é de suma importância, discutiremos sobre algumas propriedades importantes dessa transformada:

### Linearidade

Uma das mais importantes propriedades da transformada de Fourier é a linearidade.
* Demonstração da Linearidade:
Considere que: ${\mathcal {F}}\{f(t)\}={\mathcal {F}}\{ax(t)+by(t)\}$
Assim, ${\mathcal {F}}\{ax(t)+by(t)\}=\int \limits _{-\infty }^{\infty }[ax(t)+by(t)]e^{-itw}$
${\mathcal {F}}\{ax(t)+by(t)\}=a\int \limits _{-\infty }^{\infty }x(t)e^{-iwt}dt\quad +b\int \limits _{-\infty }^{\infty }y(t)e^{-iwt}dt$
utilizando a definição acima, temos que: $F[ax(t)+by(t)]=aX(w)+bY(w)$sendo X(w) e Y(w) as transformadas de Fourier.

### Transformada da Derivada

Se f(t) é uma função diferenciável, sua transformada de Fourier é
${\mathcal {F}}\{f'(t)\}=iwF(w)$
utilizando a integração por partes temos:
${\mathcal {F}}\{f'(t)\}=iw\int \limits _{-\infty }^{\infty }f(t)e^{-iwt}dt$
então ${\mathcal {F}}\{f'(t)\}=iwF(w)$

### Transformada da integral

Dada uma função integrável $f(t)$, na qual sua transformada de Fourier $F(w)$ satisfaça $F(0)=0$, e considerando uma $g(t)=\int \limits _{-\infty }^{t}f(\tau )d\tau$ e o teorema fundamental do cálculo em que $g'(t)=f(t)$, podemos ver que:
${\mathcal {F}}\{g'(t)\}={\mathcal {F}}\{f(t)\}=F(w)$
É observável que: $\lim _{t\to \infty }g(t){\displaystyle =\int \limits _{-\infty }^{\infty }f(\tau )d\tau }={\displaystyle \int \limits _{-\infty }^{\infty }f(\tau )e^{i\cdot 0\cdot \tau }d\tau =F(0)}=0$e que $\lim _{t\to \infty }g(t){\displaystyle =\int \limits _{-\infty }^{-\infty }f(\tau )d\tau }=0$Assim, ao usar a propriedade vista acima da transformada de Fourier da derivada podemos obter:
${\mathcal {F}}\{g'(t)\}=iw{\mathcal {F}}\{g(t)\}$ e assim
$F(w)=iw{\mathcal {F}}\{\int \limits _{-\infty }^{t}f(\tau )d\tau \}$
Chegando, dessa forma, à propriedade em si:
${\mathcal {F}}\{\int \limits _{-\infty }^{t}f(\tau )d\tau \}={\frac {1}{iw}}F(w)$

### Simetria ou Dualidade

Considerando uma função f(t), temos $f(-w)={\frac {1}{2\pi }}{\mathcal {F}}\{F(t)\}$
pela definição de Fourier $f(t)={\frac {1}{2\pi }}\int \limits _{-\infty }^{\infty }F(w)e^{iwt}dw$, calculamos então a f(w), para isso trocamos t e w
$f(w)={\frac {1}{2\pi }}\int \limits _{-\infty }^{\infty }F(t)e^{iwt}dt$
Então, podemos afirmar que
$f(-w)={\frac {1}{2\pi }}\int \limits _{-\infty }^{\infty }F(t)e^{-iwt}={\frac {1}{2\pi }}{\mathcal {F}}\{F(t)\}$
Essas são algumas das propriedades que auxiliam no cálculo da Transformada de Fourier

### Inversão temporal

Uma propriedade muito importante que é muito visível no cotidiano é a propriedade da inversão temporal, a qual demonstra os efeitos de quando ouvimos a música de trás para frente, por exemplo. Ao considerarmos uma $f(t)$e uma $F(w)$, temos que:
${\mathcal {F}}\{f(-t)\}=\int \limits _{-\infty }^{\infty }f(-t)e^{-itw}dt$
Ao considerarmos a mudança de variável $\tau =-t$, temos que:
${\mathcal {F}}\{f(-t)\}=\int \limits _{\infty }^{-\infty }f(\tau )e^{itw}(-d\tau )$
$=\int \limits _{-\infty }^{\infty }f(\tau )e^{itw}(d\tau )$
$=\int \limits _{-\infty }^{\infty }f(\tau )e^{-it(-w)}(d\tau )$
$=F(-w)$
E assim, temos a propriedade, em que ${\mathcal {F}}\{f(-t)\}=F(-w)$

### Deslocamento na frequência

Dada uma função f(t) e sua transformada de Fourier F($\omega$) , então
${\mathcal {F}}\{{e^{at}}f(t)\}=F(\omega +ia)$
Demonstração:
$\qquad \qquad =\int \limits _{-\infty }^{\infty }f(t)e^{at}e^{-i{\omega }t}dt$
$\qquad \qquad =\int \limits _{-\infty }^{\infty }f(t)e^{(a-i{\omega })t}dt$
$\qquad \qquad =\int \limits _{-\infty }^{\infty }f(t)e^{-i(ia+{\omega })t}dt$
$\qquad \qquad =F({\omega }+ia)$

### Deslocamento no tempo

Dada uma função f(t) e sua transformada de Fourier F($\omega$) , então
${\mathcal {F}}\{f(t-a)\}=e^{-ia\omega }F(\omega )$
Demonstração
$\qquad \qquad =\int _{-\infty }^{\infty }f(t-a)e^{-i\omega t}dt$
$\qquad \qquad =\int _{-\infty }^{\infty }f(s)e^{-i\omega (s+a)}ds$
$\qquad \qquad =\int _{-\infty }^{\infty }f(s)e^{-i\omega a}{e^{-i\omega s}}ds$
$\qquad \qquad =e^{-i\omega a}\int _{-\infty }^{\infty }f(s){e^{-i\omega s}}ds$
$\qquad \qquad =e^{-i\omega a}F(\omega )$

### Conjugação

Dada uma função real f(t) e sua transformada de Fourier $F(\omega )$, então
${\overline {F(\omega )}}=F(-\omega )$
Demonstração
${\overline {F(\omega )}}={\overline {\int _{-\infty }^{\infty }f(t)e^{-i\omega t}dt}}$
$=\int _{-\infty }^{\infty }f(t){\overline {e^{-i\omega t}}}dt,\qquad pois\ {\overline {f(t)}}=f(t)$
$=\int _{-\infty }^{\infty }f(t)e^{i\omega t}dt$
$=\int _{-\infty }^{\infty }e^{-i(-\omega )t}dt$
$=F(-\omega )$

### Modulação

Dada uma função f(t) e sua transformada de Fourier $F(w)$, então:
${\mathcal {F}}\{f(t)cos(w_{0}t)\}={\frac {1}{2}}F(w-w_{0})+{\frac {1}{2}}F(w+w_{0})$
${\mathcal {F}}\{f(t)sen(w_{0}t)\}={\frac {i}{2}}F(w+w_{0})-{\frac {i}{2}}F(w-w_{0})$
Para a função cosseno é factível que a proposição é válida porque:
${\begin{aligned}{\mathcal {F}}\{f(t)cos(w_{0}t)\}&={F}\{f(t){\frac {e^{iw_{0}t}+e^{-iw_{0}t}}{2}}\}\\{\mathcal {F}}\{f(t)cos(w_{0}t)\}&=\int \limits _{-\infty }^{\infty }f(t){\frac {e^{iat}+e^{-iat}}{2}}e^{-iwt}d(t)\\{\mathcal {F}}\{f(t)cos(w_{0}t)\}&={\frac {1}{2}}\int \limits _{-\infty }^{\infty }f(t){e^{-i(w-w_{0})t}dt+{\frac {1}{2}}\int \limits _{-\infty }^{\infty }f(t){e^{-i(w_{0}-w)t}}}dt\\{\mathcal {F}}\{f(t)cos(w_{0}t)\}&={\frac {1}{2}}F(w-w_{0})+{\frac {1}{2}}F(w+w_{0})\\\end{aligned}}$
No que se refere à função seno, pode-se provar o teorema por meio da seguinte demonstração:
${\begin{aligned}{\mathcal {F}}\{f(t)sen(w_{0}t)\}&=\int _{-\infty }^{\infty }f(t){\biggl{\frac {e^{iw_{0}t}-e^{-iw_{0}t}}{2i}}{\biggr )}e^{-iwt}dt\\{\mathcal {F}}\{f(t)sen(w_{0}t)\}&={\frac {1}{2i}}{\biggl\int _{-\infty }^{\infty }f(t)e^{i(w-w_{0})t}dt-\int _{-\infty }^{\infty }f(t)e^{i(w+w_{0})t}dt{\biggr ]}\\{\mathcal {F}}\{f(t)sen(w_{0}t)\}&={\frac {1}{2i}}{\bigglF(w-w_{0})-F(w+w_{0}){\biggr ]}\\{\mathcal {F}}\{f(t)sen(w_{0}t)\}&={\frac {F(w-w_{0})}{2i}}-{\frac {F(w+w_{0})}{2i}}\\{\mathcal {F}}\{f(t)sen(w_{0}t)\}&=-{\frac {i}{2}}F(w-w_{0})+{\frac {i}{2}}F(w+w_{0})\\\end{aligned}}$

### Mudança de Escala

Dada uma função f(t) e sua transformada de Fourier $F(w)$, então
${\mathcal {F}}\{f(at)\}={\frac {1}{|a|}}F({\frac {w}{a}})$ $\forall a\neq 0$
Demonstração:
Da definição de transformada de Fourier, temos
${\mathcal {F}}\{f(at)\}=\int \limits _{-\infty }^{\infty }f(at){e^{-iwt}}dt$
Fazendo a mudança $\tau =at$, distinguindo dois casos:a$>$0 e a$<$0.
Para = a$>$0 temos:
${\mathcal {F}}\{f(at)\}=\int \limits _{-\infty }^{\infty }f(at){e^{-iwt}}dt$
$=\int \limits _{-\infty }^{\infty }f(\tau )e^{\frac {-iw\tau }{a}}{\frac {d\tau }{a}}$
$={\frac {1}{a}}\int \limits _{-\infty }^{\infty }f(\tau )e^{\frac {-iw\tau }{a}}d\tau$
Para o caso a$<$0, temos:
${\mathcal {F}}\{f(at)\}=\int \limits _{-\infty }^{\infty }f(at){e^{-iwt}}dt$
$=\int \limits _{\infty }^{-\infty }f(\tau )e^{\frac {-iw\tau }{a}}{\frac {d\tau }{a}}$
$=-{\frac {1}{a}}\int \limits _{-\infty }^{\infty }f(\tau )e^{\frac {-iw\tau }{a}}d\tau$
Em ambos os casos, temos:
${\mathcal {F}}\{f(at)\}={\frac {1}{|a|}}\int \limits _{-\infty }^{\infty }f(\tau )e^{\frac {-iw\tau }{a}}d\tau$
$={\frac {1}{|a|}}F({\frac {w}{a}})$

## Aplicações em Acústica e Sons

Chamamos de som toda e qualquer vibração das moléculas de ar que estão imediatamente próximas de nossos tímpanos. Ao fazê-lo, tais vibrações são captadas pela membrana sensível do órgão e assimiladas pelos sensores nervoso da orelha interna, resultando na sensação de escuta em nosso cérebro. No entanto, sinais sonoros podem ser de vários tipos, desde ruídos bruscos e desagradáveis até belos acordes de um instrumento musical.
Podemos explicar a natureza sonora utilizando a teoria de Séries de Fourier e enxergando uma excitação sonora que atinge nossos ouvidos como um sinal periódico em função do tempo, com o formato de uma onda senoide muito complexa, isto é, de desenho bifurcado, repleta de amplitudes diferentes em cada ponto. Esta onda, por sua vez, nada mais é que é uma soma de inúmeras ondas senoidais simples, cada uma ponderada com sua própria amplitude (coeficientes de Fourier) e oscilando à sua própria frequência (harmônicos). Em uma expansão em Séries de Fourier, harmônicos em termos subsequentes são especialmente ordenados, de tal forma que todos acabam sendo múltiplos de um primeiro mais baixo e dominante na nossa percepção do som, o qual denotamos Frequência Fundamental $f_{0}$. Logo, as próximas devem obedecer um relação $f_{n}=nf_{0}$, onde $n$cresce conforme aumentam os harmônicos.
É aqui que reside a explicação para distinção de diferentes tipos de som, advindos de diferentes instrumentos musicais por exemplo. Instrumentos de percussão como tambores geralmente possuem expansão em Série de Fourier mais simples e com poucos termos, soando somente os primeiros harmônicos mais graves e mais ressaltados, o que resulta no som seco e rústico de uma batida de tambor. Já um violino produzirá muitas de suas vibrações em frequências mais altas, pois são estas que dão riqueza ao som, caracterizando os sons mais bonitos e elaborados do instrumento. É também por trabalhar em harmônicos altos que identificamos o som de um violino como mais fino ou agudo.
Assim, se decompormos um sinal sonoro em Série de Fourier e tocarmos apenas o primeiro termo de seno ou cosseno (nosso ouvido não é capaz de perceber a defasagem entre os dois) da série, escutaremos um som seco, grave e desagradável. À medida que adicionamos os termos subsequentes, vamos refinando este som e aos poucos tornando-o mais agradável, assimilando novas camadas sonoras vindas das novas frequências. Chamamos de timbre a singularidade que esse fenômeno dá à cada som, isto é, a capacidade de discernir entre fontes graças às ondas diferentes que decompõem o sinal resultante.
Uma nota é composta por um conjunto de frequências que são dobras umas das outras, ou seja, diferentes escalas, mas divididas em intervalos iguais de razão $2$, os quais chamamos de oitavas. Digamos, por exemplo, que tocamos a nota Dó em um piano, chamando de $f$ sua frequência. Ao tocarmos uma frequência $2f$, estaríamos produzindo novamente um Dó, uma oitava mais agudo. Em $4f$, surge mais um Dó, desta vez duas oitavas mais agudo. O mesmo ocorre para outras frequências presentes no espectro. $3f$, por exemplo, representaria uma nota diferente, sendo $6f$ esta mesma nota porém na próxima escala. Uma ferramenta poderosa e importante para este tipo estudo é construir diagramas de espectro para cada sinal. Deste modo é possível visualizar seu comportamento no domínio-frequência e identificar facilmente sob quais harmônicos este está trabalhando.

## Aplicações na Engenharia Hídrica

Conforme Ricieri(1988), embora a série de Fourier tenha sido desenvolvida como subsídio matemático ao estudo da transferência de calor, a aplicação dessa soma de senos e cossenos estendeu -se aos ramos da física , engenharia e matemática, sendo também comum encontrar o uso da série de Fourier em artigos que propõem modelos para análise da periodicidade de alguns eventos climáticos. Assim de um modo geral, pode-se dizer que a série desenvolvida por Fourier têm permitido que engenheiros e cientistas escrevam diversos tipos de funções e com elas podem administrar e prever cenários futuros nas mais diversas áreas de estudo.
Logo, vários fenômenos e grandezas físicas presentes na atmosfera terrestre, por apresentarem um caráter quase periódico, podem ser descritos com uso de séries de Fourier, a exemplo das precipitações em suas variações sazonais, radiação que chega ao topo da atmosfera, padrões de vento, temperatura, umidade do ar, etc.
Sendo a pluviosidade um elemento climático fundamental por ter influência direta nas demais variáveis que compõem o meio ambiente(temperatura do ar ,radiação solar, umidade atmosférica), é de suma importância no desenvolvimento da fauna e flora que compõem um ecossistema.
Por possuírem caráter fundamental para um ecossistema, modelos são criados a partir de modelagem atmosférica relacionando a periodicidade de chuvas e pluviometria com o fator tempo. Esses modelos possuem propósito de encontrar padrões e prever cenários futuros através de séries de Fourier além prever eventos extremos.
Outra utilização importante na análise de Fourier, foi o estudo proposto por William Thomson, que envolvia a análise e predição harmônica, para detectar o comportamento das marés ao longo do ano.

## Série de Fourier em um quadrado

Artigo principal: Fourier Series
É possível definir a série de Fourier para funções de duas variáveis, como $x$ e $y$, num intervalo $[-\pi ,\pi ]\times [-\pi ,\pi ]$
$f(x,y)=\textstyle \sum _{j,k\in \mathbb {Z} }^{N}\displaystyle c_{j,k}e^{-ijx}e^{-iky}dxdy$, onde $c_{j,k}={\frac {1}{4\pi ^{2}}}\int _{-\pi }^{\pi }\int _{-\pi }^{\pi }f(x,y)e^{-ijx}e^{-iky}dxdy$
Além dessa definição ser útil para resolver equações diferenciais parciais, uma aplicação notável da série de Fourier no quadrado é a compressão de imagens. Em particular, o padrão de compressão de imagem jpeg usa a transformação cosseno discreta bidimensional, que é uma transformação de Fourier usando as funções básicas do cosseno.