Serie Fourier

Seriile Fourier sunt o unealtă matematică folosită pentru a analiza funcțiile periodice descompunându-le într-o sumă ponderată de funcții sinusoidale componente care sunt uneori denumite armonice Fourier normale, sau pe scurt armonice. Printre generalizări se numără seriile Fourier generalizate.

Seriile Fourier au multe utilizări practice, pentru că manipularea și conceptualizarea coeficienților armonici sunt adesea mai ușoare decât lucrul cu funcția originală. Domeniile de aplicabilitate includ ingineria electrică, analiza undelor, acustică, optică, prelucrarea semnalelor și a imaginilor, și compresia datelor. Folosind uneltele și tehnicile spectroscopiei, de exemplu, astronomii pot deduce compoziția chimică a unei stele prin analizarea componentelor armonice, sau spectrului, stelei care emite lumină. Analog, inginerii pot optimiza proiectarea unui sistem de telecomunicații cu ajutorul informațiilor pe care le oferă componentele spectrale ale unui semnal de date pe care sistemul le transportă.

Seriile Fourier sunt numite după omul de știință și matematicianul francez Joseph Fourier, care le-a folosit în importanta sa lucrare despre conducția termică, Théorie Analytique de la Chaleur (Teoria analitică a căldurii), publicată în 1822.

Definiție

Forma generală

Dată fiind o funcție cu valori complexe f de argument real t, f: RC, unde f(t) este continuă și derivabilă pe porțiuni, periodică de perioadă T, și integrabilă la pătrat pe intervalul de lungime T dintre t 1 {\displaystyle t_{1}} și t 2 {\displaystyle t_{2}} , adică

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

unde

  • T = t 2 t 1 {\displaystyle T=t_{2}-t_{1}} este perioada,
  • t 1 {\displaystyle t_{1}} si t 2 {\displaystyle t_{2}} sunt limitele de integrare.

Dezvoltarea în serie Fourier a lui f este

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

unde, pentru orice întreg nenegativ n,

  • ω n = n 2 π T {\displaystyle \omega _{n}=n{\frac {2\pi }{T}}}     este armonica n (în radiani) a funcției f,
  • a 0 = 2 T t 1 t 2 f ( t ) d t {\textstyle a_{0}={\frac {2}{T}}\int _{t_{1}}^{t_{2}}f(t)dt}
  • 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}     sunt coeficienții Fourier pari ai lui f, iar
  • 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}     sunt coeficienții Fourier impari ai lui f.

Echivalent, în formă cu exponențiala complexă,

  • f ( t ) = n = + c n e i ω n t {\displaystyle f(t)=\sum _{n=-\infty }^{+\infty }c_{n}e^{i\omega _{n}t}}

unde:

  • c n = 1 T t 1 t 2 f ( t ) e i ω n t d t {\displaystyle c_{n}={\frac {1}{T}}\int _{t_{1}}^{t_{2}}f(t)e^{-i\omega _{n}t}\,dt}
  • i {\displaystyle i} este unitatea imaginară, și
  • e i ω n t = cos ( ω n t ) + i sin ( ω n t ) {\displaystyle e^{i\omega _{n}t}=\cos(\omega _{n}t)+i\sin(\omega _{n}t)}     conform cu formula lui Euler.

Forma canonică

În cazul special când perioada T = 2π, avem

ω n = n {\displaystyle \omega _{n}=n\,}

În acest caz, dezvoltarea în serie Fourier se reduce la o formă deosebit de simplă:

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

unde

  • 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}

oricare ar fi un întreg nenegativ n.

sau, echivalent:

f ( t ) = n = + c n e i n t {\displaystyle f(t)=\sum _{n=-\infty }^{+\infty }c_{n}e^{int}}

unde

  • c n = 1 2 π π π f ( t ) e i n t d t = 1 2 ( a n i b n ) . {\displaystyle c_{n}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(t)e^{-int}\,dt={\frac {1}{2}}(a_{n}-ib_{n}).}

Alegerea formei

Forma pentru perioada T poate fi ușor derivată din cea canonică cu schimbarea de variabilă definită de x = 2 π T t {\displaystyle x={\frac {2\pi }{T}}t} . De aceea, ambele formulări sunt echivalente. Totuși, forma pentru perioada T este preferată în majoritatea cazurilor practice pentru că este direct aplicabilă. Pentru teorie, este preferată forma canonică deoarece este mai elegantă și mai ușor de interpretat matematic.

Exemple

Serii Fourier simple

Fie f periodică de perioadă 2 π {\displaystyle 2\pi } , cu f ( x ) = x {\displaystyle f(x)=x} pentru x între −π și π. Se observă că această funcție este o versiune periodică a funcției identitate.

Graficul unei funcții identitate periodice - o undă dinți de fierăstrău.
Grafic animat cu primele cinci sume parțiale ale seriilor Fourier pentru funcția identitate periodică

Se vor calcula coeficienții Fourier pentr această funcție.

a n = 1 π π π f ( x ) cos ( n x ) d x = 1 π π π x cos ( n x ) d x = 0. {\displaystyle {\begin{aligned}a_{n}&{}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(x)\cos(nx)\,dx\\&{}={\frac {1}{\pi }}\int _{-\pi }^{\pi }x\cos(nx)\,dx\\&{}=0.\end{aligned}}}


b n = 1 π π π f ( x ) sin ( n x ) d x = 1 π π π x sin ( n x ) d x = 2 π 0 π x sin ( n x ) d x = 2 π ( [ x cos ( n x ) n ] 0 π + [ sin ( n x ) n 2 ] 0 π ) = 2 ( 1 ) n + 1 n . {\displaystyle {\begin{aligned}b_{n}&{}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(x)\sin(nx)\,dx\\&{}={\frac {1}{\pi }}\int _{-\pi }^{\pi }x\sin(nx)\,dx\\&{}={\frac {2}{\pi }}\int _{0}^{\pi }x\sin(nx)\,dx\\&{}={\frac {2}{\pi }}\left(\left[-{\frac {x\cos(nx)}{n}}\right]_{0}^{\pi }+\left[{\frac {\sin(nx)}{n^{2}}}\right]_{0}^{\pi }\right)\\&{}=2{\frac {(-1)^{n+1}}{n}}.\end{aligned}}}

Se observă că an sunt 0 deoarece x x cos ( n x ) {\displaystyle x\mapsto x\cos(nx)} sunt funcții pare. Deci seria Fourier pentru această funcție este:

f ( x ) = a 0 2 + n = 1 [ a n cos ( n x ) + b n sin ( n x ) ] {\displaystyle f(x)={\frac {a_{0}}{2}}+\sum _{n=1}^{\infty }\left[a_{n}\cos \left(nx\right)+b_{n}\sin \left(nx\right)\right]}
= 2 n = 1 ( 1 ) n + 1 n sin ( n x ) , x [ π , π ] . {\displaystyle =2\sum _{n=1}^{\infty }{\frac {(-1)^{n+1}}{n}}\sin(nx),\quad \forall x\in [-\pi ,\pi ].}

O aplicație a acestei serii Fourier este calculul funcției Riemann zeta la s = 2; Conform teoremei lui Parseval, avem:

1 2 π π π x 2 d x = 1 2 n > 0 [ 2 ( 1 ) n n ] 2 {\displaystyle {\frac {1}{2\pi }}\int _{-\pi }^{\pi }x^{2}dx={\frac {1}{2}}\sum _{n>0}\left[2{\frac {(-1)^{n}}{n}}\right]^{2}}

de unde rezultă: n > 0 1 n 2 = π 2 6 {\displaystyle \sum _{n>0}{\frac {1}{n^{2}}}={\frac {\pi ^{2}}{6}}} .

Ecuația undei

Ecuația undei descrie mișcarea unei coarde vibrante, care poate fi ținuta fixă de capete. Soluția acestei probleme necesită dezvoltarea în serie trigonometrică a unei funcții generale f care dispare la capetele unui interval de la x=0 la x=L. Seria Fourier pentru o asemenea funcție ia forma

f ( x ) = n = 1 b n sin ( n π L x ) {\displaystyle f(x)=\sum _{n=1}^{\infty }b_{n}\sin \left({\frac {n\pi }{L}}x\right)}

unde

b n = 2 L 0 L f ( x ) sin ( n π L x ) d x . {\displaystyle b_{n}={\frac {2}{L}}\int _{0}^{L}f(x)\sin \left({\frac {n\pi }{L}}x\right)\,dx.}

Vibrațiile aerului într-o țeavă deschisă la un capăt și închisă la celălalt sunt și ele descrise de ecuația undei. Soluția ei este dată de dezvoltarea în serie Fourier a unei funcții care dispare la x = 0 și care nu mai este derivabilă la x=L. Seria sa Fourier ia forma

f ( x ) = n = 1 b n sin ( ( 2 n + 1 ) π 2 L x ) {\displaystyle f(x)=\sum _{n=1}^{\infty }b_{n}\sin \left({\frac {(2n+1)\pi }{2L}}x\right)}

unde

b n = 2 L 0 L f ( x ) sin ( ( 2 n + 1 ) π 2 L x ) d x . {\displaystyle b_{n}={\frac {2}{L}}\int _{0}^{L}f(x)\sin \left({\frac {(2n+1)\pi }{2L}}x\right)\,dx.}

Interpretare: descompunerea unei mișcări în rotații

Seriile Fourier au o interpretare cinematică. Funcția t f ( t ) {\displaystyle t\mapsto f(t)} poate fi văzută ca mișcare a unui obiect într-un plan (t reprezentând timpul). Deoarece f ia valori complexe, se poate scrie

f ( t ) = u ( t ) + i v ( t ) . {\displaystyle f(t)=u(t)+iv(t).\,}

pentru funcții cu valori reale u și v. În această formă, putem interpreta f ca sumă de translații orizontale și verticale.

De la un moment t {\displaystyle t} la un moment t + d t {\displaystyle t+dt} , unde dt e o perioadă incrementală mică, obiectul se mișcă de la punctul A = [ u ( t ) v ( t ) ] {\displaystyle A=\left[{\begin{matrix}u(t)\\v(t)\end{matrix}}\right]} la punctul B = [ u ( t + d t ) v ( t + d t ) ] {\displaystyle B=\left[{\begin{matrix}u(t+dt)\\v(t+dt)\end{matrix}}\right]} , care corespunde unei translații infinitezimale în spațiu după vectorul A B = [ u ( t + d t ) u ( t ) v ( t + d t ) v ( t ) ] {\displaystyle {\overrightarrow {AB}}=\left[{\begin{matrix}u(t+dt)-u(t)\\v(t+dt)-v(t)\end{matrix}}\right]} . Rezultă că f se poate scrie de forma:

f ( t ) = [ u ( d t ) u ( 0 ) v ( d t ) v ( 0 ) ] + [ u ( 2 d t ) u ( d t ) v ( 2 d t ) v ( d t ) ] + + [ u ( t + d t ) u ( t ) v ( t + d t ) v ( t ) ] {\displaystyle f(t)=\left[{\begin{matrix}u(dt)-u(0)\\v(dt)-v(0)\end{matrix}}\right]+\left[{\begin{matrix}u(2dt)-u(dt)\\v(2dt)-v(dt)\end{matrix}}\right]+\cdots +\left[{\begin{matrix}u(t+dt)-u(t)\\v(t+dt)-v(t)\end{matrix}}\right]}
= 0 t 1 d x [ u ( x + d x ) u ( x ) v ( x + d x ) v ( x ) ] d x . {\displaystyle =\int _{0}^{t}{\frac {1}{dx}}\left[{\begin{matrix}u(x+dx)-u(x)\\v(x+dx)-v(x)\end{matrix}}\right]\,dx.}

În loc de a vedea pe f ca pe o sumă de translații infinitezimale, o putem vedea ca pe o sumă infinită de rotații cu diferite raze. Această interpretare este convenabilă, mai ales când mișcarea este periodică.

Fie χ n = e i n x {\displaystyle \chi _{n}=e^{inx}} rotația de n-ture pe secundă, de rază 1). Se dorește scrierea f ca f ( x ) = c n χ n {\displaystyle f(x)=\sum c_{n}\chi _{n}} . Se poate demonstra că razele de rotație (coeficienții c n {\displaystyle c_{n}} ) sunt exact cei dați în paragraful anterior.

De exemplu, graficul funcției f : t 2 cos ( t 2 ) e 3 2 i t {\displaystyle f:t\mapsto 2\cos \left({\frac {t}{2}}\right)e^{{\frac {3}{2}}it}} este închis, ceea ce înseamnă că funcția este periodică. Bucla de pe curbă sugerează că este suma a două funcții periodice, una cu o perioadă mai scurtă decât cealaltă. Într-adevăr, sepoate scrie: f ( t ) = e i t + e 2 i t = χ 1 ( t ) + χ 2 ( t ) {\displaystyle f(t)=e^{it}+e^{2it}=\chi _{1}(t)+\chi _{2}(t)} . Toți coeficienții săi Fourier sunt zero cu excepția c 1 = 1 {\displaystyle c_{1}=1} și c 2 = 1 {\displaystyle c_{2}=1} . Interpretarea grafică a unei rotații este mult mai dificil de realizat decât a translațiilor, pentru că în loc de a imagina vizual mișcarea dintr-u punct în altul, trebuie adăugată toată mișcarea pentru a avea sens descompunerea (gândirea se face în termeni de frecvență de rotație și nu de timp).

Matematic, adoptarea acestui punct de vedere reprezintă o folosire a seriilor Fourier ca pe o unealtă de înțelegere a operatorilor liniari care comută cu translația. Funcțiile χ n {\displaystyle \chi _{n}} sunt exact caracterele multiplicative ale grupului R / 2 π Z {\displaystyle \mathbb {R} /2\pi \mathbb {Z} } .

Istoric

Context

Seriile Fourier au fost denumite în onoarea lui Joseph Fourier (1768-1830), care a avut importante contribuții la studiul seriilor trigonometrice, după investigații preliminare ale lui Madhava, Nilakantha Somayaji, Jyesthadeva, Leonhard Euler, Jean le Rond d'Alembert și Daniel Bernoulli. El a aplicat această tehnică pentru a găsi soluția pentru ecuația căldurii, publicându-și rezultatele inițiale în 1807 Arhivat în , la Wayback Machine. și 1811, și publicând lucrarea Théorie analytique de la chaleur în 1822.

Dintr-un punct de vedere modern, rezultatele lui Fourier sunt puțin informale, datorită unei lipse de notație precisă a funcției și integralei la începutul secolului al XIX-lea (de exemplu, la acea vreme nu era lămurit dacă o funcție definită pe două intervale diferite cu două expresii diferite mai era tot o singură funcție). Mai târziu, Dirichlet și Riemann au exprimat rezultatele lui Fourier cu mai mare precizie și rigurozitate.

Un articol revoluționar

În lucrarea lui Fourier intitulată Mémoire sur la propagation de la chaleur dans les corps solides, la paginile 218 și 219, se pot citi următoarele:

φ ( y ) = a cos π y 2 + a cos 3 π y 2 + a cos 5 π y 2 + . {\displaystyle \varphi (y)=a\cos {\frac {\pi y}{2}}+a'\cos 3{\frac {\pi y}{2}}+a''\cos 5{\frac {\pi y}{2}}+\cdots .}
Înmulțind ambele părți cu cos ( 2 i + 1 ) π y 2 {\displaystyle \cos(2i+1){\frac {\pi y}{2}}} , și apoi integrând de la y = 1 {\displaystyle y=-1} la y = + 1 {\displaystyle y=+1} rezultă:
a i = 1 1 φ ( y ) cos ( 2 i + 1 ) π y 2 d y . {\displaystyle a_{i}=\int _{-1}^{1}\varphi (y)\cos(2i+1){\frac {\pi y}{2}}\,dy.}

În aceste câteva rânduri, surprinzător de apropiate de formalismul modern folosit în seriile Fourier, Fourier a revoluționat fără să vrea atât matematica cât și fizica. Deși serii trigonometrice similare mai fuseseră folosite anterior de Euler, d'Alembert, Daniel Bernoulli și Gauss, Fourier a fost primul care a recunoscut că astfel de serii trigonometrice pot reprezenta funcții arbitrare, chiar și funcții cu discontinuități. A trebuit să treacă mulți ani pentru a clarifica această idee, care a condus la importante teorii asupra convergenței, spațiilor de funcții, și analizei armonice.

Lucrarea lui Fourier a fost atât de originală încât, când a trimis-o în 1807, comisia (compusă din alți mari matematicieni ca Lagrange, Laplace, Malus și Legendre, printre alții) a concluzionat: ...maniera în care autorul ajunge la aceste ecuații nu este una scutită de dificultăți și analiza sa încă lasă de dorit la capitolele generalitate și chiar rigoare.

Nașterea analizei armonice

Din vremea lui Fourier până astăzi au fost descoperite multe alte abordări ale definirii și înțelegerii conceptului de serie Fourier, toate fiind corecte și echivalente matematic, dar fiecare punând accent pe alte aspecte ale subiectului. Unele dintre cele mai elegante și mai puternice abordări sunt bazate pe idei și unelte matematice care nu erau disponibile la momentul când Fourier și-a terminat lucrarea. Fourier a definit inițial seriile Fourier pentru funcții reale de argument real, și folosindu-se de funcțiile sinus și cosinus ca bază a descompunerii.

Multe alte transformări de tip Fourier au mai fost definite, extinzând cu alte noi aplicații ideea inițială de reprezentare a oricărei funcții periodice ca suprapunere de armonice. Această arie generală de studiu este uneori numită analiză armonică.

Derivări moderne ale coeficienților Fourier

Metoda lui Fourier de calcul al coeficienților seriei este foarte practică și potrivită problemei pe care o tratează (propagarea căldurii). Totuși, această metodă a fost generalizată între timp la o clasă mult mai largă de probleme: scrierea unei funcții ca sumă de funcții periodice.

Mai exact, dacă f:RC este o funcție, am vrea să scriem această funcție ca sumă de funcții trigonometrice, i.e. f ( x ) = c n e i n x {\displaystyle f(x)=\sum c_{n}e^{inx}} . Variantele de alegere a funcțiilor trebuie să fie restrânse pentru ca descompunerea să aibă sens. În primul rând, dacă f este de perioadă T, atunci prin schimbarea variabilelor, se poate studia x f ( T 2 π x ) {\displaystyle x\mapsto f\left({\frac {T}{2\pi }}x\right)} de perioadă 2π. Aceasta simplifică mult notația și permite utilizarea unei forme standard (canonice). Putem restrânge astfel studiul funcțiilor x f ( T 2 π x ) {\displaystyle x\mapsto f\left({\frac {T}{2\pi }}x\right)} la orice interval de lungime 2π, cum ar fi [-π,π].

Luând funcția f:RC din mulțimea celor continue pe porțiuni, funcții periodice de perioadă 2π cu π π | f ( x ) | 2 d x < + {\displaystyle \int _{-\pi }^{\pi }|f(x)|^{2}\,dx<+\infty } . De fapt, se iau funcții din spațiul Lp L2(μ), unde μ este măsura Lebesgue normalizată a intervalului [-π,π] (astfel încât [ π , π ] f d μ = 1 2 π π π f ( x ) d x {\displaystyle \int _{[-\pi ,\pi ]}f\,d\mu ={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x)\,dx} .

Coeficienți Fourier complecși

Putem transforma L2(μ) într-un spațiu Hilbert, ceea ce este potrivit pentru proiecții ortogonale, prin definirea produsului scalar:

f , g = [ π , π ] f g ¯ d μ = 1 2 π π π f ( x ) g ( x ) ¯ d x , {\displaystyle \langle f,g\rangle =\int _{[-\pi ,\pi ]}f{\overline {g}}\,d\mu ={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(x){\overline {g(x)}}\,dx,}

unde f ( x ) ¯ {\displaystyle {\overline {f(x)}}} reprezintă conjugata lui f(x). Vom nota cu {\displaystyle \|\cdot \|} norma asociată.

E = { t e i n t , n Z } {\displaystyle E=\{t\mapsto e^{int},n\in \mathbb {Z} \}} este o bază ortonormală din L2(μ), deci se poate scrie

f ( x ) = n Z f , e i n x e i n x . {\displaystyle f(x)=\sum _{n\in \mathbb {Z} }\left\langle f,e^{inx}\right\rangle e^{inx}.}

De regulă se definește n Z , c n = f , e i n x {\displaystyle \forall n\in \mathbb {Z} ,c_{n}=\left\langle f,e^{inx}\right\rangle } . Aceste numere se numesc coeficienți Fourier complecși. Expresia lor este

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.\,}

O formulare echivalentă este scrierea f ca sumă de funcții sinus și cosinus.

Coeficienți Fourier reali

Suma din secțiunea anterioară este simetrică în raport cu 0: într-adevăr, cu excepția lui n = 0, un coeficient cn corespunde fiecărui coeficientcn. Astfel ne amintim de formulele

cos x = e i x + e i x 2         a n d         sin x = e i x e i x 2 i . {\displaystyle \cos x={\frac {e^{ix}+e^{-ix}}{2}}{\rm {~~~~and~~~~}}\sin x={\frac {e^{ix}-e^{-ix}}{2i}}.}

Astfel se pot exprima serii Fourier cu funcții cu valori reale. Pentru a face aceasta se observă că

f ( x ) = n Z c n e i n x = c 0 + n > 0 [ c n e i n x + c n e i n x ] . {\displaystyle f(x)=\sum _{n\in \mathbb {Z} }c_{n}e^{inx}=c_{0}+\sum _{n>0}\left[c_{-n}e^{-inx}+c_{n}e^{inx}\right].}

După înlocuirea lui cn cu expresia sa și simplificarea rezultatului, obținem

f ( x ) = c 0 + n > 0 [ 1 π ( π π f ( t ) cos ( n t ) d t ) cos ( n x ) + 1 π ( π π f ( t ) sin ( n t ) d t ) sin ( n x ) ] . {\displaystyle f(x)=c_{0}+\sum _{n>0}\left[{\frac {1}{\pi }}\left(\int _{-\pi }^{\pi }f(t)\cos \left(nt\right)\,dt\right)\cos \left(nx\right)+{\frac {1}{\pi }}\left(\int _{-\pi }^{\pi }f(t)\sin \left(nt\right)\,dt\right)\sin \left(nx\right)\right].}

Dacă pentru un număr întreg nenegativ n, definim coeficienții Fourier reali an și bn prin

a n = 1 π π π f ( x ) cos ( n x ) d x , {\displaystyle a_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(x)\cos \left(nx\right)\,dx,}
b n = 1 π π π f ( x ) sin ( n x ) d x , {\displaystyle b_{n}={\frac {1}{\pi }}\int _{-\pi }^{\pi }f(x)\sin \left(nx\right)\,dx,}

obținem:

f ( x ) = a 0 2 + n > 0 [ a n cos ( n x ) + b n sin ( n x ) ] . {\displaystyle f(x)={\frac {a_{0}}{2}}+\sum _{n>0}\left[a_{n}\cos \left(nx\right)+b_{n}\sin \left(nx\right)\right].}

Proprietăți

  • Următoarele proprietăți se deduc din formula lui Euler: e i x = cos x + i sin x . {\displaystyle e^{ix}=\cos x+i\sin x.\!}
a n = c n + c n {\displaystyle a_{n}=c_{n}+c_{-n}} și b n = i ( c n c n )  for all  n  and  {\displaystyle b_{n}=i(c_{n}-c_{-n}){\mbox{ for all }}n{\mbox{ and }}\,}
c n = a n i b n 2 {\displaystyle c_{n}={\frac {a_{n}-ib_{n}}{2}}} și c n = a n + i b n 2  for all  n . {\displaystyle c_{-n}={\frac {a_{n}+ib_{n}}{2}}{\mbox{ for all }}n.}
  • Dacă f este o funcție impară, atunci a n = 0 {\displaystyle a_{n}=0} oricare ar fi n {\displaystyle n} pentru că f ( x ) cos ( n π x T ) {\displaystyle f(x)\cos \left(n\pi {\frac {x}{T}}\right)} este tot pară, deci integrala sa pe intervalul [ T , T ] {\displaystyle [-T,T]} este zero. Dacă f este o funcție pară, atunci b n = 0 {\displaystyle b_{n}=0} dintr-un motiv similar.
  • Dacă f este continuă pe porțiuni, lim n + c n ( f ) = 0 {\displaystyle \lim _{n\rightarrow +\infty }c_{n}(f)=0} , lim n + c n ( f ) = 0 {\displaystyle \lim _{n\rightarrow +\infty }c_{-n}(f)=0} , lim n + a n ( f ) = 0 {\displaystyle \lim _{n\rightarrow +\infty }a_{n}(f)=0} și lim n + b n ( f ) = 0. {\displaystyle \lim _{n\rightarrow +\infty }b_{n}(f)=0.}
  • Dacă f este continuu diferențiabilă de k ori pe porțiuni, atunci se pot calcula ușor coeficienții Fourier ai lui f ( k ) {\displaystyle f^{(k)}} din cei ai lui f:
c n ( f ( k ) ) = ( i n ) k c n ( f ) , {\displaystyle c_{n}\left(f^{(k)}\right)=(in)^{k}c_{n}(f),}

unde f ( k ) {\displaystyle f^{(k)}} este derivata de ordin k a lui f.

  • Oricare ar fi un întreg pozitiv al lui k, dacă f este Ck − 1 și Ck pe porțiuni, atunci
lim n + | n k c n ( f ) | = 0 {\displaystyle \lim _{n\rightarrow +\infty }|n^{k}c_{n}(f)|=0} pentru că n k c n ( f ) = i k c n ( f ( k ) ) 0. {\displaystyle n^{k}c_{n}(f)=i^{-k}c_{n}\left(f^{(k)}\right)\rightarrow 0.}

Aceasta înseamnă că seria c n ( f ) {\displaystyle c_{n}(f)} este rapid descrescătoare.

Cazul general

Seriile Fourier exploatează periodicitatea funcției f dar dacă f este periodică în mai multe variabile, sau chiar f neperiodică? Aceste probleme i-au condus pe matematicieni și pe fizicienii teoreticieni să încerce să definească seriile Fourier pe orice grup G. Avantajul este, de exemplu, definirea seriilor Fourier pentru funcții de mai multe variabile. Seriile Fourier și transformata Fourier folosite în prelucrarea semnalelor devin astfel cazuri speciale ale acestei teorii și sunt mai ușor de interpretat.

Dacă G grup abelian local compact și T este cercul unitate, se poate defini dualul lui G prin G ^ = { χ : G T  omomorfism } {\displaystyle {\widehat {G}}=\{\chi :G\rightarrow \mathbb {T} {\mbox{ omomorfism}}\}} . Acestea constituie mulțimea rotațiilor pe cercul unitate și elementele sale se numesc caractere. Se poate defini un produs scalar , {\displaystyle \langle \cdot ,\cdot \rangle } pe C[G] prin: χ 1 , χ 2 = G χ 1 ( g ) χ 2 ( g ) ¯ d g {\displaystyle \langle \chi _{1},\chi _{2}\rangle =\int _{G}\chi _{1}(g){\overline {\chi _{2}(g)}}\,dg} . G ^ {\displaystyle {\widehat {G}}} este atunci bază ortonormală în C[G] în raport cu acest produs scalar. Fie f :GC. Coeficienții Fourier ai lui f sunt definiți prin: f ^ ( χ ) = f , χ {\displaystyle {\widehat {f}}(\chi )=\langle f,\chi \rangle } și avem f ( g ) = G ^ f ^ ( χ ) χ ( g ) d χ {\displaystyle f(g)=\int _{\widehat {G}}{\widehat {f}}(\chi )\chi (g)\,d\chi } . Dacă grupul este discret, atunci integrala se reduce la o sumă.

De exemplu, coeficienții Fourier ai acestui articol sunt obținuți luând G = R/2πZ. Se obține

G ^ = { χ n : t e i n t , n Z } {\displaystyle {\widehat {G}}=\{\chi _{n}:t\mapsto e^{int},n\in \mathbb {Z} \}}

și

c n ( f ) = f ^ ( χ n ) = G f ( g ) χ ( g ) ¯ d g = 1 2 π π π f ( t ) e i n t d t . {\displaystyle c_{n}(f)={\widehat {f}}(\chi _{n})=\int _{G}f(g){\overline {\chi (g)}}\,dg={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(t)e^{-int}\,dt.}

Funcțiile periodice în n dimensiuni pot fi definite pe un tor n-dimensional (funcția ia valoare pe fiecare punct de pe tor). Un astfel de tor este definit prin Tn = Rn/(2πZ)n. Pentru n = 1 torul este un cerc, pentru n = 2 este produsul cartezian a două cercuri, adică un tor obișnuit. Alegând G = Tn rezultă seria Fourier corespunzătoare.

Aproximarea și convergența seriilor Fourier

Definiția unei serii Fourier

Fie χ n ( x ) = e i n π x T {\displaystyle \chi _{n}(x)=e^{in\pi {\frac {x}{T}}}} . Se numește serie Fourier a funcției f seria c n χ n {\displaystyle \sum c_{n}\chi _{n}} . Pentru orice întreg pozitiv N, f N ( x ) = n = N N c n χ n ( x ) {\displaystyle f_{N}(x)=\sum _{n=-N}^{N}c_{n}\chi _{n}(x)} se numește a N-a sumă parțială a seriei Fourier a acestei funcții.

Aproximarea cu sume parțiale

Să zicem că dorim să găsim cea mai bună aproximare a lui f folosind doar funcțiile χ n {\displaystyle \chi _{n}} pentru n de la N {\displaystyle -N} la N. Fie T N = { p = n = N N x n χ n , x n C } {\displaystyle {\mathcal {T}}_{N}=\left\{p=\sum _{n=-N}^{N}x_{n}\chi _{n},x_{n}\in \mathbb {C} \right\}} . Încercăm sa găsim coeficienții ( x N , , x N ) {\displaystyle (x_{-N},\dots ,x_{N})} astfel încât f p {\displaystyle \|f-p\|} este minim (unde {\displaystyle \|\cdot \|} este norma).

Avem f p 2 = f 2 2 Re f , p + p 2 {\displaystyle \|f-p\|^{2}=\|f\|^{2}-2{\mbox{Re}}\langle f,p\rangle +\|p\|^{2}} , unde Re(z) notează partea reală a lui z.

f , p = n = N N x n ¯ f , χ n . {\displaystyle \langle f,p\rangle =\sum _{n=-N}^{N}{\overline {x_{n}}}\langle f,\chi _{n}\rangle .}

Teorema lui Parseval (ce poate fi dedusă independent din seriile Fourier) dă

p 2 = n = N N | x n | 2 . {\displaystyle \|p\|^{2}=\sum _{n=-N}^{N}|x_{n}|^{2}.}

Prin definiție, c n = f , χ n {\displaystyle c_{n}=\langle f,\chi _{n}\rangle } ; deci

f p 2 = f 2 + n = N N [ | c n x n | 2 | c n | 2 ] . {\displaystyle \|f-p\|^{2}=\|f\|^{2}+\sum _{n=-N}^{N}\left[|c_{n}-x_{n}|^{2}-|c_{n}|^{2}\right].}

Este clar că această expresie este minimă pentru x n = c n {\displaystyle x_{n}=c_{n}} și doar pentru această valoare.

Deci există un singur f N T N {\displaystyle f_{N}\in {\mathcal {T}}_{N}} astfel încât

f f N = min p T N { f p , p T N } , {\displaystyle \|f-f_{N}\|=\min _{p\in {\mathcal {T}}_{N}}\left\{\|f-p\|,p\in {\mathcal {T}}_{N}\right\},}

este dat de

f N ( x ) = n = N N c n χ n ( x ) , {\displaystyle f_{N}(x)=\sum _{n=-N}^{N}c_{n}\chi _{n}(x),}

unde

c n = 1 2 T T T f ( t ) χ n ( t ) d t . {\displaystyle c_{n}={\frac {1}{2T}}\int _{-T}^{T}f(t)\chi _{-n}(t)\,dt.}

Deci cea mai bună aproximare a lui f ce poate fi făcută folosind doar funcțiile χ n ( x ) = e i n π x T {\displaystyle \chi _{n}(x)=e^{in\pi {\frac {x}{T}}}} pentru n de la N {\displaystyle -N} la N este exact a N-a sumă parțială a seriei Fourier.

Convergența

În timp ce coeficienții Fourier an și bn pot fi definiți formal pentru orice funcție integrabilă, dacă această funcție converge sau nu la f(x) depinde de proprietățile lui f.

Cel mai simplu răspuns este că pentru ca seria să conveargă, f trebuie să fie integrabilă la pătrat, deci

lim N π π | f ( x ) n = N N c n χ n ( x ) | 2 d x = 0. {\displaystyle \lim _{N\rightarrow \infty }\int _{-\pi }^{\pi }\left|f(x)-\sum _{n=-N}^{N}c_{n}\,\chi _{n}(x)\right|^{2}\,dx=0.}

Aceasta este convergență în norma dată de spațiul L2. Demonstrația acestui rezultat este simplă, spre deosebire de rezultatul mult mai puternic dat de Lennart Carleson conform căruia seriile converg aproape în orice caz.

Există mai multe teste care asigură că seria converge într-un punct dat, de exemplu, dacă funcția este diferențiabilă în x. Nici chiar o mică discontinuitate a derivatei nu constituie o problemă: dacă funcția are derivată la stânga și la dreapta în x, atunci seria Fourier va converge la media limitelor la stânga și la dreapta (dar vezi Fenomenul Gibbs). Totuși, fapt considerat de mulți surprinzător, seria Fourier a unei funcții continue nu trebuie neapărat să fie convergentă în fiecare punct.

Această situație neplăcută este echilibrată de o teoremă a lui Dirichlet care afirmă că dacă f este periodică de perioadă 2 T {\displaystyle 2T} și derivată cu derivata continuă, atunci seria ei Fourier converge în fiecare punct și n Z c n χ n ( x ) = f ( x + ) + f ( x ) 2 {\displaystyle \sum _{n\in \mathbb {Z} }c_{n}\chi _{n}(x)={\frac {f(x^{+})+f(x^{-})}{2}}} , unde f ( x + ) = lim t x , t > x f ( x ) {\displaystyle f(x^{+})=\lim _{t\rightarrow x,t>x}f(x)} și f ( x ) = lim t x , t < x f ( x ) {\displaystyle f(x^{-})=\lim _{t\rightarrow x,t<x}f(x)} . Dacă f este continuă și cu derivata continuă pe porțiuni, atunci seria Fourier converge uniform.

În 1922, Andrei Kolmogorov a publicat un articol intitulat Une série de Fourier-Lebesgue divergente presque partout Arhivat în , la Wayback Machine. în care a dat un exemplu de funcție integrabilă Lebesgue a cărei serie Fourier divere aproape în fiecare punct. Această funcție nu este din L 2 ( μ ) {\displaystyle L^{2}(\mu )} .

Teoremele lui Plancherel și Parseval

O altă proprietate importantă a seriilor Fourier este Teorema lui Plancherel. Fie f , g L 2 ( μ ) {\displaystyle f,g\in L^{2}(\mu )} și c n ( f ) , c n ( g ) {\displaystyle c_{n}(f),c_{n}(g)} coeficienții Fourier complecși corespunzători. Atunci

n Z c n ( f ) c n ( g ) ¯ = 1 2 T T T f ( x ) g ( x ) ¯ d x {\displaystyle \sum _{n\in \mathbb {Z} }c_{n}(f){\overline {c_{n}(g)}}={\frac {1}{2T}}\int _{-T}^{T}f(x){\overline {g(x)}}\,dx}

unde cu z ¯ {\displaystyle {\overline {z}}} s-a notat conjugatul lui z.

Teorema lui Parseval, un caz special al teoremei lui Plancherel, afirmă că:

n Z | c n ( f ) | 2 = 1 2 T T T | f ( x ) | 2 d x {\displaystyle \sum _{n\in \mathbb {Z} }|c_{n}(f)|^{2}={\frac {1}{2T}}\int _{-T}^{T}|f(x)|^{2}\,dx}

reformulabilă astfel pentru coeficienți Fourier reali:

a 0 2 4 + 1 2 n = 1 ( a n 2 + b n 2 ) = 1 2 T T T | f ( x ) | 2 d x . {\displaystyle {\frac {a_{0}^{2}}{4}}+{\frac {1}{2}}\sum _{n=1}^{\infty }\left(a_{n}^{2}+b_{n}^{2}\right)={\frac {1}{2T}}\int _{-T}^{T}|f(x)|^{2}\,dx.}

Aceste teoreme se pot demonstra folosind relațiile de ortogonalitate. Ele pot fi interpretate fizic spunând că scrierea unui semnal ca serie Fourier series nu îi modifică energia.

Bibliografie

  • Joseph Fourier, tradusă în engleză de Alexander Freeman (publicată în 1822, tradusă 1878, republicată 2003). The Analytical Theory of Heat. Dover Publications. ISBN 0-486-49531-0.  Verificați datele pentru: |date= (ajutor) republicarea completă din 2003 a traducerii în engleză din 1878 de Alexander Freeman a lucrării lui Fourier Théorie Analytique de la Chaleur, inițial publicată în 1822.
  • Yitzhak Katznelson, An introduction to harmonic analysis, Second corrected edition. Dover Publications, Inc., New York, 1976. ISBN 0-486-63331-4
  • Felix Klein, Development of mathematics in the 19th century. Mathsci Press Brookline, Mass, 1979. Translated by M. Ackerman from Vorlesungen über die Entwicklung der Mathematik im 19 Jahrhundert, Springer, Berlin, 1928.
  • Walter Rudin, Principles of mathematical analysis, Third edition. McGraw-Hill, Inc., New York, 1976. ISBN 0-07-054235-X
  • William E. Boyce and Richard C. DiPrima, Elementary Differential Equations and Boundary Value Problems, Eighth edition. John Wiley & Sons, Inc., New Jersey, 2005. ISBN 0-471-43338-1
Portal icon Portal Matematică