Вычисление фазы комплексного сигнала, Алгоритм для ПЛИС Опции

[Костян]

Ищу алгоритм вычисления фазы для комплексного сигнала.
В литературе Understanding Digital Signal Processing (2nd Edition) Richard G. Lyons дан алгоритм апроксимации atan, как
Q=I*Q/(I^2+0.28*Q^2).

Все бы ничего, но как разделить на (I^2+0.28*Q^2) , тогда как выражение очень напоминает схему АРУ (I^2+Q^2) ?

Если ли другие алгоритмы удобные для ПЛИС и с малым кол-вом ресурсов? Cordic не хочется из-за объемности. Точность нужна низкая, хватит вполне и +/- 10 градусов.