Отправлено bp 27 апреля 2006 г. 18:01

Сначала посмотрел стандартный метод Ньютона - все прекрасно в матлабе сходится, 6 итерация мне хватает. Только в этом случае используется деление. Есть ли методы его избежать ?

Посмотрел код в DSPLIB (с55) и вообще ничего не понял.
Во первых откуда они взяли такую формулу ?

Implementation Notes
The square root of a number(x) can be calculated using Newton’s method. An initial approximation is guessed and then the approximation gets recomputed using the formula,
new approximation = old - (old^2 -x)/2
The new approximation then becomes the old approximation and the process is repeated until the desired accuracy is reached.

И то что в самом коде - совсем непонятно. Если работаем с Q15 то причем здесь мантисса с экспонентой ?
MANT AC0, AC0
::NEXP AC0, T3 ;normalize the input
Потом еще почему-то появляются обращения к таблице.

• Ответ: Я тоже на них смотрел, потом мне все надоело — Alex11   (28.04.2006 00:28 195.131.139.203, 183 байт)
• sqrt(x)=x/sqrt(x). 1/sqrt(x) предпочтительнее для итерационного вычисления методом Ньютона-Рафсона, чем sqrt(x). Есть другие варианты. Аппроксимировать sqrt(x), x от 0.5 до 1, полиномом. 0.5-1 - стандартный диапазон, к которому производится нормализация. Мантисса это.Полином для 16-ти разрядной точности будет 5, кажется, порядка. Экспонента учитывается табличкой из sqrt(2^n). — -=ВН=-   (28.04.2006 00:06 213.177.101.235, пустое)
  • Сделал полиномом 4ой степени. Спасибо за совет! — bp   (28.04.2006 13:59 81.3.189.234, пустое)
  • Про скорость с полиномом. На tms54 корень из 32-х разрядного с 16-разрядным результатом занимал 30, с чем-то, циклов, если память не врет. На 55-м не делал, но из общих соображений должно быть лучше:-) — -=ВН=-   (28.04.2006 00:59 213.177.101.235, пустое)
• Это Ньютон-Рафсон для уточнения 1/sqrt(x). Формула оттуда. Мантисса с экспонентой и таблица - получение начального примерного значения (1/корень из ближайшего 2^n) — SM   (27.04.2006 18:36 195.225.131.186, пустое)

E-mail: info@telesys.ru