БПФ, "грязная" Im составляющая Опции

[AndreyVN]

Всем привет!

Отлаживаю алгоритм комплексного БПФ, формально работает корректно, функцию туда-сюда гоняет без искажений, для синуса пики приблизительно на частоте N/(2*Pi), N-N/(2*Pi).
Не нравится как выглядит аргумент результата преобразования - с какими-то шумами (первый рисунок). Для сравнения второй рисунок - результат преобразования чистого синуса.

красная линия (пики слева и справа) - модуль комплексных к-ов Фурье образа;
зеленая линия - аргумент комплексных к-ов Фурье образа;
синяя линия - действительная часть обратного преобразования Фурье (мнимая часть =0).

[Xenia]

Может быть, всё дело в шкале? Многие программы, выводя графики, стараются растянуть их до максимально возможного, т.е. увеличивают масштаб до тех пор, пока график хотя бы одной точкой не коснется рамки. В таких случаях пустые сканы оказываются растянутыми до безобразия, превращаясь в частокол, хотя на самом деле имеют микроскопически малый шум, если их изобразить в одном и том же масштабе со сканом, где присутствует полезный сигнал. Именно на этот вывод намекает то обстоятельство, что зеленая линия нигде не зашкаливает, изящно упираясь в рамку одной точкой. При таком сильном "шуме" случайное совпадение маловероятно.

Советую математически найти максимумы в числах для каждой кривой, а затем сравнить их между собой.