Derivando la Transformada Rápida de Fourier

Module by: Don Johnson Translated By Fara Meza, Erika JacksonBased on: Deriving the Fast Fourier Transform by Don Johnson

Summary: Usando el algoritmo de Cooley-Tukey para derivar la transfromada rápida.

Note: Your browser may not currently support MathML. Chrome will attempt to display MathML, but it should not be considered correct. Safari will attempt to display MathML, but it should not be considered correct. Please install the required MathPlayer plugin to view MathML correctly.Firefox requires additional mathematics fonts to display MathML correctly.Firefox requires additional mathematics fonts to display MathML correctly. See our browser support page for additional details. You can always view the correct math in the PDF version.

(Show MathML examples)

If the two examples below are mathematically the same, your browser is correctly displaying MathML, and you can dismiss this message.

Example 1:	$\frac{\sqrt{\frac{1}{2}}}{\sqrt{{f(x+y)}^2}}\frac{\sqrt{1}}{2}\langle ℝ\rangle$	Example 2:		(Hide examples)

Para derivar la FFT,asumimos que la duración de la señal es una potencia de dos: N=2l N 2 l . Considere que pasa a los elementos pares numerados y a los elemenots impares numerados de la solución de la secuencia en la DFT.

Cada término en paréntesis al cuadrado tiene la forma de una longitud N2 N 2 de DFT. La primera es una DFT de los elementos pares nuemrados,y la segunda es de los elementos impares numerados. La primer DFT es combinada con la segunda multiplicadno por el exponencial complejo ⅇ-ⅈ2πkN 2 k N . La media-longitud de la transformada es cada evaluación de los índices de frecuencia k∈0…N−1 k 0 … N 1 . Normalmente el número de indices de frecuencia del rango de una calculación de la DFT esta entre cero y la longitud de la transformación menos uno. La ventaja computacional de la FFT viene de reconocer la natulareza del periódo de la transformada discreta de Fourier. La FFT simplemente reusa las soluciones hechas en la media-longitud de la transformada y las combina a través de sumas y multiplicaciones por ⅇ-ⅈ2πkN 2 k N , que no es periódica sobre N2 N 2 , para rescribir la longitud-N DFT. La figura 1 ilustra la descomposición. Como se mantiene, ahora resolvemos dos transformadas de longitudes- N2 N 2 (complejidad 2ON24 2 O N 2 4 ), multiplicar una de ellas por el exponencial complejo (complejidad ON O N ), y sumar los resultados (complejidad ON O N ). En este punto, el total de complejidad sigue dominado por la mitad-longitud de la calculación de DFT, pero el coeficente de proporcionalidad ha sido reducido.

Ahora para la diversión. Por que N=2l N 2 l , cada una de las transformaciones de media-longitud puede ser reducida a una de dos cuartos-longitud, cada una de estas en dos octavos-longitud, etc. Esta descomposición continua hasta que nos quedamos con una transformada de longitud -2 . Esta transformada es absolutamente simple, involucrando solo sumas. Donde la primer etapa de la FFT tiene N2 N 2 transformadas de longitud-2 (vease la parte de abajo de la figura 1). Pares de estas transformadas son combinados suamando uno a otro multiplicado por el exponencial complejo. Cada par requiere 4 sumas y 4 multiplicaciones, dando un número total de computaciones igual a 8N4=N2 8 N 4 N 2 . Este número de computaciones no cambia de etapa a etapa. Ya que el número de etapas, el número de veces que la longitud puede ser dividida por dos, igual a log2N 2 N , la complejidad de la FFT es ONlogN O N N .

Figura 1: La descomposición inicial de una longitud-8 DFT en términos usando índices pares e impares entradas marca la primera fase de convertirse en el algoritmo de la FFT. Cuando esta transformada de media-longitud esta completamente descompuesta, nos quedamos con el diagrama del panel inferior que representa la solución de la FFT de longitud-8.

Descomposición de la DFT de longitud-8
(a) (b)

Haciendo un ejemplo haremos ahorros de procesos más obvios. Veamos los detalles de la DFT de longitud-8. Como se muestra en la figura 1, primero descomponemos la DFT en dos DFT de longitud-4,con las salidas sumadas y restadas en pares. Considerando la figura 1 como el índice de frecuencia que va de 0 hasta 7,reciclamos valores de la DFT de longitud-4 para los calculos finales por la periódicidad de la salida de la DFT. Examinando como los pares de salidas se recojen juntas, creamos los alementos computacionales básicos conocidos como una mariposa (figura 2).

Figura 2: Los elementos computacionales básicos de la transformada rápida de Fourier es la mariposa. Toma dos números complejos, representados por a y b, y forma las cantidades mostradas . Cada mariposa requiere una multiplicación compleja y dos sumas complejas.

Mariposa

Considerando los procesos involucrando la salida de frecuencias comunes de la DFT de dos media-longitud, vemos que las dos multiplicaciones complejas son relacionadas unas con otras, y podemos reducir nuestro proceso en un futuro. Para descomposiciones posteriores las DFT de longitud-4 en las DFT de longitud-2 y combinando sus salidas, llegamos al resumen del diagrama de la transformada rápida de Fourier de longitud-8 (figura 1). Aunque la multiplicación de los complejos es absolutamente simple (multiplicando por ⅇ-ⅈπ significa partes negativas reales y partes imaginarias), contemos ésos para los propósitos de evaluar la complejidad como el complejo se multiplica por completo. Tenemos N2=4 N 2 4 multiplicaciones complejas y 2N=16 2 N 16 sumas para cada etapa y log2N=3 2 N 3 etapas, haciendo el número de procesos básicos 3N2log2N 3 N 2 2 N como se predijo.

Otros lagoritmos "rápidos" han sido descubiertos, todos de los cuales hacen uso de cuantos factores tiene la transformada de longitud N. En teoria de números, el número de factores primos de un entero dado tiene medidas como compuesto es el. El número 16 y 81 son altamente compuestos (iguales a 24 2 4 y 34 3 4 respectivamente),el número 18 es menor asi que ( 2132 2 1 3 2 ), y el 17 no del todo (es primo). a través de treinta años de desarrollo del algoritmo de la transformada de Fourier, el algoritmo original de Cooley-Tukey esta alejado de los usados con mas frecuancia. Es tan effciente computacionalmente hablando como la transformada de longitud de potencia de dos es frecuentemente usada sin importar la longitud actual de los datos.

Comments, questions, feedback, criticisms?

Send feedback

• E-mail the author of the module, Derivando la Transformada Rápida de Fourier

Footer

This work is licensed by Erika Jackson and Fara Meza under a Creative Commons Attribution License (CC-BY 2.0), and is an Open Educational Resource.

Last edited by Fara Meza on Jul 26, 2005 7:15 pm GMT-5.