я используюscipy.fft на сигнал, с движущимся окном, чтобы построить амплитуды частот, изменяющихся со временем (здесьпример, время на X, частота на Y, а амплитуда - это цвет).

Однако меня интересуют только несколько частот (~ 3, только 4 частоты). С БПФ кажется, что я могуЯ выбираю только те частоты, которые мне нужны (потому что диапазон частот определяется алгоритмом), поэтому я вычисляю много бесполезных вещей, и моя программа даже вылетает сMemoryError если сигнал слишком длинный.

Что я должен делать ? Нужно ли использовать пользовательское преобразование Фурье - в этом случае, ссылки хороших реализаций приветствуются - или естьscipy путь ?

РЕДАКТИРОВАТЬ

После ответа @jfaller я решил (попытаться) реализовать алгоритм Гёртцела. Я придумал это:https://gist.github.com/4128537 но это нет работа (частота 440 не работаетне появляйся, не бери в голову пики, я нене пытайтесь применить правильное окно). Любая помощь !? Я плохо с DSP.