Sèrie de Fourier

De testwiki
La revisió el 20:19, 5 feb 2025 per imported>Gvc1987 (growthexperiments-addlink-summary-summary:3|0|0)
(dif.) ← Versió més antiga | Versió actual (dif.) | Versió més nova → (dif.)
Salta a la navegació Salta a la cerca

Plantilla:Transformades de Fourier

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), per exemple de període 2π, es desitja escriure-la com una combinació en la qual intervinguin únicament sinus i cosinus, que són les funcions de període 2π simples més conegudes:[1]

f(t)=a02+n=1[ancos(nt)+bnsin(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 Plantilla:Segle. 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: és periòdica de període  T>0

si és tal que

t,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 fonamentalf0=1T,ω0=2πf0).

Forma general

Sigui f(x) una funció complexa f definida en el domini dels nombres reals, on t f:, definida a trossos, periòdica en f(t) amb període T, i de quadrat integrable a l'interval de t1 a t2 de durada T:

t1t2|f(t)|2dt<+

amb T=t2t1, llavors es definex la sèrie trigonomètrica de Fourier:

f(t)=a02+n=1[ancos(ωnt)+bnsin(ωnt)]

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

  • ωn=n2πT, on n és el nombre d'harmònics de la funció f,
  • an=2Tt1t2f(t)cos(ωnt)dt, és la part parella dels coeficients de Fourier de la funció f,
  • bn=2Tt1t2f(t)sin(ωnt)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.

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

f(t)=a02+n=1[ancos(nt)+bnsin(nt)]

on

  • an=1πππf(t)cos(nt)dt
  • bn=1πππf(t)sin(nt)dt

Forma exponencial

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

Cn=12πππf(x)einxdx

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

n=1Cneinx+n=0Cneinx

De manera més compacta:

n=Cneinx

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πTt. 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ó
αxp(t)+βyp(t)αXs[k]+βYs[k]
  • Derivada
x'p(t)jk2πf0Xs[k](k0)
  • Integral
0txp(t)dtX[k]jk2πf0+C(k0)
  • Retard
xp(tα)Xs[k]ej2kπf0α
  • Escalat
xp(t) Xs[k] harmònics en f=kf0α
xp(t) Xs[k]=Xs*[k]
  • Modulació
cos(m2πf0t)xp(t)12{Xs[km]+Xs[k+m]}
12{xp(t+α)+xp(tα)}cos(2πf0α)Xs[k]
  • Convolució
xp(t)yp(t)Xs[k]*Ys[k]
xp(t)yp(t) Xs[k]Ys[k]

Exemple

Escriure la Sèrie de Fourier corresponent a la funció f definida a continuació, a l'interval π < x < π:

f(x)={0π<x0x0x<π

Els coeficients de Fourier són:

Fig.1

a0=1πππf(x)dx

=1π(0+0πxdx)

a0=π2,

an=1π0πxcos(nx)dx

=[cosnx+nxsinnxπn2]0π

an=(1)n1πn2,

bn=1π0πxsin(nx)dx

=[sinnxnxcosnxπn2]0π

bn=(1)nn


on n =1, 2,... Per tant, a l'interval (π,π),

f(x)π4+n>0[(1)n1πn2cosnx(1)nnsinnx]

Es veu així per què aquesta Sèrie de Fourier convergeix a f(x) quan (π<x<π). 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π. La teoria mostrarà que la suma de la sèrie a cadascun d'aquests punts ha de ser π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=π42πcosx+sinx12sin2x.

Convergència de la Sèrie de Fourier

Sigui f: contínua a trossos i periòdica de període T, i sigui

f(t)a02+n=1[ancos(nωt)+bnsin(nωt)]

on ω=2π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: 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] é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: una funció periòdica de període T que satisfà les condicions de Dirichlet i sigui

f(t)a02+n=1[ancos(nωt)+bnsin(nωt)]

on ω=2π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
a02+n=1[ancos(nωt)+bnsin(nω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
a02+n=1[ancos(nωt)+bnsin(nωt)]=f(t)+f(t+)2

on f(t)=limτ0,τf(tτ) indica el límit de f a t per l'esquerra i f(t+)=limτ0,τf(t+τ) 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, llavors la suma de la Sèrie de Fourier coincideix amb f(t) a cada t. 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, de període a, és diferenciable en l'interval [0,a], excepte potser en un nombre finit de punts, i la derivada f és contínua, aleshores:[2]

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

2- Els coeficients de Fourier de f satisfan:

n=+|cn(f)| < 

3- La sèrie de Fourier de f convergeix uniformement cap a f en .[3][4]

Fenomen de Gibbs

Plantilla:AP

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 12π0πsin(t)t1dt0.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)=a02+n=1[ancos(nωt)+bnsin(nωt)]

aleshores

f(t)n=1[nbnωcos(nωt)+nanωsin(nωt)]

per cada t

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)=a02+n=1[ancos(nωt)+bnsin(nωt)]

i t0,t[T2,T2], aleshores

t0tf(τ)dτ=12a0(tt0)+n=1[bn(cos(nωt0)cos(nωt))nω+an(sin(nωt)sin(nωt0))nω]

S'ha d'anar amb compte, ja que al terme 12a0t 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 [T2,T2] per t0tf(τ)dτ12a0t i repetida periòdicament.

Igualtat de Parseval

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

f(t)a02+n=1[ancos(nωt)+bnsin(nωt)]

la seva Sèrie de Fourier. Aleshores

1T0T[f(t)]2dt=14a02+12n=1[an2+bn2]

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=1T0T[f(t)]2dt 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)a02+n=1[ancos(nωt)+bnsin(nωt)]

o bé en termes l'amplitud i fase

f(t)α0+n=1αncos(nωt+ϕn)

on ω=2πT i els coeficients d'aquestes expressions estan relacionades de la següent manera: α0=12a0 i per n=1,2...

αn=an2+bn2
i
ϕn=arctan(bnan)(π,π]

o bé an=cos(ϕn)

i

bn=αnsin(ϕ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 ejnωt en lloc de cos(nωt) i sin(nωt). Per construir-les, substituïm les expressions

cos(nωt)=ejnωt+ejnωt2
i
sin(nωt)=ejnωtejnωt2

a la sèrie a02+n=1[ancos(nωt)+bnsin(nωt)], obtenint

f(t)a02+12n=1[(anjbn)ejnωt+(an+jbn)ejnωt]

Escrivint

c0=c0e0=12a0, cn=anjbn2, cn=cn=an+jbn2

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+cnejnωt

que s'anomena notació complexa de la Sèrie de Fourier de f. Els nombres cn són els coeficients de Fourier de la notació complexa i es poden calcular de forma directa

cn=1TT2T2f(t)ejnωtdt per n=0,±1,±2,

En canvi si el que tenim són els coeficients cn de la notació complexa, aleshores els coeficients de Fourier an i bn, les amplituds i les fases venen donades per

an=2Re(cn),bn=2Im(cn),αn=2|cn| i ϕn=arg(cn)(π,π]

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ω,

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=+cnejnωt i g(t)=n=+dnejnωt

on ω=2πT

  • Linealitat: si p i q son nombres complexos
pf(t)+qg(t)=+(pcn+qdn)ejnωt
  • Translació en el temps: t0 és un nombre real.
f(tt0)=n=+cnejnωt0ejnωt=n=+cnejnω(tt0)
f(t) i f(tt0) 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ω.
f(pt)=n=+cnejn(pω)t
f(t) i f(pt) tenen les mateixes amplituds y fases però corresponen a freqüències diferents
  • Convolució
hn=n=+ckdnk=n=+cnkdk
  • Multiplicació
1T0Tf(t)g(t)dt=n=+cndn
  • Forma complexa de la igualtat de Parseval
1T0T|f(t)2|dt=n=+|cn2|
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=+|cn|2. Per aquest fet, la representació dels valors |cn|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 [T2,T2] per

f(n)={k si 0|t|<d20 si d2t|<T2

són

cn=kTd2d2ejnωtdt=kdTsin(nωd2)nωd2

Els valor de cn són reals, el que correspon el fet que la funció sigui parella, així doncs el seu espectre de fases és ϕn=0 per cada n. La funció

sa(t):={1 si t=0sin(t)t si t0

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:

cn=kdTsa(nωd2)

A vegades, en lloc de la funció de mostreig es treballa amb una altra funció, anomenada sinus cardinal que es defineix com sinc(t):=sa(πt):={1 si t=0sin(πt)πt si t0

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:

Ck=1Tt1t2f(t)ejωktdt

Per tant:

f(t)=k=+Ckejωkt

Espectre d'una funció periòdica

Sigui S(x) una funció periòdica, de període T=n2πωn, que es pot escriure com:

f(t)=a02+n=1(ancos(ωnt)+bnsin(ωnt))=n=Cneiωnt

Aleshores direm que la sèrie de Fourier anterior és la descomposició espectral de la funció f(t). El conjunt de coeficients an i bn o, el conjunt de coeficients Cn s'anomena l'espectre de la funció f(t). L'espectre d'una funció, per tant, conté tota la informació d'aquesta i, si considerem la funció 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π, f(x)={0π<x0x0x<π

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

{a0}{an}{bn}={π2}{(1)n1πn2}{(1)nn},n*

o, equivalentment

{Cn}={π4}{12n[ieinπ+1πn(einπ1)]},n*

Considerem un altre exemple senzill; sigui la funció f(t)={0sit(1,12)(12,1]1sit(12,12)12sit=±12 i f(t)=f(t+2).

El seu espectre és:

{Cn}={12}{sinnπ2nπ},nZ*

Referències

Plantilla:Commonscat

Plantilla:Referències

Vegeu també

Plantilla:Trigonometria Plantilla:Autoritat

  1. Plantilla:Ref-llibre
  2. Error de citació: Etiqueta <ref> no vàlida; no s'ha proporcionat text per les refs nomenades convunif
  3. Plantilla:Ref-llibre
  4. Plantilla:Ref-llibre
  5. Plantilla:Ref-llibre Veure il·lustració a "Figure 5.8: Approximation of a square wave".
  6. Plantilla:Ref-llibre