Ищу алгоритм вычисления фазы для комплексного сигнала.
В литературе 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 градусов.

..табличный метод?

Составьте таблицу - если I такой-то, а Q такой-то то угол = такой-то , загрубите таблицу до нужной вам точности.

1. Квадрант определяется знаками квадратур. Знаки можно убрать и дальше оперировать |I| и |Q|. Если у Вас есть I^2 и Q^2, то сгодятся и они
2. Половина квадранта - |I| > |Q|, I' = max(|I|, |Q|) Q' = min(|I|, |Q|)
3. 5 делений внутри половинки квадранта можно определить 1) по минимуму евклидового расстояния (удобно, если сигнал нормирован) 2) за счет соотношений Q' < Ai * I' где Ai -сумма отрицательных степеней 2 (сколько членов в Ai и каких именно степеней - это Вам предстоит выяснить), 3) Все же сделать CORDIC на 3 уровня. Можно сделать последовательный.

..таблица определяющая угол с точностью 6 градусов будет всего из 128 значений.