Проблема точности (как повысить точность вычислений)

Проблема точности (как повысить точность вычислений) - Сообщения

#1 Опубликовано: 08.11.2010 22:36:32
hexware

hexware

0 сообщений из 10 понравились пользователям.

Группа: User

Есть формула: http://ioj.com/v/jab5d(картинка)
При попытке просчитать её же в Smath c максимально выкрученной точностью: http://ioj.com/v/sknp2
в live:
[LIVE width=329 height=156]http://smath.info/live/?file=1830[/LIVE]
Вольфрам считает корректно:
http://www.wolframalpha.com/input/?i=%28%28%2810^-14*%28cos%2810^-14%29^%281%2F3%29%29-sin%2810^-14%29%29%2F%2810^-14%29^5%29%2B1%2F45%29%2F33.5097*10^-32
#2 Опубликовано: 14.11.2010 20:11:27
Andrey Ivashov

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

Здравствуйте.

Немного потестировал Ваш пример и вижу, что SMath Studio считает верно, чего не сказать о Wolframalpha. Пример результата вычислений в маткаде:


А вот работа сайта wolframalpha.com выглядит странновато. При первом вычислении получается:


При уточнении результата ответ уже другой:


Поведение сайта крайне странное и объяснить я его никак не могу.

С уважением, Андрей Ивашов.
#3 Опубликовано: 14.11.2010 21:48:55
hexware

hexware

0 сообщений из 10 понравились пользователям.

Группа: User

Здравствуйте!
В том то и дело! Что mathcad, smath, matlab и многие другие пакеты считают этот пример неверно.
Верно вычислили только: PARI/GP (когда задаёшь точность более 80. Скрин: http://ioj.com/v/tdv4u), maple и Mathematica (wolframalpha при нажатии на more digits увеличивает точность вычислений).
Пакеты в который используется динамический precision, в вычислениях накапливают значительно меньше ошибок (если вообще накапливается ошибка). В Smath используется IEEE754 просто по умолчанию Double precision, как и другие пакеты считающие с ошибкой.
Предлагаю подумать на эту тему для повышения точности вычислений Smath. Думаю есть готовые типы для .net'a с реализацией динамических или с более широким диапазоном точности. Как вы видите ошибка накапливается очень большая.
#4 Опубликовано: 15.11.2010 04:33:16
Andrey Ivashov

Andrey Ivashov

2270 сообщений из 3734 понравились пользователям.

Группа: Super Administrator

Спасибо за пояснение! Проблему понял - что-нибудь придумаю.
#5 Опубликовано: 05.12.2010 19:13:37
hexware

hexware

0 сообщений из 10 понравились пользователям.

Группа: User

более точно по теме: http://en.wikipedia.org/wiki/Arbitrary-precision
  • Новые сообщения Новые сообщения
  • Нет новых сообщений Нет новых сообщений