Расчет и построение фазовых спектров., требуется помощь Опции

[spirit]

Впервые на таком форуме. Надеюсь, правильно выбран раздел.
Я начинаю заниматься (уже пару лет) различными алгоритмами спектрального анализа на основе преобразования Фурье. Сейчас возникла проблема построения фазового спектра. При расчете, проще всего находить мнимую и вещественную часть спектра. С амплитудами все понятно: корень из суммы квадратов мнимой и вещественной части.
А вот при нахождении фазы приходится применять функцию арктангенса. Отсюда как минимум два момента:
1) Из-за вычислительной погрешности для некоторых типовых сигналов, вместо нулевых значений мнимой и вещественной части, могут быть получены 0,000000001 или
-0,0000001. При вычислении фазы, определяется квадрант, в котором находится угол и затем арктангенс отношения мнимой части к вещественной. В таком случае знак может сыграть роль ошибки на 360 градусов.
2) Фазовый сдвиг может превышать 360 градусов. Причем, фаза, доходя до 360 может начать расти или уменьшаться.

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

Может быть, кто-то сталкивался с этим. Литература, советы, алгоритмы, опыт… Все что угодно.