Sèrie de Fourier

Transformades de Fourier
Transformada de Fourier continua
Sèrie de Fourier
Transformada Discreta de Fourier
Transformada de Fourier en Temps Discret
Transformada de Fourier sobre cossos finits
Anàlisi de Fourier
Transformades relacionades
Les primeres quatre aproximacions per a una funció periòdica esglaonada

En matemàtiques, una sèrie de Fourier descompon una funció periòdica en una suma de funcions oscil·latòries simples: el sinus i el cosinus. L'estudi de les sèries de Fourier forma part de l'anàlisi de Fourier. Les sèries de Fourier foren introduïdes per Joseph Fourier (1768–1830) amb l'objectiu de solucionar una equació de la calor en un plat metàl·lic. Va suposar una revolució dins les matemàtiques, forçant els matemàtics a reexaminar els fonaments de les matemàtiques i portant a la descoberta de teories modernes com la Integral de Lebesgue. Joseph Fourier va ser el primer que va estudiar tals sèries sistemàticament, publicant els seus resultats inicials el 1807 i 1811. Aquesta àrea d'investigació sovint s'anomena anàlisi harmònica.

Amb aquesta eina podrem analitzar un senyal periòdic en termes del seu contingut freqüencial o espectre. Ens permetrà establir la dualitat entre temps i freqüència, així, operacions realitzades en el domini temporal tindran també el seu dual en el domini freqüencial.

La sèrie de Fourier té aplicacions en moltes branques de l'enginyeria, a més de ser una eina molt útil en la teoria matemàtica abstracta. Àrees d'aplicació inclouen anàlisi de les vibracions, acústica, òptica, processament del senyal, processament d'imatge, etc. En enginyeria, per al cas dels sistemes de telecomunicacions, i a través de l'ús dels components espectrals de freqüència d'un senyal donat, es pot optimitzar el disseny d'un sistema per al senyal portador d'aquest.

Evolució històrica

Context

Donada una funció periòdica f ( t ) {\displaystyle f(t)} , per exemple de període 2 π {\displaystyle 2\pi } , es desitja escriure-la com una combinació en la qual intervinguin únicament sinus i cosinus, que són les funcions de període 2 π {\displaystyle 2\pi } simples més conegudes:[1]

f ( t ) = a 0 2 + n = 1 [ a n cos ( n t ) + b n sin ( n t ) ] {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(nt)+b_{n}\sin(nt)]}

Aquesta sèrie rep el nom de sèrie trigonomètrica o Sèrie de Fourier.

Jean Baptiste Joseph Fourier, matemàtic i físic francès, autor de la Sèrie de Fourier

El problema de la representació d'una funció mitjançant una sèrie trigonomètrica sorgeix de la resolució d'equacions en derivades parcials. Pels voltants del 1750, J. d'Alembert, D. Bernoulli i L. Euler van estudiar l'equació d'ones que governa el problema de la corba vibrant, un problema plantejat i estudiat per B. Taylor que havia obtingut solucions en forma de funcions sinusoidals. D'Alembert va donar una solució molt general i Euler va provar que si a l'instant inicial de la forma de la corda apareix una combinació finita de sinus, llavors esdevindria el mateix que en qualsevol altre instant posterior, publicant molt més tard, el 1777, les fórmules que permetien calcular els coeficients de la combinació. El 1753, Bernoulli va utilitzar aquesta representació per resoldre el problema de la corda vibrant per una posició inicial qualsevol, però la seva solució va suscitar molta controvèrsia. Va ser J. B. Fourier qui va reprendre un altre cop les idees d'Euler i Bernoulli i va obtenir resultats molt ajustats els experiments, col·locant l'estudi de les sèries trigonomètriques – que avui dia porten el seu nom - en el centre de l'escenari matemàtic del segle xix. La teoria de les sèries de Fourier va ser una forta influència per l'anàlisi matemàtica i és avui una eina bàsica de l'enginyeria de telecomunicacions.

Periodicitat

En el camp de la teoria del senyal i la comunicació, quan t és la variable temporal, es diu que f(t) és un senyal periòdic en temps continu, i quan aquesta representació en sèrie de funcions trigonomètriques és correcta, es pot dir que s'ha descompost el senyal en harmònics. A continuació s'estudia més a fons com calcular els coeficients, alguns exemples i algunes propietats, i també la notació exponencial de les sèries de Fourier. Finalment, es consideraran senyals periòdics perquè resultaran més senzills pel seu estudi. Recordem que:

f : R R {\displaystyle f:\mathbb {R} \rightarrow \mathbb {R} } és periòdica de període   T > 0 {\displaystyle \ T>0}

si és tal que

t R , f ( t + T ) = f ( t ) {\displaystyle \forall t\in \mathbb {R} ,\,f(t+T)=f(t)}

Pels casos en què les funcions no compleixen la periodicitat s'introdueix la Transformada de Fourier, que, encara que comparteix moltes característiques i propietats amb la Transformada de Laplace, s'utilitza en contextos diferents.

Definició

  • Les Sèries de Fourier descriuen els senyals periòdics com a combinació de senyals harmònics (sinusoidals).
  • Amb aquesta eina es pot analitzar un senyal periòdic en termes del seu contingut freqüencial o espectre.
  • Permet establir una dualitat entre temps i freqüència, de manera que les operacions realitzades en el domini temporal tenen el seu dual en el domini freqüencial.
  • La forma trigonomètrica de les Sèries de Fourier permet descriure una funció periòdica x(t) de període T (freqüència fonamental f 0 = 1 T , ω 0 = 2 π f 0 {\displaystyle f_{0}={\frac {1}{T}},\omega _{0}=2\pi f_{0}} ).

Forma general

Sigui f(x) una funció complexa f definida en el domini dels nombres reals, on t f : R C {\displaystyle f:\mathbb {R} \rightarrow \mathbb {C} } , definida a trossos, periòdica en f(t) amb període T, i de quadrat integrable a l'interval de t 1 {\displaystyle t_{1}} a t 2 {\displaystyle t_{2}} de durada T:

t 1 t 2 | f ( t ) | 2 d t < + {\displaystyle \int _{t_{1}}^{t_{2}}|f(t)|^{2}dt<+\infty }

amb T = t 2 t 1 {\displaystyle T=t_{2}-t_{1}} , llavors es definex la sèrie trigonomètrica de Fourier:

f ( t ) = a 0 2 + n = 1 [ a n cos ( ω n t ) + b n sin ( ω n t ) ] {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\,{\big [}a_{n}\cos(\omega _{n}t)+b_{n}\sin(\omega _{n}t){\big ]}}

on els seus coeficients venen donats per les següents expressions:

  • ω n = n 2 π T {\displaystyle \omega _{n}=n{\frac {2\pi }{T}}} , on n és el nombre d'harmònics de la funció f,
  • a n = 2 T t 1 t 2 f ( t ) cos ( ω n t ) d t {\displaystyle a_{n}={\frac {2}{T}}\int _{t_{1}}^{t_{2}}f(t)\cos(\omega _{n}t)\,dt} , és la part parella dels coeficients de Fourier de la funció f,
  • b n = 2 T t 1 t 2 f ( t ) sin ( ω n t ) d t {\displaystyle b_{n}={\frac {2}{T}}\int _{t_{1}}^{t_{2}}f(t)\sin(\omega _{n}t)\,dt} , és la part senar dels coeficients de Fourier de la funció f.

Forma canònica

En el cas especial de tenir un període T = 2π, s'obté que ω n = n {\displaystyle \omega _{n}=n\,} .

En aquest cas, els coeficients de la Sèrie de Fourier es redueixen, arribant a una expressió particularment simple:

f ( t ) = a 0 2 + n = 1 [ a n cos ( n t ) + b n sin ( n t ) ] {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(nt)+b_{n}\sin(nt)]}

on

  • a n = 1 π π π f ( t ) cos ( n t ) d t {\displaystyle a_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(t)\cos(nt)\,dt}
  • b n = 1 π π π f ( t ) sin ( n t ) d t {\displaystyle b_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(t)\sin(nt)\,dt}

Forma exponencial

Per la identitat d'Euler per l'exponencial complexa, operant adequadament, si

C n = 1 2 π π π f ( x ) e i n x d x {\displaystyle C_{n}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x)\,e^{-inx}\,dx}

La sèrie de Fourier es pot expressar com la suma de dues sèries:

n = 1 C n e i n x + n = 0 C n e i n x {\displaystyle \sum _{n=1}^{\infty }C_{-n}\,e^{-inx}+\sum _{n=0}^{\infty }C_{n}\,e^{inx}}

De manera més compacta:

n = C n e i n x {\displaystyle \sum _{n=-\infty }^{\infty }C_{n}\,e^{inx}}

Forma del període

La forma pel període T, es pot derivar fàcilment a la forma canònica amb el canvi de la variable definida per x = 2 π T t {\displaystyle x={\frac {2\pi }{T}}t} . Per tant, les dues formulacions són equivalents. Tot i això, a la pràctica la forma pel període T s'utilitza de forma majoritària degut a la seva aplicació directa. A la teoria, sempre és preferible utilitzar la forma canònica atès que és més elegant i més fàcil d'interpretar matemàticament.

Propietats

  • Superposició
α x p ( t ) + β y p ( t ) α X s [ k ] + β Y s [ k ] {\displaystyle \alpha x_{p}(t)+\beta y_{p}(t)\leftrightarrow \alpha X_{s}[k]+\beta Y_{s}[k]}
  • Derivada
x p ( t ) j k 2 π f 0 X s [ k ] ( k 0 ) {\displaystyle x'_{p}(t)\leftrightarrow jk2\pi f_{0}X_{s}[k](k\neq 0)}
  • Integral
0 t x p ( t ) d t X [ k ] j k 2 π f 0 + C ( k 0 ) {\displaystyle \int _{0}^{t}x_{p}(t)\,dt\leftrightarrow {\frac {X[k]}{jk2\pi f_{0}}}+C(k\neq 0)}
  • Retard
x p ( t α ) X s [ k ] e j 2 k π f 0 α {\displaystyle x_{p}(t-\alpha )\leftrightarrow X_{s}[k]e^{-j2k\pi f_{0}\alpha }}
  • Escalat
x p ( t )   X s [ k ] {\displaystyle x_{p}(t)\leftrightarrow \ X_{s}[k]} harmònics en f = k f 0 α {\displaystyle f=kf_{0}\alpha }
x p ( t )   X s [ k ] = X s [ k ] {\displaystyle x_{p}(-t)\leftrightarrow \ X_{s}[-k]=X_{s}^{*}[k]}
  • Modulació
cos ( m 2 π f 0 t ) x p ( t ) 1 2 { X s [ k m ] + X s [ k + m ] } {\displaystyle \cos(m2\pi f_{0}t)x_{p}(t)\leftrightarrow {\frac {1}{2}}{\big \{}X_{s}[k-m]+X_{s}[k+m]{\big \}}}
1 2 { x p ( t + α ) + x p ( t α ) } cos ( 2 π f 0 α ) X s [ k ] {\displaystyle {\frac {1}{2}}{\big \{}x_{p}(t+\alpha )+x_{p}(t-\alpha ){\big \}}\leftrightarrow \cos(2\pi f_{0}\alpha )X_{s}[k]}
  • Convolució
x p ( t ) y p ( t ) X s [ k ] Y s [ k ] {\displaystyle x_{p}(t)y_{p}(t)\leftrightarrow X_{s}[k]*Y_{s}[k]}
x p ( t ) y p ( t )   X s [ k ] Y s [ k ] {\displaystyle x_{p}(t)\bullet y_{p}(t)\leftrightarrow \ X_{s}[k]Y_{s}[k]}

Exemple

Escriure la Sèrie de Fourier corresponent a la funció f definida a continuació, a l'interval π {\displaystyle -\pi } < x < π {\displaystyle \pi } :

f ( x ) = { 0 π < x 0 x 0 x < π {\displaystyle f(x)={\begin{cases}0&-\pi <x\leq 0\\x&0\leq x<\pi \end{cases}}}

Els coeficients de Fourier són:

Fig.1

a 0 = 1 π π π f ( x ) d x {\displaystyle a_{0}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(x)\,dx}

= 1 π ( 0 + 0 π x d x ) {\displaystyle ={\frac {1}{\pi }}{\Big (}0+\int _{0}^{\pi }x\,dx{\Big )}}

a 0 = π 2 , {\displaystyle a_{0}={\frac {\pi }{2}},}

a n = 1 π 0 π x cos ( n x ) d x {\displaystyle a_{n}={\frac {1}{\pi }}\int _{0}^{\pi }x\cos(nx)\,dx}

= [ cos n x + n x sin n x π n 2 ] 0 π {\displaystyle =\left[{\frac {\cos nx+nx\sin nx}{\pi n^{2}}}\right]_{0}^{\pi }}

a n = ( 1 ) n 1 π n 2 , {\displaystyle a_{n}={\frac {(-1)^{n}-1}{\pi n^{2}}},}

b n = 1 π 0 π x sin ( n x ) d x {\displaystyle b_{n}={\frac {1}{\pi }}\int _{0}^{\pi }x\sin(nx)\,dx}

= [ sin n x n x cos n x π n 2 ] 0 π {\displaystyle =\left[{\frac {\sin nx-nx\cos nx}{\pi n^{2}}}\right]_{0}^{\pi }}

b n = ( 1 ) n n {\displaystyle b_{n}=-{\frac {(-1)^{n}}{n}}}


on n =1, 2,... Per tant, a l'interval ( π , π ) {\displaystyle (-\pi ,\pi )} ,

f ( x ) π 4 + n > 0 [ ( 1 ) n 1 π n 2 cos n x ( 1 ) n n sin n x ] {\displaystyle f(x)\sim {\frac {\pi }{4}}+\sum _{n>0}\left[{\frac {(-1)^{n}-1}{\pi n^{2}}}\cos nx-{\frac {(-1)^{n}}{n}}\sin nx\right]}

Es veu així per què aquesta Sèrie de Fourier convergeix a f(x) quan ( π < x < π ) {\displaystyle (-\pi <x<\pi )} . Així mateix aquesta sèrie representa l'extensió periòdica de la funció f indicada a les línies de punts de la Fig.1. La funció periòdica és discontínua en els punts x = ± π , ± 3 π {\displaystyle x=\pm \pi ,\pm 3\pi \ldots } . La teoria mostrarà que la suma de la sèrie a cadascun d'aquests punts ha de ser π 2 {\displaystyle {\frac {\pi }{2}}} . Una petita indicació de la convergència de la sèrie a f(x), és que es poden sumar alguns termes de la sèrie per composició d'ordenades. Es trobarà, per exemple, que la gràfica de la funció és una aproximació ondulada de la gràfica de la Fig. 1 y = π 4 2 π cos x + sin x 1 2 sin 2 x {\displaystyle y={\frac {\pi }{4}}-{\frac {2}{\pi }}\cos x+\sin x-{\frac {1}{2}}\sin 2x} .

Convergència de la Sèrie de Fourier

Sigui f : R R {\displaystyle f:\mathbb {R} \to \mathbb {R} } contínua a trossos i periòdica de període T, i sigui

f ( t ) a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)\sim {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

on ω = 2 π T {\displaystyle \omega ={\frac {2\pi }{T}}} és la seva Sèrie de Fourier. P.G.L. Dirichlet el 1829 va demostrar satisfactòriament que un grup específic de funcions, són iguals a la suma de les seves sèries de Fourier.

Condicions de Dirichlet

Sigui f : R R {\displaystyle f:\mathbb {R} \to \mathbb {R} } una funció periòdica de període T. Llavors, f satisfà les condicions de Dirichlet si cada període de la funció f : [ 0 , T ] R {\displaystyle f:[0,T]\to \mathbb {R} } és continu excepte un nombre finit de discontinuïtats de salt. En particular es pot demostrar, que si una funció periòdica és tal que ella i la seva derivada estan definides i són contínues excepte en un nombre finit de discontinuïtats de salt, llavors aquesta funció verifica les condicions de Dirichlet.

Teorema de convergència de Dirichlet

Sigui f : R R {\displaystyle f:\mathbb {R} \to \mathbb {R} } una funció periòdica de període T que satisfà les condicions de Dirichlet i sigui

f ( t ) a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)\sim {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

on ω = 2 π T {\displaystyle \omega ={\frac {2\pi }{T}}} és la seva Sèrie de Fourier.

  • Si f és contínua en un punt t, llavors la Sèrie de Fourier convergeix en aquest punt a f(t), o sigui
a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] = f ( t ) {\displaystyle {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]=f(t)}
  • Si f té una discontinuïtat de salt en un punt t, llavors la sèrie convergeix en aquest punt vers el punt mitjà del salt, o sigui
a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] = f ( t ) + f ( t + ) 2 {\displaystyle {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]={\frac {f(t^{-})+f(t^{+})}{2}}}

on f ( t ) = lim τ 0 , τ f ( t τ ) {\displaystyle f(t^{-})=\lim _{\tau \to 0,\tau }f(t-\tau )} indica el límit de f a t per l'esquerra i f ( t + ) = lim τ 0 , τ f ( t + τ ) {\displaystyle f(t^{+})=\lim _{\tau \to 0,\tau }f(t+\tau )} indica el límit de f a t per la dreta.

El teorema indica que si f satisfà les condicions de Dirichlet i es redefineix el valor de f a cada punt de la discontinuïtat com el punt mitjà del salt, és a dir f ( t ) = f ( t ) + f ( t + ) 2 {\displaystyle f(t)={\frac {f(t^{-})+f(t^{+})}{2}}} , llavors la suma de la Sèrie de Fourier coincideix amb f(t) a cada t R {\displaystyle t\in \mathbb {R} } . Suposant que això es compleixi, es pot afirmar que la condició s'acompleix.

Convergència uniforme de les sèries de Fourier

Si es compleix que la funció f {\displaystyle f} , de període a, és diferenciable en l'interval [0,a], excepte potser en un nombre finit de punts, i la derivada f {\displaystyle f'} és contínua, aleshores:[2]

1- La sèrie de Fourier de f {\displaystyle f'} s'obté derivant terme a terme els coeficients de la sèrie de Fourier de f {\displaystyle f} .

2- Els coeficients de Fourier de f {\displaystyle f} satisfan:

n = + | c n ( f ) |   <   {\displaystyle \sum _{n=-\infty }^{+\infty }|c_{n}(f)|\ <\ \infty }

3- La sèrie de Fourier de f {\displaystyle f} convergeix uniformement cap a f {\displaystyle f} en R {\displaystyle \mathbb {R} } .[3][4]

Fenomen de Gibbs

Fenomen de Gibbs

El teorema de Dirichlet diu que en els punts de discontinuïtat, la gràfica de la suma de la Sèrie de Fourier passa per un punt entremig del salt. Si es dibuixen les sumes parcials es veu que els voltants dels punts de discontinuïtat es redueix la velocitat de convergència de la sèrie i que la gràfica de la suma parcial oscil·la al voltant de la gràfica de la funció. Quan augmenta el nombre de termes, les oscil·lacions es condensen a ambdós costats del punt, però la seva amplitud no decreix.[5] Això es coneix com el fenomen de Gibbs, en honor de J.W. Gibbs, que el 1899 va provar que l'amplitud l'oscil·lació a cada costat de la gràfica d'una funció tendeix a ser 1 2 π 0 π sin ( t ) t 1 d t 0.0895 {\displaystyle {\frac {1}{2\pi }}\int _{0}^{\pi }{\frac {\sin(t)}{t}}-1\,dt\approx 0.0895} vegades la mida del salt, és a dir, aproximadament el 9%.[6]

Aplicacions Sèrie de Fourier

  • Generació de formes d'ona de corrent o tensió elèctrica per mitjà de la superposició de sinusoides generades per oscil·ladors electrònics d'amplitud variable de les quals les freqüències ja estan determinades.
  • Anàlisi en el comportament harmònic d'un senyal.
  • Reforçament de senyals.
  • Estudi de la resposta en el temps d'una variable circuital elèctrica on el senyal d'entrada no és sinusoidal o cosinusoidal, mitjançant l'ús de les transformades de Laplace i de la solució en règim permanent sinusoidal en el domini de la freqüència.
  • La resolució d'algunes equacions diferencials en derivades parcials permeten solucions particulars en sèries de Fourier fàcilment computables, i per obtenir solucions pràctiques, en la teoria de la transmissió de la calor, la teoría de plaques, etc.

Algunes propietats del desenvolupament en Sèrie de Fourier

Derivació de les sèries de Fourier

Si f és una funció contínua i periòdica de període T i la seva derivada f' verifica les condicions de Dirichlet, llavors, la sèrie de Fourier de f es pot derivar terme a terme de manera que si

f ( t ) = a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

aleshores

f ( t ) n = 1 [ n b n ω cos ( n ω t ) + n a n ω sin ( n ω t ) ] {\displaystyle f'(t)\sim \sum _{n=1}^{\infty }[nb_{n}\omega \cos(n\omega t)+na_{n}\omega \sin(n\omega t)]}

per cada t R {\displaystyle t\in \mathbb {R} }

Integració de les sèries de Fourier

Si f és una funció contínua i periòdica de període T que verifica les condicions de Dirichlet, llavors, la sèrie de Fourier de f es pot integrar terme a terme de manera que si

f ( t ) = a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

i t 0 , t [ T 2 , T 2 ] {\displaystyle t_{0},t\in \left[{\frac {-T}{2}},{\frac {T}{2}}\right]} , aleshores

t 0 t f ( τ ) d τ = 1 2 a 0 ( t t 0 ) + n = 1 [ b n ( cos ( n ω t 0 ) cos ( n ω t ) ) n ω + a n ( sin ( n ω t ) sin ( n ω t 0 ) ) n ω ] {\displaystyle \int _{t_{0}}^{t}f(\tau )\,d\tau ={\frac {1}{2}}a_{0}(t-t_{0})+\sum _{n=1}^{\infty }\left[{\frac {b_{n}(\cos(n\omega t_{0})-\cos(n\omega t))}{n\omega }}+{\frac {a_{n}(\sin(n\omega t)-\sin(n\omega t_{0}))}{n\omega }}\right]}

S'ha d'anar amb compte, ja que al terme 1 2 a 0 t {\displaystyle {\frac {1}{2}}a_{0}t} fa que el membre de la dreta no sigui una Sèrie de Fourier. El teorema ens proporciona el desenvolupament en sèrie de Fourier de la funció g definida a [ T 2 , T 2 ] {\displaystyle \left[{\frac {-T}{2}},{\frac {T}{2}}\right]} per t 0 t f ( τ ) d τ 1 2 a 0 t {\displaystyle \int _{t_{0}}^{t}f(\tau )\,d\tau -{\frac {1}{2}}a_{0}t} i repetida periòdicament.

Igualtat de Parseval

Si f és una periòdica de període T que verifica les condicions de Dirichlet, i sigui

f ( t ) a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)\sim {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

la seva Sèrie de Fourier. Aleshores

1 T 0 T [ f ( t ) ] 2 d t = 1 4 a 0 2 + 1 2 n = 1 [ a n 2 + b n 2 ] {\displaystyle {\frac {1}{T}}\int _{0}^{T}[f(t)]^{2}\,dt={\frac {1}{4}}a_{0}^{2}+{\frac {1}{2}}\sum _{n=1}^{\infty }[a_{n}^{2}+b_{n}^{2}]}

Quan f és un senyal periòdic de període fonamental T, aquesta igualtat es pot interpretar de la següent manera. La integral P = 1 T 0 T [ f ( t ) ] 2 d t {\displaystyle P={\frac {1}{T}}\int _{0}^{T}[f(t)]^{2}\,dt} s'anomena mitjana quadràtica o potència mitjana de f.

Notació complexa de la Sèrie de Fourier

Forma general

Hem vist dues maneres d'expressar la Sèrie de Fourier d'una funció periòdica de període T: en termes de sinus i cosinus

f ( t ) a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle f(t)\sim {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]}

o bé en termes l'amplitud i fase

f ( t ) α 0 + n = 1 α n cos ( n ω t + ϕ n ) {\displaystyle f(t)\sim \alpha _{0}+\sum _{n=1}^{\infty }\alpha _{n}\cos(n\omega t+\phi _{n})}

on ω = 2 π T {\displaystyle \omega ={\frac {2\pi }{T}}} i els coeficients d'aquestes expressions estan relacionades de la següent manera: α 0 = 1 2 a 0 {\displaystyle \alpha _{0}={\frac {1}{2}}a_{0}} i per n=1,2...

α n = a n 2 + b n 2 {\displaystyle \alpha _{n}={\sqrt {a_{n}^{2}+b_{n}^{2}}}}
i
ϕ n = arctan ( b n a n ) ( π , π ] {\displaystyle \phi _{n}=\arctan \left({\frac {-b_{n}}{a_{n}}}\right)\in (-\pi ,\pi ]}

o bé a n = cos ( ϕ n ) {\displaystyle a_{n}=\cos(\phi _{n})}

i

b n = α n sin ( ϕ n ) {\displaystyle b_{n}=-\alpha _{n}\sin(\phi _{n})}

A moltes aplicacions, en particular a la teoria de senyals i comunicacions, és més útil una tercera expressió, anomenada notació complexa de la sèrie de Fourier, que consisteix amb la utilització d'exponencials complexos e j n ω t {\displaystyle e^{jn\omega t}} en lloc de cos ( n ω t ) {\displaystyle \cos(n\omega t)} i sin ( n ω t ) {\displaystyle \sin(n\omega t)} . Per construir-les, substituïm les expressions

cos ( n ω t ) = e j n ω t + e j n ω t 2 {\displaystyle \cos(n\omega t)={\frac {e^{jn\omega t}+e^{-jn\omega t}}{2}}}
i
sin ( n ω t ) = e j n ω t e j n ω t 2 {\displaystyle \sin(n\omega t)={\frac {e^{jn\omega t}-e^{-jn\omega t}}{2}}}

a la sèrie a 0 2 + n = 1 [ a n cos ( n ω t ) + b n sin ( n ω t ) ] {\displaystyle {\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }[a_{n}\cos(n\omega t)+b_{n}\sin(n\omega t)]} , obtenint

f ( t ) a 0 2 + 1 2 n = 1 [ ( a n j b n ) e j n ω t + ( a n + j b n ) e j n ω t ] {\displaystyle f(t)\sim {\frac {a_{0}}{2}}+{\frac {1}{2}}\sum _{n=1}^{\infty }[(a_{n}-jb_{n})e^{jn\omega t}+(a_{n}+jb_{n})e^{-jn\omega t}]}

Escrivint

c 0 = c 0 e 0 = 1 2 a 0 {\displaystyle c_{0}=c_{0}e^{0}={\frac {1}{2}}a_{0}} , c n = a n j b n 2 {\displaystyle c_{n}={\frac {a_{n}-jb_{n}}{2}}} , c n = c n ¯ = a n + j b n 2 {\displaystyle c_{-n}={\overline {c_{n}}}={\frac {a_{n}+jb_{n}}{2}}}

separant els termes que tenen índex negatiu dels que tenen índex positiu i ordenant-los, la Sèrie de Fourier queda

f ( t ) n = 1 + c n e j n ω t {\displaystyle f(t)\sim \sum _{n=1}^{+\infty }c_{n}e^{jn\omega t}}

que s'anomena notació complexa de la Sèrie de Fourier de f. Els nombres c n {\displaystyle c_{n}} són els coeficients de Fourier de la notació complexa i es poden calcular de forma directa

c n = 1 T T 2 T 2 f ( t ) e j n ω t d t {\displaystyle c_{n}={\frac {1}{T}}\int _{\frac {-T}{2}}^{\frac {T}{2}}f(t)e^{-jn\omega t}\,dt} per n = 0 , ± 1 , ± 2 , {\displaystyle n=0,\pm 1,\pm 2,\ldots }

En canvi si el que tenim són els coeficients c n {\displaystyle c_{n}} de la notació complexa, aleshores els coeficients de Fourier a n {\displaystyle a_{n}} i b n {\displaystyle b_{n}} , les amplituds i les fases venen donades per

a n = 2 Re ( c n ) , b n = 2 Im ( c n ) , α n = 2 | c n | {\displaystyle a_{n}=2\operatorname {Re} (c_{n}),b_{n}=-2\operatorname {Im} (c_{n}),\alpha _{n}=2|c_{n}|} i ϕ n = arg ( c n ) ( π , π ] {\displaystyle \phi _{n}=\arg(c_{n})\in (-\pi ,\pi ]}

Aleshores, veiem, que una funció periòdica és par quan els seus coeficients de Fourier complexos són reals i que és imparell quan els seus coeficients complexos de Fourier són imaginaris purs. Quan f(t) és un senyal periòdic de període fonamental T, les components de la notació complexa es donen en freqüència 0 , ± ω , ± 2 ω , {\displaystyle 0,\pm \omega ,\pm 2\omega ,\ldots }

Propietats

Suposarem que f i g són funcions periòdiques de període T (reals o complexes) que verifiquen les condicions de Dirichlet, els quals desenvolupats en notació complexa de la Sèrie de Fourier, són, respectivament,

f ( t ) = n = + c n e j n ω t {\displaystyle f(t)=\sum _{n=-\infty }^{+\infty }c_{n}e^{jn\omega t}} i g ( t ) = n = + d n e j n ω t {\displaystyle g(t)=\sum _{n=-\infty }^{+\infty }d_{n}e^{jn\omega t}}

on ω = 2 π T {\displaystyle \omega ={\frac {2\pi }{T}}}

  • Linealitat: si p i q son nombres complexos
p f ( t ) + q g ( t ) = + ( p c n + q d n ) e j n ω t {\displaystyle pf(t)+qg(t)=\sum _{-\infty }^{+\infty }(pc_{n}+qd_{n})e^{jn\omega t}}
  • Translació en el temps: t 0 {\displaystyle t_{0}} és un nombre real.
f ( t t 0 ) = n = + c n e j n ω t 0 e j n ω t = n = + c n e j n ω ( t t 0 ) {\displaystyle f(t-t_{0})=\sum _{n=-\infty }^{+\infty }c_{n}e^{-jn\omega t_{0}}e^{jn\omega t}=\sum _{n=-\infty }^{+\infty }c_{n}e^{jn\omega (t-t_{0})}}
f ( t ) {\displaystyle f(t)} i f ( t t 0 ) {\displaystyle f(t-t_{0})} tenen el mateix espectre d'amplituds però fase diferent.
  • Escalat en el temps: p és un nombre real, aleshores la funció f(pt) és periòdica de període T i freqüència p ω {\displaystyle p\omega } .
f ( p t ) = n = + c n e j n ( p ω ) t {\displaystyle f(pt)=\sum _{n=-\infty }^{+\infty }c_{n}e^{jn(p\omega )t}}
f ( t ) {\displaystyle f(t)} i f ( p t ) {\displaystyle f(pt)} tenen les mateixes amplituds y fases però corresponen a freqüències diferents
  • Convolució
h n = n = + c k d n k = n = + c n k d k {\displaystyle h_{n}=\sum _{n=-\infty }^{+\infty }c_{k}d_{n-k}=\sum _{n=-\infty }^{+\infty }c_{n-k}d_{k}}
  • Multiplicació
1 T 0 T f ( t ) g ( t ) d t = n = + c n d n ¯ {\displaystyle {\frac {1}{T}}\int _{0}^{T}f(t)g(t)\,dt=\sum _{n=-\infty }^{+\infty }c_{n}{\overline {d_{n}}}}
  • Forma complexa de la igualtat de Parseval
1 T 0 T | f ( t ) 2 | d t = n = + | c n 2 | {\displaystyle {\frac {1}{T}}\int _{0}^{T}|f(t)^{2}|\,dt=\sum _{n=-\infty }^{+\infty }|{c_{n}}^{2}|}
Quan f(t) és un senyal periòdic de període fonamental T, aquesta igualtat ens diu que la potència mitjana del senyal és P = n = + | c n | 2 {\displaystyle P=\sum _{n=-\infty }^{+\infty }|{c_{n}}|^{2}} . Per aquest fet, la representació dels valors | c n | 2 {\displaystyle |{c_{n}}|^{2}} quan situem les freqüències a l'eix d'abscisses s'anomena espectre discret de potències.

Exemple

Gràfic d'una funció periódica.
Animació de les 5 primeres sèries de Fourier.

Els coeficients complexos de Fourier d'un tren de polsos rectangulars, donats com l'extensió periòdica de la funció definida a l'interval de [ T 2 , T 2 ] {\displaystyle \left[{\frac {-T}{2}},{\frac {T}{2}}\right]} per

f ( n ) = { k  si  0 | t | < d 2 0  si  d 2 t | < T 2 {\displaystyle f(n)={\begin{cases}k&{\mbox{ si }}0\leq |t|<{\frac {d}{2}}\\0&{\mbox{ si }}{\frac {d}{2}}\leq \|t|<{\frac {T}{2}}\end{cases}}}

són

c n = k T d 2 d 2 e j n ω t d t = k d T sin ( n ω d 2 ) n ω d 2 {\displaystyle c_{n}={\frac {k}{T}}\int _{\frac {-d}{2}}^{\frac {d}{2}}e^{-jn\omega t}\,dt={\frac {kd}{T}}{\frac {\sin({\frac {n\omega d}{2}})}{\frac {n\omega d}{2}}}}

Els valor de c n {\displaystyle c_{n}} són reals, el que correspon el fet que la funció sigui parella, així doncs el seu espectre de fases és ϕ n = 0 {\displaystyle \phi _{n}=0} per cada n. La funció

s a ( t ) := { 1  si  t = 0 sin ( t ) t  si  t 0 {\displaystyle sa(t):={\begin{cases}1&{\mbox{ si }}t=0\\{\frac {\sin(t)}{t}}&{\mbox{ si }}t\neq 0\end{cases}}}

s'anomena funció de mostreig (‘sa' prové de la paraula anglesa sampling) i és molt important a la teoria de senyals, precisament perquè els coeficients complexos de Fourier del tren de polsos rectangulars d'aquest exemple es poden expressar com:

c n = k d T s a ( n ω d 2 ) {\displaystyle c_{n}={\frac {kd}{T}}sa\left({\frac {n\omega d}{2}}\right)}

A vegades, en lloc de la funció de mostreig es treballa amb una altra funció, anomenada sinus cardinal que es defineix com sinc(t) := s a ( π t ) := { 1  si  t = 0 sin ( π t ) π t  si  t 0 {\displaystyle {\mbox{sinc(t)}}:=sa(\pi t):={\begin{cases}1&{\mbox{ si }}t=0\\{\frac {\sin(\pi t)}{\pi t}}&{\mbox{ si }}t\neq 0\end{cases}}}

Enginyeria

L'anàlisi de senyals en el domini de la freqüència es realitza a través de les sèries de Fourier, ja que és molt comú, reemplaçar la variable x per ωt, resultant les components:

C k = 1 T t 1 t 2 f ( t ) e j ω k t d t {\displaystyle C_{k}={\frac {1}{T}}\int _{t_{1}}^{t_{2}}f(t)e^{-j\omega kt}\,dt}

Per tant:

f ( t ) = k = + C k e j ω k t {\displaystyle f(t)=\sum _{k=-\infty }^{+\infty }C_{k}e^{j\omega kt}}

Espectre d'una funció periòdica

Sigui S ( x ) {\displaystyle S(x)} una funció periòdica, de període T = n 2 π ω n {\displaystyle T={\frac {n2\pi }{\omega _{n}}}} , que es pot escriure com:

f ( t ) = a 0 2 + n = 1 ( a n cos ( ω n t ) + b n sin ( ω n t ) ) = n = C n e i ω n t {\displaystyle f(t)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }{\bigg (}a_{n}\cos(\omega _{n}t)+b_{n}\sin(\omega _{n}t){\bigg )}=\sum _{n=-\infty }^{\infty }C_{n}e^{i\omega _{n}t}}

Aleshores direm que la sèrie de Fourier anterior és la descomposició espectral de la funció f ( t ) {\displaystyle f(t)} . El conjunt de coeficients a n {\displaystyle a_{n}} i b n {\displaystyle b_{n}} o, el conjunt de coeficients C n {\displaystyle C_{n}} s'anomena l'espectre de la funció f ( t ) {\displaystyle f(t)} . L'espectre d'una funció, per tant, conté tota la informació d'aquesta i, si considerem la funció f ( t ) {\displaystyle f(t)} com la superposició d'ones planes, l'espectre de la funció és el conjunt d'amplituds d'aquestes ones.

Exemples

Sigui la funció de període 2 π {\displaystyle 2\pi } , f ( x ) = { 0 π < x 0 x 0 x < π {\displaystyle f(x)={\begin{cases}0&-\pi <x\leq 0\\x&0\leq x<\pi \end{cases}}}

El seu espectre és, com hem vist a l'exemple d'abans,

{ a 0 } { a n } { b n } = { π 2 } { ( 1 ) n 1 π n 2 } { ( 1 ) n n } , n N {\displaystyle \{a_{0}\}\cup \{a_{n}\}\cup \{b_{n}\}={\Big \{}{\frac {\pi }{2}}{\Big \}}\cup {\Bigg \{}{\frac {(-1)^{n}-1}{\pi n^{2}}}{\Bigg \}}\cup {\Bigg \{}-{\frac {(-1)^{n}}{n}}{\Bigg \}},n\in \mathbb {N} ^{*}}

o, equivalentment

{ C n } = { π 4 } { 1 2 n [ i e i n π + 1 π n ( e i n π 1 ) ] } , n Z {\displaystyle \{C_{n}\}={\Big \{}{\frac {\pi }{4}}{\Big \}}\cup {\Bigg \{}{\frac {1}{2n}}{\bigg [}ie^{in\pi }+{\frac {1}{\pi n}}(e^{in\pi }-1){\bigg ]}{\Bigg \}},n\in \mathbb {Z} ^{*}}

Considerem un altre exemple senzill; sigui la funció f ( t ) = { 0 si t ( 1 , 1 2 ) ( 1 2 , 1 ] 1 si t ( 1 2 , 1 2 ) 1 2 si t = ± 1 2 {\displaystyle f(t)={\begin{cases}0&{\text{si}}&t\in {\Big (}-1,-{\frac {1}{2}}{\Big )}\cup {\Big (}{\frac {1}{2}},1{\Big ]}\\1&{\text{si}}&t\in {\Big (}-{\frac {1}{2}},{\frac {1}{2}}{\Big )}\\{\frac {1}{2}}&{\text{si}}&t=\pm {\frac {1}{2}}\end{cases}}} i f ( t ) = f ( t + 2 ) {\displaystyle f(t)=f(t+2)} .

El seu espectre és:

{ C n } = { 1 2 } { sin n π 2 n π } , n Z {\displaystyle \{C_{n}\}={\Big \{}{\frac {1}{2}}{\Big \}}\cup {\Bigg \{}{\frac {\sin {\frac {n\pi }{2}}}{n\pi }}{\Bigg \}},n\in Z^{*}}

Referències

  1. E. Boyce, William; C. DiPrima, Richard. Elementary Differential Equation and Boundary Value Problem (en anglès). 10a. JohnWiley & Sons, Inc, 2012, p. 597. ISBN 978-0-470-45831-0. 
  2. G. Gasquet, P.Witomski; Fourirer Analysis and Applications. Springer-Verlag, 1999; p.45.
  3. Hardy, G.H.; Rogosinski, W. «Chapter IV: Convergence of Fourier series». A: Fourier Series. Courier Dover Publications. Cambridge University Press, 1999, p. 37. ISBN 0486406814. 
  4. Jahnke, H. N.. «6.5 Convergence of Fourier series». A: A History of Analysis. American Mathematical Society, 2003, p. 178. ISBN 0821826239. 
  5. Florinsky, Igor. Digital Terrain Analysis in Soil Science and Geology. Academic Press, 2011, p. 88. ISBN 0123850371.  Veure il·lustració a "Figure 5.8: Approximation of a square wave".
  6. Jerri, Abdul J. The Gibbs Phenomenon in Fourier Analysis, Slines and Wavelet Approximations. Springer, 1998. ISBN 0792351096. 

Vegeu també

  • Vegeu aquesta plantilla
Funcions
trigonomètriques
Sinus (sin) · Cosinus (cos) · Tangent (tan) · Cotangent (cot) · Secant (sec) · Cosecant (csc) · Versinus (versin) · Coversinus (coversin) · Semiversinus (semiversin) · Vercosinus (vercos) · Exsecant (exsec) · Excosecant (excsc)
Funcions
trigonomètriques
inverses
Arcsinus (arcsin) · Arccosinus (arccos) · Arctangent (arctan) · Arccotangent (arccotan) · Arcsecant (arcsec) · Arccosecant (arccosec)
Teoremes
Fòrmules
Vegeu també
Registres d'autoritat