Модераторы: LSD, AntonSaburov
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Математические расчеты в Java, Проблемы с расчетом косинусов 
:(
    Опции темы
Гость_Владимир
Дата 1.11.2005, 11:57 (ссылка)    |    (голосов: 0) Загрузка ... Загрузка ... Быстрая цитата Цитата


Unregistered











День добрый!

Есть задачка, где нужен cos(Pi/2). Насколько я знаю из школьного курса математики, это равно 0.

Однако, получается следующая ерунда:

fi = (Math.Pi/2);
cos_fi = Math.cos(fi);

А в результате Cos угла fi оказывается равен 6.....Е-17.

Почему так?

Заранее благодарен.
С уважением.
  Вверх
LSD
Дата 1.11.2005, 12:01 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Leprechaun Software Developer
****


Профиль
Группа: Модератор
Сообщений: 15718
Регистрация: 24.3.2004
Где: Dublin

Репутация: 210
Всего: 538



Погрешность вычислений. ИМХО 17-й знак после запятой вполне приемлимая точность.


--------------------
Disclaimer: this post contains explicit depictions of personal opinion. So, if it sounds sarcastic, don't take it seriously. If it sounds dangerous, do not try this at home or at all. And if it offends you, just don't read it.
PM MAIL WWW   Вверх
Гость_Владимир
Дата 1.11.2005, 14:21 (ссылка)    |    (голосов: 0) Загрузка ... Загрузка ... Быстрая цитата Цитата


Unregistered











А если данное значение подставляется в формулу, то можно ли его считать нулем?

Заранее спасибо.
  Вверх
LSD
Дата 1.11.2005, 14:34 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Leprechaun Software Developer
****


Профиль
Группа: Модератор
Сообщений: 15718
Регистрация: 24.3.2004
Где: Dublin

Репутация: 210
Всего: 538



Если умножать некое большое число на 0.00000000000000006, то конечно ноль мы не получим. Но промежуточные результаты можно округлять например.


--------------------
Disclaimer: this post contains explicit depictions of personal opinion. So, if it sounds sarcastic, don't take it seriously. If it sounds dangerous, do not try this at home or at all. And if it offends you, just don't read it.
PM MAIL WWW   Вверх
Гость_Владимир
Дата 1.11.2005, 15:33 (ссылка)    |    (голосов: 0) Загрузка ... Загрузка ... Быстрая цитата Цитата


Unregistered











Благодарствую!
  Вверх
NotGonnaGetUs
Дата 1.11.2005, 15:40 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 92
Регистрация: 25.2.2005
Где: Москва

Репутация: 4
Всего: 12



Цитата
fi = (Math.Pi/2);
cos_fi = Math.cos(fi);

А в результате Cos угла fi оказывается равен 6.....Е-17.


Потому что число PI имеет бесконечное количество знаков после запятой smile
А ты мало того, что берёшь константу ограниченную 16 знаками, так ещё делишь её на два. А как известно:

cos(pi/2 + x) = - sin(x) = -x (при х->0)

вот ты и получаешь в результате число порядка 10^-16/17


PM MAIL   Вверх
val
Дата 1.11.2005, 15:59 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Program developer
**


Профиль
Группа: Участник Клуба
Сообщений: 992
Регистрация: 14.1.2003
Где: г. Киев

Репутация: нет
Всего: 7



Скорее всего вычисление происходит через нахождение суммы ряда, аппроксимирующего функцию косинуса, естественно, число членов этого ряда конечно и, как следствие, наличие погрешности аппроксимации, которая и отражается в семнадцатом знаке после запятой.


--------------------
Терпимость - величайшее благо человечества...
Ярчайший признак интеллекта – постоянно хорошее настроение…
PM MAIL ICQ   Вверх
NotGonnaGetUs
Дата 1.11.2005, 16:40 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 92
Регистрация: 25.2.2005
Где: Москва

Репутация: 4
Всего: 12



Цитата(val @ 1.11.2005, 15:59)
число членов этого ряда конечно

точность членов ряда тоже конечна.


как бы то ни было, что заказали, то и получили - величину отклонение fi от истинного значения pi/2 smile


вот пример:
[quot]
double fi = (Math.PI/2)+0.000001;
double cos_fi = Math.cos(fi);
System.out.println("cos: " + cos_fi);
[/quot]

Получаем "cos: -9.999999998563343E-7" - т.е. то самое отклонение.
PM MAIL   Вверх
Guest
Дата 4.11.2005, 12:03 (ссылка)    |    (голосов: 0) Загрузка ... Загрузка ... Быстрая цитата Цитата


Unregistered











Цитата(NotGonnaGetUs @ 1.11.2005, 15:40)
Цитата
fi = (Math.Pi/2);
cos_fi = Math.cos(fi);

А в результате Cos угла fi оказывается равен 6.....Е-17.


Потому что число PI имеет бесконечное количество знаков после запятой smile
А ты мало того, что берёшь константу ограниченную 16 знаками, так ещё делишь её на два. А как известно:

cos(pi/2 + x) = - sin(x) = -x (при х->0)

вот ты и получаешь в результате число порядка 10^-16/17



Хм. Почему же тогда:

fi = (Math.Pi/2);
sin_fi = Math.sin(fi);

В результате синус фи = 1. Т.е. с косинусом такие вещи не прокатывают, а с синусом со свистом?
  Вверх
maxim1000
Дата 4.11.2005, 13:55 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Участник
Сообщений: 3334
Регистрация: 11.1.2003
Где: Киев

Репутация: нет
Всего: 110



есть абсолютная погрешность, а есть относительная
при сохранении в double сохраняется мантисса и порядок числа
так что в случае с 0.00000001 сохраняется 1e-8 (если не обсчитался нулей)
а если 1.00000001, то придется сохранить 1.00000001e0
так вот точность самой мантиссы ограничена, 1 запомнить просто, а 1.00000001 - нужно достаточно бит...
Добавлено @ 13:56
т.е. в случае синуса тоже получилось 1.00...001 (например), но из-за погрешности это просто обрезалось...


--------------------
qqq
PM WWW   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Java"
LSD   AntonSaburov
powerOn   tux
javastic
  • Прежде, чем задать вопрос, прочтите это!
  • Книги по Java собираются здесь.
  • Документация и ресурсы по Java находятся здесь.
  • Используйте теги [code=java][/code] для подсветки кода. Используйтe чекбокс "транслит", если у Вас нет русских шрифтов.
  • Помечайте свой вопрос как решённый, если на него получен ответ. Ссылка "Пометить как решённый" находится над первым постом.
  • Действия модераторов можно обсудить здесь.
  • FAQ раздела лежит здесь.

Если Вам помогли, и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, LSD, AntonSaburov, powerOn, tux, javastic.

 
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Java: Общие вопросы | Следующая тема »


 




[ Время генерации скрипта: 0.0769 ]   [ Использовано запросов: 21 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.