Найти фазовую функцию, для построения комплексной огибающей

[abraziv]

Всем привет, подскажите пожалуйста как найти фазовую функцию, для построения комплексной огибающей .
Имеем сигнал: a(t)= Asin(wt+f(t))
Сопряжённый ему находим через преобразование Гильберта , а именно a1(t) = 1/pi SUM(a(x)/(t-x))  - сумма от 1 до N ,где N - длительность сигнала.
Находим амплитудную огибающую как модуль сигнала и его ортогонального дополнения  , а именно A(t) = aqrt(a(t)^2 + a1(t)^2)
Фазовая функция f(t) = arctg(a1(t)/a(t)) - wt
И тогда комплексная огибающая равна s(t) = A(t)exp(j*f(t))
Подскажите , что я делаю не так.....