Комплексное БПФ, Над вещественным сигналом Опции

[Student21]

Есть процедура БПФ, написаная на Си.
Функция принимает два массива, реальной и мнимой части. Выдаёт два массива, так же реальную и мнимую часть.

Так вот, требуется работать только с вещественными сигналами.
Ставлю опыты над функцией БПФ.

Подаю на вход только вещественный сигнал, мнимый массив заполняю нулями. Как и положено на выходе реальная и мнимая часть.
Сравнил с матлабовскими данными - всё сходится. Можно запросто вычислить хоть амплитуду хоть угол гармоник.

Вопрос: Я слышал, что используя комплексное БПФ можно за один проход алгоритма преобразовывать два вещественных сигнала.
Как это происходит? Что при этом подаётся на вход, и что снимаем с выхода. Ведь реальная и мнимая часть на выходе относится только к одному сигналу. Непонятно как можно получить тоже самое для двух сигналов...


И ещё маленький вопросик. Как происходит расчёт арктангенса в микропроцессорах с фикс. точкой. Можно конечно просто влоб на плавучке посчитать, но возможно есть более экономичные варианты? Арктангенс, как можно догадаться необходим для расчёта фазы гармоник (используя реальную и мнимую часть спектра).

Огромное спасибо всем откликнувшимся.

Сообщение отредактировал Student21 - Mar 30 2011, 16:31

[diwil]

Вопрос: Я слышал, что используя комплексное БПФ можно за один проход алгоритма преобразовывать два вещественных сигнала.
Как это происходит? Что при этом подаётся на вход, и что снимаем с выхода. Ведь реальная и мнимая часть на выходе относится только к одному сигналу. Непонятно как можно получить тоже самое для двух сигналов...


И ещё маленький вопросик. Как происходит расчёт арктангенса в микропроцессорах с фикс. точкой. Можно конечно просто влоб на плавучке посчитать, но возможно есть более экономичные варианты? Арктангенс, как можно догадаться необходим для расчёта фазы гармоник (используя реальную и мнимую часть спектра).

Огромное спасибо всем откликнувшимся.

Ответ. можно. на вход подается 2 вещественных сигнала - одни в "вещественном" массиве, второй - в "мнимом". Для разделения потом сигналов используют свойство ДПФ от вещественных сигналов - четность модуля и нечетность фазы. Т.е. X[i] = conj X[N-i] Можно почитать numerical recipes in C или самому на бумажке вывести - сложного ничего нет.

Маленький ответик - много способов начиная от использования одного деление+таблица до использования кордика. Вопорс в том, какая точность требуется. Если маленькая, то таблица, если большая, то кордик. В основном используют арифметику с фиксированной точкой.