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(incógnita){\displaystyle u(x)}yv(incógnita){\displaystyle v(x)}con transformadas de FourierU{\displaystyle U}yV{\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}}}

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

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 denota convolución, en lugar de multiplicación estándar. El símbolo del producto tensorial{\displaystyle \otimes }a veces se utiliza en su lugar.

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 )

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

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)

ConsiderarPAG{\displaystyle P}-funciones periódicasPAG{\displaystyle u_{_{P}}} y vPAG,{\displaystyle v_{_{P}},}que pueden expresarse como sumas periódicas :

PAG(incógnita) metro=(incógnitametroPAG){\displaystyle u_{_{P}}(x)\ \triangleq \sum _{m=-\infty }^{\infty }u(x-mP)}yvPAG(incógnita) metro=v(incógnitametroPAG).{\displaystyle v_{_{P}}(x)\ \triangleq \sum _{m=-\infty }^{\infty }v(x-mP).}

En la práctica, la porción no nula de los componentes{\displaystyle u}yv{\displaystyle v}a menudo están limitados a la duraciónPAG,{\displaystyle P,}pero nada en el teorema lo exige.

Los coeficientes de la serie de Fourier son:

U[k]F{PAG}[k]=1PAGPAGPAG(incógnita)mii2πkincógnita/PAGdincógnita,kZ;integración sobre cualquier intervalo de longitud PAGV[k]F{vPAG}[k]=1PAGPAGvPAG(incógnita)mii2πkincógnita/PAGdincógnita,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{integración sobre cualquier intervalo de longitud }}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}}}

dóndeF{\displaystyle {\mathcal {F}}}denota la integral de la serie de Fourier .

  • El producto:PAG(incógnita)vPAG(incógnita){\displaystyle u_{_{P}}(x)\cdot v_{_{P}}(x)}también lo esPAG{\displaystyle P}-periódica, y sus coeficientes de la serie de Fourier vienen dados por la convolución discreta de laU{\displaystyle U}yV{\displaystyle V}secuencias:F{PAGvPAG}[k]={UV}[k].{\displaystyle {\mathcal {F}}\{u_{_{P}}\cdot v_{_{P}}\}[k]=\{U*V\}[k].}
  • La convolución:{PAGv}(incógnita) PAG(incógnitaτ)v(τ) dτPAGPAG(incógnitaτ)vPAG(τ) dτ;integración sobre cualquier intervalo de longitud PAG{\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{integración sobre cualquier intervalo de longitud }}P\end{aligned}}} también esPAG{\displaystyle P}-periódica, y se denomina convolución periódica .

El teorema de convolución correspondiente es :

F{PAGv}[k]= PAGU[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 ahoraF{\displaystyle {\mathcal {F}}}denota el operador de transformada discreta de Fourier (DTFT). Consideremos dos secuencias[norte]{\displaystyle u[n]}yv[norte]{\displaystyle v[n]}con transformacionesU{\displaystyle U}yV{\displaystyle V}:

U(F)F{}(F)=norte=[norte]mii2πFnorte,FR,V(F)F{v}(F)=norte=v[norte]mii2πFnorte,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{\displaystyle u}yv{\displaystyle v}se define por:

r[norte](v)[norte]=metro=[metro]v[nortemetro]=metro=[nortemetro]v[metro].{\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{v}(F)= U(F)V(F).{\displaystyle R(f)={\mathcal {F}}\{u*v\}(f)=\ U(f)V(f).}   ( Ecuación 3 )

convolución periódica

U(F){\displaystyle U(f)}yV(F),{\displaystyle V(f),}como se definieron anteriormente, son periódicos, con un período de 1. Considerenorte{\displaystyle N}-secuencias periódicasnorte{\displaystyle u_{_{N}}}yvnorte{\displaystyle v_{_{N}}}:

norte[norte] metro=[nortemetronorte]{\displaystyle u_{_{N}}[n]\ \triangleq \sum _{m=-\infty }^{\infty }u[n-mN]}yvnorte[norte] metro=v[nortemetronorte],norteZ.{\displaystyle v_{_{N}}[n]\ \triangleq \sum _{m=-\infty }^{\infty }v[n-mN],\quad n\in \mathbb {Z} .}

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

{nortev}[norte] metro=norte[metro]v[nortemetro]metro=0norte1norte[metro]vnorte[nortemetro]{\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 lo esnorte{\displaystyle N}-periódica, y se denomina convolución periódica . Redefiniendo laF{\displaystyle {\mathcal {F}}}operador como elnorte{\displaystyle N}-longitud DFT, el teorema correspondiente es: [ 5 ] [ 4 ] : pág. 548

F{nortev}[k]= F{norte}[k]U(k/norte)F{vnorte}[k]V(k/norte),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:

{nortev}[norte]= F1{F{norte}F{vnorte}}.{\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 estonorte{\displaystyle N}-secuencia de longitud para contener un segmento sin distorsión de unav{\displaystyle u*v}convolución. Pero cuando la porción no nula de la(norte){\displaystyle u(n)}ov(norte){\displaystyle v(n)}La secuencia es igual o más larga quenorte,{\displaystyle N,}Cierta distorsión es inevitable. Tal es el caso cuando laV(k/norte){\displaystyle V(k/N)}La secuencia se obtiene muestreando directamente la DTFT de la respuesta impulsional de la transformada discreta de Hilbert infinitamente  larga . [ A ]

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

convolución circular

{nortev}[norte]= F1{F{}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í, "{\displaystyle \cdot }" representa el producto Hadamard , y "{\displaystyle *}" representa una convolución entre las dos matrices.

F{v}=F{}F{v}F{v}=F{}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

v=F1{F{}F{v}}v=F1{F{}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í,v{\displaystyle v}es una distribución temperada arbitraria:

F{v}=F{}F{v}F{v}=F{}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α=F{}{\displaystyle \alpha =F\{u\}}debe estar "disminuyendo rápidamente" hacia{\displaystyle -\infty }y+{\displaystyle +\infty }para garantizar la existencia tanto de la convolución como del producto de multiplicación. Equivalentemente, si=F1{α}{\displaystyle u=F^{-1}\{\alpha \}}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 ]

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

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 continuo y discreto de señales y sistemas (2.ª ed.). Holt, Rinehart and Winston. pág. 118 (3–102). ISBN   0-03-061703-0.
  2. 1 2 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. 1 2 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ág. 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: