Articulo de referencia

Transformadas de seno y coseno

Las transformadas de seno y coseno convierten una función en una representación en el dominio de la frecuencia como una suma de ondas sinusoidales y cosenoidales . La transforma...

Las transformadas de seno y coseno convierten una función en una representación en el dominio de la frecuencia como una suma de ondas sinusoidales y cosenoidales . La transformada inversa la convierte de nuevo al dominio del tiempo o del espacio .

En matemáticas , las transformadas de Fourier de seno y coseno son ecuaciones integrales que descomponen funciones arbitrarias en una suma de ondas sinusoidales que representan la componente impar de la función, más ondas cosenoidales que representan la componente par. La transformada de Fourier moderna, de valores complejos, incluye de forma concisa tanto la transformada de seno como la de coseno. Dado que las transformadas de seno y coseno utilizan ondas sinusoidales y cosenoidales en lugar de exponenciales complejas y no requieren números complejos ni frecuencias negativas , se corresponden más estrechamente con las ecuaciones de transformación originales de Joseph Fourier y aún se prefieren en algunas aplicaciones de procesamiento de señales y estadística , pudiendo ser más adecuadas como introducción al análisis de Fourier .

Definición

Las transformadas de Fourier relacionan una función en el dominio del tiempo (rojo) con una función en el dominio de la frecuencia (azul). Las ondas sinusoidales o cosenoidales que componen la función original aparecerán como picos en las funciones del dominio de la frecuencia generadas por la transformada sinusoidal o cosenoidal, respectivamente.

La transformada sinusoidal de Fourier deF(t){\displaystyle f(t)}es: [ nota 1 ]

transformada sinusoidal de Fourier

F^s(ξ)=F(t)pecado(2πξt)dt.{\displaystyle {\hat {f}}^{s}(\xi )=\int _{-\infty }^{\infty }f(t)\sin(2\pi \xi t)\,dt.}

Sit{\displaystyle t}significa tiempo , entoncesξ{\displaystyle \xi }es la frecuencia en ciclos por unidad de tiempo, [ nota 2 ] pero en abstracto, pueden ser cualquier par dual de variables (por ejemplo, posición y frecuencia espacial ).

La transformada sinusoidal es necesariamente una función impar de la frecuencia, es decir, para todoξ{\displaystyle \xi }:

F^s(ξ)=F^s(ξ).{\displaystyle {\hat {f}}^{s}(-\xi )=-{\hat {f}}^{s}(\xi ).}

La transformada coseno de una función rectangular simple (de altura1a{\displaystyle {\tfrac {1}{a}}}y anchoa{\displaystyle a}) es el sinc normalizado(aξ){\displaystyle (a\xi )}graficado arriba.

La transformada coseno de Fourier deF(t){\displaystyle f(t)}es: [ nota 3 ]

transformada coseno de Fourier

F^do(ξ)=F(t)porque(2πξt)dt.{\displaystyle {\hat {f}}^{c}(\xi )=\int _{-\infty }^{\infty }f(t)\cos(2\pi \xi t)\,dt.}

La transformada del coseno es necesariamente una función par de la frecuencia, es decir, para todoξ{\displaystyle \xi }:

F^do(ξ)=F^do(ξ).{\displaystyle {\hat {f}}^{c}(-\xi )={\hat {f}}^{c}(\xi ).}

Simplificación de números pares e impares

Como todas las funciones pares, la mitad izquierda de una función gaussiana es una imagen especular de su mitad derecha y su transformada seno es completamente 0. Las gaussianas tienen la formamiαt2{\displaystyle e^{-\alpha t^{2}}}y su transformada coseno: F^do(ξ)=παmi(πξ)2α{\displaystyle {\hat {f}}^{c}(\xi )={\sqrt {\tfrac {\pi }{\alpha }}}\,e^{-{\frac {(\pi \xi )^{2}}{\alpha }}}} También es una gaussiana. La gaussiana graficada usa α = π y es su propia transformada de coseno.

Las reglas de multiplicación para funciones pares e impares que se muestran entre llaves en las siguientes ecuaciones simplifican drásticamente los integrandos al transformar funciones pares e impares . Algunos autores [ 1 ] incluso definen únicamente la transformada coseno para funciones pares.Fincluso(t){\displaystyle f_{\text{even}}(t)}. Dado que el coseno es una función par y debido a que la integral de una función par de{\displaystyle {-}\infty }a{\displaystyle \infty }es el doble de su integral de0{\displaystyle 0}a{\displaystyle \infty }La transformada coseno de cualquier función par se puede simplificar para evitar valores negativos.t{\displaystyle t}:

F^do(ξ)=Fincluso(t)porque(2πξt)incluso·incluso=igualdt=20Fincluso(t)porque(2πξt)dt.{\displaystyle {\hat {f}}^{c}(\xi )=\int _{-\infty }^{\infty }\overbrace {f_{\text{even}}(t)\cdot \cos(2\pi \xi t)} ^{\text{even·even=even}}\,dt=2\int _{0}^{\infty }f_{\text{even}}(t)\cos(2\pi \xi t)\,dt.}

Y porque la integral de{\displaystyle {-}\infty }a{\displaystyle \infty }Si la transformada coseno de cualquier función impar es cero , la transformada coseno de cualquier función impar es simplemente cero:

F^do(ξ)=Fextraño(t)porque(2πξt)impar·par=impardt=0.{\displaystyle {\hat {f}}^{c}(\xi )=\int _{-\infty }^{\infty }\overbrace {f_{\text{impar}}(t)\cdot \cos(2\pi \xi t)} ^{\text{impar·par=impar}}\,dt=0.}

Las funciones impares no cambian si se rotan 180 grados alrededor del origen . Su transformada coseno es completamente cero. La función impar anterior contiene dos funciones delta de Dirac desplazadas en el tiempo de tamaño reducido . Su transformada seno es simplementepecado(aξ).{\displaystyle \sin(a\xi ).}Asimismo, la transformada sinusoidal depecado(aξ){\displaystyle \sin(a\xi )}es el gráfico anterior. Por lo tanto, la función de onda sinusoidal y la función delta de Dirac desplazada en el tiempo forman un par de transformación .

De manera similar, debido a que el seno es impar, la transformada del seno de cualquier función imparFextraño(t){\displaystyle f_{\text{impar}}(t)}También simplifica evitar lo negativot{\displaystyle t}:

F^s(ξ)=Fextraño(t)pecado(2πξt)impar·impar=pardt=20Fextraño(t)pecado(2πξt)dt{\displaystyle {\hat {f}}^{s}(\xi )=\int _{-\infty }^{\infty }\overbrace {f_{\text{impar}}(t)\cdot \sin(2\pi \xi t)} ^{\text{impar·impar=par}}\,dt=2\int _{0}^{\infty }f_{\text{impar}}(t)\sin(2\pi \xi t)\,dt}

y la transformada sinusoidal de cualquier función par es simplemente cero:

F^s(ξ)=Fincluso(t)pecado(2πξt)par·impar=impardt=0.{\displaystyle {\hat {f}}^{s}(\xi )=\int _{-\infty }^{\infty }\overbrace {f_{\text{par}}(t)\cdot \sin(2\pi \xi t)} ^{\text{par·impar=impar}}\,dt=0.}

La transformada seno representa la parte impar de una función , mientras que la transformada coseno representa la parte par de una función.

Otras convenciones

Así como la transformada de Fourier toma la forma de diferentes ecuaciones con diferentes factores constantes (ver Transformada de Fourier §  Unitaridad y definición para funciones de cuadrado integrable para discusión), otros autores también definen la transformada del coseno como [ 2 ]F^do(ξ)=2π0F(t)porque(2πξt)dt{\displaystyle {\hat {f}}^{c}(\xi )={\sqrt {\frac {2}{\pi }}}\int _{0}^{\infty }f(t)\cos(2\pi \xi t)\,dt} y la transformada del seno como F^s(ξ)=2π0F(t)pecado(2πξt)dt.{\displaystyle {\hat {f}}^{s}(\xi )={\sqrt {\frac {2}{\pi }}}\int _{0}^{\infty }f(t)\sin(2\pi \xi t)\,dt.}Otra convención define la transformada coseno como [ 3 ]Fdo(α)=2π0F(incógnita)porque(αincógnita)dincógnita{\displaystyle F_{c}(\alpha )={\frac {2}{\pi }}\int _{0}^{\infty }f(x)\cos(\alpha x)\,dx}y la transformada del seno comoFs(α)=2π0F(incógnita)pecado(αincógnita)dincógnita{\displaystyle F_{s}(\alpha )={\frac {2}{\pi }}\int _{0}^{\infty }f(x)\sin(\alpha x)\,dx}usandoα{\displaystyle \alpha }como la variable de transformación. Y mientrast{\displaystyle t}se utiliza normalmente para representar el dominio del tiempo ,incógnita{\displaystyle x}En cambio, a menudo se utiliza para representar un dominio espacial al transformarlo en frecuencias espaciales.

inversión de Fourier

La función originalF{\displaystyle f}puede recuperarse a partir de sus transformadas de seno y coseno bajo las hipótesis habituales [ nota 4 ] utilizando la fórmula de inversión: [ 4 ]

Inversión de Fourier (a partir de las transformadas de seno y coseno)

F(t)=F^s(ξ)pecado(2πξt)dξcomponente extraño de F(t)+F^do(ξ)porque(2πξt)dξincluso componente de F(t).{\displaystyle f(t)=\underbrace {\int _{-\infty }^{\infty }{\hat {f}}^{s}(\xi )\sin(2\pi \xi t)\,d\xi } _{{\text{odd component of }}f(t)}\,+\underbrace {\int _{-\infty }^{\infty }{\hat {f}}^{c}(\xi )\cos(2\pi \xi t)\,d\xi } _{{\text{even component of }}f(t)}\,.}

Simplificaciones

Tenga en cuenta que, dado que ambos integrandos son funciones pares deξ{\displaystyle \xi }El concepto de frecuencia negativa puede evitarse duplicando el resultado de la integración sobre frecuencias no negativas:

F(t)=20F^s(ξ)pecado(2πξt)dξ+20F^do(ξ)porque(2πξt)dξ.{\displaystyle f(t)=2\int _{0}^{\infty }{\hat {f}}^{s}(\xi )\sin(2\pi \xi t)\,d\xi \,+2\int _{0}^{\infty }{\hat {f}}^{c}(\xi )\cos(2\pi \xi t)\,d\xi \,.}

Además, siF{\displaystyle f}Si es una función impar , entonces la transformada del coseno es cero, por lo que su inversión se simplifica a:F(t)=F^s(ξ)pecado(2πξt)dξ, solo si F(t) es extraño.{\displaystyle f(t)=\int _{-\infty }^{\infty }{\hat {f}}^{s}(\xi )\sin(2\pi \xi t)\,d\xi ,{\text{ only if }}f(t){\text{ is odd.}}}

Asimismo, si la función originalF{\displaystyle f}Si es una función par , entonces la transformada del seno es cero, por lo que su inversión también se simplifica a:

F(t)=F^do(ξ)porque(2πξt)dξ, solo si F(t) es par.{\displaystyle f(t)=\int _{-\infty }^{\infty }{\hat {f}}^{c}(\xi )\cos(2\pi \xi t)\,d\xi ,{\text{ only if }}f(t){\text{ is even.}}}

Sorprendentemente, estas dos últimas fórmulas de inversión simplificadas se ven idénticas a las transformadas de seno y coseno originales, respectivamente, aunque cont{\displaystyle t}intercambiado conξ{\displaystyle \xi }(y conF{\displaystyle f}intercambiado conF^s{\displaystyle {\hat {f}}^{s}}oF^do{\displaystyle {\hat {f}}^{c}}Una consecuencia de esta simetría es que sus procesos de inversión y transformación siguen funcionando cuando se intercambian las dos funciones. Dos de estas funciones se denominan pares de transformación . [ nota 5 ]

Descripción general de la demostración de inversión

Utilizando la fórmula de adición para el coseno , la fórmula de inversión completa también se puede reescribir como la fórmula integral de Fourier : [ 5 ] [ 6 ]F(t)=F(incógnita)porque(2πξ(incógnitat))dincógnitadξ.{\displaystyle f(t)=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f(x)\cos(2\pi \xi (x-t))\,dx\,d\xi .} Este teorema se suele enunciar bajo diferentes hipótesis, queF{\displaystyle f}es integrable y tiene variación acotada en un intervalo abierto que contiene el puntot{\displaystyle t}, en cuyo caso 12límiteh0(F(t+h)+F(th))=20F(incógnita)porque(2πξ(incógnitat))dincógnitadξ.{\displaystyle {\tfrac {1}{2}}\lim _{h\to 0}\left(f(t+h)+f(t-h)\right)=2\int _{0}^{\infty }\int _{-\infty }^{\infty }f(x)\cos(2\pi \xi (x-t))\,dx\,d\xi .}

Esta última forma es un paso intermedio útil para demostrar las fórmulas inversas para las transformadas de coseno y . Un método para derivarla, debido a Cauchy , consiste en insertar unamiδξ{\displaystyle e^{-\delta \xi }}en la integral, dondeδ>0{\displaystyle \delta >0}Está arreglado. Entonces 20miδξporque(2πξ(incógnitat))dξF(incógnita)dincógnita=F(incógnita)2δδ2+4π2(incógnitat)2dincógnita.{\displaystyle 2\int _{-\infty }^{\infty }\int _{0}^{\infty }e^{-\delta \xi }\cos(2\pi \xi (x-t))\,d\xi \,f(x)\,dx=\int _{-\infty }^{\infty }f(x){\frac {2\delta }{\delta ^{2}+4\pi ^{2}(x-t)^{2}}}\,dx.} Ahora cuandoδ0{\displaystyle \delta \to 0}, el integrando tiende a cero excepto enincógnita=t{\displaystyle x=t}, de modo que formalmente lo anterior es F(t)2δδ2+4π2(incógnitat)2dincógnita=F(t).{\displaystyle f(t)\int _{-\infty }^{\infty }{\frac {2\delta }{\delta ^{2}+4\pi ^{2}(x-t)^{2}}}\,dx=f(t).}

Relación con exponenciales complejas

La forma exponencial compleja de la transformada de Fourier que se utiliza con más frecuencia hoy en día es [ 7 ].F^(ξ)=F(t)mi2πiξtdt{\displaystyle {\begin{aligned}{\hat {f}}(\xi )&=\int _{-\infty }^{\infty }f(t)e^{-2\pi i\xi t}\,dt\\\end{aligned}}\,}dóndei{\displaystyle i}es la raíz cuadrada de menos uno . Aplicando la fórmula de Euler(miiincógnita=porqueincógnita+ipecadoincógnita),{\textstyle (e^{ix}=\cos x+i\sin x),}Se puede demostrar (para funciones de valor real) que la componente real de la transformada de Fourier es la transformada del coseno (que representa la componente par de la función original) y la componente imaginaria de la transformada de Fourier es el negativo de la transformada del seno (que representa la componente impar de la función original): [ 8 ]F^(ξ)=F(t)(porque(2πξt)ipecado(2πξt))dtFórmula de Euler=(F(t)porque(2πξt)dt)i(F(t)pecado(2πξt)dt)=F^do(ξ)iF^s(ξ).{\displaystyle {\begin{aligned}{\hat {f}}(\xi )&=\int _{-\infty }^{\infty }f(t)\left(\cos(2\pi \xi t)-i\,\sin(2\pi \xi t)\right)dt&&{\text{Euler's Formula}}\\&=\left(\int _{-\infty }^{\infty }f(t)\cos(2\pi \xi t)\,dt\right)-i\left(\int _{-\infty }^{\infty }f(t)\sin(2\pi \xi t)\,dt\right)\\&={\hat {f}}^{c}(\xi )-i\,{\hat {f}}^{s}(\xi )\,.\end{aligned}}}Debido a esta relación, la transformada coseno de funciones cuya transformada de Fourier es conocida (por ejemplo, en la sección Tablas de transformadas de Fourier  importantes ) se puede encontrar simplemente tomando la parte real de la transformada de Fourier:F^do(ξ)=Rmi[F^(ξ)]{\displaystyle {\hat {f}}^{c}(\xi )=\mathrm {Re} {[\;{\hat {f}}(\xi )\;]}}mientras que la transformada seno es simplemente el negativo de la parte imaginaria de la transformada de Fourier:F^s(ξ)=Imetro[F^(ξ)].{\displaystyle {\hat {f}}^{s}(\xi )=-\mathrm {Im} {[\;{\hat {f}}(\xi )\;]}\,.}

Ventajas y desventajas

La suma de una onda sinusoidal (roja) y una onda cosenoidal (azul) de la misma frecuencia da como resultado una onda sinusoidal desfasada (verde) de la misma frecuencia, pero cuya amplitud y fase dependen de las amplitudes de las ondas sinusoidal y cosenoidal originales . Por lo tanto, a una frecuencia determinada, la transformada sinusoidal y la transformada cosenoidal juntas representan esencialmente una única onda sinusoidal que puede tener cualquier desfase.

Una ventaja de la transformada de Fourier moderna es que, mientras que para extraer la información de fase de una frecuencia se requieren las transformadas de seno y coseno , la transformada de Fourier moderna integra de forma compacta tanto la información de fase como la de amplitud en su resultado complejo. Sin embargo, una desventaja es que requiere conocimientos de números complejos, exponenciales complejos y frecuencias negativas.

Las transformadas de seno y coseno, por su parte, tienen la ventaja de que todas las cantidades son reales. Dado que las frecuencias positivas pueden expresarlas completamente, se puede evitar el concepto no trivial de frecuencia negativa necesario en la transformada de Fourier convencional. También pueden ser convenientes cuando la función original ya es par o impar, o puede hacerse par o impar, en cuyo caso solo se necesita la transformada de coseno o la de seno, respectivamente. Por ejemplo, aunque una entrada no sea par ni impar, una transformada discreta de coseno puede comenzar asumiendo una extensión par de su entrada, mientras que una transformada discreta de seno puede comenzar asumiendo una extensión impar de su entrada, para evitar tener que calcular toda la transformada discreta de Fourier .

Evaluación numérica

El uso de métodos estándar de evaluación numérica para integrales de Fourier, como la cuadratura gaussiana o tanh-sinh , probablemente conduzca a resultados completamente incorrectos, ya que la suma de cuadratura está (para la mayoría de los integrandos de interés) altamente mal condicionada. Se requieren métodos numéricos especiales que exploten la estructura de la oscilación, un ejemplo de la cual es el método de Ooura para integrales de Fourier [ 9 ]. Este método intenta evaluar el integrando en ubicaciones que se aproximan asintóticamente a los ceros de la oscilación (ya sea el seno o el coseno), reduciendo rápidamente la magnitud de los términos positivos y negativos que se suman.

Véase también

Notas

  1. La transformada sinusoidal a veces se denota conFs(F){\displaystyle {\mathcal {F}}_{s}(f)}en lugar deF^s{\displaystyle {\hat {f}}^{s}}.
  2. Si bien este artículo utiliza la frecuencia ordinaria paraξ{\displaystyle \xi }en ciclos por unidad de tiempo, que normalmente utiliza el Hertz y el segundo como unidades, estas transformaciones a veces se expresan utilizando la frecuencia angular en unidades angulares (por ejemplo, radianes ) por unidad de tiempo, dondeω{\displaystyle \omega } radianes por segundo es igual a2πξ{\displaystyle 2\pi \xi }.
  3. La transformada coseno a veces se denota conFdo(F){\displaystyle {\mathcal {F}}_{c}(f)}en lugar deF^do{\displaystyle {\hat {f}}^{c}}.
  4. Las hipótesis habituales son queF{\displaystyle f}y ambas transformadas deben ser absolutamente integrables. Para más detalles sobre las diferentes hipótesis, véase el teorema de inversión de Fourier .
  5. La transformada de Fourier moderna más generaltiene esta simetría incluso cuando las funciones originales no son pares ni impares. Una notación para denotar pares de transformadas de Fourier esF(t) F F^(ξ).{\displaystyle f(t)\ {\stackrel {\mathcal {F}}{\longleftrightarrow }}\ {\widehat {f}}(\xi ).}

Referencias

  • Whittaker, Edmund y James Watson, Un curso de análisis moderno , Cuarta edición, Cambridge Univ. Press, 1927, págs. 189, 211
  1. Mary L. Boas , Métodos matemáticos en las ciencias físicas , 2.ª ed., John Wiley & Sons Inc., 1983. ISBN 0-471-04409-1
  2. Nyack, Cuthbert (1996). "Transformada de Fourier, transformadas de coseno y seno" . cnyack.homestead.com . Archivado del original el 7 de junio de 2023. Consultado el 8 de octubre de 2018 .
  3. Coleman, Matthew P. (2013). Introducción a las ecuaciones diferenciales parciales con MATLAB (Segunda edición). Boca Raton. pág. 221. ISBN   978-1-4398-9846-8OCLC 822959644 {{cite book}}: CS1 mantenimiento: falta el editor de ubicación ( enlace )
  4. ^ Poincaré, Henri (1895). Teoría analítica de la propagación del calor . París: G. Carré. págs. 108 y siguientes. 
  5. Edwin Titchmarsh (1948), Introducción a la teoría de la integral de Fourier , Oxford en Clarendon Press, pág. 1 
  6. Whittaker, Edmund Taylor ; Watson, George Neville (2 de enero de 1927). Un curso de análisis moderno: Introducción a la teoría general de los procesos infinitos y de las funciones analíticas; con una descripción de las principales funciones trascendentales (4.ª ed.). Cambridge, Reino Unido: Cambridge University Press. pág. 189. ISBN   0-521-06794-4ISBN 978-0-521-06794-2.{{cite book}}: Incompatibilidad de ISBN/Fecha ( ayuda )
  7. Valentinuzzi, Max E. (25 de enero de 2016). "Aspectos destacados en la historia de la transformada de Fourier" . IEEE Pulse . Archivado del original el 15 de mayo de 2024. Consultado el 9 de septiembre de 2024 .
  8. Williams, Lance R. (2011-09-06). "Funciones pares e impares" (PDF) . www.cs.unm.edu/~williams/ . Archivado (PDF) del original el 2 de mayo de 2024. Consultado el 11 de septiembre de 2024 .
  9. Takuya Ooura, Masatake Mori, Una fórmula doble exponencial robusta para integrales de tipo Fourier , Journal of computational and applied mathematics 112.1-2 (1999): 229-241.