Точность sin, cos

@MaxCPP · Регистрация: 26.07.2013

Студворк — интернет-сервис помощи студентам

Вычисляя sin(M_PI) в C++ я получаю 1.22465e-016. Но надо 0, а не число, близкое к 0. Можно ли какими-нибудь настройками компилятора, линковщика, #define каким-нибудь сделать результат вычисления sin(M_PI) равным 0? Не хочется делать "костыль" с проверкой...

@castaway · 30.11.2013, 19:10

А проверять ты как собираешься? sin( M_PI ) == 0 ?

@MaxCPP · 30.11.2013, 19:14 **[ТС]**

дело в выводе на экран... Если проверять, то так:
to_string(sin(M_PI))=="0"
главное, чтобы вывод был читабельным и максимальной точности. Ясно дело, что 1.22465e-016 близко к 0, но это не читабельно.

@zss · 30.11.2013, 19:18

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#include <stdio.h>
#include <conio.h>
#define _USE_MATH_DEFINES
#include <math.h>
double MySin(double x)
{
    if ( x>=M_PI && x<=2.*M_PI )
        return -sin(x-M_PI);
    else
        return sin(x);
}
int main()
{
  printf("%le",MySin(M_PI));
  getch();
  return 0;
}

Т.е. аргумент x приводим к интервалу [0,пи].
Тогда сходимость разложения в ряд Тейлора наилучшая

@Nick Alte · 30.11.2013, 19:20

Точность числа типа double составляет 15-16 десятичных цифр. То есть, число "пи" представляется у нас с точностью где-то до 15 знака после запятой. Соответствующим(и даже на один знак лучшим) получается и отклонение от ожидаемого результата, если учесть, что в этом регионе sin(x) почти полностью совпадает с прямой pi - x. Так что сначала надо найти способ хранить число "пи" с идеальной точностью, потом - производить расчёты с такой же точностью...

@castaway · 30.11.2013, 19:20

http://www.cplusplus.com/refer... precision/

@programina · 30.11.2013, 19:34

C++
// g++ file.cpp -o file -std=c++11
#include <iostream>
#include <cmath>
#include <cfenv>
using namespace std;
 
int main ()
{
    fesetround(FE_TONEAREST);
    cout << nearbyint(sin(M_PI)) << endl;
}

Новые блоги и статьи Все статьи Все блоги /
Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.
Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .

@MaxCPP 0 / 0 / 0 Регистрация: 26.07.2013 Сообщений: 9

	Точность sin, cos 30.11.2013, 19:01. Показов 6585. Ответов 6 Метки нет (Все метки) Вычисляя sin(M_PI) в C++ я получаю 1.22465e-016. Но надо 0, а не число, близкое к 0. Можно ли какими-нибудь настройками компилятора, линковщика, #define каким-нибудь сделать результат вычисления sin(M_PI) равным 0? Не хочется делать "костыль" с проверкой... 0

@castaway 4986 / 3093 / 456 Регистрация: 10.11.2010 Сообщений: 11,170 Записей в блоге: 10
	30.11.2013, 19:10
	А проверять ты как собираешься? sin( M_PI ) == 0 ? 0

@MaxCPP 0 / 0 / 0 Регистрация: 26.07.2013 Сообщений: 9
	30.11.2013, 19:14 [ТС]
	дело в выводе на экран... Если проверять, то так: to_string(sin(M_PI))=="0" главное, чтобы вывод был читабельным и максимальной точности. Ясно дело, что 1.22465e-016 близко к 0, но это не читабельно. 0

@Nick Alte 1675 / 1047 / 174 Регистрация: 27.09.2009 Сообщений: 1,945
	30.11.2013, 19:20
	Точность числа типа double составляет 15-16 десятичных цифр. То есть, число "пи" представляется у нас с точностью где-то до 15 знака после запятой. Соответствующим(и даже на один знак лучшим) получается и отклонение от ожидаемого результата, если учесть, что в этом регионе sin(x) почти полностью совпадает с прямой pi - x. Так что сначала надо найти способ хранить число "пи" с идеальной точностью, потом - производить расчёты с такой же точностью... 1

@programina 2062 / 619 / 41 Регистрация: 23.10.2011 Сообщений: 4,468 Записей в блоге: 2
	30.11.2013, 19:34
	C++ // g++ file.cpp -o file -std=c++11 #include <iostream> #include <cmath> #include <cfenv> using namespace std; int main () { fesetround(FE_TONEAREST); cout << nearbyint(sin(M_PI)) << endl; } 0