Отправлено fontp 15 июля 2004 г. 12:26 В ответ на: Надо считать угол между комлексными числами (разницу фаз). А есть ли что-нибудь проще чем atan(Im1/Re1)-atan(Im2/Re2)? отправлено AntZ 15 июля 2004 г. 11:17

сложность вычислений определяется сложностью ф-и atan2. Atan2 после разборки по квадрантам вычисляется через atan(x/y). Вообще полезно найти и посмотреть эти функции в библиотеках С и ассемблерных. После этого остаётся выбрать себе аппроксимацию atan с достаточной для задачи точностью. С относительной точностью в 0.5% есть например такая апроксимация
в интервале от -pi/2 до +pi/2

atan(x) = (x<0) ? x*(1.0597+0.27131*x) : x*(1.0597-0.27131*x)

Для DSP скорее лучше вычислять без сравнений, через abs(x)

• вернее pi/4 — fontp   (15.07.2004 12:44, 1728 байт)
  • Спасибо, арктангенс есть в dsplib 55xx, можно взять за основу. Видимо без честного деления и арктангенса не обойтись. Делаю демодуляцию DPSK с мягкими решниями, 00-strong 0 , 01- soft 0, 10 - soft 1, 11- strong 1. Мне нужно проверить разницу фаз c несколькими дискретными уровненями. — AntZ   (15.07.2004 13:25, пустое)
    • А зачем для демодулятора(DPSK) с мягким решением углы и вообще арктангенсы? — alostap   (15.07.2004 15:56, пустое)
      • Вообще для декодера достаточно метрик (мягких решений) = — alostap   (15.07.2004 16:11, 158 байт)
    • тогда лучше орктангенс не вычислять, а сделать типа векторного квантования — fontp   (15.07.2004 13:37, 225 байт)
      • Очень понравилось выражение — AntZ   (15.07.2004 15:29, 118 байт)
        • у вас есть плоскость (x, y) — fontp   (15.07.2004 15:43, 808 байт)
          • конфа не ест математической нотации x меньше a[i] умножить на y ;) — fontp   (15.07.2004 15:47, пустое)
            • Никаких проблем - меньше амперсандlt - <, больше амперсандgt - > — AntZ   (15.07.2004 16:47, пустое)
            • Съест если знак меньше написать по HTML-правилам — SM   (15.07.2004 15:58, пустое)
      • Согласен, еще использовал для расчета частоты опорника при устранении частотной ошибки. По-моему на каждый бит такое считать - слишком дорого. — andy_P   (15.07.2004 14:56, пустое)