Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 31, средняя оценка - 4.74
Mirajjj
0 / 0 / 0
Регистрация: 22.10.2009
Сообщений: 5
#1

Из Sin в Cos С++ - C++

23.09.2010, 11:52. Просмотров 4149. Ответов 48
Метки нет (Все метки)

Здраствуйте, мне нужно написать програму которая за рядами Тейлора выводит косинус . Проблема заключаеться в том ,что если програма добавляет елементы сумы Тейлора то (как у меня было в разных интерпретациях цыкла) то после 30 елемента или 70 выдает -1,#IND , а значения аргумента будут вводится от 100. Порывшись на американский сайтах нашел програму которая почти идеально решает задачу для Синуса (она с циклом Тейлора только до 10, а значение начинают незначительно (0.0000001) расходиться когда вводимый аргумент больше 1 000 000. Я хотел переделать эту програму на Косинус, но увы, некоторые части кода вообще не понимаю зачем нужны, прошу помощи и обьяснений, Спасибо.

Вот сам код функций. Которую нашел, работает на ура, но это Синус.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
double SinHelp(double x)
{
    int i;
    double y=x, r=x;
        
        for (int i=0; i<10; i++)
            {
                y*=-x*x;
                r+=1.0/Factorial(1+2*(i+1))*y;
            }
    return r;
}
 
double SinCout(double number)//Функция которую буду выводить в главном файле
{
    double sign=1, x=number/3.141592653589793;
        if(x<0.0)
            {
                sign=-1;
                x=-x;
            }
        int i= static_cast<int>(x+0.5);
        double a=x-i;
            if((i-i/2*2)!=0)
                sign=-sign;
    return sign*SinHelp(a*3.141592653589793);
}
 
double Factorial(int num)
{
    double x=num, result=1,plus=1;
    for(int i=0; i!=x; ++i)
        {
            if(x!=0)
            {   result*=plus;
                ++plus;
            }
            if(x==0)
                result=1;
        }
    return result;
}
Очень интересует , что делает фрагмент кода :
C++
1
2
int i= static_cast<int>(x+0.5);
double a=x-i;
Прошу если кому не трудо ,не только сказать, а еще и проверить на своих компиляторах если значение аргумента косинуса будет больше 1000, будет ли оно соответствовать стандартной функцыи косинуса в С++.

Если кто не помнит...

Ряд Тейлора для Синуса

Из Sin в Cos С++

Ряд Тейлора для Косинуса

Из Sin в Cos С++
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.09.2010, 11:52
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Из Sin в Cos С++ (C++):

Функция sin(x+y)=sin(x)cos(y)+cos(x)sin(y) - C++
как её записать на языке с++?

y=sin(n*x)+cos(k*x)+ln(m*x) - C++
y=sin(n*x)+cos(k*x)+ln(m*x) где x= { e^z+z, при z&gt;1 z^2+1, при z меньше или равно 1 не поможете ли мне пожалуйста...

Sin() cos() - C++
Всем приветик. У меня последний Qt MinGW. Вопрос: Перед использованием функции sin(x), мне сказали что нужно градусы перевести в...

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

Выражение с cos() и sin() - C++
помагите с выражением.а то запуталась окончательно(( z=cos4x+sin y+1/4 sin22x - 1

Фунции sin и cos - C++
Здравствуйте. Начал писать игру, аналог знаменитых Tank Wars. При использовании библиотечных(math.h) функций sin и cos появилась...

48
KuKu
1558 / 1036 / 78
Регистрация: 17.04.2009
Сообщений: 2,987
23.09.2010, 13:46 #2
C++
1
static_cast<int>
это приведение типов, т.е x+0.5 переводит в тип int.

вводимый аргумент больше 1 000 000
, не знаю как считает эта американская программа, но ряд тейлора который у тебя на картинках - это разложение в районе 0, его нельзя применять для значений который "не около" нуля.
1
Mirajjj
0 / 0 / 0
Регистрация: 22.10.2009
Сообщений: 5
23.09.2010, 14:07  [ТС] #3
Цитата Сообщение от KuKu Посмотреть сообщение
C++
1
static_cast<int>
это приведение типов, т.е x+0.5 переводит в тип int.

, не знаю как считает эта американская программа, но ряд тейлора который у тебя на картинках - это разложение в районе 0, его нельзя применять для значений который "не около" нуля.
Да в том то и дело , что она работает так же как и функция в библеиотеке <cmath>, одинкойвый косинус на значении аргумента 1000 и больше , невериш сам попробуй.

Спасибо KuKu я знаю что это , самое интересное зачем оно ?

а значение косинуса всегда будет возле нуля, посмотри внимательно на код 1/... дело в том что при больших размерах вводимого значения аргумента выдает -1,#IND и я как бы незнаю что делать ...
0
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 14:30 #4
KuKu,
С чего бы это вдруг эти разложения для икс около 1? Для любого икса они работают.

Добавлено через 1 минуту
Возможно, вы спутали с арксинусом, арктангенсом и логарифмом, которые работают для икса от -1 до 1
1
Mirajjj
0 / 0 / 0
Регистрация: 22.10.2009
Сообщений: 5
23.09.2010, 14:32  [ТС] #5
Цитата Сообщение от silent_1991 Посмотреть сообщение
KuKu,
С чего бы это вдруг эти разложения для икс около 1? Для любого икса они работают.

Добавлено через 1 минуту
Возможно, вы спутали с арксинусом, арктангенсом и логарифмом, которые работают для икса от -1 до 1
The function sin returns the sine of arg, where arg is given in radians. The return value of sin will be in the range [-1,1]. If arg is infinite, sin will return NAN and raise a floating-point exception.

возвртное значение палюбому будет от -1 до 1
0
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 14:35 #6
Область определения sin(x) (cos(x)): [-1; 1]
Область значений этих же функций: [-inf; inf]
Мы говорим про область значений.
1
Mirajjj
0 / 0 / 0
Регистрация: 22.10.2009
Сообщений: 5
23.09.2010, 14:40  [ТС] #7
Цитата Сообщение от silent_1991 Посмотреть сообщение
Область определения sin(x) (cos(x)): [-1; 1]
Область значений этих же функций: [-inf; inf]
Мы говорим про область значений.
Извини. Ты прав. Просто вот надо надо както поставить вращение 2Пи. Мне кажется, что в 1 сообщении я привел код непоняный зачем , он може както вращет значение но я не понимаю принцип.
0
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 14:43 #8
Не понял, о каком "вращении" речь, до вы вообще видите разницу между двумя приведёнными разложениями? Вся разница заключается в том, что у синуса идут нечётные степени и факториал этих нечётных степеней в знаменателе, а у косинуса - чётные. Вся ваша задача - убрать +1 в коде вычисления синуча, чтобы получались чётные числа.
1
Mirajjj
0 / 0 / 0
Регистрация: 22.10.2009
Сообщений: 5
23.09.2010, 14:50  [ТС] #9
Цитата Сообщение от silent_1991 Посмотреть сообщение
Не понял, о каком "вращении" речь, до вы вообще видите разницу между двумя приведёнными разложениями? Вся разница заключается в том, что у синуса идут нечётные степени и факториал этих нечётных степеней в знаменателе, а у косинуса - чётные. Вся ваша задача - убрать +1 в коде вычисления синуча, чтобы получались чётные числа.
Прошу прощение в прошлой програме Я просто убрал +1 , а так как +1 в степени синуса не был прямо заметен была разница в вычислениях... по 10 распросматриваю програму на 11 доход что не так спасибо за советы.
0
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 15:06 #10
Вот программка, но с большими x она тоже не будет работать. Ведь скажем 100^2 - ещё нормально, 100^3 -тоже... А вот 100^20 - уже перебор, а с нормальной точностью нужно вычислить куда больше 10 членов.

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
#include <stdio.h>
#include <math.h>
 
#define epsilon 0.000000000001
 
double factorial(double);
double cosx(double);
 
int main()
{
    double x;
 
    printf("Input x: ");
    scanf("%lf", &x);
 
    printf("myCos(%lf) = %lf\nlibCos(%lf) = %lf", x, cosx(x), x, cos(x));
 
    getch();
    return 0;
}
 
double factorial(double n)
{
    double i;
 
    for (i = n - 1.0; i > 1.0; i -= 1.0)
        n *= i;
 
    return (n == 0.0) ? 1.0 : n;
}
 
double cosx(double x)
{
    double sum;
    double ak;
    double n;
 
    sum = 0.0;
    ak = 1.0;
    n = 1.0;
 
    do
    {
        sum += ak;
        ak = pow(-1.0, n) * pow(x, 2.0 * n) / factorial(2.0 * n);
        n += 1.0;
    }
    while (fabs(sum - ak) >= epsilon);
 
    return sum;
}
0
Евгений М.
23.09.2010, 15:11
  #11

Не по теме:

Цитата Сообщение от silent_1991 Посмотреть сообщение
Область определения sin(x) (cos(x)): [-1; 1]
Область значений этих же функций: [-inf; inf]
Вы перепутали. Наоборот.

1
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 15:12 #12
Евгений М.,
Разумеется, прошу прощения)))
0
KuKu
1558 / 1036 / 78
Регистрация: 17.04.2009
Сообщений: 2,987
23.09.2010, 17:31 #13
Цитата Сообщение от silent_1991 Посмотреть сообщение
KuKu,
С чего бы это вдруг эти разложения для икс около 1? Для любого икса они работают.

Добавлено через 1 минуту
Возможно, вы спутали с арксинусом, арктангенсом и логарифмом, которые работают для икса от -1 до 1
цифра 1 от меня не звучала, я сказал что это разложение Tейлора которое на картинках верно только в окрестности нуля. Возможно это вы спутали, говорил не про область значений синуса и косинуса, а о применимости ряда Тейлора.
0
silent_1991
Эксперт С++
4986 / 3043 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
23.09.2010, 17:34 #14
KuKu,
С 1 промахнулся, имел ввиду 0. Я правильно понял, вы имели ввиду, что аргумент (x) должен быть по модулю близок к нулю, иначе результат будет неверным? Если так, то это не правильно, этот ряд даёт верный результат для любого x. Если нет, то я чего-то не понимаю...
0
KuKu
1558 / 1036 / 78
Регистрация: 17.04.2009
Сообщений: 2,987
23.09.2010, 17:39 #15
Цитата Сообщение от Mirajjj Посмотреть сообщение
Да в том то и дело , что она работает так же как и функция в библеиотеке <cmath>, одинкойвый косинус на значении аргумента 1000 и больше , невериш сам попробуй.

Спасибо KuKu я знаю что это , самое интересное зачем оно ?

а значение косинуса всегда будет возле нуля, посмотри внимательно на код 1/... дело в том что при больших размерах вводимого значения аргумента выдает -1,#IND и я как бы незнаю что делать ...
процитируй пожалуйста эту функцию, которая представляет собой ряд тейлора и работает для значений много больше нуля.

Добавлено через 1 минуту
Цитата Сообщение от silent_1991 Посмотреть сообщение
KuKu,
С 1 промахнулся, имел ввиду 0. Я правильно понял, вы имели ввиду, что аргумент (x) должен быть по модулю близок к нулю, иначе результат будет неверным? Если так, то это не правильно, этот ряд даёт верный результат для любого x. Если нет, то я чего-то не понимаю...
думаю вы чего то не понимаете , ряд тейлора это разложение вокруг какой то точки, а не по всей области.

Добавлено через 1 минуту
silent_1991, чуть не засомневался Ряд Тейлора, посмтрите что такое ряд Тейлора.
0
23.09.2010, 17:39
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
23.09.2010, 17:39
Привет! Вот еще темы с ответами:

Определить sin и cos - C++
Считая, что функции sin и cos применимы только к аргументам в диапазоне , вычислить y = sin(x) для любого заданного вещественного числа x.

Символьное дифференцирование cos и sin - C++
Здравствуйте. Хочу написать символьное дифференцирование cos и sin. Производную sin определяю через cos, cos через sin. Ошибка понятна:...

Вычисление sin, cos, tg, ctg - C++
хотел написать прогу для вычисления sin, cos, tg, ctg по градусной мере. но никак не получается выдаёт не та значения. проверял по таблице...

Не хочет посчитать cos и sin - C++
Пишу вот такой цикл: int A=5; int alpha=1; for(int i=0; i&lt;4; i++){ int re=A*std::cos(i*alpha); int im=A*std::sin(i*alpha);}...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
15
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru