Anàlisi de Fourier

branca de les matemàtiques que estudia les senyals periòdiques i contínues

En matemàtiques, l'anàlisi de Fourier (/ˈfʊrieɪ, -iər/) és l'estudi de la forma com una funció general es pot representar o aproximar a partir de sumes o de funcions trigonomètriques més simples. L'anàlisi de Fourier va sorgir de l'estudi de les sèries de Fourier, i du el nom de Joseph Fourier, que va demostrar que representar una funció com a sumatori de funcions trigonomètriques simplifica en gran manera l'estudi de la transmissió tèrmica.

Transformades de Fourier
Transformada de Fourier continua
Sèrie de Fourier
Transformada Discreta de Fourier
Transformada de Fourier en Temps Discret
Transformada de Fourier sobre cossos finits
Anàlisi de Fourier
Transformades relacionades
Senyal temporal d'un baix elèctric de la corda La (55 Hz).
Transformada de Fourier de la senyal temporal d'un ubaix elèctric de la corda A (55 Hz).L'anàlisi de Fourier revela els components oscil·latoris de senyals i funcions.

Avui en dia, el tema de l'anàlisi de Fourier engloba una vast espectre de les matemàtiques. En les ciències i en enginyeria, s'anomena sovint anàlisi de Fourier al procés de descompondre una funció en components oscil·latoris, mentre que l'operació de reconstruir una funció a partir d'aquestes peces és conegut com síntesi de Fourier. Per exemple, determinar quin component freqüencials són presents en una nota musical implica calcular la transformada de Fourier d'una nota musical. Es pot, doncs, resintetitzar el mateix so incloent els components freqüencials com es revelen en l'anàlisi de Fourier. En matemàtiques, el terme anàlisi de Fourier sovint fa referència a l'estudi de totes dues operacions.

El procés de descomposició en si s'anomena transformació de Fourier. Al seu resultat, la transformada de Fourier, sovint se li dona un nom més específic, que depèn del domini i d'altres propietats de la funció que es transforma. A més, el concepte original de l'anàlisi de Fourier s'ha estès al llarg del temps per ser aplicat a situacions cada vegada més abstractes i generals, i el camp general és sovint anomenat anàlisi harmònica. Cada transformació utilitzada en l'ànalisi té la seva corresponent transformada inversa que pot ser utilitzada en la síntesi.

Aplicacions modifica

L'anàlisi de Fourier té moltes aplicacions científiques – en física, en equacions diferencials en derivades parcials, en teoria de nombres, en combinatòria, en processament de senyals, en processament digital d'imatges, en teoria de la probabilitat, en estadística, en ciències forenses, en valoració d'opcions, en criptografia, en anàlisi numèrica, en acústica, en oceanografia, en sonars, en òptica, en difracció, en geometria, en anàlisi de l'estructura de les proteïnes, entre d'altres.

Aquesta àmplia aplciabilitat rau en moltes propietats útils de les transformades:

En ciències forenses, els espectrofotòmetres infra-rojos dels laboratoris utilitzen l'anàlisi de la transformada de Fourier per mesurar les longituds d'ona de la llum dins de l'espectre infra-roig que un material absorbeix. S'utilitza el mètode FT per descodificar les senyals mesurades i obtenir les dades de longitud d'ona. I, utilitzant un ordinar, aquèsts càlculs de Fourier es fan ràpidament; així, en qüestió de segons, un instrument operat amb un ordinador pot produir un patró d'absorció infra-roja comparable al de l'instrument.[5]

La transformada de Fourier també és útil com a representació compacta d'una senyal. Per exemple, la compressió JPEG utilitza una variant de la transformada de Fourier (transformada cosinus discreta) de trossos quadrats petits de la imatge digital. Els components de Fourier de cada quadrat són arrodonits a una precisió més baixa, i els components baixos s'eliminen completament, de tal manera que els components restants es poden emmagatzemar més compactament. En la reconstrucció de la imatge, cada imatge és reagrupada a partir de les components aproximades restants de la transformada de Fourier, que se'ls aplica la transformada inversa per produir l'aprosimació de la imatge original.

Aplicacions en processament de senyals modifica

En el processament de senyals, com ara d'àudio, d'ones de ràdio, ones de llum, ones sísmiques, o fins i tot imatges, l'ànalisi de Fourier pot aïllar components d'amplada de banda estreta d'una ona composta, i detectar-les o eliminar-les més fàcilment. Una àmplia família de tècniques de processament de senyals estan basades en fer la transformada de Fourier de la senyal, manipular-la posteriorment de forma simple i finalment invertint la transformada.[6]

Alguns d'aquests exemples són:

Variants de l'anàlisi de Fourier modifica

 
Una transformada de Fourier i 3 variacions obtingudes amb mostratge periòdic (amb interval T) i/o sumatori periòdic (amb interval P) d'una funció en el domini temporal. La facilitat relativa de la computació de la DFT (transformada discreta de Fourier) així com la informació que proporciona en S( f ) fan que sigui una eina popular en l'anàlisi.

Transformada (contínua) de Fourier modifica

Molt sovint, el terme transformada de Fourier fa referència a la transformada de funcions en el domini continu dels nombres reals, que produeix una funció contínua en freqüència, coneguda com la distribució freqüencial. Una funció és transformada en una altra, i l'operador és reversible. Quan el domini de la funció d'entrada (input) és el temps (t), i el domini de la funció de sortida (output) és la freqüència ordinària, la transformada de al funció s(t) a la freqüència f ve donada pel nombre complex:

 

Avaluant aquesta quantitat pels valor de f produeix la funció en el domini freqüencial. Llavors es pot representar s(t) com una recombinació d'exponencials complexes de totes les freqüències possibles:

 

que és la fórmula de la transformada inversa. El nombre complex, S( f ), transmet tant l'amplitud com la fase a la freqüència f.

Vegeu transformada de Fourier per més informació.

Sèrie de Fourier modifica

La transformada de Fourier d'una funció periòdica, sP(t), amb període P, esdevé una funció pinta de Dirac, modulada per una seqüència de coeficients complexes:

      (on P és la integral en l'interval de longitud P).

La transformada inversa, anomenada sèrie de Fourier, és una representació de sP(t) en termes d'un sumatori d'un nombre potencialment infinit de sinusoïdes harmònicament relacionades o de funcions exponencials complexes, cadascuna d'elles amb la seva amplitud i fase, especificades per un dels coeficients:

 

Es pot expressar qualsevol sP(t) com un sumatori periòdic d'una altra funció, s(t):

 

i els coeficients són proporcionals a mostres de S( f ) en intervals discrets de 1/P:

 [A]

Noti's que en el sumatori periòdic es pot utilitzar qualsevol s(t) la transformada de la qual tingui els mateixos valors discrets de mostreig. Una condició suficient per recuperar s(t) (i, per tant S( f )) a partir d'aquestes mostres només (és a dir, a partir de la sèrie de Fourier) és que la porció no-zero de s(t) estigui confinada en un interval conegut de duració P, que és el domini freqüencial dual del teorema de mostratge de Nyquist-Shannon.

Transformada de Fourier de senyal discret modifica

La DTFT és el dual matemàtic de la sèrie de Fourier en el domini temporal. Així, es pot representar com a sèrie de Fourier un sumatori periòdic convergent en el domini freqüencial, que tindrà com a coeficients mostres d'una funció relacionada contínua en el temps:

 

que és coneguda com DTFT. Llavors la DTFT de la seqüència s[n] és també la transformada de Fourier de la funció pinta de Dirac modulada.[B]

Els coeficients de la sèrie de Fourier (i la transformada inversa) estan definits com:ç

 

El paràmetre T correspon a l'interval de mostratge, i ara es pot reconèixer aquesta sèrie de Fourier com una forma de fórmula del sumatori de Poisson. El resultat que segueix és que quan una seqüència discreta de dades, s[n], és proporcional a mostres d'una funció contínua subjacent, s(t), es pot observar un sumatori periòdic de la transformada contínua de Fourier, S( f ). Noti's que qualsevol s(t) amb els mateixos valors de mostratge discret produeix la mateixa DTFT  Però sota una certes condicions es pot, teòricament, recuperar S( f ) i s(t) de forma exacta. Una condició suficient per a la recuperació perfecta és que la porció (no-nul·la) de S( f ) estigui confinada en un interval freqüencial d'amplada 1/T. Quan l'interval és [−1/2T, 1/2T], la fórmula aplicable de reconstrucció és la fórmula d'interpolació de Whittaker–Shannon. Aquesta és la pedra angular del processament de senyals digitals.

Una altra raó d'interès de S1/T( f ) és que sovint proporciona informació sobre la quantitat d'aliàsing causat pel procés de mostratge.

Transformada discreta de Fourier modifica

De forma similar a la sèrie de Fourier, la DTFT d'una seqüència periòdica, sN[n], amb període N, es converteix en una funció pinta de Dirac, modulada per una seqüència de coeficients complexos:

      (on n és el sumatori al llarg de qualsevol seqüència de longitud N).

La seqüència S[k] és el que es coneix normalment com la DFT d'un cicle de sN. És també N-periòdica, així doncs mai no cal calcula més de N coeficients (1 cicle). La transformada inversa, també coneguda com a sèrie discreta de Fourier, ve donada per:

    on k és el sumatori de qualsevol seqüència de longitud N.

Quan s'expressa sN[n] com a sumatori periòdic d'una altra funció:

    and    [C]

els coeficients són proporcionals a les mostres de S1/T( f ) en intervals discrets de 1/P = 1/NT:

 [D]

En canvi, quan es vol calcular un nombre arbitrari (N) de mostres discretes d'un cicle d'una DTFT contínua, S1/T( f ), es pot fer calculant la (relativament simple) DFT de sN[n], com s'ha definit més amunut. En la majoria dels casos, es tria N igual a la longitud de la porció (no nul·la) de s[n]. Augmentar N, mètode conegut com zero-padding o interpolació, fa que les mostres estiguin més espaiades en un cicle de S1/T( f ). Disminuir N, causa un solapament (suma) en el domini temporal (anàleg a l'aliàsing, solapament en freqüència), que correspon al submostratge en el domini freqüencial. En la majoria de casos que tenen interès pràctic, la seqüència s[n] representa una seqüència més llarga que ha estat truncada per l'aplicació d'una funció finestra de longitud finita o d'un filtre FIR.

Es pot calcular la DFT utilitzant un algorisme de transformada ràpida de Fourier (FFT, per les seves sigles en anglès), que fa que sigui una transformació pràctica i important en computació.

Resum modifica

Per a funcions periòdiques, tant la transformada de Fourier com la DTFT consten d'un conjunt discret de components freqüencials (sèrie de Fourier), i les transformades divergeixen en aquestes freqüències. Una pràctica habitual és resoldre aquesta divergència amb funcions delta de Dirac i pinta de Dirac. No obstant això, a partir de només un cicle de la funció periòdica es pot extreure la mateixa informació espectral, ja que en tots els cicles la funció és idèntica. Similarment, es poden representar en forma de sèrie de Fourier funcions que tenen una duració finita, sense pèrdua real d'informació més enllà del fet que la periodicitat de la transformada inversa és un mer constructe.

En la pràctica, és habitual que la durada de s(•) estigui limitada al període, P o N. Però per a les fórmules que es mostren a continuació no es requereix aquesta condició.

Transformades de s(t) (temps continu)
Freqüència contínua Freqüència discreta
Transforma    
Inversa    
Transformades de s(nT) (temps discret)
Freqüència contínua Freqüència discreta
Transformada  

 

Inversa  

 

 

Història modifica

Una forma primerenca de sèrie harmònica va aparèixer en les matemàtiques de Babilònia, en què eren usades per calcular efemèrides (taules de posicions astronòmiques).[7][8][9][10]

El conceptes grecs clàssics de deferents i epicicles en el sistema ptolemaic d'astronomia estaen relacionades amb les sèries de Fourier.

En l'època moderna, Alexis Clairaut va utilitzar variants de la transformada discreta de Fourier l'any 1754 per calcular una òrbita,[11] en el que s'ha considerat l'ha primera fórmula de la DFT,[12] i l'any 1759 la va utilitzar Joseph Louis Lagrange per calcular els coeficients de sèries trigonomètriques d'una corda en vibració.[12] Tècnicament, l'obra de Clairaut només es basava en sèries de cosinus (una espècie de transformada cosinus discreta), mentre que l'obra de Lagrange es basava en sèries de només sinus (una forma de transformada sinus discreta). L'any 1805, Gauß va utilitzar una sèrie de sinus i cosinus en una interpolació trigonomètrica d'òrbites d'asteroides.[13] Euler i Lagrange van discretitzar el domini en el problema de la corda en vibració, usant el que avui s'anomenaria mostratge.[12]

Un primer desenvolupament modern cap a l'anàlisi de Fourier va ser l'articles de 1770 Réflexions sur la résolution algébrique des équations (Reflexions sobre la resolució algebraica d'equacions) de Lagrange, que en el mètodes dels resolvents de Lagrange utilitzava una descomposició complexa de Fourier per estudiar la solució d'una cúbica.[14] Lagrange va transformada les arrels x1, x₂, x en els resolvents:

 

on ζ és una arrel cúbica de la unitat, que és la DFT d'ordre 3.

Uns quants autors, notablement Jean le Rond d'Alembert i Carl Friedrich Gauß van utilitzar sèries trigonomètriques per estudiar l'equació de la calor,[15] però el desenvolupament més revolucionari va ser l'article de 1807 Mémoire sur la propagation de la chaleur dans les corps solides (Disseració sobre la propagació de la calor en cossos sòlids) de Joseph Fourier, que va tenir com a principal aportació el modelatge de totes les funcions en sèries trigonomètriques i va introduir les sèries de Fourier.

Els historiadors de les matemàtiques discuteixen fins a quin punt Lagrange i altres matemàics van fer aportacions clau en el desenvolupament de la teoria de Fourier: Daniel Bernoulli i Leonhard Euler havien introduït representacions trigonomètriques de funcions, i Lagrange va donar la oslució en sèries de Fourier de l'equació d'ones. Així doncs, la contribució de Fourier va ser principalment afirmar que es pot representar una funció arbitrària en forma de sèrie de Fourier.[12]

El desenvolupament posterior del camp és conegut com anàlisi harmònica, i és també un dels primers exemples de teoria de la representació.

El primer algorisme de transformada ràpida de Fourier (FFT per les seves sigles en anglès) per la DFT va ser creat al voltant de l'any 1805 per Carl Friedrich Gauß per interpolar mesures de les òrbites dels asteroides (3) Juno i (2) Pal·les, tot i que aquest algorisme d'FFT en particular és sovint atribuït als seus redescobridors moderns Cooley i Tukey.[13][11]

Vegeu també modifica

Notes modifica

  1.  
  2. Noti's també que:
     
    En conseqüència, una pràctica habitual és modelar el "mostratge" com la multiplació per una funció pinta de Dirac, que evidentment només és "possible" en un sentit purament matemàtic.
  3. Noti's que aquesta definició difereix -intncionadament- de la de la secció de la DTFT en un factor T. Això facilita la taula de les "Transformades  ". Alternativament, es pot definir   com  . Llavors  
  4.  

Referències modifica

  1. Rudin, Walter. Fourier Analysis on Groups. Wiley-Interscience, 1990. ISBN 978-0-471-52364-2. 
  2. Evans, L. Partial Differential Equations. American Mathematical Society, 1998. ISBN 978-3-540-76124-2. 
  3. Knuth, Donald E. The Art of Computer Programming Volume 2: Seminumerical Algorithms. 3rd. Addison-Wesley Professional, 1997. ISBN 978-0-201-89684-8. 
  4. Conte, S. D.; de Boor, Carl. Elementary Numerical Analysis. Third. Nova York: McGraw Hill, Inc., 1980. ISBN 978-0-07-066228-5. 
  5. Saferstein, Richard. Criminalistics: An Introduction to Forensic Science, 2013. 
  6. Rabiner, Lawrence R.; Gold, Bernard. Theory and Application of Digital Signal Processing, 1975. 
  7. Prestini, Elena. The Evolution of Applied Harmonic Analysis: Models of the Real World. Birkhäuser, 2004, p. 62. ISBN 978-0-8176-4125-2. 
  8. Rota, Gian-Carlo; Palombi, Fabrizio. Indiscrete Thoughts. Birkhäuser, 1997, p. 11. ISBN 978-0-8176-3866-5. 
  9. Neugebauer, Otto. The Exact Sciences in Antiquity. 9. 2nd. Dover Publications, 1969, p. 1–191. ISBN 978-0-486-22332-2. 
  10. Brack-Bernsen, Lis; Brack, Matthias «Analyzing shell structure from Babylonian and modern times». International Journal of Modern Physics E, 13, 1, 2004, pàg. 247. arXiv: physics/0310126. Bibcode: 2004IJMPE..13..247B. DOI: 10.1142/S0218301304002028.
  11. 11,0 11,1 Terras, Audrey. Fourier Analysis on Finite Groups and Applications. Cambridge University Press, 1999, p. 30-32. ISBN 978-0-521-45718-7. 
  12. 12,0 12,1 12,2 12,3 Briggs, William L.; Henson, Van Emden. The DFT: An Owner's Manual for the Discrete Fourier Transform. SIAM, 1995, p. 2–4. ISBN 978-0-89871-342-8. 
  13. 13,0 13,1 Heideman, M.T.; Johnson, D. H.; Burrus, C. S. «Gauss and the history of the fast Fourier transform». IEEE ASSP Magazine, 1, 4, 1984, pàg. 14–21. DOI: 10.1109/MASSP.1984.1162257.
  14. Knapp, Anthony W. Basic Algebra. Springer, 2006, p. 501. ISBN 978-0-8176-3248-9. 
  15. Narasimhan, T.N. «Fourier's heat conduction equation: History, influence, and connections». Reviews of Geophysics, 37, 1, febrer 1999, pàg. 151–172. Bibcode: 1999RvGeo..37..151N. DOI: 10.1029/1998RG900006. ISSN: 1944-9208. OCLC: 5156426043.

Bibliografia complementària modifica

Enllaços externs modifica