Отправлено -=ВН=- 12 апреля 2004 г. 13:40 В ответ на: Я сделал комплексный FFT. Делать действительный времени нет, проект горит, по своему опыту знаю, что оптимизировать FFT можно бесконечно. отправлено AntZ 12 апреля 2004 г. 13:17

Нужно после комплексного FFT добавить еще одну ступень. FFT комплексное в 2 раза меньшей длины. Входной действительный массив длиной N интерпретируется следующим образом. Четные отсчеты - реальная часть некоторого комплексного массива длиной N/2, нечетные - мнимая часть. Над этим "комплексным" массвом FFT комплексное, но размером N/2. Из результата достаются спектры действительной и мнимой частей комплексного входного массива. По другому - спектры последовательностей, составленных из четных отсчетов (действит) и из нечетных (мнимая) действительного входного массива.
Т.е. спектры последов. x(2n),x(2n+1)
При доставании руководствуются правилом, что для действит. сигнала X(k)=X'(N-k). X - спектр. отсчет, апостроф - компл. сопряжение. Затем выделенные спектры сшиваются, используется бабочка с прореж. по времени.
X(k)=XE(k)+exp(-j*2pi*k/N)*XO(k).
XE(k), XO(k) - спектры последов. x(2n),x(2n+1) соответственно , X(k) -результир. спектр действительного входного массива, т.е. спектр последов. x(n). Точнее его половина, частоты 0-N/2.
А имеющееся компл. БПФ нет никакой нужды трогать. К нему просто добавка нужна, можно как отдельную функцию оформить. И вызывать ее после FFT.
Но - дело Ваше.

• Спасибо за идеи — AntZ   (12.04.2004 13:56, пустое)
  • Пожалуйста, нет там идей, одна математика. Кучу раз опробованная. — -=ВН=-   (12.04.2004 13:59, пустое)