En matemáticas, la transformada de Fourier de tiempo discreto, del inglés discrete-time Fourier transform (DTFT) es una forma de análisis de Fourier que es aplicable a una secuencia de valores.

La DTFT se usa frecuentemente para analizar un muestreo discreto de una función continua. El término tiempo discreto se refiere al hecho de que se trabaja sobre segmentos discretos de una función, comúnmente se trata de segmentos de intervalos de tiempo. De segmentos de función de intervalos fijos de tiempo se produce una función con dominio en la frecuencia que es la sumatoria periódica de la transformada de Fourier continua. Bajo ciertas condiciones teoréticas, descritas en el teorema de muestreo de Nyquist-Shannon, la función continua original puede ser reconcstruida perfectamente desde la DTFT y por lo tanto desde las muestras discretas originales. La DTFT en sí misma es una función continua con dominio en la frecuencia, pero las muestras discretas de la misma pueden ser calculadas con la transformada de Fourier discreta (DFT), que es, por mucho, el más usual de lo métodos en el anális de Fourier moderno.

Ambas transformadas son invertibles. La inversa de la DTFT es la función muestreada original. La inversa de la DFT es la suma periódica de la función original. La Transformada rápida de Fourier (FFT) es un algoritmo para calcular un ciclo de la DFT, y su inversa es un ciclo de la inversa de la DFT.

Definición

editar

La transformada de Fourier en tiempo discreto de secuencias de números reales o complejos x[n], para n miembro de los Enteros, es una serie de Fourier que produce una función periódica de frecuencia variable. Cuando la frecuencia variable, ω, está normalizada a radianes/muestra, la periodicidad es , y la serie de Fourier es:[1]: p.147 

 

 

 

 

 

(Eq.1)

La utilidad de esta función en el dominio de la frecuencia está basada en la en:Poisson summation formula. 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], ej. Tx(nT) = x[n].[2]​  Entonces la función periódica representada por las series de Fourier es una sumatoria periódica de X(f) en términos de su frecuencia f en hertz (ciclos/segundo):Oppenheim y Schafer,[1]​ p 147 (4.20), p 694 (10.1), y Prandoni y Vetterli,[3]​ p 255, (9.33), donde:         y   

Cuando la dependencia en T es importante, es práctica común reemplazarla con    Entonces  f  tiene las unidades de (ciclos/muestra), llamadas normalized frequency.

 

 

 

 

 

(Eq.2)

 
Fig 1. Ilustración de una transformada de Fourier (arriba a la izquierda) y su suma periódica (DTFT) abajo a la izquierda. La esquina inferior izquierda nos enseña muestras de DTFT que fueron calculadas por la transformada discreta de Fourier (DFT).

El entero k tiene las unidades de ciclos/muestra, y 1/T es la rapidez del muestreo, fs (muestras/segundo). Así X1/T(f) contiene copias exactas de X(f) que están desplazadas por múltiplos de fs hertz y combinadas por suma. Para un fs lo suficientemente grande el término k = 0 puede ser observado en la región [−fs/2, fs/2] con poca o ninguna distorsión por( solapamiento sincrónico) con otros términos. En Fig.1, los extremos de la distribución en la esquina superior izquierda están enmascaradas por solapamiento sincrónico en la sumatoria periódica (abajo a la izquierda).

También nótese que ei2πfTn es la transformada de Fourier de δ(tnT). Por lo tanto, una definición alternativa al DTFT es: De hecho la Eq.2 seguido se justifica como sigue:[1]: p.143 

 


 

 

 

 

 

(Eq.3)

La función modulada (un peine de Dirac) es una abstracción matemática algunas veces referida como tren de impulsos o función de muestreo.[4]

Transformada Inversa

editar

Una operación que reconstruye la secuencia de datos discreta desde la función DTFT se llama una DTFT inverso. Por ejemplo, la transformada de Fourier inversa contínua de ambos lados de {EquationNote|Eq.3}} produce una secuencia con la forma de función de peine de Dirac modulado:

 

Sin embargo, notando que X1/T(f) es periódica, toda la información necesaria está contenida dentro de cualquier intervalo de longitud 1/T. En Eq.1 y Eq.2, las sumatorias sobre n son series de Fourier, con coeficientes x[n]. Las fórmulas estándar para los coeficientes de Fourier son también las transformadas inversas:

 

 

 

 

 

(Eq.4)

Funciones Periódicas

editar

Cuando la secuencia de entrada x[n] es N-periodica, Eq.2 puede ser computacionalmente reducida a una transformada de Fourier discreta (DFT), porque:

  • Toda la información puede ser contenida dentro N muestras.
  • X1/T(f) converge a cero en todas partes excepto en los múltiplos enteros 1/(NT), conocidos como frecuencias armónicas (armónicos). A esas frecuencias, la DTFT diverge a diferentes velocidades dependientes de la frecuencia. Y esas frecuencias están dadas por la DFT de un ciclo de la secuencia x[n].
  • La DTFT es periódic, por lo que el máximo número de amplitudes armónicas es (1/T) / (1/(NT)) = N

Los coeficientes de la DFT están dados por:

      y la DTFT es:
       Oppenheim y Schafer,[1]​ p 551 (8.35), y Prandoni y Vetterli,[3]​ p 82, (4.43), donde:           y   


Substituyendo esta expresión en la fórmula de la transformada inversa confirma:

  (all integers)

como se esperaba. La inversa de la DFT en la línea de arriba a veces es referida como en:Discrete Fourier series (DFS).[1]: p 542 

Muestreando la 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 X1/T: [1]​pp 557–559 & 703

 

donde   es una en:periodic summation:

      (ver en:Discrete Fourier series)

La secuencia   es la inversa de la DFT. Por lo que nuestro muestreo de la DTFT causa que la transformada inversa se vuelva periódica. El arreglo de valores   conocidos como en:periodogram, y el parámetro N es llamado NFFT en la función de Matlab del mismo nombre.[5]

A fin de evaluar un ciclo de   numéricamente, se requiere la secuencia de longitud finita x[n]. Por ejemplo, una secuencia larga puede ser truncada por una función ventana de longitud L resultando los 3 casos siguientes, (para simplificar la notación, consideraremos que los valores x[n] abajo representan los valores modificados por la función de ventana):

Caso 1: Agotamiento de la frecuencia. L = NI, para algún entero I (tipicamente 6 u 8)

Un ciclo de   se reduce a una sumatoria de I segmentos de longitud N.  Entonces la DFT se denomina por varios nombres como:

Ejemplo WOLA: File:WOLA channelizer example.png


  • DFT de polifase[10][11]
  • Batería de filtrado polifase[13]
  • Ventaneo de bloque múltiple y reforzado con sobremuestreo.[14]

Hay que recordar que el agotamiento de los datos muestreados en un dominio (del tiempo o frecuencia) produce solapamiento (a veces conocido como aliasing) en el otro y viceversa. Comparado a una DFT de L-largo, la sumatoria/traslape   causa el agotamiento en frecuencia,[1]: p.558  dejando solo las muestras de la DTFT menos afectadas por manchado espectral o leakage. Esto es prioritario cuando se implementa un banco de filtros (canalizador) de FFT. Con una función de ventana de longitud L, el manchado espectral sería inaceptable. Así que multibloques de ventanas son creados usando herramientas de diseño en:FIR filter.[15][16]​  Su perfil de frecuencia es plano en el punto más alto y cae rápidamente a la mitad de las muestras DTFT restantes. Mientras más grande sea el valor del parámetro I, más probabilidades de un mejor comportamiento.

Caso: L = N+1.

Cuando es simétrica la ventana (función) de longitud L, el coeficiente ( ) truncado por 1 se llama periódica o DFT-suave. El truncamiento afecta a la DTFT.  Una DFT de muestras de secuencia truncada de la DTFT a intervalos de frecuencia de 1/N. To sample   at the same frequencies, for comparison, the DFT is computed for one cycle of the periodic summation,  An example is figure Sampling the DTFT. La DFT de muestras con reales son el resultado de la simetría de DFT-suave symmetry[17]: p.52 

 
Fig 2. DFT of ei2πn/8 for L = 64 and N = 256
 
Fig 3. DFT of ei2πn/8 for L = 64 and N = 64

Caso: Interpolación de frecuencia. LN

In this case, the DFT simplifies to a more familiar form:

 

In order to take advantage of a fast Fourier transform algorithm for computing the DFT, the summation is usually performed over all N terms, even though NL of them are zeros. Therefore, the case L < N is often referred to as zero-padding.

Spectral leakage, which increases as L decreases, is detrimental to certain important performance metrics, such as resolution of multiple frequency components and the amount of noise measured by each DTFT sample. But those things don't always matter, for instance when the x[n] sequence is a noiseless sinusoid (or a constant), shaped by a window function. Then it is a common practice to use zero-padding to graphically display and compare the detailed leakage patterns of window functions. To illustrate that for a rectangular window, consider the sequence:

  and  

Figures 2 and 3 are plots of the magnitude of two different sized DFTs, as indicated in their labels. In both cases, the dominant component is at the signal frequency: f = 1/8 = 0.125. Also visible in Fig 2 is the spectral leakage pattern of the L = 64 rectangular window. The illusion in Fig 3 is a result of sampling the DTFT at just its zero-crossings. Rather than the DTFT of a finite-length sequence, it gives the impression of an infinitely long sinusoidal sequence. Contributing factors to the illusion are the use of a rectangular window, and the choice of a frequency (1/8 = 8/64) with exactly 8 (an integer) cycles per 64 samples. A Hann window would produce a similar result, except the peak would be widened to 3 samples (see DFT-even Hann window).

Convolución

editar

Art principal VVV Teorema de convolución#Functions of a discrete variable (sequences) El teorema de convolución para secuencias es:

 [18]: p.297 , Oppenheim and Schafer,[1]​ p 60, (2.169), y Prandoni and Vetterli,[3]​ p 122, (5.21)

An important special case is the en:circular convolution of sequences x and y defined by   where   is a periodic summation. The discrete-frequency nature of   means that the product with the continuous function   is also discrete, which results in considerable simplification of the inverse transform:

 [19][1]: p.548 

For x and y sequences whose non-zero duration is less than or equal to N, a final simplification is:

 

The significance of this result is explained at Circular convolution and Fast convolution algorithms.

Symmetry properties

editar

When the real and imaginary parts of a complex function are decomposed into their even and odd parts, there are four components, denoted below by the subscripts RE, RO, IE, and IO. And there is a one-to-one mapping between the four components of a complex time function and the four components of its complex frequency transform:[18]: p.291 

 

From this, various relationships are apparent, for example: *The transform of a real-valued function ({{math|''x''{{sub|{{sub|RE}}}}+ ''x''{{sub|{{sub|RO}}}}}}) is the [[Even and odd functions#Complex-valued functions|even symmetric]] function {{math|X{{sub|{{sub|RE}}}}+ i X{{sub|{{sub|IO}}}}}}. Conversely, an even-symmetric transform implies a real-valued time-domain. *The transform of an imaginary-valued function ({{math|i ''x''{{sub|{{sub|IE}}}}+ i ''x''{{sub|{{sub|IO}}}}}}) is the [[Even and odd functions#Complex-valued functions|odd symmetric]] function {{math|X{{sub|{{sub|RO}}}}+ i X{{sub|{{sub|IE}}}}}}, and the converse is true. *The transform of an even-symmetric function ({{math|''x''{{sub|{{sub|RE}}}}+ i ''x''{{sub|{{sub|IO}}}}}}) is the real-valued function {{math|X{{sub|{{sub|RE}}}}+ X{{sub|{{sub|RO}}}}}}, and the converse is true. *The transform of an odd-symmetric function ({{math|''x''{{sub|{{sub|RO}}}}+ i ''x''{{sub|{{sub|IE}}}}}}) is the imaginary-valued function {{math|i X{{sub|{{sub|IE}}}}+ i X{{sub|{{sub|IO}}}}}}, and the converse is true.

Relationship to the Z-transform

editar

  is a Fourier series that can also be expressed in terms of the bilateral Z-transform.  I.e.:

 

where the   notation distinguishes the Z-transform from the Fourier transform. Therefore, we can also express a portion of the Z-transform in terms of the Fourier transform:

 

Note that when parameter T changes, the terms of   remain a constant separation   apart, and their width scales up or down. The terms of X1/T(f) remain a constant width and their separation 1/T scales up or down.

Table of discrete-time Fourier transforms

editar

Some common transform pairs are shown in the table below. The following notation applies:

  •   is a real number representing continuous angular frequency (in radians per sample). (  is in cycles/sec, and   is in sec/sample.) In all cases in the table, the DTFT is 2π-periodic (in  ).
  •   designates a function defined on  .
  •   designates a function defined on  , and zero elsewhere. Then:
 
  •   is the Dirac delta function
  •   is the normalized sinc function
  •   is the rectangle function
  •   is the triangle function
  • n is an integer representing the discrete-time domain (in samples)
  •   is the discrete-time unit step function
  •   is the Kronecker delta  
Time domain
x[n]
Frequency domain
X(ω)
Remarks Reference
    [18]: p.305 
    integer  
   

      odd M
      even M

integer  
   

 

The   term must be interpreted as a distribution in the sense of a Cauchy principal value around its poles at  .
      [18]: p.305 
        -π < a < π

 

real number  
   

 

real number   with  
    real number   with  
    integers   and  
    real numbers   with  
    real number  ,  
    it works as a differentiator filter
    real numbers   with  
   
    Hilbert transform
     real numbers  
complex  

Properties

editar

This table shows some mathematical operations in the time domain and the corresponding effects in the frequency domain.

Property Time domain
x[n]
Frequency domain
 
Remarks Reference
Linearity     complex numbers   [18]: p.294 
Time reversal / Frequency reversal     [18]: p.297 
Time conjugation     [18]: p.291 
Time reversal & conjugation     [18]: p.291 
Real part in time     [18]: p.291 
Imaginary part in time     [18]: p.291 
Real part in frequency     [18]: p.291 
Imaginary part in frequency     [18]: p.291 
Shift in time / Modulation in frequency     integer k [18]: p.296 
Shift in frequency / Modulation in time     real number   [18]: p.300 
Decimation      La plantilla {{Fn}} ya no debe usarse. En su lugar usa el nuevo sistema de referencias.. integer  
Time Expansion     integer   [1]: p.172 
Derivative in frequency     [18]: p.303 
Integration in frequency    
Differencing in time    
Summation in time    
Convolution in time / Multiplication in frequency     [18]: p.297 
Multiplication in time / Convolution in frequency     Periodic convolution [18]: p.302 
Cross correlation    
Parseval's theorem     [18]: p.302 

Ver también

editar
  1. a b c d e f g h i j Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). «4.2, 8.4». Discrete-time signal processing (2nd edición). Upper Saddle River, N.J.: Prentice Hall. ISBN 0-13-754920-2. (requiere registro). «samples of the Fourier transform of an aperiodic sequence x[n] can be thought of as DFS coefficients of a periodic sequence obtained through summing periodic replicas of x[n].»   url=https://d1.amobbs.com/bbs_upload782111/files_24/ourdev_523225.pdf
  2. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Ahmed
  3. a b c Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Prandoni
  4. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Rao
  5. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Matlab
  6. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Gumas
  7. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Crochiere
  8. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Wang
  9. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Lyons
  10. a b Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Lillington
  11. a b Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Lillington2
  12. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Hochgürtel
  13. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Chennamangalam
  14. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Dahl
  15. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Lin
  16. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Harrisbook
  17. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Harris
  18. a b c d e f g h i j k l m n ñ o p q Proakis, John G.; Manolakis, Dimitri G. (1996). Digital Signal Processing: Principles, Algorithms and Applications (en inglés) (3 edición). New Jersey: Prentice-Hall International. Bibcode:1996dspp.book.....P. ISBN 9780133942897. sAcfAQAAIAAJ. (requiere registro). 
  19. Error en la cita: Etiqueta <ref> no válida; no se ha definido el contenido de las referencias llamadas Rabiner

Referencias

editar

Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Proakis» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Rabiner» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Oppenheim» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Crochiere» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Rao» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Matlab» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Gumas» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Lillington» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Lillington2» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Chennamangalam» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Lyons» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Dahl» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Wang» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Lin» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Harrisbook» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Harris» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Hochgürtel» no se utiliza en el texto anterior.
Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Ahmed» no se utiliza en el texto anterior.

Error en la cita: La etiqueta <ref> definida en las <references> con nombre «Prandoni» no se utiliza en el texto anterior.

Further reading

editar
  • Porat, Boaz (1996). A Course in Digital Signal Processing. John Wiley and Sons. pp. 27-29 and 104-105. ISBN 0-471-14961-6. 
  • Siebert, William M. (1986). Circuits, Signals, and Systems. MIT Electrical Engineering and Computer Science Series. Cambridge, MA: MIT Press. ISBN 0262690950. 
  • Lyons, Richard G. (2010). Understanding Digital Signal Processing (3rd edición). Prentice Hall. ISBN 978-0137027415.