Articulo de referencia

Teorema de convolución

En matemáticas , el teorema de convolución establece que, bajo ciertas condiciones, la transformada de Fourier de la convolución de dos funciones (o señales ) es el producto de ...

En matemáticas , el teorema de convolución establece que, bajo ciertas condiciones, la transformada de Fourier de la convolución de dos funciones (o señales ) es el producto de sus transformadas de Fourier. De forma más general, la convolución en un dominio (por ejemplo, el dominio del tiempo ) equivale a la multiplicación punto a punto en el otro dominio (por ejemplo, el dominio de la frecuencia ). Existen otras versiones del teorema de convolución aplicables a diversas transformadas relacionadas con Fourier .

Funciones de una variable continua

Consideremos dos funciones y con transformadas de Fourier y : (incógnita){\displaystyle u(x)}v(incógnita){\displaystyle v(x)}U{\displaystyle U}V{\displaystyle V}

U(F)F{}(F)=(incógnita)mii2πFincógnitadincógnita,FRV(F)F{v}(F)=v(incógnita)mii2πFincógnitadincógnita,FR{\displaystyle {\begin{aligned}U(f)&\triangleq {\mathcal {F}}\{u\}(f)=\int _{-\infty }^{\infty }u(x)e^{-i2\pi fx}\,dx,\quad f\in \mathbb {R} \\[1ex]V(f)&\triangleq {\mathcal {F}}\{v\}(f)=\int _{-\infty }^{\infty }v(x)e^{-i2\pi fx}\,dx,\quad f\in \mathbb {R} \end{aligned}}}

donde denota el operador de transformada de Fourier . La transformada puede normalizarse de otras maneras, en cuyo caso aparecerán factores de escala constantes (típicamente o ) en el teorema de convolución que se presenta a continuación. La convolución de y se define por: F{\displaystyle {\mathcal {F}}}2π{\displaystyle 2\pi }2π{\displaystyle {\sqrt {2\pi }}}{\displaystyle u}v{\displaystyle v}

r(incógnita)={v}(incógnita)(τ)v(incógnitaτ)dτ=(incógnitaτ)v(τ)dτ.{\displaystyle r(x)=\{u*v\}(x)\triangleq \int _{-\infty }^{\infty }u(\tau )v(x-\tau )\,d\tau =\int _{-\infty }^{\infty }u(x-\tau )v(\tau )\,d\tau .}

En este contexto, el asterisco indica convolución, en lugar de la multiplicación estándar. En ocasiones, se utiliza el símbolo de producto tensorial . {\displaystyle \otimes }

El teorema de convolución establece que: [ 1 ] [ 2 ] : ecuación 8

R(F)F{r}(F)=U(F)V(F).FR{\displaystyle R(f)\triangleq {\mathcal {F}}\{r\}(f)=U(f)V(f).\quad f\in \mathbb {R} }   ( Ecuación 1a )

La aplicación de la transformada inversa de Fourier produce el corolario : [ 2 ] : ecuaciones 7, 10 F1,{\displaystyle {\mathcal {F}}^{-1},}

Teorema de convolución

r(incógnita)={v}(incógnita)=F1{UV}.{\displaystyle r(x)=\{u*v\}(x)={\mathcal {F}}^{-1}\{U\cdot V\}.}   ( Ecuación 1b )

El teorema también se aplica generalmente a funciones multidimensionales.

Este teorema también es válido para la transformada de Laplace , la transformada de Laplace bilateral y, con las modificaciones adecuadas, para la transformada de Mellin y la transformada de Hartley (véase el teorema de inversión de Mellin ). Puede extenderse a la transformada de Fourier del análisis armónico abstracto definido sobre grupos abelianos localmente compactos .

Convolución periódica (coeficientes de la serie de Fourier)

Consideremos funciones -periódicas  que   pueden expresarse como sumas periódicas : P{\displaystyle P}uP{\displaystyle u_{_{P}}}vP,{\displaystyle v_{_{P}},}

uP(x) m=u(xmP){\displaystyle u_{_{P}}(x)\ \triangleq \sum _{m=-\infty }^{\infty }u(x-mP)}yvP(x) m=v(xmP).{\displaystyle v_{_{P}}(x)\ \triangleq \sum _{m=-\infty }^{\infty }v(x-mP).}

En la práctica, la parte no nula de los componentes y a menudo se limita a la duración, pero nada en el teorema exige eso. u{\displaystyle u}v{\displaystyle v}P,{\displaystyle P,}

Los coeficientes de la serie de Fourier son:

U[k]F{uP}[k]=1PPuP(x)ei2πkx/Pdx,kZ;integration over any interval of length PV[k]F{vP}[k]=1PPvP(x)ei2πkx/Pdx,kZ{\displaystyle {\begin{aligned}U[k]&\triangleq {\mathcal {F}}\{u_{_{P}}\}[k]={\frac {1}{P}}\int _{P}u_{_{P}}(x)e^{-i2\pi kx/P}\,dx,\quad k\in \mathbb {Z} ;\quad \quad \scriptstyle {\text{integration over any interval of length }}P\\V[k]&\triangleq {\mathcal {F}}\{v_{_{P}}\}[k]={\frac {1}{P}}\int _{P}v_{_{P}}(x)e^{-i2\pi kx/P}\,dx,\quad k\in \mathbb {Z} \end{aligned}}}

donde denota la integral de la serie de Fourier . F{\displaystyle {\mathcal {F}}}

  • El producto: también es -periódico, y sus coeficientes de la serie de Fourier vienen dados por la convolución discreta de las secuencias y :uP(x)vP(x){\displaystyle u_{_{P}}(x)\cdot v_{_{P}}(x)}P{\displaystyle P}U{\displaystyle U}V{\displaystyle V}F{uPvP}[k]={UV}[k].{\displaystyle {\mathcal {F}}\{u_{_{P}}\cdot v_{_{P}}\}[k]=\{U*V\}[k].}
  • La convolución: también es -periódica, y se denomina convolución periódica .{uPv}(x) uP(xτ)v(τ) dτPuP(xτ)vP(τ) dτ;integration over any interval of length P{\displaystyle {\begin{aligned}\{u_{_{P}}*v\}(x)\ &\triangleq \int _{-\infty }^{\infty }u_{_{P}}(x-\tau )\cdot v(\tau )\ d\tau \\&\equiv \int _{P}u_{_{P}}(x-\tau )\cdot v_{_{P}}(\tau )\ d\tau ;\quad \quad \scriptstyle {\text{integration over any interval of length }}P\end{aligned}}}P{\displaystyle P}

El teorema de convolución correspondiente es :

F{uPv}[k]= PU[k] V[k].{\displaystyle {\mathcal {F}}\{u_{_{P}}*v\}[k]=\ P\cdot U[k]\ V[k].}   ( Ecuación 2 )

Funciones de una variable discreta (secuencias)

Mediante una derivación similar a la ecuación 1, existe un teorema análogo para secuencias, como muestras de dos funciones continuas, donde ahora denota el operador de la transformada discreta de Fourier (DTFT). Consideremos dos secuencias y con transformadas y : F{\displaystyle {\mathcal {F}}}u[n]{\displaystyle u[n]}v[n]{\displaystyle v[n]}U{\displaystyle U}V{\displaystyle V}

U(f)F{u}(f)=n=u[n]ei2πfn,fR,V(f)F{v}(f)=n=v[n]ei2πfn,fR.{\displaystyle {\begin{aligned}U(f)&\triangleq {\mathcal {F}}\{u\}(f)=\sum _{n=-\infty }^{\infty }u[n]\cdot e^{-i2\pi fn}\;,\quad f\in \mathbb {R} ,\\V(f)&\triangleq {\mathcal {F}}\{v\}(f)=\sum _{n=-\infty }^{\infty }v[n]\cdot e^{-i2\pi fn}\;,\quad f\in \mathbb {R} .\end{aligned}}}

La convolución discreta de y se define por: u{\displaystyle u}v{\displaystyle v}

r[n](uv)[n]=m=u[m]v[nm]=m=u[nm]v[m].{\displaystyle r[n]\triangleq (u*v)[n]=\sum _{m=-\infty }^{\infty }u[m]\cdot v[n-m]=\sum _{m=-\infty }^{\infty }u[n-m]\cdot v[m].}

El teorema de convolución para secuencias discretas es: [ 3 ] [ 4 ] : p.60 (2.169)

R(f)=F{uv}(f)= U(f)V(f).{\displaystyle R(f)={\mathcal {F}}\{u*v\}(f)=\ U(f)V(f).}   ( Ec. 3 )

convolución periódica

U(f){\displaystyle U(f)}y como se definió anteriormente, son periódicas, con un período de 1. Consideremos secuencias -periódicas y : V(f),{\displaystyle V(f),}N{\displaystyle N}uN{\displaystyle u_{_{N}}}vN{\displaystyle v_{_{N}}}

uN[n] m=u[nmN]{\displaystyle u_{_{N}}[n]\ \triangleq \sum _{m=-\infty }^{\infty }u[n-mN]}yvN[n] m=v[nmN],nZ.{\displaystyle v_{_{N}}[n]\ \triangleq \sum _{m=-\infty }^{\infty }v[n-mN],\quad n\in \mathbb {Z} .}

Estas funciones aparecen como resultado del muestreo y a intervalos de y realizando una transformada discreta de Fourier inversa (DFT) en las muestras (véase § Muestreo de la DTFT ). La convolución discreta :U{\displaystyle U}V{\displaystyle V}1/N{\displaystyle 1/N}N{\displaystyle N}

{uNv}[n] m=uN[m]v[nm]m=0N1uN[m]vN[nm]{\displaystyle \{u_{_{N}}*v\}[n]\ \triangleq \sum _{m=-\infty }^{\infty }u_{_{N}}[m]\cdot v[n-m]\equiv \sum _{m=0}^{N-1}u_{_{N}}[m]\cdot v_{_{N}}[n-m]}

También es -periódica y se denomina convolución periódica . Redefiniendo el operador como la DFT de longitud -, el teorema correspondiente es: [ 5 ] [ 4 ] : pág. 548 N{\displaystyle N}F{\displaystyle {\mathcal {F}}}N{\displaystyle N}

F{uNv}[k]= F{uN}[k]U(k/N)F{vN}[k]V(k/N),kZ.{\displaystyle {\mathcal {F}}\{u_{_{N}}*v\}[k]=\ \underbrace {{\mathcal {F}}\{u_{_{N}}\}[k]} _{U(k/N)}\cdot \underbrace {{\mathcal {F}}\{v_{_{N}}\}[k]} _{V(k/N)},\quad k\in \mathbb {Z} .}   ( Ecuación 4a )

Y por lo tanto:

{uNv}[n]= F1{F{uN}F{vN}}.{\displaystyle \{u_{_{N}}*v\}[n]=\ {\mathcal {F}}^{-1}\{{\mathcal {F}}\{u_{_{N}}\}\cdot {\mathcal {F}}\{v_{_{N}}\}\}.}   ( Ecuación 4b )

En las condiciones adecuadas, es posible que esta secuencia de longitud contenga un segmento de convolución sin distorsión. Pero cuando la porción no nula de la secuencia es igual o mayor que cierta distorsión es inevitable. Tal es el caso cuando la secuencia se obtiene muestreando directamente la DTFT de la respuesta impulsional de la transformada discreta de Hilbert infinitamente larga . [ A ]N{\displaystyle N}uv{\displaystyle u*v}u(n){\displaystyle u(n)}v(n){\displaystyle v(n)}N,{\displaystyle N,}V(k/N){\displaystyle V(k/N)}

Para secuencias cuya duración no nula es menor o igual a una simplificación final es: u{\displaystyle u}v{\displaystyle v}N,{\displaystyle N,}

convolución circular

{uNv}[n]= F1{F{u}F{v}}.{\displaystyle \{u_{_{N}}*v\}[n]=\ {\mathcal {F}}^{-1}\{{\mathcal {F}}\{u\}\cdot {\mathcal {F}}\{v\}\}.}   ( Ecuación 4c )

Esta forma se utiliza a menudo para implementar de forma eficiente la convolución numérica mediante ordenador . (Véase § Algoritmos de convolución rápidos y § Ejemplo ).

Como recíproco parcial, se ha demostrado [ 6 ] que cualquier transformada lineal que convierte la convolución en un producto es la DFT (salvo una permutación de coeficientes).

Teorema de convolución para la transformada inversa de Fourier

También existe un teorema de convolución para la transformada inversa de Fourier:

Aquí, " " representa el producto de Hadamard , y " " representa una convolución entre las dos matrices. {\displaystyle \cdot }{\displaystyle *}

F{uv}=F{u}F{v}F{uv}=F{u}F{v}{\displaystyle {\begin{aligned}&{\mathcal {F}}\{u*v\}={\mathcal {F}}\{u\}\cdot {\mathcal {F}}\{v\}\\&{\mathcal {F}}\{u\cdot v\}={\mathcal {F}}\{u\}*{\mathcal {F}}\{v\}\end{aligned}}}

de modo que

uv=F1{F{u}F{v}}uv=F1{F{u}F{v}}{\displaystyle {\begin{aligned}&u*v={\mathcal {F}}^{-1}\left\{{\mathcal {F}}\{u\}\cdot {\mathcal {F}}\{v\}\right\}\\&u\cdot v={\mathcal {F}}^{-1}\left\{{\mathcal {F}}\{u\}*{\mathcal {F}}\{v\}\right\}\end{aligned}}}

Teorema de convolución para distribuciones temperadas

El teorema de convolución se extiende a las distribuciones temperadas . Aquí, es una distribución temperada arbitraria: v{\displaystyle v}

F{uv}=F{u}F{v}F{uv}=F{u}F{v}.{\displaystyle {\begin{aligned}&{\mathcal {F}}\{u*v\}={\mathcal {F}}\{u\}\cdot {\mathcal {F}}\{v\}\\&{\mathcal {F}}\{u\cdot v\}={\mathcal {F}}\{u\}*{\mathcal {F}}\{v\}.\end{aligned}}}

Pero debe ser "rápidamente decreciente" hacia y para garantizar la existencia tanto del producto de convolución como del producto de multiplicación. De forma equivalente, si es una función ordinaria suave de "crecimiento lento", garantiza la existencia tanto del producto de multiplicación como del producto de convolución. [ 7 ] [ 8 ] [ 9 ]α=F{u}{\displaystyle \alpha =F\{u\}}{\displaystyle -\infty }+{\displaystyle +\infty }u=F1{α}{\displaystyle u=F^{-1}\{\alpha \}}

En particular, toda distribución temperada con soporte compacto, como la delta de Dirac , es "rápidamente decreciente". De forma equivalente, las funciones de banda limitada , como la función que es constantemente, son funciones ordinarias suaves de "crecimiento lento". Si, por ejemplo, es el peine de Dirac , ambas ecuaciones producen la fórmula de suma de Poisson y si, además, es la delta de Dirac, entonces es constantemente uno y estas ecuaciones producen la identidad del peine de Dirac . 1{\displaystyle 1}vШ{\displaystyle v\equiv \operatorname {\text{Ш}} }uδ{\displaystyle u\equiv \delta }α1{\displaystyle \alpha \equiv 1}

Véase también

Notas

  1. ^ Un ejemplo es la función de MATLAB , hilbert(u,N) .

Referencias

  1. ^ McGillem, Clare D.; Cooper, George R. (1984). Análisis de señales y sistemas continuos y discretos (2.ª ed.). Holt, Rinehart and Winston. pág. 118 (3–102). ISBN 0-03-061703-0.
  2. ^ a b Weisstein, Eric W. "Teorema de convolución" . De MathWorld--A Wolfram Web Resource . Recuperado el 8 de febrero de 2021 .
  3. ^ Proakis, John G.; Manolakis, Dimitri G. (1996), Procesamiento digital de señales: principios, algoritmos y aplicaciones (3.ª ed.), Nueva Jersey: Prentice-Hall International, pág. 297, Bibcode : 1996dspp.book.....P , ISBN 9780133942897, sAcfAQAAIAAJ
  4. ^ a b Oppenheim, Alan V. ; Schafer, Ronald W. ; Buck, John R. (1999). Procesamiento de señales en tiempo discreto (2.ª ed.). Upper Saddle River, NJ: Prentice Hall. ISBN 0-13-754920-2.
  5. ^ Rabiner, Lawrence R. ; Gold, Bernard (1975). Teoría y aplicación del procesamiento digital de señales . Englewood Cliffs, NJ: Prentice-Hall, Inc. pág. 59 (2.163). ISBN 978-0139141010.
  6. ^ Amiot, Emmanuel (2016). Música a través del espacio de Fourier . Ciencia musical computacional. Zúrich: Springer. p. 8. doi : 10.1007/978-3-319-45581-5 . ISBN 978-3-319-45581-5. S2CID  6224021 .
  7. ^ Horváth, John (1966). Espacios vectoriales topológicos y distribuciones . Reading, MA: Addison-Wesley Publishing Company.
  8. ^ Barros-Neto, José (1973). Introducción a la teoría de las distribuciones . Nueva York, Nueva York: Dekker.
  9. ^ Petersen, Bent E. (1983). Introducción a la transformada de Fourier y operadores pseudodiferenciales . Boston, MA: Pitman Publishing.

Lecturas adicionales

  • Katznelson, Yitzhak (1976), Introducción al análisis armónico , Dover, ISBN 0-486-63331-4
  • Li, Bing; Babu, G. Jogesh (2019), "Teorema de convolución y eficiencia asintótica", A Graduate Course on Statistical Inference , Nueva York: Springer, pp.  295–327 , ISBN 978-1-4939-9759-6
  • Crutchfield, Steve (9 de octubre de 2010), "The Joy of Convolution" , Universidad Johns Hopkins , consultado el 19 de noviembre de 2010.

Recursos adicionales

Para una representación visual del uso del teorema de convolución en el procesamiento de señales , consulte:

Obtenido de " https://en.wikipedia.org/w/index.php?title=Convolution_theorem&oldid=1359305879 "