En las matemáticas, el que Fourier discreto transforma (DFT) es una de las formas específicas del análisis de Fourier . Como tal, transforma una función en otra, que se llama la representación del dominio de frecuencia, o simplemente el DFT, de la función original (que es a menudo una función en el Tiempo-dominio ). Pero el DFT requiere una función de entrada que sea el '' discreto '' y cuyos valores diferentes a cero tengan ( finito) una duración limitada. Tales entradas son creadas a menudo por el muestreo una función continua, como la voz de una persona. Y desemejante DTFT, evalúa solamente bastantes componentes de la frecuencia para reconstruir el segmento finito que era analizado. Su lo contrario transforma no puede reproducir el tiempo-dominio entero, a menos que la entrada suceda ser periódica (por siempre). Por lo tanto se dice a menudo que el DFT es una transformación para el análisis de Fourier del finito-dominio, tiempo discreto funciona. Las funciones de base sinusoidales de la descomposición tienen las mismas características.

Puesto que la función de entrada es una secuencia finita verdadero o los números complejos el DFT es ideal para la información de proceso almacenada en las computadoras particularmente, el DFT se emplea extensamente en el tratamiento de señales y los campos relacionados para analizar las frecuencias contenidas en un muestreado señalan, para solucionar las ecuaciones diferenciales parciales, y realizar otras operaciones tales como circunvoluciones El DFT se puede computar eficientemente en la práctica usar un que Fourier rápido transforma algoritmo de (FFT).

Puesto que los algoritmos de FFT se emplean tan comúnmente para computar el DFT, los dos términos son de uso frecuente alternativamente en ajustes familiares, aunque haya una distinción clara: " DFT" refiere a una transformación matemática, sin importar cómo se computa, mientras que " FFT" refiere de varios algoritmos eficientes para el DFT. Esta distinción se empaña más a fondo, sin embargo, por el Fourier finito del sinónimo transformar para el DFT, que precede al parecer el " del término; ayuna el transform" de Fourier; (Cooley y otros, 1969) pero tiene el mismo Initialism .

Definición

La secuencia del x 0 de los números complejos N ,…, el N −1 del del x es transformada en la secuencia del X 0 de los números complejos del N ,…, el N −1 del del X por el DFT según la fórmula:

X_k = \ el sum_ {n=0} el ^ {N-1} x_n. ^ {- \} \ patio \ patio del frac {2 \ pi i} {N} k n k = 0, \ puntea, N-1 donde está la base el e del logaritmo natural, i \, es la unidad imaginaria (i^2=-1), y el π es el pi . Transformación es a veces denotado por símbolo \ mathcal {F} , como en \ mathbf {X} = \ mathcal {F} \ a la izquierda \ {\ mathbf {x} \ derecho \} o \ mathcal {F} \ a la izquierda (\ mathbf {x} \ derecho) o \ mathcal {} \ mathbf {x} de F.

El que Fourier discreto inverso transforma (IDFT) se da cerca

x_n = \ el frac {1} {N} \ el sum_ {k=0} el ^ {N-1} X_k e^ {\} \ patio \ patio del frac {2 \ pi i} {N} k n n = 0, \ puntea, N-1.

Una descripción simple de estas ecuaciones es que los números complejos X_k representan la amplitud y la fase de los diversos componentes sinusoidales del " de la entrada; signal" x_n. El DFT computa el X_k del x_n, mientras que el IDFT demuestra cómo computar el x_n como suma del X_k sinusoidal \ exp de los componentes (2 \ pi i k n/N)/N con la frecuencia k/N completa un ciclo por muestra. Escribiendo las ecuaciones en esta forma, estamos haciendo el uso extenso de la fórmula de Euler de expresar sinusoids en términos de exponentials complejos, que son mucho más fáciles de manipular. (De la misma manera, escribiendo X_k en la forma polar, obtenemos inmediatamente la amplitud sinusoide del |X_k| y la fase de la discusión compleja .) Una delicadeza importante de esta representación, alias, se discute abajo.

Observar que el factor de la normalización que multiplica el DFT e IDFT (aquí 1 y 1 N ) y las muestras de los exponentes son simplemente convenciones, y diferenciar en algunos tratamientos. Los únicos requisitos de estas convenciones son que los DFT y los IDFT tienen exponentes de la opuesto-muestra y que el producto de sus factores de la normalización sea 1 N . Una normalización de 1/\ de la raíz cuadrada {N} para el DFT e IDFT hace transforma el unitario, que tiene algunas ventajas teóricas, pero es a menudo más práctico en el cómputo numérico realizar el escalamiento de una vez como arriba (y un escalamiento de la unidad puede ser conveniente de otras maneras).

(La convención de un negativo firma adentro el exponente es a menudo conveniente porque significa que X_k es la amplitud de un " frequency" positivo; 2 \ pi k/N. Equivalente, el DFT se piensa a menudo en como filtro emparejado : al buscar una frecuencia de +1, una correlaciona la señal entrante con una frecuencia de −1.)

En la discusión siguiente el " de los términos; sequence" y " vector" será considerado permutable.

Características

Lo completo

El Fourier discreto transforma es un inversible, transformación linear del

l \ mathcal {F}: \ ^N del mathbb {C} \ \ mathbb {C} ^N

con el \ el mathbb {C} denotando el sistema de los números complejos es decir para cualquie   del N ; >  0, un N - vector complejo dimensional tiene un DFT y un IDFT que son alternadamente el N - vectores complejos dimensionales.

Ortogonalidad

La forma del e^ de los vectores {\ kn del frac {2 \ pi i} {N}} una base ortogonal sobre el sistema de N - vectores complejos dimensionales: del

l \ ^ del sum_ {n=0} {N-1} \ ido (e^ {\ kn del frac {2 \ pi i} {N}} \ derecho) \

dejado (e^ {- \ k'n del frac {2 \ pi i} {N}} \ derecho) N~ \ delta_ {kk'}

donde ~ \ delta_ {kk'} es el delta de Kronecker. Esta condición de la ortogonalidad se puede utilizar para derivar la fórmula para el IDFT de la definición del DFT.

El teorema de Plancherel y teorema de Parseval

Si el k
del del k y del Y del del X es el DFTs de los estados del teorema de Plancherel n del del n y del y del del x respectivamente entonces: = \ frac {1} {N} \ ^ del sum_ {k=0} {N-1} X_k Y^*_k del y^*_n del x_n del ^ del \ del sum_ del

l {n=0} {N-1}

donde la estrella denota la conjugación compleja. El teorema de Parseval es una caja especial del teorema y de los estados de Plancherel: del

l \ ^ del sum_ {n=0} {N-1} |x_n|^2 = \ frac {1} {N} \ ^ del sum_ {k=0} {N-1} |X_k|^2.

Periodicidad

Si la expresión que define el DFT se evalúa para todos los números enteros k en vez de apenas para el k = 0, \ puntea, N-1 , después la secuencia infinita resultante es una extensión periódica del DFT, periódica con el N del período.

La periodicidad se puede demostrar directo de la definición: X_ {k+N} = \ ^ del x_n. del ^ del sum_ {n=0} {N-1} {- \ frac {2 \ pi i} {N} (k+N) n} = \ e^ del ^ del x_n. del ^ del sum_ {n=0} {N-1} {- \ frac {2 \ pi i} {N} k n} {- 2 \ pi i n} = \ ^ del x_n. del ^ del sum_ {n=0} {N-1} {- \ frac {2 \ pi i} {N} k n} = X_k

donde hemos utilizado el hecho ese e^ del {- 2 \ pi i} = 1. Puede ser demostrado de la misma manera que la fórmula de IDFT lleva a una extensión periódica.

El teorema del cambio

Multiplicando x_n por un linear \ exp de la fase del (\ frac {2 \ pi i} {N} n m) para un cierto número entero m corresponde a un cambio circular del de la salida X_k: X_k es substituido por el X_ {kilómetro} , donde está el Modulo el subíndice interpretado N (es decir, periódico). Semejantemente, un cambio circular de la entrada x_n corresponde a multiplicar la salida X_k por una fase linear. Matemáticamente, si el \ {x_n \} entonces representa el x del vector

l si \ {F} (\ {x_n \}) _k=X_k mathcal

entonces \ mathcal {F} (\ {^ del x_n. {\} \} del frac {2 \ pi i} {N} n m) _k=X_ {kilómetro}

l y \ {F} (\ {x_ {} \} del nanómetro) e^ mathcal del _k=X_k {- \ frac {2 \ pi i} {N} k m}

Teorema circular de la circunvolución y teorema de la correlación cruzada

El circular x * y de la circunvolución del cíclico o del x de dos vectores =   del xk ; y y =   del yn ; es el x * y del vector con los componentes

(\ mathbf {x*y}) _n = \ el sum_ {m=0} el ^ {N-1} x_m y_ {} \ patio \ patio del nanómetro n = 0, \ puntea, N-1

donde continuamos el y cíclico de modo que

y_ {- m} = y_ {} \ patio \ quad~~~~~~~~~~ del nanómetro m = 0, \ puntea, N-1

El Fourier discreto transforma circunvoluciones cíclicas de las vueltas en la multiplicación componente-sabia. Es decir, si z_n del

l = (\ mathbf {x*y}) _n

entonces el Z_k=X_k Y_k \ patio \ quad~~~~~~~~~~ del

l k = 0, \ puntea, N-1

donde las mayúsculas ( X, Y, Z ) representan el DFTs de las secuencias representadas por las pequeñas letras ( x, y, z ). Observar que si adoptan a una diversa convención de la normalización para el DFT (e., la normalización unitaria), después en general habrá un factor constante que multiplica la relación antedicha.

La evaluación directa de la adición de la circunvolución, arriba, requeriría el O (las operaciones de N^2), solamente el DFT (vía un FFT) proporcionan un O (N \ método del registro N) para computar la misma cosa. Inversamente, las circunvoluciones se pueden utilizar para computar eficientemente DFTs vía el algoritmo de FFT de Rader y el algoritmo de FFT de Bluestein. El método se puede ampliar a las señales non-circular usar Traslapar-agrega el método .

El considera también: teorema de la circunvolución de

En una manera análoga, puede ser demostrado que si z_n es la correlación cruzada de x_n y de y_n: _n del z_n= del

l (\ mathbf {x * y}) = \ x_m^* del ^ del sum_ {m=0} {N-1} \, y_ {m+n}

donde está otra vez cíclico la suma en el m, después el Fourier discreto transforma de z_n es: Z_k = X_k^* del

l \, Y_k

donde las mayúsculas se utilizan otra vez para significar el Fourier discreto transformar.

Polinomio trigonométrico de la interpolación

El p polinómico del de la interpolación trigonométrica (t) = \ frac {X_0} {N} + \ e^ del frac {X_1} {N} {él} + \ + \ frac {X_ {N/2}} {N} de los cdots \ e^ de lechuga romana (Nt/2) + \ frac {X_ {N/2+1}} {N} {(- N/2+1) él} + \ cdots + \ e^ del frac {X_ {N-1}} {N} {- él} para el incluso, p de N del
(t) = \ + \ frac {X_1} {N del frac {X_0} {N}} e^ {él} + \ cdots + \ e^ del frac {X_ {\ lfloor N/2 \ rfloor}} {N} {\ lfloor N/2 \ rfloor él} + \ e^ del frac {X_ {\ lfloor N/2 \ rfloor+1}} {N} {- \ lfloor N/2 \ rfloor él} + \ cdots + \ e^ del frac {X_ {N-1}} {N} {- él} para N impar, donde el k
/ N del del X de los coeficientes es dado por el DFT del n del del x arriba, satisface el p de la característica de la interpolación (2 \ pi n/N) = x_n para n=0, \ los ldots, N-1.

Para incluso N, notan que Nyquist componente \ frac {X_ {N/2}} {} \ lechuga romana de N (Nt/2) se dirige especialmente.

Esta interpolación es el no único: el alias implica que uno podría agregar el N a las frecuencias complejo-sinusoides unas de los (e. e^ cambiante {- él} al e^ {i (N-1) t} ) sin el cambio de la característica de la interpolación, pero dando a diverso valora entre los puntos de x_n. La opción arriba, sin embargo, es típica porque tiene dos características útiles. Primero, consiste en los sinusoids cuyas frecuencias tienen las magnitudes posibles más pequeñas, y por lo tanto reduce al mínimo el \ internacional de la cuesta de la media cuadrada |p'(t)|^2 dt de la función de interpolación. En segundo lugar, si los x_n son números verdaderos, después p (t) es verdadero también.

En cambio, el polinomio trigonométrico más obvio de la interpolación es el en el cual la gama de frecuencias a partir de la 0 a N-1 (en vez áspero de -N/2 a +N/2 como arriba), similar a la fórmula inversa de DFT. Esta interpolación hace el no reduce al mínimo la cuesta, y es el no generalmente con valores reales para x_n verdadero; su uso es un error común.

El DFT unitario

Otra manera de mirar el DFT es observar que en la discusión antedicha, el DFT se puede expresar como matriz de Vandermonde: del

l \ mathbf {F} = \ comenzar {el bmatrix} \ y \ omega_N^ {0 \ cdot 1} del omega_N^ {0 \ cdot 0} y \ y \ omega_N^ {0 \ cdot (N-1) de los ldots} \ \ \ y \ omega_N^ {1 \ cdot 1} del omega_N^ {1 \ cdot 0} y \ y \ omega_N^ {1 \ cdot (N-1) de los ldots} \ \ \ vdots y \ vdots y \ ddots y \ de los vdots \ \ \ omega_N^ {(N-1) \ cdot 0} y \ omega_N^ {(N-1) \ cdot 1} y \ ldots y \ omega_N^ {(N-1) \ cdot (N-1)} \ \ \ extremo {bmatrix}

donde

\ omega_N = e^ {-} \, de 2 \ pi i/N

es una raíz primitiva Nth de la unidad . Lo contrario transforma entonces es dado por lo contrario de la matriz antedicha:

\ mathbf {F} ^ {- 1} = \ frac {1} {} \ mathbf {F} ^* de N

Con los constantes unitarios 1/\ raíz cuadrada {N} de la normalización, el DFT se convierte en una transformación unitaria, definida por una matriz unitaria: del

l \ mathbf {U} = \ del
/\ raíz cuadrada {N} del mathbf {F} \ ^ del mathbf {U} {- 1} = \ mathbf {U} ^*|\ det (\ mathbf {U})|=1

donde   del det del () ; es la función determinante . El determinante es el producto de los valores propios, y por lo tanto (véase abajo) está siempre el \ P. 1 o el \ P. En un espacio de vector verdadero, una transformación unitaria se puede pensar en como simplemente rotación rígida del sistema coordinado, y todas las características de una rotación rígida se pueden encontrar en el DFT unitario.

La ortogonalidad del DFT ahora se expresa como una condición de Orthonormality (que se presenta en muchas áreas de las matemáticas según lo descrito en la raíz de la unidad ): del

l \ ^*= de U_ del ^ del sum_ {m=0} {N-1} {kilómetro} U_ {manganeso} \ delta_ {kn}

Si el \ el mathbf {X} se define como el DFT unitario del \ del mathbf {x} del vector entonces X_k= del

l \ ^ del sum_ {n=0} {N-1} U_ {kn} x_n

y se expresa el teorema de Plancherel como: del

l \ y_n^* del x_n del ^ del sum_ {n=0} {N-1} = \ ^ del sum_ {k=0} {N-1} X_k Y_k^*

Si vemos el DFT como apenas transformación coordinada que especifique simplemente los componentes de un vector en un nuevo sistema coordinado, después el antedicho es apenas la declaración que el producto de punto de dos vectores está preservado bajo transformación unitaria de DFT. Para el del caso especial \ = \ mathbf {y} del mathbf {x}, esto implica que la longitud de un vector es también &mdash preservado; éste es apenas el teorema de Parseval: del

l \ ^ del sum_ {n=0} {N-1}|x_n|^2 = \ ^ del sum_ {k=0} {N-1}|X_k|^2

Expresión del DFT inverso en términos de DFT

Una característica útil del DFT es que el DFT inverso se puede expresar fácilmente en términos de DFT (delantero), vía varios " bien conocido; tricks". (Por ejemplo, en cómputos, es a menudo conveniente ejecutar solamente un Fourier rápido transforma la correspondencia a uno transforma la dirección y entonces conseguir el otro transformar la dirección del primera.)

Primero, podemos computar el DFT inverso invirtiendo las entradas: del

l \ ^ mathcal {F} {- 1} (\ {x_n \}) = \ {F} (\ {x_ {N -} \} de n) mathcal/N

(Como de costumbre, los subíndices son el Modulo interpretado N; así, para n=0, tenemos x_ {N-0} =x_0.)

En segundo lugar, uno puede también conjugar las entradas y las salidas: del

l \ ^ mathcal {F} {- 1} (\ mathbf {x}) = \ {F} (\ ^* del mathbf {x}) ^* mathcal/N

Tercero, una variante de este truco de la conjugación, que es a veces preferible porque no requiere ninguna modificación de los valores de datos, implica el intercambiar de las piezas verdaderas e imaginarias (que se pueden hacer en una computadora simplemente modificando los indicadores . Definir el intercambio (x_n) como x_n con su swapped&mdash de las piezas verdaderas e imaginarias; es decir, si el x_n = a + b i entonces intercambian (x_n) está el b + un i. Equivalente, el intercambio (x_n) iguala el i x_n^*. Entonces del

l \ = mathcal \ textrm {intercambio} del ^ {F} {- 1} (\ mathbf {x}) (\ mathcal {F} (\ textrm {intercambio} (\ mathbf {x}))) /N

Es decir, lo contrario transforma es igual que el delanteros transforman con las piezas verdaderas e imaginarias intercambiadas por entrada y salida, hasta una normalización ( y otros, 1988 de Duhamel).

El truco de la conjugación se puede también utilizar para definir un nuevo transforma, estrechamente vinculado al DFT, que es &mdash de Involutary ; es decir, que es su propio lo contrario. Particularmente, el T (\ mathbf {x}) = \ {F} (\ ^* del mathbf {x})/mathcal \ raíz cuadrada {N} es claramente su propio lo contrario: T (T (\ mathbf {x})) = \ mathbf {x} . Una transformación involutary estrechamente vinculada (por un factor de (1+ el i )/√2) está el H (\ mathbf {x}) = \ {F} ((1+i) \ ^* del mathbf {x})/mathcal \ raíz cuadrada {2N} , desde los factores del (1+i) en el H (H (\ mathbf {x}))cancelación de los 2. Para el \ el mathbf verdaderos {x} de las entradas, la parte real del H (\ mathbf {x}) no es ninguna con excepción que Hartley discreto transforma, que es también involutary.

Valores propios y vectores propios

Los valores propios de la matriz de DFT son simples y bien sabido, mientras que los vectores propios son complicados, no únicos, y son el tema de la investigación en curso.

Considerar el \ el mathbf unitarios {U} de la forma definidos arriba para el DFT de la longitud N, donde el _ del \ del mathbf {U} {m, n} = \ = \ exp/\ raíz cuadrada {N} del omega_N^ {manganeso} (- 2 \ pi i mn/N)/\ raíz cuadrada {N} . Esta matriz satisface la ecuación: = \ mathbf {I} del \ del mathbf del {U} ^4. Esto se puede ver de las características inversas arriba: el del funcionamiento \ el mathbf {U} da dos veces los datos originales en orden reversa, así que el del funcionamiento \ el mathbf {U} cuatro veces de da detrás los datos originales y es así la matriz de identidad . Esto significa que el \ lambda de los valores propios satisface una ecuación característica : \ lambda^4 = 1. del Por lo tanto, los valores propios del \ del mathbf {U} son las cuartas raíces de la unidad : el \ lambda es +1, −1, + el i, o el i del −.

Puesto que hay solamente cuatro valores propios distintos para este N \ matriz de las épocas N, tienen cierta multiplicidad . La multiplicidad da el número de vectores propios independientes linear que corresponden a cada valor propio. (Nota que hay vectores propios de la independiente del N ; la matriz no es el defectuoso.)

El problema de su multiplicidad fue solucionado por McClellan y Parks (1972), aunque fuera más adelante demostrado para haber sido equivalente a un problema solucionado por el gauss (Dickinson y Steiglitz, 1982). La multiplicidad depende del valor del Modulo 4 de N, y es dada por la tabla siguiente:

La verdadero-entrada DFT

Si x_0, \ los ldots, el x_ {N-1} son los números verdaderos pues están a menudo en usos prácticos, después el DFT obedece la simetría: X_k del

l = ^* de X_ {N-k},

donde la estrella denota la conjugación compleja y los subíndices son el interpretado N del modulo.

Por lo tanto, el DFT hecho salir para las entradas verdaderas es a medias redundante, y una obtiene la información completa solamente mirando áspero la mitad de las salidas X_0, \ ldots, X_ {N-1} . En este caso, el " DC" el elemento X_0 es puramente verdadero, y para incluso el N el " Nyquist" el X_ del elemento {N/2} es también verdadero, tan allí es exactamente números verdaderos non-redundant del N en la primera mitad + el elemento de Nyquist del complejo X de la salida.

Usar la fórmula de Euler, el polinomio trigonométrico de interpolación se puede entonces interpretar como suma de funciones del seno y de coseno.

Generalizado/cambió de puesto DFT

Es posible cambiar de puesto el muestreo de la transformación a tiempo y/o el dominio de frecuencia por un cierto verdadero de los cambios un y el b, respectivamente. Esto se conoce a veces como DFT generalizado el (o GDFT ), también llamado el DFT cambiado de puesto o DFT compensado, y tiene características análogas al DFT ordinario:

X_k = \ sum_ {n=0} ^ {N-1} x_n. ^ {- \ el frac {2 \ pi i} {N} (k+b) (n+a)} \ el patio \ el patio k = 0, \ puntea, N-1

Lo más a menudo posible, los cambios de 1/2 (mitad de una muestra) se utilizan. Mientras que el DFT ordinario corresponde a una señal periódica en dominios del tiempo y de frecuencia, a=1/2 produce una señal que sea anti-periódica en el dominio de frecuencia (X_ {k+N} = - X_k) y viceversa para b=1/2. Así, el caso específico del a = de b = de 1/2 se conoce como impar-frecuencia del impar-tiempo del que Fourier discreto transforma (u O2 DFT). Tal cambiado de puesto transforma es el más de uso frecuente para los datos simétricos, representar diversas simetrías del límite, y para los datos verdadero-simétricos corresponden a diversas formas del coseno discreto y el seno transforma.

Otra opción interesante es el a=b=- (N-1) /2, que se llama el DFT centrado el (o CDFT ). El DFT centrado tiene la característica útil que, cuando N es un múltiplo de cuatro, los cuatro de sus valores propios (véase arriba) tener multiplicities iguales (Rubio y Santhanam, 2005).

El Fourier discreto transforma se puede ver como caso especial del Z-transform, evaluado en el círculo de unidad en el plano complejo; más general z-transforma corresponden al complejo de los cambios del un y al b arriba.

¡DFT< multidimensional! -- Esta sección se liga de que Fourier rápido transforma -->

El DFT ordinario computa la transformación de un " one-dimensional" grupo de datos: una secuencia (o arsenal ) x_n que es una función de una variable discreta n. Más generalmente, uno puede definir el multidimensional DFT de un x_ de la matriz multidimensional {n_1, n_2, \ puntos, n_d} que sea una función del n_ \ de la ana discretos de las variables de d = 0, 1, \ puntea, N_ \ ell-1 para el \ ell en 1, 2, \ puntea, d:

X_ {k_1, k_2, \ puntea, k_d} = \ sum_ {n_1=0} ^ {N_1-1} \ a la izquierda (\ omega_ {N_1} ^ {~k_1 n_1} \ sum_ {n_2=0} ^ {N_2-1} \ a la izquierda (\ ^ del omega_ {N_2} {~k_2 n_2} \ cdots \ ^ del sum_ {n_d=0} {N_d-1} \ ^ del omega_ {N_d} {n_d del ~k_d} \ x_ del cdot {n_1, n_2, \ puntea, n_d} \ derecho) \ cdots \) derecho \,

donde = \ exp (- 2 \ pi i/N_ \ ana) como arriba y los índices del \ del omega_ {N_ \ ana} de la salida de d funciona de k_ \ de ana = 0, 1, \ puntea, N_ \ ell-1. Esto más compacto se expresa en la notación del vector, donde definimos el \ el mathbf {n} = (n_1, n_2, \ puntea, n_d) y el \ el mathbf {k} = (k_1, k_2, \ puntea, k_d) como vectores de d-dimensional de índices a partir de la 0 al \ al mathbf {N} - 1, que definimos como el \ mathbf {N} - 1 = (N_1 - 1, N_2 - 1, \ puntea, N_d - 1):

X_ \ mathbf {k} = \ sum_ {\ mathbf {n} =0} ^ {\ mathbf {N} - 1} e^ {- 2 \ pi i \ mathbf {} \ cdot (\/\ mathbf {N} de k del mathbf {n})} x_ \ mathbf {} \, de n,

donde el de la división \/\ mathbf {N} del mathbf {n} se define como/\ mathbf {N} del \ del mathbf {n} = (n_1/N_1, \ puntea, n_d/N_d) para ser elemento-sabio realizado, y la suma denota el sistema de adiciones jerarquizadas arriba.

Lo contrario del DFT multidimensional es, análogo al caso unidimensional, dado cerca:

x_ \ mathbf {n} = \ frac {1} {\ ^d N_ \ ana del prod_ {\ ell=1}} \ sum_ {\ mathbf {k} =0} ^ {\ mathbf {N} - 1} e^ {2 \ pi i \ mathbf {} \ cdot (\/\ mathbf {N} de n del mathbf {k})} X_ \ mathbf {} \, de k.

El DFT multidimensional tiene una interpretación simple. Apenas mientras que el DFT unidimensional expresa la entrada x_n como superposición de sinusoids, el DFT multidimensional expresa la entrada como superposición de las ondas planas o de los sinusoids que oscilan a lo largo el de la dirección \ de/\ mathbf {N} del mathbf {k} en espacio y que tienen el X_ \ mathbf {k} de la amplitud. Tal descomposición es de gran importancia para todo del tratamiento de la imagen de Digitaces ( d = 2) a solucionar las ecuaciones diferenciales parciales en tres dimensiones ( d = 3) rompiendo la solución para arriba en ondas planas.

De cómputo, el DFT multidimensional es simplemente la composición de una secuencia de DFTs unidimensional a lo largo de cada dimensión. Por ejemplo, en el x_ de dos dimensiones del caso {n_1, n_2} uno puede primero computar el N_1 DFTs independiente de las filas (es decir, a lo largo de n_2) para formar un nuevo y_ del arsenal {n_1, k_2} , y después computa el N_2 DFTs independiente de y a lo largo de las columnas (a lo largo de n_1) para formar el X_ del resultado final {k_1, k_2} . O, uno puede transformar las columnas y entonces el rows— la orden es inmaterial porque las adiciones jerarquizadas sobre conmutan .

Debido a esto, dado una manera de computar un DFT unidimensional (e. un algoritmo unidimensional ordinario de FFT), uno tiene inmediatamente una manera de computar eficientemente el DFT multidimensional. Esto se conoce como algoritmo de la fila-columna del, aunque haya también intrínseco los algoritmos multidimensionales FFT.

Usos

El DFT ha visto uso amplio a través de una gran cantidad de campos; bosquejamos solamente algunos ejemplos abajo (véase también las referencias en el extremo). Todos los usos del DFT dependen crucial de la disponibilidad de un algoritmo rápido para computar Fourier discreto transforman y sus lo contrario, un Fourier rápido transforman .

Análisis espectral

Cuando DFT es utilizado para espectral análisis, \ {x_n \} \, la secuencia de representa generalmente un sistema finito de tiempo-muestras uniformly-spaced de un cierto x de la señal (t) \, , donde el t representa tiempo. La conversión a partir del tiempo continuo a las muestras (tiempo discreto) cambia el subyacente Fourier transforma de x (t) en un tiempo discreto Fourier transforma (DTFT), que exige generalmente un tipo del alias llamado distorsión . La opción de una muestra-tarifa apropiada (véase la frecuencia de Nyquist) es la llave a reducir al mínimo esa distorsión. Semejantemente, la conversión de una secuencia muy larga (o infinita) a un tamaño manejable exige un tipo de llamado distorsión '' salida '', que se manifiesta como pérdida de detalle (resolución del aka) en el DTFT. La opción de una longitud apropiada de la secundario-secuencia es la llave primaria a reducir al mínimo ese efecto. Cuando los datos disponibles (y la hora de procesarlos) es más que la cantidad necesaria para lograr la resolución deseada de la frecuencia, una técnica estándar es realizar DFTs múltiple, por ejemplo para crear un espectrograma . Si el resultado deseado es un espectro de energía y el ruido o la aleatoriedad está presente en los datos, hacer un promedio de los componentes de la magnitud del DFTs múltiple es un procedimiento útil para reducir la variación del espectro (también llamado un periodograma en este contexto); dos ejemplos de tales técnicas son el método galés y el método de Bartlett.

Una fuente final de distorsión (o quizás de la ilusión del ) es el DFT sí mismo, porque es apenas un muestreo discreto del DTFT, que es una función de un dominio de frecuencia continuo. Eso puede ser atenuada aumentando la resolución del DFT. Ese procedimiento se ilustra en el tiempo discreto Fourier transforma el artículo de .
El procedimiento se refiere a veces como cero-acolchado del, que es una puesta en práctica particular usada conjuntamente con el que Fourier rápido transforma algoritmo de (FFT). La ineficacia de realizar multiplicaciones y adiciones con el " cero-valorado; samples" es más que compensó por la eficacia inherente del FFT.
Según lo observado ya, la salida impone un límite ante la resolución inherente del DTFT. Tan hay un límite práctico a la ventaja que se puede obtener de un DFT de grano fino.

Compresión de datos

El campo del proceso de señal numérica confía pesadamente en operaciones en el dominio de frecuencia (es decir en el Fourier transformar). Por ejemplo, varios la imagen del lossy y los métodos de la compresión del sonido emplean el Fourier discreto transforman: la señal se corta en segmentos cortos, se transforma cada uno, y entonces los coeficientes de Fourier de los de alta frecuencia, que se asumen para ser unnoticeable, se desechan. El descompresor computa lo contrario transforma basado en este número reducido de coeficientes de Fourier. (Los usos de la compresión utilizan a menudo una forma especializada del DFT, el coseno discreto transforma o el coseno discreto modificado transforma a veces ).

Ecuaciones diferenciales parciales

Fourier discreto transforma es de uso frecuente solucionar las ecuaciones diferenciales parciales, donde el DFT se utiliza otra vez como aproximación para la serie de Fourier Del (que se recupera en el límite del infinito N ). La ventaja de este acercamiento es que amplía la señal en el inx complejo del del del e de los exponentials, que son funciones propias de la diferenciación: d /inx del e del dx del = en el inx del del del e de . Así, en la representación de Fourier, la diferenciación es simple-nosotros apenas se multiplica por el i n . Una ecuación diferencial linear con coeficientes constantes se transforma en una ecuación algebraica fácilmente soluble. Uno entonces utiliza el DFT inverso para transformar el resultado nuevamente dentro de la representación espacial ordinaria. Tal acercamiento se llama un método espectral .

Multiplicación polinómica

Suponer que deseamos computar el polinómico c ( x ) del producto = un ( x ) · b ( x ). La expresión ordinaria del producto para los coeficientes del c implica una circunvolución (acíclica) linear, donde no lo hacen los índices " abrigo around." Esto se puede reescribir como circunvolución cíclica tomando a los vectores del coeficiente para el un ( x ) y el b ( x ) con término constante primero, después añadiendo ceros de modo que el de los vectores del coeficiente del resultado un y el b tengan d de la dimensión > grado ( un ( x )) + grado ( b ( x )). Entonces, del

l \ mathbf {c} = \ * \ mathbf {b} del mathbf {a}

Donde está el vector de los coeficientes para el c ( x ), y el * el c del operador de la circunvolución \, se define tan

c_n = \ sum_ {m=0} ^ {d-1} a_m b_ {nanómetro \ \ mathrm {} \ qquad \ qquad \ qquad n=0,1} \ d de la MOD,…, d-1

Pero la circunvolución se convierte en multiplicación debajo del DFT: del

l \ mathcal {F} (\ mathbf {c}) = \ {F} (\ mathbf {b}) mathcal {F} (\ mathbf {a}) \ mathcal

Aquí el producto de vector se toma el elementwise. Así los coeficientes del polinómico c ( x ) del producto son apenas los términos 0,…, grado ( un ( x )) + grado ( b ( x )) del vector del coeficiente del

l \ mathbf {c} = \ ^ {F} {- 1} (\ mathcal {F} (\ mathbf {a}) \ mathcal mathcal {F} (\ mathbf {b})).

Con un que Fourier rápido transforma, el algoritmo resultante toma operaciones aritméticas de O ( N del registro del N ). Debido a su simplicidad y velocidad, el algoritmo de Cooley-Tukey FFT, que se limita a los tamaños compuestos, se elige a menudo para la operación de la transformación. En este caso, el d se debe elegir como el número entero más pequeño mayor que la suma de los grados polinómicos de la entrada que es factorizable en los pequeños factores primeros (e. 2, 3, y 5, dependiendo de la puesta en práctica de FFT).

Multiplicación de números enteros grandes

Los algoritmos sabidos más rápidos para la multiplicación del uso muy grande de los números enteros que el método polinómico de la multiplicación contorneó arriba. Los números enteros se pueden tratar como el valor de un polinomio evaluado específicamente en la base de número, con los coeficientes de la correspondencia polinómica a los dígitos en esa base. Después de la multiplicación polinómica, un paso de la llevar-propagación de la bajo-complejidad termina relativamente la multiplicación.

Algún Fourier discreto transforma pares

Derivación como serie de Fourier

considera también: La derivación A del Fourier discreto transforma el

El DFT se puede derivar como truncamiento de la serie de Fourier Del de una secuencia periódica de los impulsos

Ver también


matriz DFT del


la lista Fourier-relacionado transforma

.

  • Zenithic
  • Friedrich Wieck
    Random links:Guillermo Harcourt | Sannan, Hyōgo | Ron del espolón de gallo | RRS James Clark Ross | Primera república del este de Turkestan

  • © 2007-2008 enciclopediaespana.com; article text available under the terms of GFDL, from en.wikipedia.org
    ="http://pagead2.googlesyndication.com/pagead/show_ads.js">