De Wikipedia, la enciclopedia libre
Saltar a navegación Saltar a búsqueda

En matemáticas , la transformada de Fourier de tiempo discreto ( DTFT ) es una forma de análisis de Fourier que se aplica a una secuencia de valores.

El DTFT se utiliza a menudo para analizar muestras de una función continua. El término tiempo discreto se refiere al hecho de que la transformada opera sobre datos discretos, a menudo muestras cuyo intervalo tiene unidades de tiempo. A partir de muestras uniformemente espaciadas, produce una función de frecuencia que es una suma periódica de la transformada de Fourier continua de la función continua original. En determinadas condiciones teóricas, descritas por el teorema de muestreo , la función continua original se puede recuperar perfectamente de la DTFT y, por tanto, de las muestras discretas originales. La DTFT en sí es una función continua de la frecuencia, pero las muestras discretas se pueden calcular fácilmente mediante la transformada discreta de Fourier.(DFT) (ver § Muestreo de DTFT ), que es, con mucho, el método más común de análisis de Fourier moderno.

Ambas transformaciones son invertibles. La DTFT inversa es la secuencia de datos muestreada original. La DFT inversa es una suma periódica de la secuencia original. La transformada rápida de Fourier (FFT) es un algoritmo para calcular un ciclo de la DFT, y su inversa produce un ciclo de la DFT inversa.

Definición [ editar ]

La transformada de Fourier de tiempo discreto de una secuencia discreta de números reales o complejos x [ n ] , para todos los enteros n , es una serie de Fourier , que produce una función periódica de una variable de frecuencia. Cuando la variable de frecuencia, ω, tiene unidades normalizadas de radianes / muestra , la periodicidad es y la serie de Fourier es : [1] : p.147

La utilidad de esta función en el dominio de la frecuencia se basa en la fórmula de suma de Poisson . Sea X ( f ) la transformada de Fourier de cualquier función, x ( t ) , cuyas muestras en algún intervalo T ( segundos ) son iguales (o proporcionales) a la secuencia x [ n ] , es decir, Tx ( nT ) = x [ n ] . [2]   Entonces, la función periódica representada por la serie de Fourier es una suma periódica de X ( f) en términos de frecuencia f en hercios ( ciclos / seg ) : [a] [A]

Fig 1. Representación de una transformada de Fourier (superior izquierda) y su suma periódica (DTFT) en la esquina inferior izquierda. La esquina inferior derecha muestra muestras de la DTFT que se calculan mediante una transformada discreta de Fourier (DFT).

El número entero k tiene unidades de ciclos / muestra , y 1 / T es la frecuencia de muestreo, f s ( muestras / seg ). Por tanto, X 1 / T ( f ) comprende copias exactas de X ( f ) que se desplazan en múltiplos de f s hertz y se combinan mediante la suma. Para f s suficientemente grande, el término k = 0 se puede observar en la región [- f s / 2, f s / 2 ] con poca o ninguna distorsión (aliasing ) de los otros términos. En la figura 1, los extremos de la distribución en la esquina superior izquierda están enmascarados por aliasing en la suma periódica (inferior izquierda).

También notamos que e - i2πfTn es la transformada de Fourier de δ ( t - nT ) . Por lo tanto, una definición alternativa de DTFT es: [B]

La función de peine de Dirac modulada es una abstracción matemática que a veces se denomina muestreo de impulsos . [4]

Transformación inversa [ editar ]

Una operación que recupera la secuencia de datos discretos de la función DTFT se denomina DTFT inversa . Por ejemplo, la transformada de Fourier continua inversa de ambos lados de la ecuación 3 produce la secuencia en la forma de una función de peine de Dirac modulada:

Sin embargo, observando que X 1 / T ( f ) es periódica, toda la información necesaria está contenida dentro de cualquier intervalo de longitud 1 / T . Tanto en la Ec. 1 como en la Ec. 2 , las sumas sobre n son una serie de Fourier , con coeficientes x [ n ] . Las fórmulas estándar para los coeficientes de Fourier también son las transformadas inversas:

Datos periódicos [ editar ]

Cuando la secuencia de datos de entrada x [ n ] es N -periódica , la ecuación 2 se puede reducir computacionalmente a una transformada discreta de Fourier (DFT), porque :

  • Toda la información disponible está contenida en N muestras.
  • X 1 / T ( f ) converge a cero en todas partes excepto en múltiplos enteros de 1 / ( NT ) , conocidos comofrecuencias armónicas . En esas frecuencias, la DTFT diverge a diferentes tasas dependientes de la frecuencia. Y esas tasas vienen dadas por la DFT de un ciclo de lasecuencia x [ n ] .
  • La DTFT es periódica, por lo que el número máximo de amplitudes armónicas únicas es (1 / T ) / (1 / ( NT )) = N

Los coeficientes DFT vienen dados por :

    y la DTFT es :
      [B]

La sustitución de esta expresión en la fórmula de transformación inversa confirma :

( todos los enteros )

como se esperaba. La DFT inversa en la línea anterior a veces se denomina serie discreta de Fourier (DFS). [1] : pág. 542

Muestreo de DTFT [ editar ]

Cuando la DTFT es continua, una práctica común es calcular un número arbitrario de muestras ( N ) de un ciclo de la función periódica X 1 / T :  [1] : págs. 557–559 y 703

donde es una suma periódica :

    (ver la serie discreta de Fourier )

La secuencia es la DFT inversa. Por lo tanto, nuestro muestreo de la DTFT hace que la transformada inversa se vuelva periódica. La matriz de | X k | 2 valores se conoce como periodograma , y el parámetro N se llama NFFT en la función Matlab del mismo nombre. [5]

Para evaluar un ciclo numéricamente, necesitamos una secuencia x [ n ] de longitud finita . Por ejemplo, una secuencia larga puede ser truncada por una función de ventana de longitud L, dando como resultado tres casos dignos de mención especial. Para simplificar la notación, considere los valores x [ n ] siguientes para representar los valores modificados por la función de ventana.

Caso: Decimación de frecuencia. L = NI , para algún número entero I (típicamente 6 u 8)

Un ciclo de reduce a una suma de I segmentos de longitud N . La DFT luego tiene varios nombres, tales como :

  • ventana-presum FFT [6]
  • Peso, superposición, agregar (WOLA) [7] [8] [9] [10] [11] [12] [C] [D]
  • polifase DFT [10] [11]
  • banco de filtros polifásicos [13]
  • ventanas de bloques múltiples y alias de tiempo . [14]

Recuerde que la diezmación de los datos muestreados en un dominio (tiempo o frecuencia) produce superposición (a veces conocida como aliasing ) en el otro y viceversa. En comparación con una DFT de longitud L , la suma / superposición causa una reducción en la frecuencia, [1] : p.558 dejando solo las muestras DTFT menos afectadas por la fuga espectral . Eso suele ser una prioridad cuando se implementa un banco de filtros FFT (canalizador). Con una función de ventana convencional de longitud L , la pérdida de festón sería inaceptable. Por lo tanto, las ventanas de bloques múltiples se crean utilizando herramientas de diseño de filtros FIR . [15] [16]  Su perfil de frecuencia es plano en el punto más alto y cae rápidamente en el punto medio entre las muestras restantes de DTFT. Cuanto mayor sea el valor del parámetro I , mejor será el rendimiento potencial.

Caso: L = N +1 .

Cuando una función de ventana simétrica de longitud L ( ) se trunca por 1 coeficiente, se denomina periódica o DFT-par . El truncamiento afecta a la DTFT. Una DFT de las secuencia truncada muestras de la DTFT a intervalos de frecuencia de 1 / N . Para muestrear a las mismas frecuencias, a modo de comparación, la DFT se calcula para un ciclo de la suma periódica, [E]

Fig 2. DFT de e i2πn / 8 para L = 64 y N = 256
Fig 3. DFT de e i2πn / 8 para L = 64 y N = 64

Caso: interpolación de frecuencia. LN

En este caso, el DFT se simplifica a una forma más familiar:

Para aprovechar un algoritmo de transformada de Fourier rápido para calcular la DFT, la suma se realiza generalmente sobre todos los N términos, aunque N - L de ellos son ceros. Por lo tanto, el caso L < N a menudo se denomina relleno de ceros .

La fuga espectral, que aumenta a medida que disminuye L , es perjudicial para ciertas métricas de rendimiento importantes, como la resolución de múltiples componentes de frecuencia y la cantidad de ruido medido por cada muestra de DTFT. Pero esas cosas no siempre importan, por ejemplo, cuando la secuencia x [ n ] es una sinusoide silenciosa (o una constante), formada por una función de ventana. Entonces, es una práctica común utilizar el relleno de ceros para mostrar y comparar gráficamente los patrones de fuga detallados de las funciones de la ventana. Para ilustrar eso para una ventana rectangular, considere la secuencia:

y

Las Figuras 2 y 3 son gráficos de la magnitud de dos DFT de diferentes tamaños, como se indica en sus etiquetas. En ambos casos, el componente dominante está en la frecuencia de la señal: f = 1/8 = 0,125 . También visible en la Fig. 2 es el patrón de fuga espectral de la ventana rectangular L = 64 . La ilusión de la figura 3 es el resultado de muestrear la DTFT solo en sus cruces por cero. En lugar de la DTFT de una secuencia de longitud finita, da la impresión de una secuencia sinusoidal infinitamente larga. Los factores que contribuyen a la ilusión son el uso de una ventana rectangular y la elección de una frecuencia (1/8 = 8/64) con exactamente 8 (un número entero) ciclos por 64 muestras. Una ventana de Hannproduciría un resultado similar, excepto que el pico se ampliaría a 3 muestras (consulte la ventana DFT-even Hann ).

Convolución [ editar ]

El teorema de convolución para sucesiones es :

[18] : p . 297 [c]

Un caso especial importante es la convolución circular de secuencias de x y y definidos por donde es una suma periódica. La naturaleza de frecuencia discreta de significa que el producto con la función continua también es discreto, lo que resulta en una simplificación considerable de la transformada inversa :

[19] [1] : pág.548

Para x y Y secuencias cuya duración no nula es menor que o igual a N , una simplificación final es :

La importancia de este resultado se explica en los algoritmos de convolución circular y convolución rápida .

Propiedades de simetría [ editar ]

Cuando las partes real e imaginaria de una función compleja se descomponen en sus partes pares e impares , hay cuatro componentes, indicados a continuación por los subíndices RE, RO, IE e IO. Y hay un mapeo uno a uno entre los cuatro componentes de una función de tiempo compleja y los cuatro componentes de su transformada de frecuencia compleja : [18] : p.291

A partir de esto, son evidentes varias relaciones, por ejemplo :

  • La transformada de una función de valor real ( x RE + x RO ) es la función simétrica par X RE + i X IO . Por el contrario, una transformación par simétrica implica un dominio de tiempo de valor real.
  • La transformada de una función con valores imaginarios ( i x IE + i x IO ) es la función simétrica impar X RO + i X IE , y lo contrario es cierto.
  • La transformada de una función par-simétrica ( x RE + i x IO ) es la función de valor real X RE + X RO , y lo contrario es cierto.
  • La transformada de una función de simetría impar ( x RO + i x IE ) es la función de valor imaginario i X IE + i X IO , y lo contrario es cierto.

Relación con la transformada Z [ editar ]

es una serie de Fourier que también se puede expresar en términos de la transformada Z bilateral . Es decir:

donde la notación distingue la transformada Z de la transformada de Fourier. Por lo tanto, también podemos expresar una parte de la transformada Z en términos de la transformada de Fourier:

Tenga en cuenta que cuando el parámetro T cambia, los términos de permanecen separados por una separación constante y su ancho aumenta o disminuye. Los términos de X 1 / T ( f ) permanecen con un ancho constante y su separación 1 / T aumenta o disminuye.

Tabla de transformadas de Fourier de tiempo discreto [ editar ]

En la siguiente tabla se muestran algunos pares de transformadas comunes. Se aplica la siguiente notación:

  • es un número real que representa la frecuencia angular continua (en radianes por muestra). ( está en ciclos / seg, y está en seg / muestra). En todos los casos de la tabla, la DTFT es 2π-periódica (pulgadas ).
  • designa una función definida en .
  • designa una función definida en , y cero en otra parte. Luego:
  • es la función delta de Dirac
  • es la función sinc normalizada
  • es la función del rectángulo
  • es la función triangular
  • n es un número entero que representa el dominio de tiempo discreto (en muestras)
  • es la función de paso de la unidad de tiempo discreto
  • es el delta de Kronecker

Propiedades [ editar ]

Esta tabla muestra algunas operaciones matemáticas en el dominio del tiempo y los efectos correspondientes en el dominio de la frecuencia.

  • es la convolución discreta de dos secuencias
  • es el conjugado complejo de x [ n ] .

Ver también [ editar ]

  • Análisis espectral de mínimos cuadrados
  • Transformación multidimensional
  • Transformación Zak

Notas [ editar ]

  1. ^ Cuando la dependencia de T no es importante, una práctica común es reemplazarla con  Entonces   f   tiene unidades de ( ciclos / muestra ), llamada frecuencia normalizada .
  2. ^ De hecho, la ecuación 2 a menudo se justifica de la siguiente manera: [1] : p.143
  3. ^ WOLA no debe confundirse con el método Overlap-add de convolución por partes.
  4. ^ Ejemplo de WOLA: Archivo: ejemplo de canalizador WOLA.png
  5. ^ Un ejemplo es la figura Muestreo de DTFT . Las muestras DFT con valores reales son el resultado de la simetría DFT-par [17] : p.52
  6. ^ Esta expresión se deriva de la siguiente manera: [1] : p.168

Citas de página [ editar ]

  1. Oppenheim y Schafer, [1] p 147 (4.20), p 694 (10.1), y Prandoni y Vetterli, [3] p 255, (9.33), donde:    y    
  2. Oppenheim y Schafer, [1] p 551 (8.35), y Prandoni y Vetterli, [3] p 82, (4.43), donde :     y     
  3. Oppenheim y Schafer, [1] p 60, (2.169) y Prandoni y Vetterli, [3] p 122, (5.21)

Referencias [ editar ]

  1. ↑ a b c d e f g h i j k Oppenheim, Alan V .; Schafer, Ronald W .; Buck, John R. (1999). "4,2, 8,4". Procesamiento de señales en tiempo discreto (2ª ed.). Upper Saddle River, Nueva Jersey: Prentice Hall. ISBN 0-13-754920-2. Las muestras de la transformada de Fourier de una secuencia aperiódica x [n] pueden considerarse como coeficientes DFS de una secuencia periódica obtenida mediante la suma de réplicas periódicas de x [n].  url = https://d1.amobbs.com/bbs_upload782111/files_24/ourdev_523225.pdf
  2. ^ Ahmed, N .; Rao, KR (10 de julio de 1975). Transformaciones ortogonales para el procesamiento de señales digitales (1 ed.). Berlín Heidelberg Nueva York: Springer-Verlag. doi : 10.1007 / 978-3-642-45450-9 . ISBN 9783540065562.
  3. ^ a b c Prandoni, Paolo; Vetterli, Martin (2008). Procesamiento de señales para comunicaciones (PDF) (1 ed.). Boca Raton, FL: CRC Press. págs. 72, 76. ISBN  978-1-4200-7046-0. Consultado el 4 de octubre de 2020 . los coeficientes DFS para la señal periodizada son un conjunto discreto de valores para su DTFT
  4. ^ Rao, R. (2008). Señales y Sistemas . Prentice-Hall Of India Pvt. Limitado. ISBN 9788120338593.
  5. ^ "Estimación de la densidad espectral de potencia del periodograma - periodograma de MATLAB" .
  6. ^ Gumas, Charles Constantine (julio de 1997). "Window-presum FFT logra un alto rango dinámico, resolución" . Noticias de ingeniería personal e instrumentación : 58–64. Archivado desde el original el 10 de febrero de 2001.CS1 maint: bot: original URL status unknown (link)
  7. ^ Crochiere, RE; Rabiner, LR (1983). "7,2". Procesamiento de señales digitales multivelocidad . Englewood Cliffs, Nueva Jersey: Prentice-Hall. págs. 313–326. ISBN 0136051626.
  8. ^ Wang, Hong; Lu, Youxin; Wang, Xuegang (16 de octubre de 2006). "Receptor canalizado con WOLA Filterbank". 2006 CIE International Conference on Radar . Shanghái, China: IEEE: 1–3. doi : 10.1109 / ICR.2006.343463 . ISBN 0-7803-9582-4.
  9. ^ Lyons, Richard G. (junio de 2008). "Trucos DSP: construcción de un analizador de espectro práctico" . EE Times . Consultado el 20 de febrero de 2020 .   Sin embargo, tenga en cuenta que contiene un enlace etiquetado como estructura de superposición-adición ponderada que va incorrectamente al método Overlap-add .
  10. ↑ a b Lillington, John (marzo de 2003). "Comparación de arquitecturas de canalización de banda ancha" (PDF) . Dallas: Conferencia internacional de procesamiento de señales. pag. 4 (figura 7). Archivado desde el original (PDF) el 8 de marzo de 2019 . Consultado el 6 de septiembre de 2020 . El "Weight Overlap and Add" o WOLA o su subconjunto el "Polyphase DFT", está cada vez más establecido y es ciertamente muy eficiente donde se requieren bancos de filtros grandes y de alta calidad.
  11. ^ a b Lillington, John. "Una revisión de las técnicas del banco de filtros - RF y digital" (PDF) . armms.org . Isla de Wight, Reino Unido: Libra Design Associates Ltd. p. 11 . Consultado el 6 de septiembre de 2020 . Afortunadamente, existe una solución mucho más elegante, como se muestra en la Figura 20 a continuación, conocida como Polyphase o WOLA (Weight, Overlap and Add) FFT.
  12. ^ Hochgürtel, Stefan (2013). "Implementaciones eficientes de espectrómetros FFT de banda ancha de alta resolución y su aplicación a un levantamiento de línea del Centro Galáctico APEX" (PDF) . hss.ulb.uni-bonn.de . Bonn: Universidad Rhenish Friedrich Wilhelms de Bonn. págs. 26-27 . Consultado el 6 de septiembre de 2020 . Para realizar M-fold WOLA para una DFT de N puntos, M · N muestras de entrada reales a j primero multiplicadas por una función de ventana w j del mismo tamaño
  13. Chennamangalam, Jayanth (18 de octubre de 2016). "La técnica del banco de filtros polifásicos" . Grupo CASPER . Consultado el 30 de octubre de 2016 .
  14. Dahl, Jason F. (6 de febrero de 2003). Métodos de estimación de espectro de alias en el tiempo (Ph.D.). Universidad Brigham Young . Consultado el 31 de octubre de 2016 .
  15. ^ Lin, Yuan-Pei; Vaidyanathan, PP (junio de 1998). "Un enfoque de ventana de Kaiser para el diseño de filtros prototipo de bancos de filtros modulados con coseno" (PDF) . Cartas de procesamiento de señales IEEE . 5 (6): 132-134. Código bibliográfico : 1998ISPL .... 5..132L . doi : 10.1109 / 97.681427 . Consultado el 16 de marzo de 2017 .
  16. Harris, Frederic J. (24 de mayo de 2004). "9". Procesamiento de señales multivelocidad para sistemas de comunicación . Upper Saddle River, Nueva Jersey: Prentice Hall PTR. págs. 226-253. ISBN 0131465112.
  17. ^ Harris, Fredric J. (enero de 1978). "Sobre el uso de Windows para el análisis armónico con la transformada discreta de Fourier" (PDF) . Actas del IEEE . 66 (1): 51–83. Código bibliográfico : 1978IEEEP..66 ... 51H . CiteSeerX 10.1.1.649.9880 . doi : 10.1109 / PROC.1978.10837 .  
  18. ^ a b c d e f g h i j k l m n o p q r Proakis, John G .; Manolakis, Dimitri G. (1996). Procesamiento de señales digitales: principios, algoritmos y aplicaciones (3 ed.). Nueva Jersey: Prentice-Hall International. Bibcode : 1996dspp.book ..... P . ISBN 9780133942897. sAcfAQAAIAAJ.
  19. ^ Rabiner, Lawrence R .; Oro, Bernard (1975). Teoría y aplicación del procesamiento de señales digitales . Englewood Cliffs, Nueva Jersey: Prentice-Hall, Inc. p. 59 (2,163). ISBN 978-0139141010.

Lectura adicional [ editar ]

  • Porat, Boaz (1996). Curso de Procesamiento de Señales Digitales . John Wiley e hijos. págs. 27-29 y 104-105. ISBN 0-471-14961-6.
  • Siebert, William M. (1986). Circuitos, señales y sistemas . Serie de Ingeniería Eléctrica e Informática del MIT. Cambridge, MA: MIT Press. ISBN 0262690950.
  • Lyons, Richard G. (2010). Comprensión del procesamiento de señales digitales (3ª ed.). Prentice Hall. ISBN 978-0137027415.