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.
Sk=s0+s2ⅇ-ⅈ2π2kN+…+sN−2ⅇ-ⅈ2πN−2kN+s1ⅇ-ⅈ2πkN+s3ⅇ-ⅈ2π2+1kN+…+sN−1ⅇ-ⅈ2πN−2+1kN=s0+s2ⅇ-ⅈ2πkN2+…+sN−2ⅇ-ⅈ2πN2−1kN2+s1+s3ⅇ-ⅈ2πkN2+…+sN−1ⅇ-ⅈ2πN2−1kN2ⅇ-ⅈ2πkN
S
k
s
0
s
2
2
2
k
N
…
s
N
2
2
N
2
k
N
s
1
2
k
N
s
3
2
2
1
k
N
…
s
N
1
2
N
2
1
k
N
s
0
s
2
2
k
N
2
…
s
N
2
2
N
2
1
k
N
2
s
1
s
3
2
k
N
2
…
s
N
1
2
N
2
1
k
N
2
2
k
N
(1)
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
.
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).
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.
Notese que el orden de la secuencia de entrada en dos partes de
figura 1 no estan absolutamente iguales. ¿por que no?, ¿cómo esta determinado este orden?
El panel de arriba no usa el algoritmo de FFT para procesar la DFT de longitud-4 mientras que el abajo si.El orden es determinado por el algoritmo.
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.
"Señales y Sistemas is a Spanish translation of Dr. Rich Baraniuk's collection Signals and Systems (col10064). The translation was coordinated by an an assistant electrical engineering professor […]"