¿Alguien sabe cuál es la mejor manera de hacer frente a los errores de coma flotante en el proceso de cálculo de FFT? El problema es que estos errores se acumulan por todas las sumatorias. Son especialmente molestos cuando se calcula la fase para un componente de frecuencia dado:

atan(imaginario/real) puede dar un cambio de fase bastante significativo si tanto el real como el imaginario son pequeños debido a un truncamiento de coma flotante, pero del mismo orden.

La única solución que he encontrado hasta ahora es observar la magnitud máxima en la señal e ignorar cualquier frecuencia con una magnitud inferior al 1% del valor máximo y establecer su parte real en 0 para que la fase salga como NaN. Sin embargo, esto me parece un poco torpe. Cualquier idea mejor sería bienvenida.