Частота и амплитуда FFT

[WorldMaster]

Здравствуйте, возникла необходимость провести спектральный анализ сигнала.
Взял библиотеку AlgLib с реализацией дискретного преобразования Фурье.
Сгенерировал тестовый сигнал Sin(x) при x от 0 до 1024. Подал его на вход функции FFT. Получил преобразование.
Дальше нужно получить частоту гармоник и их амплитуду.
для получения частоты нашел 2 формулы но не могу определится какая правильная.
вот формула 1: Frq = i * (samples / nFFT / 2.0f);
вот формула 2: Frq = samples * fd / f.Length;

где i - номер элемента массива, полученного после FFT
samples - действительное значение i-го элемента
nFFT - число элементов FFT
fd - частота дискретизации в моем случае 1024

но график получается непонятный ... все значения стремятся к 0.
Менял функцию Sin(x) на Sin(x/2) результат один и тотже.

Вопрос:
Подскажите пожалуйста как правильно получить из массива, полученного после FFT, частоты всех гармоник и их амплитуды.
И как учитывать при этом мнимые части. Или они не нужны вовсе.
Спасибо.

[raxp]

см. в блоге.

[WorldMaster]

Сообщение от raxp

см. в блоге.

а ссылочку?

[raxp]

[WorldMaster]

антивирь сказал что там опасное содержимое ...

[WorldMaster]

Правильно ли я понял что частота расчитывается по формуле arctg ( Imk / Rek ); ???

с используемой формулой получилось следующее.
Почему максимальная частота отрицательная получилась???
частота - нижний график ось Х
Y - амплитуда в попугаях

[raxp]

антивирь сказал что там опасное содержимое ...

NOD?

arctg ( Imk / Rek );

это фаза.

амплитуду сигнала для каждой гармоники получаем складывая квадраты действительной и мнимой части и извлекая из суммы корень квадратный

[WorldMaster]

Сообщение от raxp

NOD?

это фаза.

мне частота нужна

[raxp]

cм. внимательнее. Невнимательные даже в комментариях переспрашивали.

F(i) = i * Fд / N;

[UA3GDW]

Есть у меня спектроанализатор сигнала со звуковой карты. Частота, амплитуда отображаются на графике. Писал для себя на Delphi. Могу исходником поделиться.