Разложение сигнала на гармоники.

[Marabou]

Сразу к сути.
Имеется дискретный сигнал, состоящий из 2001 отсчета. В состав сигнала входят <= 40 гармоник.
Задача. Получение амплитуды, фазы и синусной и/или косинусной составляющей каждой гармоники для последующего расчета и графического представления каждой гармоники.
Прошу обратить внимание, что нужно получить не амплитудно-частотную характеристику, а именно параметры каждой гармоники, входящей в состав сигнала. Просмотрел немалое количество алгоритмов реализующих преобразование Фурье(как ДПФ, так и БПФ). Не во всех смог разобраться. Те алгоритмы, в которых удалось разобрался, дают АЧХ с ошибками (неверная амплитуда, "зеркальный эффект").
Буду рад любой помощи в решении этой задачи. Спасибо!

[Marabou]

Забыл добавить. Среда разработки - Delphi. С другими, к сожалению, не знаком совершенно.

[raxp]

среда разработки не имеет значения, изучите сначала саму теорию ...Марпл, "Спектроанализ" как раз вам подойдет.

дают АЧХ с ошибками (неверная амплитуда, "зеркальный эффект")

отнюдь, все имеет под собой физическое основание и граничные условия, которые ВЫ не учитываете ...неверная амплитуда - необходима нормировка, зеркальный спектр - образуется при определенных отношениях частот сигнала и частоты дискретизации (Зоны Найквиста – зеркальные отображения спектра при использовании частот выше половины частоты дискретизации)

http://www.programmersclub.ru/Быстро...-Фурье-Практи/
http://pblog.ru/?p=658


Добавлено
"Marabou", чей то знакомое сочетание с одной лишней буквой, есть одно подозреньице

[Marabou]

http://pblog.ru/?p=658

Мне попадалась Ваша статья в самом начале моих поисков, написано хорошо, но код показался мне слишком сложным. Что ж, спасибо за разъяснения, буду разбираться.

Не знаю кого Вам напомнил мой ник. Разница в одной букве - это просто совпадение.