Sari la conținut

Transformata Fourier

De la Wikipedia, enciclopedia liberă
Reprezentare timp–frecvență a unui acord de Do major interpretat la pian, obținută printr-o transformare cu factor Q constant. Frecvențele sunt reprezentate pe axa orizontală, pe o scară logaritmică marcată prin note muzicale, iar timpul pe axa verticală, crescând de jos în sus. Cele trei componente fundamentale din stânga corespund notelor do, mi și sol. Celelalte componente reprezintă armonicele superioare ale acestor sunete.

În matematică, transformata Fourier este o transformată integrală care asociază unei funcții o altă funcție, descriind componentele sale de frecvență. Atunci când variabila inițială reprezintă timpul, transformata oferă o reprezentare în domeniul frecvențelor, în locul reprezentării în domeniul timpului. Transformata are, în general, valori complexe, care codifică atât amplitudinea, cât și faza componentelor. Termenul desemnează atât operația matematică, numită și transformare Fourier, cât și rezultatul aplicării ei unei funcții.[notă 1]

Transformata Fourier pune în corespondență reprezentarea în timp, cu roșu, și reprezentarea în frecvență, cu albastru. Componentele sinusoidale ale semnalului corespund unor maxime ale reprezentării spectrale.

Localizarea unei funcții și a transformatei sale sunt legate prin principiul de incertitudine. O funcție concentrată într-un interval scurt are, de regulă, un spectru mai larg, iar concentrarea spectrului limitează localizarea funcției. Funcțiile de tip Gauss realizează cazul de egalitate în forma clasică a acestui principiu. Transformata unei funcții gaussiene este tot una gaussiană. Aceste funcții apar în teoria probabilităților, în statistică și în descrierea unor fenomene precum difuzia. Joseph Fourier a introdus transformările prin sinus și cosinus în studiul propagării căldurii. Acestea sunt legate de părțile imaginară și reală ale transformatei complexe moderne.

În forma de bază, transformata este exprimată printr-o integrală improprie. Pentru funcții absolut integrabile, definiția se formulează cu integrala Lebesgue. Alte aplicații necesită extinderi la funcții de pătrat integrabil, la măsuri sau la distribuții.[1] De exemplu, delta lui Dirac este folosită frecvent ca un impuls ideal, însă nu este o funcție obișnuită și trebuie tratată în cadrul măsurilor sau al distribuțiilor.

Transformata se definește și pentru funcții de mai multe variabile pe spații euclidiene. În fizică, ea permite trecerea de la o reprezentare în spațiul pozițiilor la una în spațiul impulsurilor sau de la variabilele spațiu–timp la cele ale cuadrimpulsului. Această legătură este importantă în studiul undelor și în mecanica cuantică. Funcțiile considerate pot avea valori complexe, vectoriale sau, în anumite aplicații, valori operatoriale.[notă 2]

Generalizarea la funcții definite pe grupuri reunește mai multe construcții înrudite. Dreapta reală și spațiile ℝn conduc la transformata integrală, grupul ℤ la transformata Fourier în timp discret, grupul ciclic ℤ/Nℤ la transformata Fourier discretă, iar grupul cercului la seriile Fourier. Cercul poate fi reprezentat printr-un interval ale cărui extremități sunt identificate, ceea ce îl face potrivit pentru funcții periodice. Transformata Fourier rapidă (abreviată FFT), este o familie de algoritmi pentru calcularea transformatei discrete.

Definiție și inversare

[modificare | modificare sursă]

Fie f:ℝ→ℂ o funcție măsurabilă și absolut integrabilă, adică ‖f‖1=∫ℝ|f(x)|dx<∞. În convenția cu frecvență obișnuită, transformata Fourier a lui f este funcția[2] f^(ξ)=∫ℝf(x)e−2πixξdx,ξ∈ℝ. Integrala converge absolut pentru orice ξ. Funcția f^ este mărginită și uniform continuă, iar f^(ξ)→0 când |ξ|→∞. Ultima proprietate constituie lema Riemann–Lebesgue.[3]

Notația L1(ℝ) desemnează spațiul acestor funcții, cu identificarea funcțiilor egale aproape peste tot. Transformata este injectivă pe acest spațiu, însă imaginea sa nu are o caracterizare elementară. În particular, transformata unei funcții din L1 nu este neapărat integrabilă. Dacă f și f^ sunt ambele integrabile, formula de inversare este:

f(x)=∫ℝf^(ξ)e2πixξdξ.

Egalitatea este valabilă aproape peste tot și în fiecare punct în care f este continuă. Individual, ea este valabilă peste tot pentru funcții din spațiul Schwartz. Alte forme ale teoremei de inversare utilizează ipoteze și sensuri de convergență diferite, care trebuie precizate separat.[4]

Funcțiile f și f^ formează o pereche de transformate Fourier.[5] Această corespondență se notează frecvent:[6] f(x) ⟷ℱ f^(ξ). De exemplu, în sensul distribuțiilor, δ(x)⟷ℱ1. Pentru asemenea obiecte, integrala inițială este înlocuită de o definiție prin dualitate. Transformata se extinde astfel la toate distribuțiile temperate.[7] Pentru unele funcții local integrabile care nu aparțin lui L1, ea poate fi definită și prin integrale improprii ori printr-o procedură de regularizare, atunci când limitele respective există în sensul indicat.[8]

Transformata integrală nu este definită direct prin aceeași integrală Lebesgue pentru orice funcție din L2(ℝ). De exemplu, funcția f(x)=11+x2 aparține spațiului L2(ℝ), dar nu și spațiului L1(ℝ). Transformata Fourier, definită inițial pe subspațiul dens L1(ℝ)∩L2(ℝ), admite o unică prelungire continuă la L2(ℝ). Această prelungire este un operator unitar, adică un automorfism liniar care păstrează produsul scalar și norma. Spre deosebire de transformata Fourier definită pe L1(ℝ), operatorul astfel obținut este bijectiv pe L2(ℝ). Transformata unei funcții din acest spațiu poate fi definită prin aproximarea cu funcții din L1(ℝ)∩L2(ℝ) și trecerea la limită în norma L2, fără a fi necesară convergența punctuală a integralei Fourier. În anumite condiții, transformata poate fi obținută și prin alte procedee de trecere la limită, inclusiv prin convergență slabă sau prin integrale în valoare principală.[9][10]

Gaussienele au un rol special în justificarea acestei extinderi. Combinații liniare finite de gaussiene, admițând translații și variații de scară, pot aproxima funcții din L1∩L2, iar proprietăți precum unitaritatea se verifică mai întâi pe aceste funcții. Două relații fundamentale sunt ℱ(e−πx2)(ξ)=e−πξ2 și ∫ℝe−πx2dx=1.[11]

Convenția cu exponentul −2πixξ are avantajul că transformarea este unitară pe L2 și transformă convoluția în produs, fără factori suplimentari. Mai precis, ea este un homomorfism al algebrei Banach L1, dotată cu convoluția, în algebra funcțiilor continue care se anulează la infinit, dotată cu produsul punctual și normă uniformă.[12]

Notații și convenții

[modificare | modificare sursă]

Transformata poate fi notată f^, f~, F, ℱ(f) sau ℱf. Valoarea sa în punctul ξ se scrie f^(ξ) ori (ℱf)(ξ). Parantezele fac explicit faptul că operatorul se aplică întâi funcției, după care rezultatul este evaluat în punctul indicat. Notațiile ℱ{f} și ℱ{f(t)} sunt de asemenea frecvent folosite, ultima precizând variabila în raport cu care se efectuează transformarea.

În științe și inginerie, variabila inițială este frecvent timpul t, semnalul este notat x(t), iar frecvența este notată f. Perechea se scrie atunci x(t)⟷ℱX(f). Această notație cu majuscule poate deveni incomodă pentru produse sau derivate, pentru care expresiile ℱ{fg} și ℱ{f′} sunt mai explicite. Uneori, în fizica particulelor elementare, funcția și transformata ei sunt notate cu același simbol și sunt deosebite numai prin argument. Astfel, un argument de tip impuls, ca în f(k1+k2), indică reprezentarea transformată, iar un argument de tip poziție, ca în f(x0+π𝒓), indică reprezentarea inițială. Tilda poate desemna și alte modificări ale unei mărimi, precum măsura Lorentz-invariantă d3k~=d3k/((2π)32ωk), unde ωk=|k|2+m2>0, în unități c=ℏ=1, pe ramura energiei pozitive.[13]

Expresiile ℱ(f(x)) sau ℱ(rect⁡(x)) sunt uzuale, dar trebuie înțelese ca transformate ale unor funcții, nu ale valorilor numerice ale acestora într-un punct fix. De exemplu, relațiile ℱ(rect⁡)(ξ)=sinc⁡(ξ),ℱ{f(x+x0)}(ξ)=e2πix0ξf^(ξ) presupun că integrarea se face în raport cu x, iar x0 este un parametru. Prima formulă utilizează funcția sinc normalizată, sinc⁡(ξ)=sin⁡(πξ)/(πξ), cu valoarea 1 la origine, și funcția dreptunghiulară asociată intervalului [−1/2,1/2].

Frecvență obișnuită și frecvență unghiulară

[modificare | modificare sursă]

Dacă variabila x reprezintă timpul, exprimat în secunde, variabila ξ reprezintă frecvența, măsurată în hertzi (Hz). Frecvența unghiulară, denumită și pulsație, este dată de relația ω=2πξ și se exprimă în radiani pe secundă (rad/s). Prin schimbarea de variabilă ξ=ω2π, se obțin două convenții curente de definire a transformatei Fourier. În prima convenție, transformata directă nu conține un factor de normalizare, iar transformata inversă conține factorul 12π. În a doua convenție, normalizarea este simetrică, ambele transformate având factorul 12π.[14]

Convenții frecvente în cazul unei variabile
Variabilă spectralăNormalizareTransformata directă, relațiile dintre convenții și inversa
Frecvența obișnuită ξUnitarăf^1(ξ)=∫ℝf(x)e−2πixξdx=2πf^2(2πξ)=f^3(2πξ)f(x)=∫ℝf^1(ξ)e2πixξdξ
Frecvența unghiulară ωUnitarăf^2(ω)=12π∫ℝf(x)e−iωxdx=12πf^1(ω2π)=12πf^3(ω)f(x)=12π∫ℝf^2(ω)eiωxdω
Neunitarăf^3(ω)=∫ℝf(x)e−iωxdx=f^1(ω2π)=2πf^2(ω)f(x)=12π∫ℝf^3(ω)eiωxdω

Caracterul unitar se referă la prelungirea pe L2, cu măsura Lebesgue obișnuită pentru fiecare variabilă. Trecerea la pulsație fără normalizare simetrică nu păstrează norma L2, deoarece schimbarea de variabilă modifică elementul de integrare.

În cazul funcțiilor definite pe ℝn, factorii de normalizare corespunzători celor două convenții sunt 1(2π)n/2 și 1(2π)n. În tabelul următor, x, ξ și ω reprezintă vectori din ℝn, iar punctul dintre doi vectori indică produsul scalar euclidian. Dacă variabilele inițiale reprezintă coordonate spațiale, variabilele corespunzătoare din domeniul spectral se exprimă în unități inverse de lungime, spre deosebire de frecvența temporală, exprimată în hertzi, sau de pulsație care este exprimată în radiani pe secundă.

Convenții frecvente în dimensiune n
Variabilă spectralăNormalizareTransformata directă, relațiile dintre convenții și inversa
Variabila ξUnitarăf^1(ξ)=∫ℝnf(x)e−2πiξ⋅xdx=(2π)n/2f^2(2πξ)=f^3(2πξ)f(x)=∫ℝnf^1(ξ)e2πiξ⋅xdξ
Variabila unghiulară ωUnitarăf^2(ω)=1(2π)n/2∫ℝnf(x)e−iω⋅xdx=1(2π)n/2f^1(ω2π)=1(2π)n/2f^3(ω)f(x)=1(2π)n/2∫ℝnf^2(ω)eiω⋅xdω
Neunitarăf^3(ω)=∫ℝnf(x)e−iω⋅xdx=f^1(ω2π)=(2π)n/2f^2(ω)f(x)=1(2π)n∫ℝnf^3(ω)eiω⋅xdω

În toate aceste convenții, semnul exponentului din transformata directă poate fi inversat dacă se inversează simultan semnul din transformata inversă. Semnele celor două nuclee exponențiale trebuie să fie opuse. Atunci când sunt comparate formule din surse diferite, trebuie verificate împreună semnul, variabila spectrală și factorii de normalizare.

Unități și variabile duale

[modificare | modificare sursă]

Variabila spectrală are unități inverse celor ale variabilei inițiale, astfel încât exponentul din nucleul transformării să fie adimensional. Un timp exprimat în secunde corespunde unei frecvențe exprimate în cicluri pe secundă sau unei pulsații în radiani pe secundă, în funcție de convenție. O coordonată de lungime corespunde unei variabile de tip număr de undă. Spațiul variabilei inițiale și spațiul variabilei spectrale au astfel semnificații și unități diferite.

În fizică este frecventă convenția cu semnul pozitiv în transformata directă și semnul negativ în inversă.[15] Aceasta schimbă alegerea undelor elementare ale reprezentării, fără a schimba fenomenul descris. De exemplu, coeficientul spectral poate fi asociat undei e−2πiξx, în locul undei e2πiξx. Exponențiale cu semnul negativ apar în dependența temporală a undelor electromagnetice plane și a funcțiilor de undă cuantice. Identitățile transformării trebuie rescrise consecvent, înlocuind semnele factorilor care conțin explicit unitatea imaginară. În electrotehnică, aceasta este notată de obicei j, deoarece litera i este rezervată intensității curentului electric.

În teoria probabilităților, funcția caracteristică a unei variabile aleatorii X cu densitatea f este definită de obicei prin φ(λ)=𝔼(eiλX)=∫ℝf(x)eiλxdx. Această convenție are semnul pozitiv și nu conține factorul 2π.

Interpretarea spectrală

[modificare | modificare sursă]

Amplitudine, fază și sinusoide complexe

[modificare | modificare sursă]
O sinusoidă, cu amplitudinea de vârf indicată prin 1, valoarea vârf la vârf prin 2, valoarea efectivă prin 3 și perioada prin 4.
Sinusoidele roșie și albastră au un defazaj θ.

Un coeficient spectral complex se poate exprima în formă polară sau algebrică, prin formula lui Euler, ca f^(ξ)=A(ξ)eiθ(ξ)=A(ξ)cos⁡θ(ξ)+iA(ξ)sin⁡θ(ξ). Aici A(ξ)=|f^(ξ)| reprezintă modulul, iar θ(ξ)=arg⁡f^(ξ) este faza, definită modulo 2π atunci când coeficientul nu este nul.[16]

În integrala inversă, contribuția corespunzătoare frecvenței ξ are expresia f^(ξ)e2πiξx=A(ξ)ei(2πξx+θ(ξ))=A(ξ)cos⁡(2πξx+θ(ξ))+iA(ξ)sin⁡(2πξx+θ(ξ)). Formula de inversare recompune funcția din aceste sinusoide complexe, cu amplitudinile și fazele lor,[17] f(x)=∫ℝA(ξ)ei(2πξx+θ(ξ))dξ.

Transformata directă poate fi interpretată ca măsurarea răspunsului funcției la o oscilație de probă. Înmulțirea cu e−2πiξx deplasează frecvența urmărită la zero. Integrarea însumează această contribuție, în timp ce oscilațiile la alte frecvențe tind să se compenseze. Interpretarea este euristică și trebuie adaptată cadrului funcțional utilizat.[notă 3]

Frecvențe negative

[modificare | modificare sursă]

Exponențiala complexă permite distingerea celor două sensuri de rotație în planul complex, corespunzătoare frecvențelor pozitive și negative. Transformata este definită pentru orice ξ∈ℝ. Unele semnale complexe au spectrul nul pentru ξ<0, cum este cazul semnalelor analitice. Pentru alte funcții complexe sunt necesare ambele semiaxe spectrale. Astfel de reprezentări apar în prelucrarea semnalelor, radar, optică neliniară, ecuații cu derivate parțiale și mecanică cuantică.

Pentru o funcție reală, transformata are operator hermitic, exprimat prin f^(−ξ)=f^(ξ)‾. Un cosinus real cos⁡(2πξ0x) conține componente la ambele frecvențe ±ξ0, spre deosebire de exponențiala complexă e2πiξ0x. Totuși, semnul frecvenței unui cosinus nu poate fi determinat din funcția reală, deoarece cos⁡(2πξ0x)=cos⁡(2π(−ξ0)x).[18]

Funcții periodice și eșantionarea spectrului

[modificare | modificare sursă]

O funcție periodică integrabilă pe o perioadă, care nu este nulă aproape peste tot, nu este absolut integrabilă pe întreaga dreaptă. Transformata sa nu se definește, în general, prin integrala Lebesgue folosită pentru L1(ℝ). Seria Fourier oferă însă o reprezentare adecvată, iar transformata integrală se poate extinde la funcția periodică în sensul distribuțiilor temperate.

Pentru o funcție P-periodică, cu P>0 și coeficienții Fourier cn=1P∫0Pf(x)e−2πinx/Pdx, corespondența distribuțională este f^(ξ)=∑n∈ℤcnδ(ξ−nP). Spectrul este astfel o sumă de impulsuri Dirac situate la multiplii frecvenței fundamentale, ponderate cu coeficienții seriei Fourier. Forma și sensul convergenței seriei funcției inițiale depind de regularitatea acesteia. Pentru o funcție de pătrat integrabil pe o perioadă, seria converge în norma L2, ceea ce nu trebuie confundat cu o afirmație de convergență punctuală sau uniformă.[19]

Periodizare și formula de sumare Poisson

[modificare | modificare sursă]

Eșantionarea transformatei la frecvențe echidistante este legată de periodizarea funcției inițiale. Pentru o funcție Schwartz f, periodizarea de perioadă P>0 este fP(x)=∑n∈ℤf(x+nP). În acest cadru, seriile și operațiile necesare sunt convergente, iar formula de sumare Poisson dă[20] fP(x)=1P∑k∈ℤf^(kP)e2πikx/P. Coeficientul Fourier de ordin k al funcției periodizate este deci P−1f^(k/P). Invers, fiecare eșantion spectral se determină prin integrarea pe o perioadă, f^(kP)=∫0PfP(x)e−2πikx/Pdx.

Pentru f∈L1(ℝ), seria care definește periodizarea converge absolut aproape peste tot, iar funcția rezultată este integrabilă pe orice interval de lungime egală cu perioada. Coeficienții Fourier ai acestei funcții se exprimă prin valorile f^(kP). Totuși, simpla integrabilitate a funcției f nu asigură convergența punctuală a seriei Fourier către funcția periodizată. Dacă f aparține spațiului Schwartz, condițiile suplimentare de regularitate și descreștere rapidă permit justificarea identității dintre seria Fourier și periodizarea funcției.

Exemplu de analiză a unui semnal gaussian modulat

[modificare | modificare sursă]

Funcția f(t)=cos⁡(2π⋅3t)e−πt2 descrie un cosinus de 3 Hz modulat printr-o anvelopă gaussiană. Anvelopa face ca oscilația să crească și să scadă treptat. Transformata este suma a două gaussiene deplasate, f^(ξ)=12(e−π(ξ−3)2+e−π(ξ+3)2). Pentru frecvența de 3Hz, partea reală a funcției de sub semnul integralei Fourier este e−πt2(1+cos⁡(12πt))2. Aceasta nu este negativă și are o integrală strict pozitivă. Partea imaginară este o funcție impară, astfel încât integrala sa pe întreaga axă reală este nulă. Acest rezultat se explică prin faptul că factorii cos⁡(6πt) au aceeași fază, iar componenta sinusoidală este defazată cu un sfert de perioadă față de componenta cosinusoidală.

Funcția inițială, cu o componentă dominantă la 3 Hz, și părțile reală și imaginară ale integrandului transformatei la frecvența de 3 Hz. Imaginea utilizează factorul e+2πi⋅3t, conjugat față de convenția principală a articolului. Partea imaginară are astfel semnul opus, dar integrala ei rămâne zero.

Valoarea exactă este f^(3)=12(1+e−36π), foarte apropiată de 0,5. Aceeași valoare apare la −3 Hz datorită simetriei funcției reale. Cele două componente au împreună contribuția dominantă asociată cosinusului de amplitudine 1. Întrucât anvelopa are un spectru gaussian, acesta nu este concentrat într-o singură frecvență, iar amplitudinea cosinusului nu trebuie identificată exact cu suma valorilor transformatei în numai două puncte.

La 5 Hz, oscilațiile integrandului produc o compensare mult mai puternică. Valoarea transformatei este f^(5)=12(e−4π+e−64π), pozitivă, dar mică. Exemplul arată cum transformata diferențiază o regiune spectrală dominantă de una în care contribuția este redusă, fără ca aceasta din urmă să fie identic nulă.

În schimb, semnalul complex e2πi⋅3te−πt2 are transformata e−π(ξ−3)2, o singură gaussiană centrată la 3 Hz. Valoarea ei maximă este 1, egală cu integrala anvelopei. Apariția celor două maxime pentru cosinusul real rezultă din descompunerea lui în două exponențiale complexe de frecvențe opuse.

Proprietăți

[modificare | modificare sursă]

În convenția cu frecvență obișnuită, nucleul transformării este e−2πixξ. Proprietățile elementare de mai jos sunt valabile pentru funcții din L1(ℝ), dacă nu se precizează alte ipoteze. Ele se extind și la alte spații de funcții sau de distribuții, cu interpretarea corespunzătoare a transformatei și a egalităților.[21]

Liniaritate, translație și schimbare de scară

[modificare | modificare sursă]

Transformarea Fourier este liniară. Pentru a,b∈ℂ și f,g∈L1(ℝ), are loc ℱ(af+bg)=af^+bg^. Translația argumentului funcției produce o schimbare a fazei spectrale, iar înmulțirea cu o exponențială complexă produce o translație în frecvență. Mai precis, ℱ{f(x−x0)}(ξ)=e−2πix0ξf^(ξ),x0∈ℝ,ℱ{e2πiξ0xf(x)}(ξ)=f^(ξ−ξ0),ξ0∈ℝ. Aceste relații explică folosirea transformatei în studiul semnalelor decalate în timp și al modulației.[21][22]

Pentru un factor real a≠0, schimbarea de scară este descrisă de ℱ{f(ax)}(ξ)=1|a|f^(ξa). Astfel, comprimarea funcției în variabila x dilată transformata în variabila ξ. Cazul a=−1 dă relația de reflexie ℱ{f(−x)}(ξ)=f^(−ξ).[21]

În stânga sunt reprezentate pulsurile, iar în dreapta părțile reală, violet, și imaginară, verde, ale transformatelor. Sus, pulsul este centrat și transformata sa este reală și pară. Jos, translația introduce o componentă imaginară impară.

Simetrie și conjugare

[modificare | modificare sursă]

Părțile reală și imaginară ale unei funcții complexe se pot descompune, fiecare, într-o parte pară și una impară. Dacă indicii RE, RO, IE și IO desemnează, în această ordine, partea reală pară, reală impară, imaginară pară și imaginară impară, corespondențele sunt Componentă în domeniul inițialComponentă a transformateifREf^REfROif^IOifIEif^IEifIOf^RO În acest tabel, funcțiile notate prin componente sunt reale. Factorul i este scris separat.[23]

În consecință, transformata unei funcții reale este hermitiană, adică f^(−ξ)=f^(ξ)‾. Reciproc, o transformată adjunctă corespunde unei funcții reale, în sensul egalității adoptat în spațiul respectiv. Transformata unei funcții pur imaginare este anti-adjunctă, satisfăcând f^(−ξ)=−f^(ξ)‾. Transformata unei funcții adjuncte este reală, iar a unei funcții anti-adjuncte este pur imaginară.[23]

Relația generală pentru conjugarea complexă este ℱ{f(x)‾}(ξ)=f^(−ξ)‾. Din aceasta rezultă formulele pentru părțile reală și imaginară, ℱ{Re⁡f}(ξ)=f^(ξ)+f^(−ξ)‾2,ℱ{Im⁡f}(ξ)=f^(ξ)−f^(−ξ)‾2i. Simetria adjunctă nu trebuie confundată cu paritatea obișnuită a unei funcții complexe. Ea pune în relație valorile în puncte opuse prin conjugare.[21]

Valoarea în origine și comportamentul la infinit

[modificare | modificare sursă]

Valoarea transformatei la frecvența zero este f^(0)=∫ℝf(x)dx. Pentru o funcție integrabilă, aceasta reprezintă integrala totală sau aria algebrică a semnalului. Interpretarea ei trebuie deosebită de media temporală a unui semnal periodic sau staționar.[21]

Funcția dreptunghiulară este integrabilă Lebesgue.
Transformata funcției dreptunghiulare este o funcție sinc, continuă și mărginită, care nu este integrabilă în modul.

Pentru orice f∈L1(ℝ), transformata este uniform continuă și satisface estimarea ‖f^‖∞≤‖f‖1. Lema Riemann–Lebesgue afirmă că f^(ξ)→0 când |ξ|→∞.[24][25]

Transformata nu este însă neapărat integrabilă. Un exemplu este perechea formată din funcția dreptunghiulară și funcția sinc normalizată. Integrala improprie a funcției sinc converge, dar integrala modulului său diverge. În consecință, pentru o funcție arbitrară din L1, inversarea nu se poate interpreta automat ca o integrală Lebesgue. Dacă atât f, cât și f^ aparțin lui L1(ℝ), atunci f(x)=∫ℝf^(ξ)e2πixξdξ aproape peste tot. Într-un punct de continuitate al reprezentantului ales pentru f, egalitatea este punctuală. Transformarea Fourier este injectivă pe L1, considerat ca spațiu de clase de funcții egale aproape peste tot.[26]

Relația lui Parseval și teorema lui Plancherel

[modificare | modificare sursă]

Pentru f,g∈L1(ℝ)∩L2(ℝ), relația lui Parseval are forma ⟨f,g⟩L2=∫ℝf(x)g(x)‾dx=∫ℝf^(ξ)g^(ξ)‾dξ. În special, ‖f‖22=∫ℝ|f(x)|2dx=∫ℝ|f^(ξ)|2dξ. Această egalitate a normelor permite prelungirea unică prin continuitate a transformării Fourier la un operator unitar pe întregul spațiu L2(ℝ). Pe intersecția L1∩L2, prelungirea coincide cu transformarea definită prin integrală.[27] Definițiile sunt astfel compatibile pe suma de spații L1(ℝ)+L2(ℝ), care conține și Lp(ℝ) pentru 1≤p≤2. În aplicații, egalitatea normelor pătratice exprimă conservarea energiei semnalului.[27][28]

Convoluție și corelație

[modificare | modificare sursă]

Transformarea Fourier transformă convoluția în produs. Dacă f,g∈L1(ℝ) și h(x)=(f∗g)(x)=∫ℝf(y)g(x−y)dy, atunci h^(ξ)=f^(ξ)g^(ξ). Sub alte convenții de normalizare poate apărea un factor constant.[26][29]

În teoria sistemelor liniare invariante în timp, g poate reprezenta răspunsul la impuls al unui sistem, f semnalul de intrare, iar h semnalul de ieșire. Relația δ∗g=g explică denumirea răspunsului la impuls, iar g^ este răspunsul sistemului în frecvență. Produsul spectral permite calcularea ieșirii fără evaluarea directă a integralei de convoluție.[30]

În sens dual, dacă p,q∈L2(ℝ), produsul pq este integrabil și ℱ(pq)(ξ)=(p^∗q^)(ξ). Integrala din membrul drept este definită prin inegalitatea Cauchy–Schwarz. Ipoteza că ambii factori ar fi numai integrabili nu ar garanta, în general, integrabilitatea produsului.[26]

Pentru convenția de corelație încrucișată (f⋆g)(x)=∫ℝf(y)‾g(x+y)dy, relația corespunzătoare este ℱ(f⋆g)(ξ)=f^(ξ)‾g^(ξ). Ea este valabilă, de exemplu, pentru f,g∈L1∩L2. Pentru date doar din L2, corelația este bine definită, iar identitatea se interpretează prin inversarea produsului spectral integrabil sau în sensul distribuțiilor. Cazul particular g=f dă autocorelația și relația ℱ(f⋆f)(ξ)=|f^(ξ)|2. Această identitate pentru semnale de energie este înrudită cu teorema Wiener–Hincin pentru spectrele proceselor staționare, unde ipotezele și normalizarea corelației sunt diferite.[21][30]

Derivare și regularitate

[modificare | modificare sursă]

Dacă f este local absolut continuă și f,f′∈L1(ℝ), integrarea prin părți conduce la f′^(ξ)=2πiξf^(ξ). Continuitatea absolută este esențială pentru această formulare cu derivată clasică aproape peste tot. Simpla existență aproape peste tot a derivatei și integrabilitatea ei nu sunt suficiente. O formulare echivalentă folosește derivata generalizată și condiția f∈W1,1(ℝ).[31][32]

Mai general, pentru f∈Wn,1(ℝ), f(n)^(ξ)=(2πiξ)nf^(ξ). În cealaltă direcție, dacă f și xnf(x) sunt integrabile, transformata admite derivatele corespunzătoare și ℱ{xnf(x)}(ξ)=(i2π)ndnf^dξn(ξ). În particular, pe spațiul Schwartz, toate aceste operații sunt definite și se poate scrie și ℱ−1{dnf^dξn}(x)=(−2πix)nf(x).[33][34]

Derivarea unei funcții devine astfel înmulțire cu un polinom în variabila spectrală, proprietate care transformă anumite ecuații diferențiale în ecuații algebrice. Formulele exprimă și legătura dintre regularitate și descreșterea transformatei. Derivatele integrabile produc factori suplimentari de descreștere în frecvență, iar existența momentelor integrabile produce regularitatea transformatei. Pe spațiul Schwartz, relația se exprimă exact prin invarianța acestui spațiu la transformarea Fourier.[33]

Funcții proprii și funcții Hermite

[modificare | modificare sursă]

O funcție proprie nenulă a transformării Fourier satisface ℱψ=λψ, pentru un scalar λ∈ℂ. Funcțiile proprii pot fi obținute din operatori diferențiali compatibili cu transformarea. Dacă o ecuație omogenă definește un spațiu de soluții invariant la ℱ și acest spațiu este unidimensional, soluțiile sale nenule sunt funcții proprii. Invarianța singură nu implică această concluzie pentru fiecare soluție a unei ecuații cu un spațiu de soluții de dimensiune mai mare.[35][36]

Un procedeu formal pornește de la ecuația [U(12πddx)+U(x)]ψ(x)=0. Pentru un polinom U, notația operatorială se obține înlocuind variabila cu operatorul de derivare. Dacă toate monomiile sale au grade congruente modulo 4, transformarea ecuației introduce un singur factor comun de forma ik, care poate fi simplificat. Spațiul soluțiilor potrivite este astfel invariant. Pentru o funcție U dată printr-o serie de puteri, aceeași scriere cere în plus justificarea convergenței seriei de operatori și precizarea domeniului. Alegerea U(x)=x conduce la ψ′(x)+2πxψ(x)=0, ale cărei soluții sunt multipli ai gaussianei e−πx2.[35]

O altă construcție folosește operatorul A=W(i2πddx)+W(x). Pentru un polinom par W, acesta comută cu transformarea Fourier pe spațiul Schwartz. Ecuația Aψ=Cψ permite alegerea unor funcții proprii comune atunci când realizarea operatorului și proprietățile spațiilor sale proprii o permit. Cazul W(x)=x2 este cel al oscilatorului armonic cuantic și furnizează baza funcțiilor Hermite.[36]

În convenția transformatei Fourier cu nucleul e−2πixξ, funcțiile Hermite pot fi exprimate cu ajutorul polinoamelor Hermite definite prin convenția utilizată în teoria probabilităților. Aceste funcții au forma:

ψn(x)=21/4n!e−πx2Hen(2πx),n=0,1,2,…,

unde polinoamele Hen sunt definite prin relația:

Hen(x)=(−1)nex2/2(ddx)ne−x2/2

Funcțiile ψn formează o bază ortonormată completă a spațiului L2(ℝ) și sunt funcții proprii ale transformatei Fourier, cu valorile proprii (−i)n. Mai precis, ele satisfac relația ψ^n(ξ)=(−i)nψn(ξ). Dacă în definiția transformatei Fourier se utilizează nucleul cu exponent pozitiv, valorile proprii corespunzătoare sunt in.[21][37][38]

Transformarea are numai patru valori proprii, 1,−1,i,−i, iar combinațiile liniare ale funcțiilor proprii cu aceeași valoare proprie rămân funcții proprii. Spațiul L2(ℝ) este suma directă ortogonală a celor patru subspații proprii. De exemplu, subspațiul corespunzător valorii ik este generat, după închidere, de funcțiile ψn pentru care (−i)n=ik.[38]

Dezvoltarea unei funcții în această bază diagonalizează operatorul Fourier, ℱf=∑n=0∞(−i)n⟨f,ψn⟩ψn. Seria converge în norma L2. Scrierea echivalentă a nucleului, e−2πixξ=∑n=0∞(−i)nψn(x)ψn(ξ), se interpretează în sensul distribuțiilor temperate, fără a presupune convergență absolută punctuală. Această reprezentare spectrală, legată de o construcție propusă de Wiener, explică și extinderea la transformata Fourier fracționară.[39][40]

Operatorul număr al oscilatorului armonic are forma:

N=14π(2πx−ddx)(2πx+ddx)=−14πd2dx2+πx2−12.

Operatorul N, considerat în realizarea sa autoadjunctă, are funcțiile Hermite ψn drept funcții proprii, corespunzătoare valorilor proprii n, potrivit relației Nψn=nψn. Prin calculul spectral se obține familia de operatori unitari e−itN, definită pentru orice t∈ℝ. Transformata Fourier, în convenția cu nucleul e−2πixξ, este un caz particular al acestei familii, corespunzător valorii t=π2, astfel încât ℱ=e−iπN/2. Funcțiile Hermite rămân funcții proprii pentru fiecare operator al familiei, având valorile proprii e−itn. Reprezentarea integrală a acestor operatori este legată de formula lui Mehler, care exprimă nucleul integral prin intermediul funcțiilor Hermite.[41][42]

Funcțiile Hermite au descreștere gaussiană atât în variabila inițială, cât și în cea spectrală. Această localizare explică rolul lor în analiza timp–frecvență și în transformata Fourier fracționară. În fizică, o construcție a acestei familii de transformări a fost dată de Condon. Prin extinderea la distribuții, și pieptenele Dirac cu pasul 1 este invariant la transformarea Fourier.[43][44]

Inversare și periodicitatea operatorului

[modificare | modificare sursă]

Pe spațiul Schwartz, pe L2 și pe spațiul distribuțiilor temperate, aplicarea de două ori a transformării Fourier produce reflexia argumentului. Dacă (𝒫f)(x)=f(−x), atunci:

ℱ0=I,ℱ1=ℱ,ℱ2=𝒫,ℱ3=ℱ−1=𝒫ℱ=ℱ𝒫,ℱ4=I.

În L2, egalitățile se referă la clase de funcții egale aproape peste tot. În 𝒮′, sunt egalități de distribuții. Pentru o funcție arbitrară din L1, iterațiile nu trebuie confundate cu aplicări succesive ale integralei Lebesgue, deoarece prima transformată poate să nu mai aparțină lui L1.[45][38]

Periodicitatea de ordin patru poate fi pusă în legătură cu o rotație de 90∘ în planul timp–frecvență. Iterarea de două ori corespunde rotației de 180∘ și reflexiei ambelor variabile. Transformata fracționară permite unghiuri intermediare. Transformările canonice liniare generalizează această interpretare la acțiunea grupului SL2(ℝ) asupra planului fazelor, cu păstrarea formei simplectice. Realizarea lor prin operatori unitari implică alegerea convențiilor de fază.[46][47]

Legătura cu grupul Heisenberg

[modificare | modificare sursă]

Pe L2(ℝ), translațiile și modulațiile sunt operatori unitari, (Tyf)(x)=f(x+y),(Mξf)(x)=e2πiξxf(x). Ei nu comută cu aceste semne, TyMξ=e2πiξyMξTy, astfel încât Mξ−1Ty−1MξTy=e−2πiξyI. Comutatorul este o înmulțire cu un scalar de modul 1, independent de variabila funcției.[11][48]

O realizare coerentă a grupului generat de acești operatori folosește triple (x,ξ,z)∈ℝ2×U(1), cu reprezentarea ρ(x,ξ,z)=zMξTx. Legea de compoziție este (x1,ξ1,z1)(x2,ξ2,z2)=(x1+x2,ξ1+ξ2,z1z2e2πix1ξ2). Aceasta este o versiune a grupului Heisenberg cu centru circular. Formulele legii de grup depind de ordinea aleasă pentru translație și modulație, precum și de semnul din definiția translației.[11][48]

Aplicația liniară J(xξ)=(−ξx),J2=−I, admite, în aceste coordonate, ridicarea la automorfismul j(x,ξ,z)=(−ξ,x,ze−2πixξ). Transformarea Fourier leagă cele două reprezentări prin ρ∘j=ℱρℱ∗. Relația exprimă schimbarea dintre translațiile în variabila inițială și modulațiile în variabila duală.[11][47]

Teorema Stone–von Neumann afirmă echivalența unitară a reprezentărilor unitare ireductibile, continue în sens tare, ale grupului Heisenberg care au aceeași acțiune netrivială a centrului. Operatorul unitar care realizează echivalența este unic până la un factor scalar de modul 1, conform lemei lui Schur. În cazul automorfismului j, alegerea fazei corespunzătoare normalizării Fourier dă operatorul ℱ. În acest cadru, ℱ2=𝒫 reflectă relația J2=−I.[49][11]

Prelungirea în domeniul complex și transformata Laplace

[modificare | modificare sursă]

Integrala f^(ζ)=∫ℝe−2πitζf(t)dt poate fi studiată și pentru ζ=σ+iτ∈ℂ. Domeniul în care converge depinde de descreșterea lui f. Uneori nu admite o prelungire dincolo de axa reală, alteori definește o funcție olomorfă într-o bandă sau chiar o funcție întreagă.[50]

Suport compact și cauzalitate

[modificare | modificare sursă]

O formă a teoremei Paley–Wiener caracterizează funcțiile netede cu suport compact. O funcție f∈Cc∞(ℝ) are suportul în [−a,a] dacă și numai dacă transformata ei se prelungește la o funcție întreagă pentru care, pentru fiecare întreg n≥0, există o constantă Cn astfel încât |ζnf^(ζ)|≤Cne2πa|Im⁡ζ|. Transformata descrește rapid pe fiecare dreaptă orizontală, iar creșterea în direcția imaginară este controlată de suportul funcției inițiale. Există și versiuni pentru funcții de pătrat integrabil, în care condiția de descreștere rapidă este înlocuită de o condiție potrivită de integrabilitate pătratică. Spațiile de funcții întregi astfel obținute sunt numite spații Paley–Wiener. Teoreme de acest tip au fost extinse și la grupuri Lie reductive.[51][52][53]

O funcție care se anulează pentru t<0 este numită cauzală. Pentru răspunsul la impuls al unui sistem cauzal, această condiție exprimă absența unui răspuns anterior impulsului. Dacă f∈L1(ℝ) este cauzală, aceeași integrală definește o funcție holomorfă în semiplanul inferior τ<0, iar valorile ei tind la zero când τ→−∞. Rezultatele reciproce cer ipoteze suplimentare de creștere sau apartenența la spații Hardy. Simpla holomorfie într-un semiplan nu caracterizează, singură, toate transformatele funcțiilor cauzale.[50][54][55]

Relația cu transformata Laplace

[modificare | modificare sursă]

Pentru o funcție cauzală local integrabilă care satisface |f(t)|≤Ceat pentru t≥0, integrala f^(iτ)=∫0∞e2πτtf(t)dt, converge atunci când 2πτ<−a. Dacă F(s)=∫0∞f(t)e−stdt, este transformata Laplace unilaterală, atunci f^(ζ)=F(2πiζ) în domeniul comun de convergență și, în particular, f^(iτ)=F(−2πτ).[56]

Factorul exponențial suplimentar permite transformatei Laplace să trateze funcții pentru care integrala Fourier nu converge pe axa reală. În acest sens, o transformată Laplace este o transformată Fourier a unei funcții ponderate exponențial, evaluată pe o dreaptă adecvată din planul complex. Această legătură este folosită în rezolvarea ecuațiilor diferențiale și în analiza filtrelor și a sistemelor liniare, inclusiv în situațiile în care comportarea temporală nu permite o descriere Fourier clasică.[56]

Deplasarea conturului și inversarea

[modificare | modificare sursă]

Dacă f^ este olomorfă într-o bandă orizontală care conține dreptele Im⁡ζ=a și Im⁡ζ=b, iar descreșterea sa permite anularea contribuțiilor laturilor verticale ale conturului, teorema integrală a lui Cauchy dă ∫ℝf^(σ+ia)e2πi(σ+ia)tdσ=∫ℝf^(σ+ib)e2πi(σ+ib)tdσ. Formula de inversare poate fi astfel exprimată prin integrare pe drepte paralele cu axa reală. Descreșterea în bandă și sensul convergenței integralei fac parte din ipotezele rezultatului.[57]

Pentru o funcție cauzală f, a cărei creștere este de ordin exponențial, cel mult eat, se alege un număr real τ<−a2π. În condițiile care asigură existența și inversarea transformatei, funcția poate fi reconstituită prin integrarea transformatei Fourier complexe de-a lungul dreptei orizontale din planul complex, definită de partea imaginară τ

f(t)=∫ℝf^(σ+iτ)e2πi(σ+iτ)tdσ

Prin schimbarea de variabilă s=2πi(σ+iτ), această relație conduce la formula de inversare a transformatei Laplace, cunoscută sub denumirea de formula Bromwich–Mellin

f(t)=12πi∫b−i∞b+i∞F(s)estds,

unde F(s) reprezintă transformata Laplace a funcției f, iar b=−2πτ>a. Integrarea se efectuează de-a lungul dreptei verticale Re⁡(s)=b, orientată de jos în sus, situată în semiplanul de convergență al transformatei Laplace.[56]

O condiție suficientă mai generală cere ca e−btf(t) să fie integrabilă și ca f să aibă variație mărginită într-o vecinătate închisă a punctului analizat. Integrala de inversare este atunci luată ca limită simetrică. La un salt, ea restituie media limitelor laterale, iar la un punct de continuitate restituie valoarea funcției. Variante în norma L2 se obțin prin aplicarea teoremei lui Plancherel funcției ponderate corespunzătoare.[58]

Transformata Fourier în spațiul euclidian

[modificare | modificare sursă]

Pentru o funcție integrabilă f:ℝn→ℂ, transformata este f^(𝝃)=∫ℝnf(𝐱)e−2πi𝝃⋅𝐱d𝐱. Aici 𝝃⋅𝐱 este produsul scalar. Într-o formulare independentă de coordonate, variabila spectrală aparține spațiului vectorial dual, iar exponentul folosește împerecherea dintre vectori și forme liniare. Proprietățile elementare, lema Riemann–Lebesgue și teoremele Parseval și Plancherel au versiuni analoage în orice dimensiune finită.[26]

Principiul incertitudinii

[modificare | modificare sursă]

Localizarea unei funcții și cea a transformatei sale Fourier nu pot fi oricât de precise simultan. Această proprietate este sugerată de comportamentul transformatei Fourier la schimbarea scării. Comprimarea unei funcții în variabila inițială determină o dilatare corespunzătoare în domeniul frecvențelor. Principiul incertitudinii exprimă cantitativ această relație prin momentele de ordinul al doilea ale densităților asociate funcției și transformatei sale.[59]

În cazul unidimensional, pentru o funcție f∈L2(ℝ) normalizată prin ‖f‖2=1, teorema lui Plancherel asigură egalitatea ‖f^‖2=1. Funcțiile pozitive |f|2 și |f^|2 au integrala egală cu 1 și pot fi interpretate ca densități de probabilitate. Momentul de ordinul al doilea față de origine este definit prin

D0(f)=∫ℝx2|f(x)|2dx.

Dacă f este local absolut continuă, xf∈L2(ℝ) și f′∈L2(ℝ), principiul incertitudinii se exprimă prin inegalitatea

D0(f)D0(f^)≥116π2.

Această relație arată că momentele de ordinul al doilea ale celor două densități nu pot fi simultan arbitrar de mici.[59][60]

Egalitatea este atinsă de funcțiile gaussiene centrate, fără modulație, de forma

f(x)=Ce−πx2σ2,f^(ξ)=σCe−πσ2ξ2,

unde σ>0 și |C|=21/4σ. Argumentul complex al constantei C poate fi ales arbitrar. Pentru aceste funcții, momentele de ordinul al doilea sunt

D0(f)=σ24π,D0(f^)=14πσ2.

Deoarece densitățile |f|2 și |f^|2 sunt centrate în origine, aceste momente coincid cu varianțele lor. Trebuie făcută distincția între funcțiile gaussiene care reprezintă amplitudinile și densitățile obținute prin ridicarea la pătrat a modulelor acestora. Dacă funcțiile e−πx2σ2 și e−πσ2ξ2 sunt normalizate separat astfel încât integralele lor să fie egale cu 1, varianțele corespunzătoare sunt σ22π și 12πσ2. Această diferență de normalizare explică apariția unor constante diferite în formulele pentru varianțe.[59][61]

Prin translație și modulație se obține forma centrată a principiului incertitudinii, în care momentele sunt calculate față de punctele arbitrare x0 și ξ0

(∫ℝ(x−x0)2|f(x)|2dx)(∫ℝ(ξ−ξ0)2|f^(ξ)|2dξ)≥116π2.

Egalitatea este atinsă de funcțiile de forma f(x)=Ce2πiξ0xe−s(x−x0)2, unde s>0, iar constanta C este aleasă astfel încât ‖f‖2=1. În mecanica cuantică, relația dintre impuls și frecvența spațială permite exprimarea acestei inegalități cu ajutorul constantei lui Planck, obținându-se relația de incertitudine dintre poziție și impuls.[61][62]

Principiul incertitudinii admite și o formulare în termenii entropiei diferențiale. Pentru o funcție f normalizată prin ‖f‖2=1, inegalitatea entropică este

H(|f|2)+H(|f^|2)≥log⁡(e2), unde entropia diferențială a unei densități de probabilitate p este definită prin
H(p)=−∫ℝp(x)log⁡p(x)dx.

Relația presupune existența entropiilor respective, de exemplu ca numere reale finite. Logaritmii trebuie calculați în aceeași bază, mai mare decât 1. Inegalitatea entropică a fost studiată de Hirschman, iar constanta optimă a fost stabilită prin rezultatele lui Białynicki-Birula și Mycielski, respectiv Beckner. Această formulare implică inegalitatea bazată pe varianțe, folosind proprietatea conform căreia distribuția gaussiană are entropie maximă dintre distribuțiile cu aceeași varianță finită. Egalitatea este atinsă și în cazul entropic de funcțiile gaussiene corespunzătoare.[63]

Transformatele sinus și cosinus

[modificare | modificare sursă]

Formularea inițială a lui Fourier folosea funcții sinus și cosinus, fără exponențiale complexe. Pentru o funcție integrabilă pentru care inversarea este valabilă, descompunerea poate fi scrisă numai cu frecvențe pozitive, f(t)=∫0∞[a(λ)cos⁡(2πλt)+b(λ)sin⁡(2πλt)]dλ, cu a(λ)=2∫ℝf(t)cos⁡(2πλt)dt,b(λ)=2∫ℝf(t)sin⁡(2πλt)dt. Aceasta este o dezvoltare în integrală trigonometrică. Normalizarea trebuie precizată, deoarece denumirile de transformată sinus și transformată cosinus sunt folosite și pentru integrale pe semiaxă cu alți factori. Evitarea frecvențelor negative poate fi convenabilă în unele interpretări fizice și statistice.[64][56]

Identitățile trigonometrice conduc și la formula integrală Fourier f(t)=2∫0∞[∫ℝf(τ)cos⁡(2πλ(τ−t))dτ]dλ. Integrala exterioară trebuie interpretată în sensul de convergență corespunzător teoremei de inversare utilizate. Integrabilitatea lui f singură nu justifică schimbarea arbitrară a ordinii celor două integrale.[65][66][67]

Armonice sferice și funcții radiale

[modificare | modificare sursă]

În ℝn, polinoamele omogene armonice de grad k formează un spațiu notat 𝒜k. Ele sunt numite și armonice sferice solide. Dacă P∈𝒜k și f(𝐱)=e−π|𝐱|2P(𝐱), atunci f^(𝝃)=i−kf(𝝃). Ele oferă astfel funcții proprii în mai multe dimensiuni, într-un rol analog celui al funcțiilor Hermite.[26]

Fie ℋk închiderea în L2(ℝn) a combinațiilor liniare de funcții f0(|𝐱|)P(𝐱), cu P∈𝒜k, care aparțin acestui spațiu. Atunci L2(ℝn) este suma directă ortogonală a spațiilor ℋk, iar transformarea Fourier lasă fiecare dintre ele invariant. Pentru o funcție integrabilă de forma f(𝐱)=f0(|𝐱|)P(𝐱), are loc

f^(𝝃)=F0(|𝝃|)P(𝝃), unde, pentru r>0,

F0(r)=2πi−kr−(n+2k−2)/2∫0∞f0(s)J(n+2k−2)/2(2πrs)s(n+2k)/2ds. Funcția Jν este funcția Bessel de speța întâi și ordinul ν. Cazul k=0 dă transformata unei funcții radiale și o exprimă printr-o transformată Hankel.[26][68]

Există și o recurență între dimensiunile n și n+2. Dacă același profil f0 definește funcții radiale integrabile în ambele dimensiuni, iar ℱn(f0)(r) desemnează profilul radial al transformatei în dimensiunea n, atunci ℱn+2(f0)(r)=−12πrddrℱn(f0)(r),r>0. Aceasta permite, de exemplu, calcularea transformatei tridimensionale din cea unidimensională. Rezultatul are și versiuni pentru distribuții radiale temperate.[69]

Restricția transformatei și trunchierea spectrală

[modificare | modificare sursă]

Transformata unei funcții integrabile este continuă și poate fi restrânsă la orice submulțime a spațiului frecvențelor. Pentru o funcție arbitrară din L2(ℝn), transformata este însă o clasă de echivalență aproape peste tot, astfel încât valorile sale pe o mulțime de măsură Lebesgue nulă nu sunt determinate de acea clasă. Problema restricției studiază condițiile în care operatorul definit inițial pentru funcții suficient de regulate admite o prelungire continuă la alte clase Lp. Curbura mulțimii pe care se face restricția are un rol esențial.[39]

Pentru sfera unitate Sn−1⊂ℝn, cu n≥2, teorema Tomas–Stein afirmă că operatorul Rf=f^|Sn−1 se prelungește la un operator mărginit

R:Lp(ℝn)⟶L2(Sn−1,dσ),1≤p≤2(n+1)n+3, unde dσ este măsura de suprafață.

Afirmația nu presupune alegerea unor valori punctuale arbitrare ale unei clase Lp pe sferă.[39]

O altă diferență între dimensiunea 1 și dimensiunile superioare apare la inversarea prin trunchierea domeniului spectral. Pentru o familie crescătoare de mulțimi mărginite ER, se poate defini, inițial pentru f∈L1∩Lp, fR(𝐱)=∫ERf^(𝝃)e2πi𝐱⋅𝝃d𝝃. În dimensiunea 1, pentru ER=(−R,R) și 1<p<∞, funcțiile fR converg la f în norma Lp. Rezultatul este legat de mărginirea transformării Hilbert pe aceste spații. În mai multe dimensiuni, convergența în Lp se păstrează pentru trunchieri prin cuburi.[39]

Pentru bile euclidiene, situația este diferită. Teorema lui Fefferman afirmă că, în dimensiune n≥2, multiplicatorul Fourier definit de funcția indicatoare a bilei unitate este mărginit pe Lp(ℝn) numai în cazul p=2. Pentru p≠2, aceasta împiedică o teoremă generală de convergență în norma Lp pentru toate datele. Concluzia privește mărginirea operatorilor și convergența familiei de trunchieri, fără a afirma că fiecare trunchiere a unei funcții integrabile ar ieși din Lp.[70][71]

Transformata Fourier pe spații de funcții și de distribuții

[modificare | modificare sursă]

Spațiile L¹ și L²

[modificare | modificare sursă]

În orice dimensiune finită, definiția integrală dă un operator mărginit ℱ:L1(ℝn)→L∞(ℝn), deoarece sup𝝃∈ℝn|f^(𝝃)|≤∫ℝn|f(𝐱)|d𝐱. Lema Riemann–Lebesgue precizează că imaginea sa se află în C0(ℝn), spațiul funcțiilor continue care tind la zero la infinit. Imaginea lui L1 este totuși o submulțime proprie a acestui spațiu.[25][72][73]

Pentru f∈L2(ℝn), integrala Fourier nu este neapărat o integrală Lebesgue convergentă. Transformata se definește prin prelungirea unitară a operatorului de pe subspațiul dens L1∩L2. O realizare explicită este f^=limR→∞∫|𝐱|≤Rf(𝐱)e−2πi𝝃⋅𝐱d𝐱, unde limita este luată în norma L2, ca limită de funcții ale variabilei 𝝃. Mai general, orice șir din L1∩L2 care converge la f în L2 produce același rezultat prin trecerea la limită a transformatelor.[28][74]

Unitaritatea înseamnă bijectivitate și conservarea produsului scalar. Identitatea biliniară ∫ℝnf(𝐱)ℱg(𝐱)d𝐱=∫ℝnℱf(𝐱)g(𝐱)d𝐱 este valabilă pentru f,g∈L2. Împreună cu inversarea și relația ℱf‾=ℱ−1f‾, ea conduce la ∫ℝnℱf‾ℱg=∫ℝn(ℱ−1f‾)ℱg=∫ℝnℱ(ℱ−1f‾)g=∫ℝnf‾g. În această scriere, toate integralele sunt în raport cu măsura Lebesgue. Formula exprimă conservarea produsului scalar, distinctă de identitatea biliniară fără conjugare.[28]

Alte spații Lᵖ

[modificare | modificare sursă]

Pentru 1<p<2, o funcție din Lp poate fi descompusă într-o parte din L2, unde |f|≤1, și o parte din L1, unde |f|>1. Compatibilitatea definițiilor permite astfel definirea transformatei. Teorema de interpolare Riesz–Thorin și inegalitatea Hausdorff–Young dau estimarea ‖f^‖q≤‖f‖p,1p+1q=1,1≤p≤2. Pentru p=1, se ia q=∞. În afara cazului p=2, imaginea nu admite o caracterizare la fel de simplă ca întregul spațiu țintă.[75][26]

Funcțiile din Lp, pentru 2<p<∞, definesc distribuții temperate și au transformate în acest sens generalizat. Transformata lor nu este întotdeauna reprezentabilă printr-o funcție. Extinderea la distribuții nu constituie, așadar, numai o convenție de notație pentru o integrală existentă.[75][26]

Spațiul Schwartz și distribuțiile temperate

[modificare | modificare sursă]

O distribuție pe ℝn este o funcțională liniară continuă definită pe spațiul Cc∞(ℝn) al funcțiilor indefinit derivabile cu suport compact, înzestrat cu topologia corespunzătoare. Deși acest spațiu este dens în L2(ℝn), el nu este păstrat de transformata Fourier. Transformata unei funcții nenule din Cc∞(ℝn) admite o prelungire la o funcție întreagă pe ℂn și nu poate avea suport compact. Singura funcție pentru care ambele proprietăți sunt îndeplinite este funcția nulă.[76][77]

Pentru extinderea transformatei Fourier la distribuții se utilizează spațiul Schwartz 𝒮(ℝn), alcătuit din funcțiile indefinit derivabile care, împreună cu toate derivatele lor, scad la infinit mai repede decât orice putere inversă a normei variabilei. Mai precis, o funcție f aparține acestui spațiu dacă, pentru orice pereche de multiindici α și β, este finită seminorma

pα,β(f)=sup𝐱∈ℝn|𝐱α∂βf(𝐱)|.

Spațiul funcțiilor netede cu suport compact este inclus în spațiul Schwartz, Cc∞(ℝn)⊂𝒮(ℝn). Transformata Fourier este un automorfism liniar continuu al spațiului 𝒮(ℝn), iar transformata inversă are aceeași proprietate.[78] Dualul topologic al spațiului Schwartz, notat cu 𝒮′(ℝn), constituie spațiul distribuțiilor temperate. Acesta cuprinde distribuțiile cu suport compact și distribuțiile asociate funcțiilor local integrabile cu creștere cel mult polinomială. De asemenea, funcțiile din spațiile Lp(ℝn), pentru 1≤p≤∞, definesc distribuții temperate. Pentru o funcție integrabilă f, distribuția asociată Tf acționează asupra unei funcții-test φ∈𝒮(ℝn) prin relația

⟨Tf,φ⟩=∫ℝnf(𝐱)φ(𝐱)d𝐱.

Această construcție permite considerarea funcțiilor integrabile ca elemente ale spațiului distribuțiilor temperate.[26][45]

Definirea transformatei Fourier a unei distribuții temperate se bazează pe identitatea

∫ℝnf^(𝐱)g(𝐱)d𝐱=∫ℝnf(𝐱)g^(𝐱)d𝐱, valabilă pentru funcții integrabile f și g.

Pornind de la această relație, transformata Fourier a unei distribuții T∈𝒮′(ℝn) este definită prin dualitate, potrivit formulei

⟨ℱT,φ⟩=⟨T,ℱφ⟩, pentru orice φ∈𝒮(ℝn).

Dualitatea considerată este biliniară, fără conjugarea complexă a funcției-test. Pentru distribuțiile asociate funcțiilor integrabile, această definiție coincide cu transformata Fourier clasică. Transformarea astfel definită este un automorfism liniar continuu al spațiului 𝒮′(ℝn), având drept inversă transformata Fourier inversă extinsă prin dualitate.[26][79]

Transformata Fourier a distribuțiilor păstrează relațiile fundamentale dintre derivare și înmulțirea cu variabila spectrală. Pentru derivata distribuțională în raport cu coordonata xj, se obține

ℱ(∂jT)=2πiξjℱT.

Proprietățile referitoare la convoluție se extind, de asemenea, la distribuții, cu condiția ca operațiile respective să fie bine definite. Pentru T∈𝒮′(ℝn) și φ∈𝒮(ℝn), sunt valabile relațiile

ℱ(T∗φ)=(ℱT)φ^

și

ℱ(Tφ)=(ℱT)∗φ^.

Aceste identități exprimă corespondența dintre convoluție și înmulțire în cadrul teoriei distribuțiilor temperate. Totuși, convoluția nu este definită în general pentru orice pereche de distribuții temperate, fiind necesare condiții suplimentare asupra distribuțiilor considerate.[45][80][81]

Generalizări

[modificare | modificare sursă]

Transformata Fourier–Stieltjes și măsurile

[modificare | modificare sursă]

Pentru o măsură Borel complexă finită μ pe ℝn, transformata Fourier–Stieltjes este definită prin

μ^(ξ)=∫ℝne−2πix⋅ξdμ(x).

Integrala există pentru orice ξ∈ℝn, deoarece funcția de sub semnul integralei are modulul egal cu 1, iar măsura are variație totală finită. Transformata este mărginită și uniform continuă și satisface ‖μ^‖∞≤‖μ‖, unde ‖μ‖ reprezintă variația totală a măsurii. Denumirea reflectă legătura dintre integrarea în raport cu o măsură și integrala de tip Stieltjes.[82][83][84]

Dacă măsura μ este absolut continuă în raport cu măsura Lebesgue și are densitatea f∈L1(ℝn), adică dμ(x)=f(x)dx, atunci μ^=f^. Definiția se aplică și măsurilor concentrate în puncte. Pentru măsura Dirac δa, concentrată în punctul a, transformata este δa^(ξ)=e−2πia⋅ξ. În cazul a=0, se obține δ0^=1. Concluzia lemei Riemann–Lebesgue nu este valabilă pentru toate măsurile finite, deoarece transformatele acestora nu tind neapărat la zero la infinit.[85][86]

În teoria probabilităților, funcția caracteristică a unei variabile aleatoare X cu valori în ℝn se definește prin

φX(u)=𝔼(eiu⋅X)=∫ℝneiu⋅xdμX(x)=μX^(−u2π), unde μX este legea de probabilitate a lui X.

Această reprezentare nu presupune existența unei densități în raport cu măsura Lebesgue și se aplică inclusiv legilor discrete sau mixte. Dacă legea admite o densitate, funcția caracteristică se exprimă prin transformata Fourier a acesteia, cu schimbarea de semn și de scară indicată în formulă.[87][88]

Teorema lui Bochner caracterizează transformatele Fourier–Stieltjes ale măsurilor pozitive finite. O funcție continuă pe ℝn este transformata unei asemenea măsuri dacă și numai dacă este pozitiv definită. Condiția suplimentară φ(0)=1 corespunde unei măsuri de probabilitate. Teorema are și o formulare pentru grupurile abeliene local compacte, în care măsura este definită pe grupul dual.[89][90]

Grupuri abeliene local compacte

[modificare | modificare sursă]

Pe un grup topologic abelian local compact și Hausdorff G, transformata Fourier se definește folosind o măsură Haar μ. Aceasta este invariantă la translații și este determinată până la înmulțirea cu o constantă pozitivă. Un caracter al grupului G este un morfism continuu de grupuri χ:G→𝕋, unde 𝕋={z∈ℂ∣|z|=1}. Mulțimea caracterelor, înzestrată cu înmulțirea punctuală și cu topologia convergenței uniforme pe mulțimi compacte, formează un grup abelian local compact.[89][75]

Pentru f∈L1(G,μ), transformata Fourier este definită prin

f^(χ)=∫Gf(x)χ(x)‾dμ(x), pentru orice χ∈G^.

Transformata aparține spațiului C0(G^) al funcțiilor continue care tind la zero la infinit. Acest rezultat reprezintă forma lemei Riemann–Lebesgue pentru grupuri abeliene local compacte. În cazul G=ℝn, caracterele sunt χξ(x)=e2πix⋅ξ, iar definiția coincide cu cea folosită în articol. Formulele de inversiune și Plancherel necesită alegerea unei măsuri Haar pe G^ compatibile cu măsura fixată pe G.[89][91]

Pentru grupul cercului 𝕋≃ℝ/ℤ, caracterele sunt funcțiile ek(x)=e2πikx, cu k∈ℤ. Măsura Haar normalizată corespunde integrării pe [0,1), iar transformata Fourier este șirul coeficienților

f^(k)=∫01f(x)e−2πikxdx.

Funcțiile ek formează o bază ortonormată a spațiului L2(𝕋). Pentru f∈L2(𝕋), dezvoltarea

f=∑k∈ℤf^(k)ek converge în norma L2.

Grupul dual al cercului este grupul discret ℤ, iar dualul lui ℤ este cercul. Pentru un grup ciclic finit, aceeași construcție conduce la transformata Fourier discretă. Seriile Fourier, transformata Fourier a șirurilor și transformata Fourier discretă se încadrează astfel într-o definiție comună.[89][75]

Descrierea prin caractere face legătura cu teoria reprezentărilor. Pentru un grup finit G, caracterele reprezentărilor ireductibile formează o bază ortonormată a spațiului funcțiilor constante pe clasele de conjugare, în raport cu produsul scalar

⟨α,β⟩=1|G|∑g∈Gα(g)β(g)‾.

Caracterul unei reprezentări de dimensiune finită este o combinație liniară a caracterelor ireductibile, iar coeficienții sunt multiplicitățile reprezentărilor corespunzătoare. Pentru un grup abelian, fiecare clasă de conjugare conține un singur element, iar reprezentările ireductibile sunt unidimensionale. În cazul cercului, descompunerea după caractere ia forma unei serii Fourier, cu convergență în norma spațiului Hilbert, în locul unei sume finite.[92][93]

Transformata Gelfand

[modificare | modificare sursă]

Pentru un grup abelian local compact G, spațiul L1(G), înzestrat cu produsul de convoluție, este o algebră Banach comutativă. Involuția este definită prin f∗(x)=f(x−1)‾, grupul fiind unimodular. Prin completarea în norma universală de tip C∗ se obține algebra C∗(G). Această normă este definită cu ajutorul reprezentărilor unitare și, pe L1(G), este mărginită superior de norma L1.[89]

Pentru o algebră C∗ comutativă A, spectrul Gelfand Δ(A) este spațiul funcționalelor liniare multiplicative nenule, înzestrat cu topologia slabă-*. Transformata Gelfand asociază fiecărui element a∈A funcția definită pe spectru prin

Γ(a)(φ)=φ(a), pentru φ∈Δ(A).

Această transformare este un izomorfism izometric de algebre cu involuție între A și C0(Δ(A)). Dacă algebra are unitate, spectrul este compact, iar algebra funcțiilor corespunzătoare este C(Δ(A)).[94]

În cazul A=C∗(G), spectrul Gelfand se identifică cu grupul dual G^. Restricția transformatei Gelfand la subalgebra densă L1(G) coincide cu transformata Fourier, dacă identificarea caracterelor respectă conjugarea complexă din definiția acesteia.[89]

Grupuri compacte neabeliene

[modificare | modificare sursă]

Transformata Fourier se poate extinde și la grupuri topologice compacte Hausdorff care nu sunt abeliene. Reprezentările unitare ireductibile pot avea în acest caz dimensiune mai mare decât 1. Valorile transformatei sunt operatori, care pot fi reprezentați prin matrice, în locul valorilor scalare din cazul abelian. Această extindere aparține teoriei reprezentărilor și analizei armonice necomutative.[95][93]

Fie λ măsura Haar normalizată a grupului compact G. Pentru fiecare clasă de echivalență σ a reprezentărilor unitare ireductibile se alege o reprezentare U(σ) pe un spațiu Hilbert Hσ de dimensiune finită dσ. Transformata unei funcții f∈L1(G) este definită prin

f^(σ)=∫Gf(g)Ug−1(σ)dλ(g).

Caracterul unitar al reprezentării implică relația Ug−1(σ)=(Ug(σ))∗. Integrala poate fi definită prin coeficienții matriceali. Pentru ξ,η∈Hσ, se are

⟨f^(σ)ξ,η⟩=∫Gf(g)⟨Ug−1(σ)ξ,η⟩dλ(g).

Produsul scalar este considerat liniar în prima variabilă.[95]

Pentru o măsură Borel complexă finită μ pe G, definiția corespunzătoare este

μ^(σ)=∫GUg−1(σ)dμ(g), iar transformata satisface supσ‖μ^(σ)‖≤‖μ‖.

Dacă dμ=fdλ, se obține transformata funcției f. Coeficienții matriceali determină în mod unic măsura, ceea ce asigură injectivitatea aplicației μ↦μ^.[95]

Ordinea factorilor în formula convoluției depinde de convenția aleasă pentru transformare. Dacă

(f∗g)(x)=∫Gf(y)g(y−1x)dλ(y), atunci definiția cu Ug−1(σ) conduce la relația
f∗g^(σ)=g^(σ)f^(σ).

Aceeași inversare a ordinii este valabilă pentru convoluția măsurilor finite. Dacă transformata este definită folosind Ug(σ) în locul lui Ug−1(σ), ordinea factorilor se păstrează. Involuția funcțiilor este dată de f∗(g)=f(g−1)‾, deoarece grupurile compacte sunt unimodulare.[95]

Teorema Peter–Weyl stabilește completitudinea coeficienților matriceali ai reprezentărilor ireductibile. Pentru f∈L2(G), formula de inversiune este

f(g)=∑σdσtr⁡(f^(σ)Ug(σ)), unde suma converge în norma L2, iar egalitatea este înțeleasă în acest spațiu.

Identitatea Plancherel corespunzătoare are forma

‖f‖22=∑σdσ‖f^(σ)‖HS2, unde ‖⋅‖HS desemnează norma Hilbert–Schmidt.

Pentru un grup abelian, toate dimensiunile dσ sunt egale cu 1, iar relațiile se reduc la formulele scalare.[95][93] Dualitatea Tannaka–Krein oferă o altă descriere a grupurilor compacte neabeliene, folosind întreaga categorie a reprezentărilor acestora, nu doar caracterele scalare.[96]

Aplicarea transformatei Fourier poate simplifica rezolvarea anumitor probleme matematice, inclusiv a ecuațiilor diferențiale. Soluția problemei inițiale se obține apoi prin transformata inversă.

Transformata Fourier permite exprimarea anumitor operații efectuate asupra funcțiilor în domeniul timpului sau în cel spațial prin operații mai simple în domeniul frecvențelor. În convenția cu nucleul e−2πiξx, derivarea în raport cu x corespunde înmulțirii transformatei cu 2πiξ, în condițiile de regularitate necesare. Această proprietate permite reducerea unor ecuații diferențiale cu coeficienți constanți la ecuații algebrice sau, în cazul ecuațiilor cu derivate parțiale, la ecuații diferențiale cu mai puține variabile. O altă proprietate fundamentală este corespondența dintre convoluția funcțiilor și înmulțirea transformatelor lor Fourier. După rezolvarea problemei în domeniul frecvențelor, aplicarea transformatei inverse permite determinarea soluției în variabilele inițiale. Aceste proprietăți sunt studiate în cadrul analizei armonice și stau la baza numeroaselor aplicații ale transformatei Fourier în matematică, fizică și prelucrarea semnalelor.

Ecuații cu derivate parțiale

[modificare | modificare sursă]

Una dintre aplicațiile importante ale transformatei Fourier este rezolvarea ecuațiilor cu derivate parțiale, în special a celor cu coeficienți constanți. Metoda permite transformarea derivatelor în raport cu anumite variabile în expresii algebrice, reducând astfel complexitatea ecuației inițiale. Numeroase ecuații fundamentale ale fizicii matematice, studiate încă din secolul al XIX-lea, pot fi abordate prin această metodă.

Un exemplu îl constituie ecuația propagării căldurii în cazul unidimensional, scrisă în variabile adimensionale sub forma

∂u∂t=∂2u∂x2.

Aplicarea transformatei Fourier în raport cu variabila spațială x, păstrând timpul t ca parametru, conduce la ecuația diferențială ordinară

∂u^∂t(ξ,t)=−4π2ξ2u^(ξ,t).

Pentru condiția inițială u(x,0)=f(x), soluția în domeniul frecvențelor este

u^(ξ,t)=e−4π2ξ2tf^(ξ), unde t>0.

Aplicarea transformatei Fourier inverse permite exprimarea soluției ca o convoluție între funcția inițială f și nucleul căldurii.[97]

Problema Cauchy pentru ecuația undelor

[modificare | modificare sursă]

Pentru ecuația undelor unidimensională, cu viteza de propagare egală cu unitatea,

∂2u∂t2=∂2u∂x2, unde x,t∈ℝ, problema Cauchy constă în determinarea unei soluții care satisface condițiile inițiale
u(x,0)=f(x),∂u∂t(x,0)=g(x).

Funcțiile f și g reprezintă deplasarea inițială, respectiv viteza inițială. Spre deosebire de ecuația propagării căldurii, care este de ordinul întâi în raport cu timpul și necesită o singură condiție inițială, ecuația undelor este de ordinul al doilea în raport cu această variabilă. Cunoașterea deplasării inițiale nu este suficientă pentru determinarea soluției, fiind necesară și precizarea vitezei inițiale. Dacă datele inițiale îndeplinesc condițiile de regularitate corespunzătoare, cele două funcții determină în mod unic soluția clasică.[97][98]

Rezolvarea prin transformata Fourier poate fi realizată mai întâi pentru f,g∈𝒮(ℝ), unde 𝒮(ℝ) este spațiul Schwartz. Se notează cu

v(ξ,t)=u(⋅,t)^(ξ), transformata Fourier a funcției u în raport cu variabila spațială x.

Proprietatea transformatei Fourier privind derivarea conduce, pentru fiecare ξ∈ℝ, la ecuația diferențială ordinară

∂2v∂t2+(2πξ)2v=0, însoțită de condițiile inițiale v(ξ,0)=f^(ξ) și ∂v∂t(ξ,0)=g^(ξ).

Soluția acestei probleme este

v(ξ,t)=f^(ξ)cos⁡(2πξt)+g^(ξ)sin⁡(2πξt)2πξ.

Pentru ξ=0, raportul sin⁡(2πξt)2πξ se prelungește prin continuitate, valoarea sa fiind t. În acest punct, soluția ia forma v(0,t)=f^(0)+tg^(0).

Pentru fiecare valoare fixată a timpului t, funcția v(⋅,t) aparține spațiului Schwartz. Se justifică astfel aplicarea transformatei Fourier inverse și derivarea sub semnul integralei, obținându-se reprezentarea

u(x,t)=∫−∞∞[f^(ξ)cos⁡(2πξt)+g^(ξ)sin⁡(2πξt)2πξ]e2πiξxdξ.

Transformarea Fourier reduce astfel ecuația cu derivate parțiale la o familie de ecuații diferențiale ordinare de tipul oscilatorului armonic, parametrizate prin frecvența ξ. Soluția problemei inițiale se reconstituie prin transformarea inversă, care reunește contribuțiile tuturor frecvențelor.[98]

Reprezentarea obținută poate fi exprimată direct prin funcțiile inițiale f și g. Aplicarea transformatei inverse primului termen conduce la relația

ℱ−1[f^(ξ)cos⁡(2πξt)](x)=f(x+t)+f(x−t)2.

Pentru t≥0, se consideră funcția

Kt(x)=12𝟏[−t,t](x), unde 𝟏[−t,t] este funcția caracteristică a intervalului [−t,t]. Transformata sa Fourier este
K^t(ξ)=sin⁡(2πξt)2πξ, cu valoarea t în punctul ξ=0, obținută prin continuitate.

Potrivit teoremei convoluției, transformata inversă a celui de-al doilea termen este:

(g∗Kt)(x)=12∫x−tx+tg(s)ds.

Prin reunirea celor două expresii se obține formula lui d'Alembert

u(x,t)=f(x+t)+f(x−t)2+12∫x−tx+tg(s)ds.

Pentru t<0, formula rămâne valabilă dacă integrala este considerată orientată, cu limitele de integrare în ordinea indicată. Prin derivare se verifică faptul că funcția astfel definită satisface atât ecuația undelor, cât și cele două condiții inițiale. Formula lui d'Alembert este valabilă și pentru date inițiale f∈C2(ℝ) și g∈C1(ℝ), caz în care determină unica soluție clasică a problemei Cauchy. Această generalizare nu necesită existența transformatelor Fourier ale datelor inițiale ca integrale absolut convergente, deoarece expresia soluției poate fi definită și verificată direct.[97][98]

Suprapunerea undelor elementare

[modificare | modificare sursă]

Soluția ecuației undelor poate fi reprezentată prin suprapunerea unor unde elementare. Funcțiile cos⁡(2πξ(x±t)) și sin⁡(2πξ(x±t)) satisfac fiecare ecuația undelor, semnele corespunzând celor două sensuri de propagare. În condițiile care justifică convergența integralelor și derivarea sub semnul integralei, o suprapunere continuă a acestor funcții constituie, la rândul ei, o soluție

u(x,t)=∫0∞[a+(ξ)cos⁡(2πξ(x+t))+a−(ξ)cos⁡(2πξ(x−t))+b+(ξ)sin⁡(2πξ(x+t))+b−(ξ)sin⁡(2πξ(x−t))]dξ.

Pentru date inițiale reale, coeficienții acestei reprezentări pot fi determinați cu ajutorul transformatelor în cosinus și sinus. Dacă, pentru o funcție h, se definesc

Ah(ξ)=2∫−∞∞h(x)cos⁡(2πξx)dx,Bh(ξ)=2∫−∞∞h(x)sin⁡(2πξx)dx, cele două condiții inițiale ale problemei Cauchy duc la sistemul
a++a−=Af,b++b−=Bf,2πξ(−a++a−)=Bg,2πξ(b+−b−)=Ag.

Pentru fiecare ξ≠0, sistemul conține patru ecuații liniare care determină coeficienții a+, a−, b+ și b−. Datele inițiale și soluția sunt astfel exprimate prin aceeași familie de unde elementare, iar transformata Fourier inversă permite reconstituirea funcției u din amplitudinile componentelor sale.

Această descompunere admite și o reprezentare cu exponențiale complexe. Pentru ξ≠0, soluția transformată poate fi scrisă sub forma

v(ξ,t)=s+(ξ)e2πiξt+s−(ξ)e−2πiξt, unde coeficienții sunt
s±(ξ)=12[f^(ξ)±g^(ξ)2πiξ].

Chiar dacă datele inițiale aparțin spațiului Schwartz, amplitudinile s+ și s− pot prezenta separat o singularitate de ordinul 1|ξ| în vecinătatea originii. Din acest motiv, ele nu aparțin în mod necesar spațiilor L1(ℝ) sau L2(ℝ). Reprezentarea soluției prin sinus și cosinus evită această dificultate, deoarece raportul sin⁡(2πξt)2πξ admite o prelungire continuă în ξ=0, cu valoarea t. Frecvența nulă este astfel inclusă fără introducerea unei singularități în expresia soluției.

O interpretare suplimentară se obține prin aplicarea transformatei Fourier atât în variabila spațială x, cât și în variabila temporală t. Pentru date inițiale din spațiul Schwartz, formula lui d'Alembert definește o funcție mărginită pe ℝ2, care poate fi considerată o distribuție temperată. Dacă u~(ξ,ν) reprezintă transformata Fourier a acestei distribuții în ambele variabile, ecuația undelor conduce la relația

(ξ2−ν2)u~(ξ,ν)=0.

Această egalitate implică faptul că suportul distribuției u~ este inclus în reuniunea dreptelor ν=ξ și ν=−ξ din planul frecvențelor. Cele două drepte corespund familiilor de unde care se propagă în sensuri opuse, iar contribuțiile lor pot fi reprezentate prin distribuții Dirac concentrate pe aceste drepte. Condiția privind suportul este necesară, dar nu este suficientă, în general, pentru satisfacerea ecuației distribuționale de mai sus. Relațiile trebuie interpretate prin dualitatea cu funcțiile-test din spațiul Schwartz, deoarece exponențialele complexe utilizate în reprezentarea undelor nu aparțin ele însele acestui spațiu.[99] Din punct de vedere al calculului, metoda transformatei Fourier presupune transformarea datelor inițiale, rezolvarea ecuațiilor obținute în domeniul frecvențelor și aplicarea transformatei inverse.

Transformata Fourier neliniară

[modificare | modificare sursă]

Metodele bazate pe transformata Fourier au fost extinse, în cursul secolului al XX-lea, de la ecuațiile liniare la anumite clase de ecuații diferențiale neliniare integrabile. În cazul ecuațiilor liniare cu coeficienți polinomiali, proprietățile transformatei Fourier permit înlocuirea derivării și a înmulțirii cu variabila prin operații corespunzătoare în domeniul frecvențelor. Ecuația obținută poate conține în continuare derivate, astfel încât aplicarea transformatei nu conduce întotdeauna la o ecuație algebrică. Pentru anumite ecuații neliniare de evoluție, rezolvarea poate fi raportată la o problemă spectrală liniară, ale cărei valori proprii se păstrează în timpul evoluției. Această proprietate este exprimată prin formularea lui Peter Lax, în care evoluția neliniară este asociată unei familii de operatori cu același spectru.[100][101]

Procedeul presupune determinarea datelor spectrale corespunzătoare condiției inițiale, calcularea evoluției acestora și reconstituirea soluției prin rezolvarea problemei spectrale inverse. Această metodă, denumită metoda transformatei inverse de împrăștiere, stă la baza unor transformări cunoscute sub numele de transformate Fourier neliniare. Spre deosebire de transformata Fourier clasică, acestea nu constituie o transformare liniară aplicabilă oricărei ecuații neliniare, ci sunt definite în raport cu structura spectrală a anumitor sisteme integrabile. În telecomunicațiile optice, metodele respective sunt utilizate pentru analiza unor modele de propagare descrise de ecuația Schrödinger neliniară și pentru dezvoltarea unor tehnici de transmitere a informației prin reprezentarea spectrală neliniară a semnalelor.[102]

Spectroscopie prin transformare Fourier

[modificare | modificare sursă]

În rezonanța magnetică nucleară, transformata Fourier permite obținerea spectrului din semnalul înregistrat după excitarea probei printr-un impuls de radiofrecvență. Acest semnal, numit semnal de inducție liberă, conține contribuții oscilante ale căror amplitudini scad în timp. Transformarea lor separă componentele de frecvență și permite identificarea liniilor de rezonanță. Pentru o componentă cu descreștere exponențială, profilul de absorbție este lorentzian, după corectarea corespunzătoare a fazei. Un timp de descreștere mai scurt corespunde unei linii spectrale mai largi.[103] În spectroscopia în infraroșu prin transformare Fourier, informația spectrală este obținută dintr-o interferogramă, care descrie variația semnalului detectat în funcție de diferența de drum optic dintre brațele interferometrului. Variabila transformată este aici diferența de drum optic, nu timpul de evoluție a probei.[104] Spectrul de transmitanță se obține prin raportarea spectrului măsurat cu proba la un spectru de referință, iar din acesta se calculează absorbanta. Benzile de absorbție oferă informații despre vibrațiile moleculare.[105]

Transformata Fourier este utilizată și în anumite metode de spectrometrie de masă, între care spectrometria prin rezonanță ciclotronică ionică cu transformare Fourier, abreviată FT-ICR. Mișcarea coerentă a ionilor într-un câmp magnetic induce un semnal electric în electrozii de detecție. Transformata Fourier separă frecvențele acestui semnal, iar acestea sunt convertite, prin relațiile de calibrare ale instrumentului, în valori ale raportului masă–sarcină.[106] O aplicație înrudită este reconstrucția imaginilor în imagistica prin rezonanță magnetică. Gradienții câmpului magnetic introduc o dependență a frecvenței și a fazei semnalului de poziția spațială. Transformata Fourier inversă permite folosirea acestei informații pentru reconstruirea imaginii. În acest caz, scopul principal este localizarea spațială a surselor semnalului, nu separarea componentelor după mediul lor chimic.[107]

Mecanică cuantică

[modificare | modificare sursă]

În mecanica cuantică, transformata Fourier face legătura dintre reprezentarea unei stări în spațiul pozițiilor și reprezentarea sa în spațiul impulsurilor. Pentru o particulă care se deplasează într-o singură dimensiune, fără a considera gradele de libertate asociate spinului, o stare pură poate fi descrisă prin funcția de undă ψ(q) în reprezentarea poziției sau prin funcția φ(p) în reprezentarea impulsului. Cele două funcții descriu aceeași stare, dar exprimă diferit informațiile despre rezultatele măsurării poziției și impulsului.[108]

Notând cu ℏ=h2π constanta Planck redusă, trecerea în reprezentarea impulsului este dată de

φ(p)=12πℏ∫−∞∞ψ(q)e−ipqℏdq, iar revenirea în reprezentarea poziției se face prin transformata inversă
ψ(q)=12πℏ∫−∞∞φ(p)eipqℏdp.

Produsul pq are dimensiunea unei acțiuni, la fel ca ℏ, iar raportul pqℏ este adimensional. Aceste formule corespund convenției Fourier cu nucleul e−2πiqξ prin relația p=2πℏξ. Schimbarea variabilei spectrale este însoțită de factorul de normalizare care asigură caracterul unitar al transformării.[109]

Pentru o stare normalizată, transformarea păstrează norma și se obține

∫−∞∞|ψ(q)|2dq=∫−∞∞|φ(p)|2dp=1.

Funcțiile |ψ(q)|2 și |φ(p)|2 sunt densitățile de probabilitate pentru poziție, respectiv impuls. Pentru funcții de undă din L2(ℝ), transformarea este înțeleasă prin prelungirea unitară la acest spațiu, chiar dacă integralele nu converg absolut.[110][111]

Poziția și impulsul sunt observabile canonic conjugate. Legătura Fourier dintre cele două reprezentări conduce la relația de incertitudine

ΔqΔp≥ℏ2, unde Δq și Δp sunt abaterile standard ale distribuțiilor de probabilitate corespunzătoare, presupuse finite.

Relația limitează precizia cu care poziția și impulsul pot fi determinate simultan într-o aceeași stare.[112] Pentru o particulă nerelativistă de masă m, aflată într-un potențial real V(q) independent de timp, ecuația lui Schrödinger unidimensională este

iℏ∂ψ∂t=−ℏ22m∂2ψ∂q2+V(q)ψ.

În cazul particulei libere, pentru care V=0, aplicarea transformatei Fourier în variabila spațială conduce la

iℏ∂φ∂t(p,t)=p22mφ(p,t).

Pentru fiecare valoare a impulsului p, aceasta este o ecuație diferențială ordinară în raport cu timpul. Soluția care corespunde stării inițiale φ(p,0) este

φ(p,t)=e−ip2t2mℏφ(p,0).

Transformata inversă determină funcția de undă ψ(q,t). Factorul exponențial are modulul egal cu 1 și modifică numai faza fiecărei componente de impuls. Evoluția liberă păstrează astfel densitatea de probabilitate în spațiul impulsurilor și norma stării.[113]

Pentru un potențial care depinde de poziție, termenul V(q)ψ(q,t) introduce, în general, o legătură între componentele de impuls, iar ecuațiile acestora nu mai pot fi rezolvate separat. Dacă potențialul este independent de timp, stările staționare au forma ψ(q,t)=e−iEtℏψE(q), unde ψE satisface ecuația de valori proprii H^ψE=EψE pentru operatorul hamiltonian. Dezvoltarea spectrală a stării inițiale permite studierea evoluției sale, ținând seama de eventualele componente discrete și continue ale spectrului energetic.[114][115]

În teoria cuantică a câmpurilor, transformata Fourier permite descompunerea unui câmp liber în moduri care evoluează independent. Pentru un câmp scalar Φ, într-o singură dimensiune spațială și cu variabile adimensionale alese astfel încât termenul de masă să aibă coeficientul egal cu 1, ecuația Klein–Gordon se scrie

(∂2∂t2−∂2∂x2+1)Φ(x,t)=0.

Aplicând transformata Fourier numai în variabila x, cu nucleul e−2πixξ, se obține

∂2Φ^∂t2(ξ,t)+[1+(2πξ)2]Φ^(ξ,t)=0.

Fiecare mod Fourier satisface o ecuație de oscilator armonic cu pulsația ωξ=1+(2πξ)2, exprimată în variabilele adimensionale considerate. Cuantizarea acestor moduri permite descrierea câmpului liber prin operatori de creare și anihilare. În prezența interacțiunilor, modurile sunt, în general, cuplate și nu mai evoluează independent.[116]

Legătura cu oscilatorul armonic apare și în proprietățile transformatei Fourier ca operator. După alegerea unor variabile adimensionale și a normalizării corespunzătoare, transformata poate fi exprimată printr-un operator unitar generat de operatorul numărului de excitații. Funcțiile Hermite sunt funcții proprii comune celor doi operatori, iar formula lui Mehler permite exprimarea nucleului integral al familiei de transformări asociate.[117]

Prelucrarea semnalelor

[modificare | modificare sursă]

Transformata Fourier are un rol important în prelucrarea semnalelor, fiind utilizată pentru analiza spectrală a semnalelor și a seriilor temporale. În cazul semnalelor de energie finită, distribuția energiei în domeniul frecvențelor poate fi determinată direct din transformata Fourier. Pentru semnalele persistente, care pot avea energie totală infinită, dar putere medie finită, transformata definită printr-o integrală pe întreaga axă reală nu există întotdeauna în sens clasic. Analiza spectrală se poate realiza atunci prin intermediul autocorelației, care descrie dependența dintre valorile unui semnal corespunzătoare unor momente diferite.

Pentru un semnal complex f∈L1(ℝ)∩L2(ℝ), funcția de autocorelație de energie este definită prin

Rf(E)(τ)=∫−∞∞f(t)‾f(t+τ)dt.

Variabila τ reprezintă decalajul temporal dintre valorile comparate, iar conjugarea complexă permite extinderea definiției la semnale cu valori complexe. Din proprietățile convoluției și ale conjugării rezultă relația

ℱ[Rf(E)](ξ)=|f^(ξ)|2.

Aceasta exprimă faptul că transformata Fourier a autocorelației de energie este egală cu pătratul modulului transformatei semnalului, care reprezintă densitatea sa spectrală de energie. Pentru semnalele reale, conjugarea complexă nu modifică expresia autocorelației.[118]

În cazul unui semnal persistent, autocorelația poate fi definită prin mediere temporală, dacă limita corespunzătoare există

Rf(τ)=limT→∞12T∫−TTf(t)‾f(t+τ)dt.

Factorul 12T reprezintă inversul duratei intervalului de observație. Spre deosebire de autocorelația de energie, această definiție se bazează pe puterea medie a semnalului și este potrivită pentru analiza semnalelor persistente cu putere finită. Valoarea Rf(0) corespunde puterii medii, atunci când limita există.

Pentru un proces aleatoriu complex X(t), staționar în sens larg și cu moment de ordinul al doilea finit, autocorelația se definește prin media statistică

RX(τ)=𝔼[X(t)‾X(t+τ)].

Staționaritatea în sens larg implică independența acestei funcții față de momentul inițial t, astfel încât autocorelația depinde numai de decalajul temporal τ. Dacă procesul are media constantă m=𝔼[X(t)], funcția de autocovarianță este

CX(τ)=𝔼[X(t)−m‾(X(t+τ)−m)]=RX(τ)−|m|2.

Pentru un proces centrat, caracterizat prin m=0, autocorelația coincide cu autocovarianța. Centrarea procesului, realizată prin scăderea mediei, trebuie deosebită de normalizarea autocorelației prin valoarea sa la origine, atunci când aceasta este nenulă. De asemenea, egalitatea dintre media temporală și media statistică nu rezultă numai din staționaritate, fiind necesare condiții suplimentare de ergodicitate.[119]

Funcția de autocorelație a unui proces staționar în sens larg satisface proprietatea de simetrie hermitiană și inegalitatea

RX(−τ)=RX(τ)‾,|RX(τ)|≤RX(0).

Pentru un proces real, autocorelația este o funcție pară și are valoarea maximă în origine. Dacă procesul este continuu în medie pătratică, autocorelația este continuă și, datorită caracterului său pozitiv definit, este uniform continuă. Aceste proprietăți permit studierea dependențelor temporale și a fenomenelor periodice. De exemplu, în analiza unui șir de măsurători ale temperaturii, valorile autocorelației pentru decalaje apropiate de 24 de ore pot evidenția existența unui ciclu zilnic.

Pentru un proces staționar în sens larg care este continuu în medie pătratică, teorema lui Bochner asigură existența unei măsuri Borel pozitive finite μX, numită măsură spectrală, astfel încât autocorelația permite reprezentarea

RX(τ)=∫ℝe2πiξτdμX(ξ).

Masa totală a acestei măsuri este μX(ℝ)=RX(0), adică puterea medie a procesului. Reprezentarea spectrală permite descrierea distribuției puterii între diferitele frecvențe, inclusiv în situațiile în care spectrul nu poate fi exprimat printr-o funcție de densitate.[120]

Dacă autocorelația aparține și spațiului L1(ℝ), măsura spectrală este absolut continuă în raport cu măsura Lebesgue și admite densitatea

PX(ξ)=∫−∞∞RX(τ)e−2πiξτdτ, astfel încât dμX(ξ)=PX(ξ)dξ.

Funcția PX este continuă și nenegativă și se numește densitate spectrală de putere. Relația dintre funcția de autocorelație și spectrul de putere este exprimată de teorema Wiener–Hincin. Această corespondență rămâne valabilă în formularea generală prin măsuri spectrale, fără a impune existența unei densități.[121]

Componentele periodice ale unui proces pot determina apariția unor linii spectrale, reprezentate prin mase punctuale ale măsurii spectrale. Dacă procesul are o medie nenulă m, spectrul său conține, în plus față de cel al procesului centrat, o masă spectrală |m|2 concentrată la frecvența zero. Spectrul procesului centrat descrie repartizarea varianței între frecvențe, în timp ce spectrul obținut din autocorelația procesului necentrat include și contribuția mediei.[122] În aplicațiile fizice și tehnice, densitatea spectrală de putere caracterizează contribuția diferitelor benzi de frecvență la puterea medie a semnalului.

Metode de calcul

[modificare | modificare sursă]

Metoda de calcul depinde de reprezentarea datelor de intrare și de forma rezultatului cerut. O funcție de variabilă continuă poate fi dată printr-o expresie analitică, în timp ce un semnal măsurat este adesea cunoscut numai prin perechi de valori (xj,fj). O transformare poate fi urmărită ca funcție de frecvență pe un întreg interval sau numai la un număr finit de frecvențe. Trecerea de la integrală la o sumă trebuie să țină seama dacă datele descriu o măsură discretă, o secvență sau eșantioane folosite pentru aproximarea unei funcții continue.

Transformata Fourier discretă și algoritmii rapizi

[modificare | modificare sursă]

Pentru un șir (an)n∈ℤ de valori corespunzătoare unor momente echidistante nΔ, unde Δ>0 este intervalul de eșantionare, transformata Fourier în timp discret este definită, în condițiile de convergență corespunzătoare, prin

A(ξ)=∑n∈ℤane−2πiξnΔ.

Deși semnalul este discret în timp, variabila de frecvență ξ rămâne continuă, iar transformata este periodică, având perioada 1Δ. Dacă șirul conține un număr finit de termeni nenuli, transformata sa este un polinom trigonometric.

Evaluarea transformatei unui șir finit în puncte echidistante din domeniul frecvențelor conduce la transformata Fourier discretă (TFD). Pentru un șir format din N valori, aceasta este definită prin

Ak=∑n=0N−1ane−2πiknN, unde k=0,1,…,N−1.

Valorile Ak corespund frecvențelor ξk=kNΔ, considerate modulo perioada spectrală 1Δ. Transformata Fourier discretă asociază astfel unui șir finit un alt șir finit, ale cărui elemente descriu componentele sale spectrale.

Calculul direct al celor N coeficienți necesită un număr de operații de ordinul N2. Algoritmii de transformată Fourier rapidă valorifică structura sumei pentru a reduce volumul calculelor, în numeroase implementări, la un ordin de complexitate Nlog⁡N. Acești algoritmi calculează aceeași transformată discretă, diferențele față de rezultatul matematic exact provenind din erorile de rotunjire ale calculului numeric.

Lungimea secvenței și intervalul de eșantionare determină distribuția punctelor de evaluare în domeniul frecvențelor, distanța dintre frecvențele consecutive fiind Δξ=1NΔ. Această distanță nu trebuie confundată cu posibilitatea efectivă de a separa două componente spectrale apropiate, care depinde și de durata observației și de metoda de analiză. Transformata Fourier discretă reprezintă calculul unei sume finite, nu evaluarea directă a integralei Fourier a unui semnal continuu. Interpretarea rezultatelor ca aproximații ale transformatei continue necesită luarea în considerare a eșantionării, a trunchierii semnalului și a metodei de aproximare a integralei.[123]

Integrare simbolică

[modificare | modificare sursă]

Pentru anumite funcții exprimate analitic, transformata Fourier poate fi calculată exact prin evaluarea integralei corespunzătoare. Rezultatele unor asemenea calcule sunt reunite în tabele de transformate Fourier, care oferă expresii valabile pentru diferite clase de funcții, în condițiile de existență și pentru valorile parametrilor precizate. Aceste relații permit determinarea transformatei fără repetarea calculelor de integrare.[124]

Calculul simbolic al transformatelor Fourier poate fi realizat și cu ajutorul programelor specializate, precum MATLAB, prin intermediul pachetului Symbolic Math Toolbox. Pentru interpretarea corectă a rezultatelor trebuie cunoscute convenția de semn, factorii de normalizare și variabila spectrală utilizată. De exemplu, pentru funcția f(t)=cos⁡(6πt)e−πt2, transformata Fourier în convenția cu nucleul e−2πiξt este

f^(ξ)=12[e−π(ξ−3)2+e−π(ξ+3)2].

Rezultatul poate fi obținut în MATLAB prin următorul cod, care utilizează funcția simbolică fourier

syms t w xi real

f = cos(6*sym(pi)*t)*exp(-sym(pi)*t^2);

F_w = fourier(f, t, w);
F_xi = simplify(subs(F_w, w, 2*sym(pi)*xi));

disp(F_xi)

Funcția fourier din MATLAB utilizează implicit transformarea cu nucleul e−iωt. În exemplul prezentat, substituția ω=2πξ exprimă rezultatul în funcție de frecvența ξ, potrivit convenției adoptate în articol. Transformata obținută este suma a două funcții gaussiene centrate la frecvențele ξ=−3 și ξ=3. Această proprietate rezultă din descompunerea funcției cosinus în exponențiale complexe și din proprietatea de modulație a transformatei Fourier.

MATLAB permite și modificarea directă a convenției de transformare prin funcția sympref. Parametrul FourierParameters stabilește factorul de normalizare și coeficientul exponentului complex. Valorile implicite sunt [1, -1], iar alegerea [1, -2*pi] corespunde convenției cu nucleul e−2πiξt. La efectuarea calculelor trebuie verificată convenția activă, deoarece aceasta poate modifica forma expresiilor rezultate.[125]

Integrare numerică a funcțiilor continue

[modificare | modificare sursă]

Atunci când transformata Fourier nu poate fi determinată analitic sau expresia sa exactă este dificil de utilizat, integrala care o definește poate fi aproximată prin metode numerice. Calculul se efectuează pentru frecvențele considerate, folosind reguli de integrare numerică adaptate proprietăților funcției. În cazul integralelor pe întreaga axă reală, aproximarea presupune, de regulă, restrângerea domeniului de integrare la un interval finit și aplicarea unei formule de cuadratură. Precizia rezultatului depinde de eroarea introdusă prin trunchiere, de regularitatea funcției, de alegerea punctelor de integrare și de comportamentul oscilant al nucleului Fourier.[123] Dacă transformata trebuie evaluată pentru un număr mare de frecvențe, calcularea separată a fiecărei integrale poate necesita un volum considerabil de operații.[126][127]

Integrare numerică din date discrete

[modificare | modificare sursă]

Atunci când valorile unei funcții sunt cunoscute numai într-un număr finit de puncte, transformata Fourier poate fi aproximată prin aplicarea unei formule de integrare numerică. Integrala este înlocuită cu o sumă ponderată de forma

f^(ξ)≈∑jwjf(xj)e−2πiξxj.

Punctele xj reprezintă pozițiile de eșantionare, iar ponderile wj sunt stabilite prin metoda de integrare aleasă. Pentru puncte echidistante, o aproximare simplă utilizează ponderi egale cu pasul de eșantionare, cu eventualele ajustări la capetele intervalului. Formulele de cuadratură de ordin superior, precum regula lui Simpson, folosesc ponderi diferite și pot oferi o precizie mai bună pentru funcțiile care îndeplinesc condițiile de regularitate necesare. Erorile aproximării depind atât de formula de integrare, cât și de distribuția punctelor și de comportamentul oscilant al integrandului.[128]

O sumă de aceeași formă poate avea și o interpretare exactă, dacă datele reprezintă o măsură discretă finită. Pentru măsura μ=∑jcjδxj, unde δxj este măsura Dirac concentrată în punctul xj, transformata Fourier–Stieltjes este μ^(ξ)=∑jcje−2πiξxj.

În acest caz, egalitatea este exactă, iar coeficienții cj reprezintă masele atribuite punctelor, nu ponderi introduse pentru aproximarea integralei unei funcții continue. Cele două interpretări trebuie deosebite, chiar dacă expresiile matematice au aceeași structură.

Tabele de transformate Fourier

[modificare | modificare sursă]

Relațiile de mai jos folosesc cele trei convenții definite anterior, în aceeași ordine. Numerele rândurilor permit compararea proprietăților duale. Identitățile între funcții de pătrat integrabil sunt înțelese în sensul lui Plancherel, iar cele care implică distribuții sunt egalități în 𝒮′. Definițiile funcțiilor auxiliare și ipotezele asupra parametrilor fac parte din fiecare pereche.[129][130]

Relații operaționale într-o variabilă

[modificare | modificare sursă]
Nr.FuncțieTransformată unitară în frecvență obișnuită ξTransformată unitară în frecvență unghiulară ωTransformată neunitară în frecvență unghiulară ωCondiții și observații
101 af(x)+bg(x) af^(ξ)+bg^(ξ) af^(ω)+bg^(ω) af^(ω)+bg^(ω) Liniaritatea, pentru constante complexe a,b.
102 f(x−a) e−i2πξaf^(ξ) e−iaωf^(ω) e−iaωf^(ω) Translație în variabila inițială, pentru a∈ℝ.
103 f(x)eiax f^(ξ−a2π) f^(ω−a) f^(ω−a) Modulație, echivalentă cu o translație în frecvență.
104 f(ax) 1|a|f^(ξa) 1|a|f^(ωa) 1|a|f^(ωa) Schimbare de scară, pentru a∈ℝ∖{0}. O comprimare în variabila inițială produce o dilatare a spectrului.
105 f^n(x) f^1(x) ⟷ℱ1 f(−ξ) f^2(x) ⟷ℱ2 f(−ω) f^3(x) ⟷ℱ3 2πf(−ω) Aplicarea de două ori a aceleiași transformări. Variabila frecvență a primei aplicări este redenumită x.
106 dnf(x)dxn (i2πξ)nf^(ξ) (iω)nf^(ω) (iω)nf^(ω) Derivata de ordin n. Formula este valabilă, în particular, pentru f∈𝒮.
106.5 ∫−∞xf(τ)dτ pv⁡f^(ξ)i2πξ+Cδ(ξ) pv⁡f^(ω)iω+2πCδ(ω) pv⁡f^(ω)iω+2πCδ(ω) Integrare, pentru f∈𝒮. Câturile se interpretează prin valoare principală și C=12∫ℝf(t)dt. Această constantă reprezintă media limitelor primitivei la cele două capete ale axei.[131]
107 xnf(x) (i2π)ndnf^(ξ)dξn indnf^(ω)dωn indnf^(ω)dωn Înmulțire cu un monom. Pentru formula clasică este suficient ca momentele necesare să fie integrabile. Pe 𝒮 toate ordinele sunt admise.
108 (f∗g)(x) f^(ξ)g^(ξ) 2π f^(ω)g^(ω) f^(ω)g^(ω) Convoluție, pentru f,g∈L1.
109 f(x)g(x) (f^∗g^)(ξ) 12π(f^∗g^)(ω) 12π(f^∗g^)(ω) Produs punctual. Este valabil, de exemplu, pentru f,g∈𝒮, iar formula se extinde în condițiile de integrabilitate corespunzătoare.
110 Pentru f cu valori reale f^(−ξ)=f^(ξ)‾ f^(−ω)=f^(ω)‾ f^(−ω)=f^(ω)‾ Simetrie hermitiană pentru o funcție reală. Bara reprezintă conjugarea complexă.
113 Pentru f cu valori pur imaginare f^(−ξ)=−f^(ξ)‾ f^(−ω)=−f^(ω)‾ f^(−ω)=−f^(ω)‾ Simetrie antihermitiană pentru o funcție cu valori pur imaginare.
114 f(x)‾ f^(−ξ)‾ f^(−ω)‾ f^(−ω)‾ Conjugare complexă, generalizând cazurile reale și pur imaginare.
115 f(x)cos⁡(ax) f^(ξ−a2π)+f^(ξ+a2π)2 f^(ω−a)+f^(ω+a)2 f^(ω−a)+f^(ω+a)2 Rezultă din liniaritate, modulație și cos⁡(ax)=(eiax+e−iax)/2.
116 f(x)sin⁡(ax) f^(ξ−a2π)−f^(ξ+a2π)2i f^(ω−a)−f^(ω+a)2i f^(ω−a)−f^(ω+a)2i Rezultă din liniaritate, modulație și sin⁡(ax)=(eiax−e−iax)/(2i).

Funcții cu pătratul integrabil

[modificare | modificare sursă]

Perechile următoare sunt întâlnite și în tabelele de integrale Fourier pentru aplicații.[132]

Nr.FuncțieTransformată unitară în frecvență obișnuită ξTransformată unitară în frecvență unghiulară ωTransformată neunitară în frecvență unghiulară ωCondiții și observații
201 rect⁡(ax) 1|a|sinc⁡(ξa) 12πa2sinc⁡(ω2πa) 1|a|sinc⁡(ω2πa) Impuls dreptunghiular și sinus cardinal normalizat. Parametrul a este real și nenul.
202 sinc⁡(ax) 1|a|rect⁡(ξa) 12πa2rect⁡(ω2πa) 1|a|rect⁡(ω2πa) Perechea duală a dreptunghiului, pentru a real și nenul. Funcția sinc aparține lui L2, dar nu lui L1. În interpretarea ca filtru trece-jos ideal, răspunsul la impuls este necauzal.
203 sinc2(ax) 1|a|tri⁡(ξa) 12πa2tri⁡(ω2πa) 1|a|tri⁡(ω2πa) Funcția tri⁡(x)=max⁡(1−|x|,0) este triunghiulară. Parametrul a este real și nenul.
204 tri⁡(ax) 1|a|sinc2(ξa) 12πa2sinc2(ω2πa) 1|a|sinc2(ω2πa) Perechea duală a pătratului funcției sinc, pentru a real și nenul.
205 e−axu(x) 1a+i2πξ 12π(a+iω) 1a+iω u este funcția treaptă Heaviside și a>0.
206 e−αx2 παe−(πξ)2α 12αe−ω24α παe−ω24α Re⁡α>0, cu rădăcina pătrată aleasă prin continuare din semiaxa reală pozitivă. Pentru o anumită scară, gaussiana este invariantă la transformarea unitară.
208 e−a|x| 2aa2+4π2ξ2 2πaa2+ω2 2aa2+ω2 Re⁡a>0. Exponențiala bilateral descrescătoare se transformă într-un profil de tip Lorentz.
209 sech⁡(ax) πasech⁡(π2aξ) 1aπ2sech⁡(π2aω) πasech⁡(π2aω) a>0. Secanta hiperbolică este invariantă la transformare pentru scara potrivită.
210 e−a2x22Hn(ax) 2π(−i)nae−2π2ξ2a2Hn(2πξa) (−i)nae−ω22a2Hn(ωa) (−i)n2πae−ω22a2Hn(ωa) a>0. Hn este polinomul Hermite în convenția fizicienilor. Pentru transformata unitară în frecvență unghiulară și a=1, funcțiile Gauss–Hermite sunt funcții proprii. Cazul n=0 este gaussian.

Distribuții într-o variabilă

[modificare | modificare sursă]
Nr.FuncțieTransformată unitară în frecvență obișnuită ξTransformată unitară în frecvență unghiulară ωTransformată neunitară în frecvență unghiulară ωCondiții și observații
301 1 δ(ξ) 2πδ(ω) 2πδ(ω) δ este distribuția Dirac.
302 δ(x) 1 12π 1 Perechea duală a constantei.
303 eiax δ(ξ−a2π) 2πδ(ω−a) 2πδ(ω−a) Exponențială complexă, cu a∈ℝ.
304 cos⁡(ax) δ(ξ−a2π)+δ(ξ+a2π)2 2πδ(ω−a)+δ(ω+a)2 π(δ(ω−a)+δ(ω+a)) Descompunerea cosinusului în două exponențiale complexe.
305 sin⁡(ax) δ(ξ−a2π)−δ(ξ+a2π)2i 2πδ(ω−a)−δ(ω+a)2i iπ(δ(ω+a)−δ(ω−a)) Descompunerea sinusului în două exponențiale complexe.
306 cos⁡(ax2) πacos⁡(π2ξ2a−π4) 12acos⁡(ω24a−π4) πacos⁡(ω24a−π4) a>0. Formula rezultă din combinarea gaussianelor oscilante, interpretate distribuțional.
307 sin⁡(ax2) −πasin⁡(π2ξ2a−π4) −12asin⁡(ω24a−π4) −πasin⁡(ω24a−π4) a>0. Formula rezultă din partea imaginară a transformatei unei gaussiene oscilante.
308 e−πiαx2 1|α|e−iπ4sgn⁡(α)eiπξ2α 12π|α|e−iπ4sgn⁡(α)eiω24πα 1|α|e−iπ4sgn⁡(α)eiω24πα α∈ℝ∖{0}. Factorul de fază depinde de semnul lui α. Cazul cu amortizare complexă se obține din perechea gaussiană. Limita oscilantă este înțeleasă în 𝒮′.
309 xn (i2π)nδ(n)(ξ) in2πδ(n)(ω) 2πinδ(n)(ω) n∈ℕ0, iar δ(n) este derivata distribuțională de ordin n. Prin liniaritate se transformă orice polinom.
310 δ(n)(x) (i2πξ)n (iω)n2π (iω)n Perechea duală a monomului, cu aceeași semnificație pentru δ(n).
311 pv⁡1x −iπsgn⁡(ξ) −iπ2sgn⁡(ω) −iπsgn⁡(ω) pv⁡(1/x) este valoarea principală Cauchy. Formula intervine în transformarea Hilbert.
312 Pf⁡1xn:=(−1)n−1(n−1)!dndxnlog⁡|x| −iπ(−i2πξ)n−1(n−1)!sgn⁡(ξ) −iπ2(−iω)n−1(n−1)!sgn⁡(ω) −iπ(−iω)n−1(n−1)!sgn⁡(ω) Partea finită este definită prin derivarea distribuțională a lui log⁡|x|, conform expresiei din prima coloană. Pentru n=1 se regăsește valoarea principală.
313 |x|α −2sin⁡(πα2)Γ(α+1)|2πξ|α+1 −22πsin⁡(πα2)Γ(α+1)|ω|α+1 −2sin⁡(πα2)Γ(α+1)|ω|α+1 Formula ca funcție local integrabilă în variabila frecvență este dată pentru −1<α<0. Familia |x|α, inițial definită pentru Re⁡α>−1, are o continuare meromorfă distribuțională, cu poli la −1,−3,…. Pentru alți exponenți este necesară regularizarea și, în anumite cazuri, apar termeni concentrați în origine.
313a 1|x| 1|ξ| 1|ω| 2π|ω| Cazul α=−1/2 al perechii precedente.
314 sgn⁡(x) pv⁡1iπξ 2πpv⁡1iω pv⁡2iω sgn este funcția semn. Câturile din transformate sunt valori principale.
315 u(x) 12(pv⁡1iπξ+δ(ξ)) π2(pv⁡1iπω+δ(ω)) π(pv⁡1iπω+δ(ω)) Funcția treaptă Heaviside. Câturile se interpretează prin valoare principală. Valoarea izolată a lui u(0) nu modifică distribuția.
316 ∑n=−∞∞δ(x−nT) 1T∑k=−∞∞δ(ξ−kT) 2πT∑k=−∞∞δ(ω−2πkT) 2πT∑k=−∞∞δ(ω−2πkT) T>0. Pieptenele Dirac. Identitatea este distribuțională, ca în formula de sumare Poisson.
317 J0(x) 2rect⁡(πξ)1−4π2ξ2 2πrect⁡(ω2)1−ω2 2rect⁡(ω2)1−ω2 J0 este funcția Bessel de prima speță de ordin zero. Expresiile sunt nule în afara intervalului spectral indicat și au singularități integrabile la extremități.
318 Jn(x) 2(−i)nTn(2πξ)rect⁡(πξ)1−4π2ξ2 2π(−i)nTn(ω)rect⁡(ω2)1−ω2 2(−i)nTn(ω)rect⁡(ω2)1−ω2 Jn este funcția Bessel de prima speță, iar Tn este polinomul Cebîșev de prima speță, cu n∈ℕ0. Aceeași interpretare la extremități ca în rândul precedent.
319 log⁡|x| −12P(ξ)−(γ+log⁡(2π))δ(ξ) −π2P(ω)−2πγδ(ω) −πP(ω)−2πγδ(ω) γ este constanta Euler–Mascheroni. Partea finită P este definită după tabel, cu aceeași scară de regularizare în toate coloanele. Schimbarea scării modifică termenul Dirac.
320 (∓ix)−α (2π)αΓ(α)u(±ξ)(±ξ)α−1 2πΓ(α)u(±ω)(±ω)α−1 2πΓ(α)u(±ω)(±ω)α−1 0<α<1, cu puterea complexă în ramura principală. u este treapta Heaviside. Derivarea distribuțională furnizează extensii la alți exponenți, cu tratarea termenilor singulari.

În rândul 319, partea finită Hadamard, definită cu scara de referință egală cu 1, este distribuția temperată P, a cărei acțiune asupra unei funcții-test φ∈𝒮(ℝ) este dată de

⟨P,φ⟩=∫|s|<1φ(s)−φ(0)|s|ds+∫|s|≥1φ(s)|s|ds.

Scăderea valorii φ(0) în prima integrală înlătură singularitatea neintegrabilă din origine, permițând definirea distribuției asociate formal expresiei 1|s|. În convenția transformatei Fourier exprimată prin frecvența ξ, factorul log⁡(2π) apare ca urmare a schimbării scării în definiția părții finite. Alegerea unei alte scări de regularizare modifică distribuția printr-un multiplu al distribuției Dirac concentrate în origine.[133]

Funcții cu pătratul integrabil de două variabile

[modificare | modificare sursă]
Nr.FuncțieTransformată unitară în frecvență obișnuită ξTransformată unitară în frecvență unghiulară ωTransformată neunitară în frecvență unghiulară ωCondiții și observații
400 f(x,y) f^(ξx,ξy)≜∬f(x,y)e−i2π(ξxx+ξyy)dxdy f^(ωx,ωy)≜12π∬f(x,y)e−i(ωxx+ωyy)dxdy f^(ωx,ωy)≜∬f(x,y)e−i(ωxx+ωyy)dxdy Integralele sunt pe ℝ2. Toate variabilele de frecvență sunt reale.
401 e−π(a2x2+b2y2) 1|ab|e−π(ξx2a2+ξy2b2) 12π|ab|e−14π(ωx2a2+ωy2b2) 1|ab|e−14π(ωx2a2+ωy2b2) a,b∈ℝ∖{0}. Pereche gaussiană bidimensională, fără impunerea unei integrale egale cu unu.
402 circ⁡(x2+y2) J1(2πξx2+ξy2)ξx2+ξy2 J1(ωx2+ωy2)ωx2+ωy2 2πJ1(ωx2+ωy2)ωx2+ωy2 circ⁡(r)=1 pentru 0≤r≤1 și zero în rest. Formula descrie amplitudinea de difracție asociată discului Airy. J1 este funcția Bessel de prima speță de ordin unu. În origine se ia limita.[134]
403 1x2+y2 1ξx2+ξy2 1ωx2+ωy2 2πωx2+ωy2 Pereche radială în sensul distribuțiilor temperate, legată de transformata Hankel a lui r−1.[135]
404 ix+iy 1ξx+iξy 1ωx+iωy 2πωx+iωy Identitate în sensul distribuțiilor temperate. Singularitățile de ordin 1/r sunt local integrabile în două dimensiuni.

Funcții pe spațiul euclidian de dimensiune arbitrară

[modificare | modificare sursă]
Nr.FuncțieTransformată unitară în frecvență obișnuită ξTransformată unitară în frecvență unghiulară ωTransformată neunitară în frecvență unghiulară ωCondiții și observații
500 f(𝐱) f^1(𝝃)≜∫ℝnf(𝐱)e−i2π𝝃⋅𝐱d𝐱 f^2(𝝎)≜1(2π)n2∫ℝnf(𝐱)e−i𝝎⋅𝐱d𝐱 f^3(𝝎)≜∫ℝnf(𝐱)e−i𝝎⋅𝐱d𝐱 𝐱,𝝃,𝝎∈ℝn.
501 χ[0,1](|𝐱|)(1−|𝐱|2)δ Γ(δ+1)πδ|𝝃|γJγ(2π|𝝃|) 2δΓ(δ+1)|𝝎|γJγ(|𝝎|) Γ(δ+1)πδ|𝝎2π|−γJγ(|𝝎|) δ>−1 și γ=n/2+δ. χ[0,1] este funcția indicatoare, Γ funcția gamma și Jγ funcția Bessel de prima speță. Cazul n=2,δ=0 dă perechea discului.[136]
502 |𝐱|−α,0<Re⁡α<n (2π)αcn,α|𝝃|−(n−α) (2π)n2cn,α|𝝎|−(n−α) (2π)ncn,α|𝝎|−(n−α) 0<Re⁡α<n. Constanta cn,α este dată după tabel. Relația este legată de potențialele Riesz. Continuarea analitică la α∉{n,n+2,n+4,…} necesită distribuții regularizate.[137]
503 1|𝝈|(2π)n2e−12𝐱T𝝈−T𝝈−1𝐱 e−2π2𝝃T𝝈𝝈T𝝃 (2π)−n2e−12𝝎T𝝈𝝈T𝝎 e−12𝝎T𝝈𝝈T𝝎 𝝈 este o matrice reală inversabilă, iar |𝝈|=|det⁡𝝈|. Densitate normală multivariată, de medie zero și covarianță 𝜮=𝝈𝝈T, normalizată la integrală unu. Inversa covarianței este 𝝈−T𝝈−1.
504 e−2πα|𝐱| cnα(α2+|𝝃|2)n+12 cn(2π)n+22α(4π2α2+|𝝎|2)n+12 cn(2π)n+1α(4π2α2+|𝝎|2)n+12 Re⁡α>0 și cn=Γ((n+1)/2)/π(n+1)/2.[138]

Constanta cn,α din formula corespunzătoare rândului 502 este dată de

cn,α=2απn/2Γ(α2)Γ(n−α2), unde Γ este funcția gamma a lui Euler.

Pentru 0<Re⁡α<n, funcțiile omogene |x|−α și |ξ|α−n sunt local integrabile și definesc distribuții temperate. În această situație, identitatea dintre transformatele lor Fourier este valabilă în sensul distribuțiilor, fără ca integrala Fourier să fie neapărat absolut convergentă. În afara intervalului indicat, extinderea formulei necesită interpretarea corespunzătoare în teoria distribuțiilor, prin continuare analitică sau regularizare, după caz.[137]

Funcția dreptunghiulară considerată în tabel este egală cu 1 pentru |x|<12 și cu 0 pentru |x|>12. Valorile atribuite în punctele x=±12 nu influențează transformata Fourier, deoarece modificarea unei funcții într-un număr finit de puncte nu schimbă clasa sa de echivalență în L2(ℝ) și nici valoarea integralei Lebesgue. Funcția sinus cardinal normalizată este definită prin sinc⁡(x)=sin⁡(πx)πx pentru x≠0, iar în origine se prelungește prin continuitate, cu sinc⁡(0)=1. Această definiție trebuie deosebită de forma nenormalizată sin⁡xx, întâlnită în unele lucrări de analiză matematică. Cele două convenții diferă prin schimbarea scării argumentului și conduc la formule Fourier corespunzătoare.

  1. ↑ Distincția rezultă adesea din context. Aplicarea transformării unei funcții desemnează operația, iar transformata unei funcții desemnează rezultatul. În literatura românească se întâlnește și distincția explicită dintre transformarea Fourier ca operator și transformata Fourier ca imagine a unei funcții. Conform Olteanu, p. 17.
  2. ↑ În mecanica cuantică relativistă apar transformate ale funcțiilor de undă cu mai multe componente. În teoria cuantică a câmpurilor sunt utilizate transformate ale unor funcții operatoriale de spațiu–timp. Greiner & Reinhardt 1996.
  3. ↑ Deplasarea frecvenței are forma ℱ{f(x)e−2πiξ0x}(ξ)=f^(ξ+ξ0). Valoarea la ξ=0 este f^(ξ0), ceea ce exprimă deplasarea componentei ξ0 la zero. Semnul din argumentul transformatei nu indică o deplasare fizică în sens opus.
  1. ↑ Vretblad 2003.
  2. ↑ Pinsky 2002, p. 91.
  3. ↑ Stade 2005, pp. 298–299.
  4. ↑ Holhoș, §10.1, teorema 10.18.
  5. ↑ Rahman 2011, p. 10.
  6. ↑ Oppenheim, Schafer & Buck 1999, p. 58.
  7. ↑ Gelfand & Shilov 1968, p. 128.
  8. ↑ Lieb & Loss 2001, pp. 123–125.
  9. ↑ Titchmarsh 1986.
  10. ↑ Dym & McKean 1985.
  11. 1 2 3 4 5 Howe 1980.
  12. ↑ Folland 1989.
  13. ↑ Tong, §2.4.1, ecuațiile 2.134–2.141.
  14. ↑ Holhoș, p. 1.
  15. ↑ Arfken 1985.
  16. ↑ Grosu, p. 4.
  17. ↑ Lupșa 2008, p. 64.
  18. ↑ Smith 2007.
  19. ↑ Olteanu, pp. 16–18.
  20. ↑ Miller & Piatetski-Shapiro 1998, §1.4, propoziția 1.2.
  21. 1 2 3 4 5 6 7 Pinsky 2002.
  22. ↑ Holhoș, § 10.1.1.
  23. 1 2 Proakis & Manolakis 1996, p. 291.
  24. ↑ Katznelson 2004, p. 153.
  25. 1 2 Stein & Weiss 1971, pp. 1–2.
  26. 1 2 3 4 5 6 7 8 9 10 Stein & Weiss 1971.
  27. 1 2 Rudin 1987, pp. 186–187.
  28. 1 2 3 Stein & Weiss 1971, teorema 2.3.
  29. ↑ Holhoș, teorema 10.21.
  30. 1 2 Proakis & Manolakis 1996.
  31. ↑ Hunter 2014, pp. 78–80.
  32. ↑ Holhoș, teorema 10.15.
  33. 1 2 Hunter 2014, teorema 5.66.
  34. ↑ Holhoș, teoremele 10.13–10.14.
  35. 1 2 Folland 1992, p. 216.
  36. 1 2 Wolf 1979, pp. 307–312.
  37. ↑ Folland 1989, p. 53.
  38. 1 2 3 Celeghini, Gadella & del Olmo 2021, § 2.
  39. 1 2 3 4 Duoandikoetxea 2001.
  40. ↑ Celeghini, Gadella & del Olmo 2021, § 2, ecuațiile 2.15–2.20.
  41. ↑ Wolf 1979, pp. 312, 320.
  42. ↑ Folland 1989, p. 52.
  43. ↑ Boashash 2003.
  44. ↑ Condon 1937.
  45. 1 2 3 Hunter 2014, apendicele 5.B.
  46. ↑ Wolf 1979, p. 320.
  47. 1 2 Woit 2023, § 4.1 și § 4.3.
  48. 1 2 Woit 2023, § 3.1–3.3.
  49. ↑ Woit 2023, § 3.3 și § 4.1.
  50. 1 2 Paley & Wiener 1934.
  51. ↑ Gelfand & Vilenkin 1964.
  52. ↑ Kirillov & Gvishiani 1982.
  53. ↑ Clozel & Delorme 1985, pp. 331–333.
  54. ↑ de Groot & Mazur 1984, p. 146.
  55. ↑ Champeney 1987, p. 80.
  56. 1 2 3 4 Kolmogorov & Fomin 1999.
  57. ↑ Wiener 1949.
  58. ↑ Champeney 1987, p. 63.
  59. 1 2 3 Pinsky 2002, § 2.4.3.
  60. ↑ Wigderson & Wigderson 2021, § 4.3, teorema 4.9.
  61. 1 2 Strohmer.
  62. ↑ Stein & Shakarchi 2003, § 5.4.
  63. ↑ Białynicki-Birula & Mycielski 1975, pp. 129–131.
  64. ↑ Chatfield 2004, p. 113.
  65. ↑ Fourier 1822, p. 441.
  66. ↑ Poincaré 1895, p. 102.
  67. ↑ Whittaker & Watson 1927, p. 188.
  68. ↑ Grafakos 2004.
  69. ↑ Grafakos & Teschl 2013, teorema 1.1.
  70. ↑ Duoandikoetxea 2001, teorema 8.3.
  71. ↑ Fefferman 1971, pp. 330–336.
  72. ↑ Rudin 1987, pp. 182–183.
  73. ↑ Chandrasekharan 1989, pp. 7–8, 84.
  74. ↑ Candès 2016, prelegerea 3.
  75. 1 2 3 4 Katznelson 2004.
  76. ↑ Mallat 2009, p. 45.
  77. ↑ Strichartz 1994, p. 150.
  78. ↑ Hunter 2014, apendicele 5.A și teorema 5.66.
  79. ↑ Hunter 2014, definiția 5.68.
  80. ↑ Strichartz 1994.
  81. ↑ Salo 2024, teoremele 3.9.1 și 3.9.9.
  82. ↑ Pinsky 2002, p. 256.
  83. ↑ Rudin 1991, p. 15.
  84. ↑ Edwards 1982, pp. 53, 67, 72–73.
  85. ↑ Edwards 1982, p. 53.
  86. ↑ Katznelson 2004, pp. 40, 155, 164.
  87. ↑ Katznelson 2004, p. 173.
  88. ↑ Billingsley 1995, p. 345.
  89. 1 2 3 4 5 6 Rudin 1991.
  90. ↑ Bell 2014, secțiunile despre funcții pozitiv definite și teorema lui Bochner.
  91. ↑ Bell 2014, secțiunile despre transformata Fourier și măsura Haar pe grupul dual.
  92. ↑ Hewitt & Ross 1970, capitolul despre reprezentări și grupuri compacte.
  93. 1 2 3 Knapp 2001.
  94. ↑ Shapiro 2024, teorema 8.36.
  95. 1 2 3 4 5 Hewitt & Ross 1970, capitolul VIII.
  96. ↑ Hewitt & Ross 1970.
  97. 1 2 3 Hunter 2014, § 7.1.
  98. 1 2 3 Precup 2004, § 5.3.
  99. ↑ Hunter 2014, §§ 5.A.2–5.A.3.
  100. ↑ Ablowitz et al. 1974, pp. 249–315.
  101. ↑ Lax 1968, pp. 467–490.
  102. ↑ Yousefi & Kschischang 2014, pp. 4312–4328.
  103. ↑ Keeler 2010, §§ 5.1–5.3.
  104. ↑ Griffiths & de Haseth 2007, §§ 2.1–2.2.
  105. ↑ Griffiths & de Haseth 2007, pp. 5, 12–13.
  106. ↑ Gross 2017, § 4.7.10.
  107. ↑ Brown et al. 2014, pp. 7–8.
  108. ↑ Beu 2020, „Mecanica ondulatorie”, diapozitivele 9–16.
  109. ↑ Beu 2020, „Mecanica ondulatorie”, diapozitivul 32.
  110. ↑ Beu 2020, „Mecanica ondulatorie”, diapozitivele 13–16.
  111. ↑ Beu 2020, „Formalismul mecanicii ondulatorii”, diapozitivele 5–6.
  112. ↑ Beu 2020, „Formalismul mecanicii ondulatorii”, diapozitivele 32–34.
  113. ↑ Beu 2020, „Mecanica ondulatorie”, diapozitivele 18–21 și 32.
  114. ↑ Beu 2020, „Mecanica ondulatorie”, diapozitivele 25–28.
  115. ↑ Beu 2020, „Formalismul mecanicii ondulatorii”, diapozitivele 25–26.
  116. ↑ Tong, §§ 2.1–2.2.
  117. ↑ Wolf 1979, p. 312.
  118. ↑ Grosu, § 9.1, definiția 9.1.5 și teorema 9.1.14.
  119. ↑ Shalizi 2006, §§ 21.1–21.2.
  120. ↑ Shalizi 2006, § 21.1, teorema 240.
  121. ↑ Shalizi 2006, § 21.1, teoremele 240 și 244.
  122. ↑ Shalizi 2006, § 21.1.
  123. 1 2 Press et al. 1992, § 13.9.
  124. ↑ Gradshteyn et al. 2015.
  125. ↑ „fourier — Fourier transform of symbolic expression or function”. MATLAB Documentation. MathWorks. Accesat în .
  126. ↑ Bailey & Swarztrauber 1994.
  127. ↑ Lado 1971.
  128. ↑ Simonen & Olkkonen 1985.
  129. ↑ Erdélyi 1954.
  130. ↑ Kammler 2000, anexa cu tabele.
  131. ↑ NIST, § 1.16(viii).
  132. ↑ Campbell & Foster 1948.
  133. ↑ Gelfand & Shilov 1964.
  134. ↑ Stein & Weiss 1971, teorema IV.3.3.
  135. ↑ Easton 2010.
  136. ↑ Stein & Weiss 1971, teorema 4.15.
  137. 1 2 Gelfand & Shilov 1964, p. 363.
  138. ↑ Stein & Weiss 1971, p. 6.

Legături externe

[modificare | modificare sursă]
Commons
Commons
Wikimedia Commons conține materiale multimedia legate de transformate Fourier