Отправлено ssassh 21 апреля 2006 г. 11:58 В ответ на: Ответ: смлтри CORDIC -> преобразование декартовых координат в полярные отправлено vitus_strom 21 апреля 2006 г. 10:23

А в реальности, делал ли кто по данному алгоритму arctg(x/y)?
И ещё вопросик реально ли это всё провернуть без применения математики с плавающей точкой.

• Я делал "почти" всю плавучку на С64ХХ. не доделал только некоторые библиотечные функции. — KPAH   (21.04.2006 15:32 213.85.191.245, пустое)
  • И как по затратам времени получилось? — ssassh   (21.04.2006 16:01 213.184.233.178, 58 байт)
    • Времени потратил достаточно много — KPAH   (21.04.2006 16:47 213.85.191.245, 1147 байт)
      • Просто к сведению - библиотека эмуляции плавающей точки от TI (та что в RTS входит) сравнения не по стандарту обрабатывает. У них два NAN равны между собой :-) — andy_P   (21.04.2006 16:59 80.82.63.185, пустое)
        • Я NaNы убрал. Да у меня и представление другое... А вот насчёт сравнений - помню-помню, замечал на тестах! — KPAH   (21.04.2006 17:26 213.85.191.245, пустое)
          • У них еще где-то в недрах fprintf/sprintf/snprintf и тп глюки - при выводе в формате %f - откуда-то дробная часть лезет, когда ее в оригинальном числе нет + — andy_P   (21.04.2006 17:43 80.82.63.185, 87 байт)
            • Дробной части нет на самом деле только в двух случаях: X=2^n, n>=0 или X<2^24 — KPAH   (21.04.2006 17:55 213.85.191.245, пустое)
              • Глюк касается только float. C double все в порядке. Фишка в том что x = 123457e17 (больше чем 2^24) выводится с дробной частью :-) — andy_P   (21.04.2006 18:00 80.82.63.185, пустое)
• Ответ: Да кстати на этом методе все сопры работают — vitus_strom   (21.04.2006 14:08 80.95.102.226, пустое)
• Ответ: Я делал, на ФПГА только, работает и корень тоже извлекает, и без плавучки, сейчас думаю это дело в ДСПю загнть — vitus_strom   (21.04.2006 14:06 80.95.102.226, пустое)
  • А алгоритмом поделиться не можешь? — ssassh   (21.04.2006 16:09 213.184.233.178, 191 байт)
    • Ответ: Смысл такой крутишь вектор до тех пор пока не совпадет с входным или ошибка будет в пределах, при этом для каждой итерации подсуммируется конечный результат в конце имеешь значение арктангенса для первого квадранта и корень квадратный из суммы квадратов(конечно если надо), только значение корня получается маштабированное на множитель зависящий от количества итераций, ну а из первого квадранта это перевести в реальный арктангенс больших проблем нет — vitus_strom   (21.04.2006 17:12 80.95.102.226, пустое)
      • А сами формулы не подкинешь, или где почитать. А то уже находил несколько разных интерпритаций, которые сильно отличаются... — ssassh   (21.04.2006 17:24 213.184.233.178, пустое)
        • Ответ: — ux   (23.04.2006 21:56 212.92.163.192, 59 байт, ссылка)
      • Ответ: Может в пределах месяца переложу на си, к сожалению пока нет вемени... — vitus_strom   (21.04.2006 17:13 80.95.102.226, пустое)
        • А сейчас на чём делал? — ssassh   (21.04.2006 17:20 213.184.233.178, пустое)
          • На маткаде удобно... — ux   (23.04.2006 21:58 212.92.163.192, пустое)
  • IMHO кордик хорош там где умножителей нет. В DSP часто умножение почти так же дешево как сложение. А еще человеку не atan2 a atan надо... — andy_P   (21.04.2006 14:37 80.82.63.185, пустое)
    • Ответ: Получить из atan atan2 - задачка для ученика средней школы, не знаю сколько занимает умножение в 6414 но точно никак не медленней сдвига.... — vitus_strom   (21.04.2006 17:06 80.95.102.226, пустое)
      • Из atan atan2 никак не получишь. Имея atan2 подсчитать atan легко. Вопрос в другом - если есть быстрое умножение то зачем кордик? Есть более быстро сходящиеся к искомому результату варианты. — andy_P   (21.04.2006 17:09 80.82.63.185, пустое)
        • А реално ли иих сделать не переходя на плавающую арифметику. Вот по кордику реально... — ssassh   (21.04.2006 17:25 213.184.233.178, пустое)
          • Умножение двух интов есть? . Тогда читайте мой совет про полином, если не считаете что вычислить значение полинома по схеме Горнера - это не реально. — andy_P   (21.04.2006 17:31 80.82.63.185, пустое)
            • Умножение есть, но судя по формуле x там дробный и привести его к целому виду не получится — ssassh   (21.04.2006 18:22 213.184.233.178, пустое)
              • Просто умножьте коэффициенты на 2^N и округлите до целых. — andy_P   (21.04.2006 18:28 80.82.63.185, пустое)
                • А как сами рассчитать коэффициенты? — ssassh   (21.04.2006 18:29 213.184.233.178, пустое)
                  • Ответ: на малабе так (xi значения аргумента функции, yi -значения функции в точках xi): — andy_P   (21.04.2006 18:39 80.82.63.185, 356 байт)
• Можно просто затабулировать, если ресурсы по памяти позволяют — rmn   (21.04.2006 12:17 195.19.43.20, пустое)

E-mail: info@telesys.ru