¿Por qué el DFT supone que la señal transformada es periódica?

10

En muchos libros de procesamiento de señales, se afirma que el DFT supone que la señal transformada es periódica (y que esta es la razón por la cual puede ocurrir una fuga espectral, por ejemplo).

Ahora, si nos fijamos en la definición de DFT, simplemente no hay ese tipo de suposición. Sin embargo, en el artículo de Wikipedia sobre la transformada de Fourier de tiempo discreto (DTFT), se afirma que

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

  • Entonces, ¿esta suposición se deriva de la DTFT?
  • En realidad, al calcular el DFT, ¿estoy calculando el DTFT con el supuesto de que la señal es periódica?
usuario10839
fuente
Debido a que DFT X [k] de x [n] es el primer período de la serie discreta de Fourier (DFS) de la señal periódica xp [n] cuyo primer período se toma como x [n]
Fat32
1
Parece que tendré que escribir una respuesta disidente a esto. el DFT supone que la señal transformada es periódica porque está ajustando un conjunto de funciones básicas a la señal transformada, todas las cuales son periódicas.
robert bristow-johnson
1
El DFT es solo la expresión simplificada del DFS, por lo tanto, el supuesto periódico existe inherentemente.
lxg

Respuestas:

12

Ya hay algunas buenas respuestas, pero todavía tengo ganas de agregar otra explicación, porque considero que este tema es extremadamente importante para comprender muchos aspectos del procesamiento de señales digitales.

En primer lugar, es importante comprender que el DFT no "asume" la periodicidad de la señal a transformar. El DFT se aplica simplemente a una señal finita de longitud N y los coeficientes DFT correspondientes se definen por

(1)X[k]=n=0N1x[n]ej2πnk/N,k=0,1,,N1

A partir de (1) es obvio que solo se consideran muestras de x[n] en el intervalo [0,N1] , por lo que no se supone periodicidad. Por otro lado, los coeficientes X[k] pueden interpretarse como coeficientes de Fourier de la continuación periódica de la señal x[n] . Esto se puede ver desde la transformación inversa

(2)x[n]=k=0N1X[k]ej2πnk/N

que calcula correctamente en el intervalo [ 0 , N - 1 ] , pero también calcula su continuación periódica fuera de este intervalo debido a que el lado derecho de (2) es periódica con periodo N . Esta propiedad es inherente a la definición del DFT, pero no tiene por qué molestarnos porque normalmente solo nos interesa el intervalo [ 0 , N - 1 ] .x[n][0,N1]N[0,N1]

Considerando la DTFT de x[n]

(3)X(ω)=n=x[n]ejnω

podemos ver al comparar (3) con (1), que si es una secuencia finita en el intervalo [ 0 , N - 1 ] , los coeficientes DFT X [ k ] son muestras de DTFT X ( ω ) :x[n][0,N1]X[k]X(ω)

(4)X[k]=X(2πk/N)

Entonces, un uso del DFT (pero ciertamente no el único) es calcular muestras del DTFT. Pero esto solo funciona si la señal a analizar es de longitud finita . Por lo general, esta señal de longitud finita se construye mediante una ventana más larga. Y es esta ventana la que causa la fuga espectral.

Como último comentario, tenga en cuenta que la DTFT de la continuación periódica de la secuencia finita x [ n ] puede expresarse en términos de los coeficientes DFT de x [ n ] :x~[n]x[n]x[n]

˜ X (ω)=2π

(5)x~[n]=k=x[nkN]
(6)X~(ω)=2πNk=X[k]δ(ω2πk/N)

EDITAR: El hecho de que y ˜ X ( ω ) dados anteriormente son un par de transformada DTFT se puede mostrar de la siguiente manera. Primera nota que el DTFT de un peine de impulso de tiempo discreto es un peine de Dirac:x~[n]X~(ω)

(7)k=δ[nkN]2πNk=δ(ω2πk/N)

La secuencia se puede escribir como la convolución de x [ n ] con un peine de impulsos:x~[n]x[n]

(8)x~[n]=x[n]k=δ[nkN]

Como la convolución corresponde a la multiplicación en el dominio DTFT, la DTFT de ˜ x [ n ] viene dada por la multiplicación de X ( ω ) con un peine Dirac:X~(ω)x~[n]X(ω)

(9)X~(ω)=X(ω)2πNk=δ(ω2πk/N)=2πNk=X(2πk/N)δ(ω2πk/N)

La combinación de con ( 4 ) establece el resultado ( 6 ) .(9)(4)(6)

Matt L.
fuente
Bajé la flecha con esta respuesta por la misma razón por la que tengo la respuesta más reciente de @ hotpaw2. en esta declaración: "De (1) es obvio que solo se consideran muestras de en el intervalo [ 0 , N - 1 ] , por lo que no se asume ninguna periodicidad". x[n][0,N1]La conclusión no se sigue de la premisa.
robert bristow-johnson
44
@ robertbristow-johnson: lo hace. Dame muestras consecutivas y te doy el DFT. No necesito asumir nada sobre la señal fuera del rango [ 0 , N - 1 ] , ni siquiera su existencia. Esto es lo único que reclamo en esa oración, y obviamente es cierto. Para calcular el DFT no necesito saber nada excepto los valores en el intervalo [ 0 , N - 1 ] . No estoy seguro de cómo podría malinterpretar o leer mal mi declaración. Si se trata de un problema de formulación, me alegraría aclarar mi oración, pero en términos de contenido es realmente trivial. N[0,N1][0,N1]
Matt L.
44
lea la otra respuesta a continuación y mi respuesta en el otro hilo. no se trata de lo que supone sobre fuera de 0 n N - 1 . se trata de lo que "asume" la transformación (si se nos permite antropomorfizar un poco) acerca de x [ n ] fuera de 0 n N - 1 . podemos averiguar qué supone la transformación cuando invocamos una operación en un dominio que desplaza al otro dominio en una cantidad entera. x[n]0nN1x[n]0nN1
robert bristow-johnson
@MattL. (9) debería leer lugar de=2π
=2πNk=X[k]δ(ω2πk/N)
=2πNk=X(2πk/N)δ(ω2πk/N)
jomegaA
@jomegaA: No en ambos casos. Como se indicó en la última oración de mi respuesta, el resultado final (6) se concluye de la combinación de (9) con (4), por lo que, por supuesto, , pero en (9) se deriva de la DTFT X ( ω ) . Y con respecto al factor de escala 2 π / N , definitivamente necesita estar allí. No confunda las expresiones con ω y f , tienen diferentes factores de escala. X[k]=X(2πk/N)X(ω)2π/Nωf
Matt L.
8

Proviene de la definición de la señal del dominio del tiempo:

x[n]=k=0N1X[k]e2πinkN

Puede ver por definición que x[n]=x[n+N] .
Por otro lado, el DFT reconstruye perfectamente las N muestras de la señal.
Por lo tanto, puede concluir que supone una continuación periódica de la misma.

Otro punto de vista sería mirar el DFT como una serie de Fourier discreta finita (en realidad es, eche un vistazo a la serie de Fourier discreta - DFS ), que por supuesto señala que la señal es periódica (la suma finita de señales con período T es una señal que tiene un período T ).

Royi
fuente
2
No veo cómo proviene de la definición.
user10839
1
@ user10839: simplemente evalúe y verá que es igual a x [ n ] . Como se señaló en la respuesta, el DFT es solo una serie de Fourier de la señal del dominio del tiempo. La longitud finita de la señal en el dominio del tiempo se considera el período fundamental. x[n+N]x[n]
Matt L.
@ user10839, solo conéctelo a la ecuación. El exponente se puede definir con las funciones coseno y seno, que como se puede ver tienen un período de . nkN
Royi
1
DFT no es el DFS. Esto es pedante, pero DFT le ofrece los coeficientes de la serie de Fourier. Es importante tener en cuenta que DFT es como cualquier otra transformación lineal. Es una multiplicación matricial. La matriz es ortonormal, lo que la hace agradable. También se puede demostrar que los coeficientes de salida equivalen a la correspondiente expansión de la serie de Fourier de los datos, pero la transformada de Fourier no es la serie de Fourier (desajuste de tipo: p).
thang
@thang, no tengo idea de lo que quieres decir. El DFT es el DFS. Ellos son iguales. Es fácil ver eso. Preste atención, esta es la serie discreta de Fourier y no la serie de Fourier (con integrales). Eche un vistazo aquí en.wikipedia.org/wiki/Discrete_Fourier_series y vea que es DFT.
Royi
5

Es una suposición innecesaria (y a menudo falsa). El DFT es solo una transformación base de un vector finito.

Los vectores básicos de la DFT son fragmentos de funciones periódicas infinitamente extensibles. Pero no hay nada intrínsecamente periódico sobre la entrada o resultados DFT a menos que extienda los vectores base fuera de la apertura DFT. Muchas formas de análisis de señal no requieren ninguna extensión o suposición fuera de la ventana muestreada o del vector de datos finitos.

También se puede suponer que cualquier artefacto de "fuga" proviene de una convolución de la ventana rectangular predeterminada con una señal que no es periódica o es de periodicidad o estacionariedad desconocida. Esto tiene mucho más sentido al analizar ventanas FFT superpuestas, donde cualquier supuesto de periodicidad fuera de cualquier ventana DFT o FFT puede ser inconsistente con los datos en otras ventanas.

La periodicidad puede hacer que las matemáticas que relacionan el DFT con el DTFT sean más manejables. Pero cualquier relación con la DTFT puede o no ser necesaria cuando realmente se utiliza una FFT para el procesamiento de la señal (dependiendo exactamente de qué propiedades de transformación de Fourier se necesitan para un análisis posterior del método de procesamiento).

hotpaw2
fuente
Flecha hacia abajo por la misma razón por la que bajé su respuesta más reciente sobre esto.
robert bristow-johnson
5

Ok, mi respuesta será algo diferente a las otras respuestas. mi respuesta acepta la premisa de la pregunta en lugar de negar la premisa de la pregunta.

La razón por la que el DFT "asume" la señal de entrada (la señal a transformar, lo que supongo que OP significa "señal transformada") es periódica es porque el DFT ajusta una colección de funciones básicas a esa señal de entrada, todo lo cual son periódicos

considere un conjunto diferente de funciones básicas:

gk(u)uk0k<N

N

x[n]0n<N

gk(n)

x[n]=k=0N1X[k]gk(n)=k=0N1X[k]nk

with judicious selection of the coefficients X[k]. calculating all X[k] requires solving N linear equations with N unknowns. you can use Gaussian elimination to do it.

with the N correct values for X[k] for 0kN1, we can make sure that the sum of these power functions (which is an (N1)th-order polynomial) will evaluate exactly to x[n] for each n such that 0nN1.

now what if you use that summation to go beyond the interval of 0nN1? you can evaluate it for any n. you will notice that the behavior of that function will be that of an (N1)th-order polynomial because that is what it is. for n large enough, only the highest power with a non-zero coefficient will set the trend for the extrapolated x[n].

so now, with the DFT we are fitting a different set of basis functions to our input sequence:

gk(u)1Ne+j2πku/N0k<N

x[n]=k=0N1X[k]gk(n)=1Nk=0N1X[k]e+j2πnk/N

and the coefficients, X[k], can be solved for and are:

X[k]=n=0N1x[n] ej2πnk/N

the placement of that 1N is a matter of convention. i am putting it where most of the literature puts the 1N factor. it could be removed from the x[n] equation and put in the X[k] equation, instead. or "half" of it (1N) could be placed with both equations. it's just a matter of convention.

but here we are fitting a set of basis functions that are all periodic with period N to the original x[n]. so even if x[n] came from a longer sequence was not periodic, the DFT is considering that x[n] is the sum of a bunch of basis functions each that are periodic with period N. if you add up a bunch of periodic functions, all with the same period, the sum must also be periodic with the same period.

robert bristow-johnson
fuente
para un poco más polémico, donde disputo la noción de que el DFT no necesariamente extiende periódicamente los datos que se le pasan, por favor, mire esta respuesta anterior mía . Prefiero no repetirlo aquí.
robert bristow-johnson
1

DFT is discrete. DTFT is continuous. We can get DFT from DTFT by sampling it with the pulse train of the right period, which is actually equal to multiplying it with the pulse train. Multiplication in the transform domain is equal to convolution in discrete-time domain, this implies periodicity of signal.

learner
fuente
DTFT is continuous? How come?
jojek
2
The result of the DTFT is continous (in frequency).
Deve
Indeed - thus you should state it clearly to avoid any misunderstanding and supply adequate equations.
jojek
@jojek Eso es cierto, también creo que esta respuesta podría mejorarse con algunas ecuaciones
Deve
1
Agregaré más detalles muy pronto.
aprendiz
0

Solo DFT es práctico en el mundo digital discreto debido a suposiciones periódicas en ambos dominios. (Si lo llama así.) Debido a que la señal no periódica en un dominio causa una señal continua en el otro y solo puede almacenar una señal discreta en la memoria digital. Por lo tanto, debe asumir que las señales son periódicas en ambos dominios para que sea discreto en ambos dominios.

Cuando calcula DTFT, obtiene una señal continua en el dominio de frecuencia como salida.
No creo que utilice el mismo procedimiento cuando calcule DFT en la práctica. Cuando realmente calculó DTFT y DFT, comprenderá que ambos cálculos de transformación son historias diferentes.

Sam
fuente
0

Since the signal is periodic, the time shifted signal doesn't change the absolute magnitude of the frequency domain.

X[k]=k=0N1x[n]e2πinkN

e2πiDkNX[k]=k=0N1x[nD]e2πinkNe2πiDkN

By the way, there is nothing stopping you from taking the FFT of a non-periodic signal, but there it little practical use if none of the transformations work.

FreedomToWin
fuente