Форум программистов, компьютерный форум CyberForum.ru

Найти ошибку (Вычисление функции с помощью разложения в ряд) - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 13, средняя оценка - 4.62
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
23.03.2013, 18:23     Найти ошибку (Вычисление функции с помощью разложения в ряд) #1
Итак,проблема с разницей в показаниях

Задание
Вычисление функции с помощью разложения в ряд
Вычислить и вывести на экран в виде таблицы значения функции, заданной с помощью ряда Тейлора, на интервале от хнач до хкон с шагом dx с точностью e . Таблицу снабдить заголовком и шапкой. Каждая строка таблицы должна содержать значение аргумента,значение функции, значение функции через ряд тейлора, количество просуммированных членов ряда,точность e.
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
#include <iostream>
#include <conio.h>
#include <clocale>
#include <cmath>
#include <iomanip>
 
 
using namespace std;
 
 
inline double my_abs(double val)
{
        double ret = val;
        if(val < 0)
                ret *= -1;
        return ret;
}
 
inline double Taylor(double x, double e, long &n)
{
        double sum = 0;
        double an = 1;
        n = 1;
        do
        {
                sum += an;
                n   += 1;
                an *= (-1)*(x*x/(2.0*n*2*(n - 1)));
        }
        while(e <= my_abs(an));
        return sum;
        
}
 
 
int main()
{
    setlocale(LC_ALL, "Russian");
 
    long n = 0;
    double xn, xk, dx, e;
 
    cout << "Программа для вычисления значения функции\ny = cos(x) с помощью ряда Тейлора\n";
 
    cout<<" xn = ";
    cin>>xn;
    cout<<"Конечное значение xk = ";
    cin>>xk;
    cout<<"Шаг dx = ";
    cin>>dx;
    cout<<"Погрешность e = ";
    cin>>e;
 
    cout << "-----------------------------------------------------------\n";
    cout << "|   x  | f(x) | f(x) в ряде Тейлора  |   n   | точность e | \n";
    cout << "-----------------------------------------------------------\n";
 
    while(xn <= xk)
    {
        
        cout<< "|" << setw(10)<< setprecision(8) << xn << "  |    "
            << setw(10) << setprecision(8) << cos(xn) <<"| "
            << setw(10) << setprecision(8) << Taylor(xn, e, n) <<"  | "
            << setw(10) << n << "|" 
            << setw(10) << e << "|" << endl;
        xn = xn + dx;
    }
    cout << "--------------------------------------------------------------------------\n";
 
    _getch();
    return 0;
}

Буду благодарен
Миниатюры
Найти ошибку (Вычисление функции с помощью разложения в ряд)  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Дмитрий 93
74 / 59 / 5
Регистрация: 03.06.2011
Сообщений: 309
23.03.2013, 19:06     Найти ошибку (Вычисление функции с помощью разложения в ряд) #2
зависит от вводимой точности
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
24.03.2013, 11:36  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #3
нет,не из-за точности. даже если e=0.000000000001 , разница между значениями огромна

Добавлено через 16 часов 4 минуты
Актуально
пытался разобраться с программой из Написать программу для нахождения значения функции с помощью ряда Тейлора темы,но безрезультатно
-=ЮрА=-
Заблокирован
Автор FAQ
24.03.2013, 22:53     Найти ошибку (Вычисление функции с помощью разложения в ряд) #4
valenti2, что от меня требуется? в заголовке темы мой код, толкьо не помню для какого разложения(+его кто то "довёл" до внестандартного, но это не важно), линки тоже на мои коды. Ошибка потмоу что тебе надо найти скажем разложение для тангенса а взял мой код для косинуса или ещё чего то там, неужели думаем что итератор
Цитата Сообщение от valenti2 Посмотреть сообщение
an *= (-1)*(x*x/(2.0*n*2*(n - 1)));
один и тоже для всех разложений. Уфф его выводить надо. Какую функцию раскладываешь?
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
25.03.2013, 11:10  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #5
Код переписал немного я
Мне нужно для косинуса, но вывод немного других данных.Но ошибки приличные,и не пойму в чём проблема
-=ЮрА=-
Заблокирован
Автор FAQ
25.03.2013, 13:24     Найти ошибку (Вычисление функции с помощью разложения в ряд) #6
Цитата Сообщение от valenti2 Посмотреть сообщение
Мне нужно для косинуса, но вывод немного других данных.Но ошибки приличные,и не пойму в чём проблема
тогда напиши конкретно задание

Не по теме:

PS:На скрине 123 это что значение в градусах?

valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
25.03.2013, 14:41  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #7
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
тогда напиши конкретно задание

Не по теме:

PS:На скрине 123 это что значение в градусах?

задание
Вычислить и вывести на экран в виде таблицы значения функции, заданной с помощью ряда Тейлора, на интервале от хнач до хкон с шагом dx с точностью e . Таблицу снабдить заголовком и шапкой. Каждая строка таблицы должна содержать значение аргумента,значение функции(сделал как cos(xn)), значение функции через ряд тейлора(твой код), количество просуммированных членов ряда,точность e.



123 это число
-=ЮрА=-
Заблокирован
Автор FAQ
25.03.2013, 15:09     Найти ошибку (Вычисление функции с помощью разложения в ряд) #8
Цитата Сообщение от valenti2 Посмотреть сообщение
123 это число
- это должна быть радианная мера угла вывести надо косинус от xn радиан до xk радиан 123 это 123 РАДИАНЫ не ГРАДУСА!(ЭТО ОЧЕНЬ МНОГО И НЕЛЬЗЯ БРАТЬ ТАКОЕ ЧИСЛО!)
Во-первых при вводе надо переводить угол в радианы, либо перводить градусную меру програмно а во вторых разложение в ряд Тэйлора валидно вблизи малых значений х-а т.е -2Пи туда 2Пи обратно и больше 7-ми х вообще не должен быть. Конечно с 123 радианами разложение даст огромнейшую погрешность.
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
25.03.2013, 15:47  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #9
Понял.На самом деле не обязательно переводить,главное знать про интервал
Но,опять же. функцию через cos(xn) считает нормально,а через ряд число другое
Скрин при интервале [-6;6]
Миниатюры
Найти ошибку (Вычисление функции с помощью разложения в ряд)  
-=ЮрА=-
Заблокирован
Автор FAQ
25.03.2013, 16:26     Найти ошибку (Вычисление функции с помощью разложения в ряд) #10
valenti2, Я то сделал задание, но мне интересно посмотреть на твою работу! Жду твоих соображений, а не вопросов почему не работает, тогда реально помогуПока посмотрим, что сделаешь самостоятельно.

Не по теме:

Не работает потому что был слепо сделан копипаст даже не вдумываясь, что в нём да как

Миниатюры
Найти ошибку (Вычисление функции с помощью разложения в ряд)  
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
25.03.2013, 17:11  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #11
при числе радиан,заметно меньших единицы, всё более-менее стабильно работает
модуль точно верный

Место,где может быть ошибка-работа ф-ции Taylor,куда потом подставляется xn
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
inline double Taylor(double x, double e, long &n)
{
        double sum = 0;
        double an = 1;
        n = 1;
        do
        {
                sum += an;
                n   += 1;
                an *= (-1)*(x*x/(2.0*n*2*(n - 1)));
        }
        while(e <= my_abs(an));
        return sum;
        
}
На этом мысли иссякают
-=ЮрА=-
Заблокирован
Автор FAQ
25.03.2013, 17:56     Найти ошибку (Вычисление функции с помощью разложения в ряд) #12
Цитата Сообщение от valenti2 Посмотреть сообщение
На этом мысли иссякают
Ладно лови код, да в разложении біла помарка, я коода выводил итератор стратил, никто меня не поправил и только твой вопрос заставил меня перепроверить свои же выкладки (получается все остальные тупо копипастили)
Вот код с корректным итератором (руссификацию заменяй как тебе угодно мне просто system больше по душе)
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
 
double Radian(double x);
double Period(double x);
double Taylor(double x, double e, long &n, double &an);
 
 
int main()
{
    long n;
    double x, y, xn, xk, dx, e, err;
 
    system("chcp 1251");//это нестандартная функция потому заменяйте на свою ф-цию руссификации
    cout<<"Программа для вычисления значения функции"<<endl
        <<"y = cos(x) с помощью ряда Тейлора"<<endl;
    cout<<"\tВвод значений углов в градусах"<<endl;
    cout<<"Начальное значение xn = ";cin>>xn;
    cout<<"Конечное значение  xk = ";cin>>xk;
    cout<<"Шаг аргумента      dx = ";cin>>dx;
    cout<<"Погрешность e = ";cin>>e;
    cout<<"Радианные меры углов : "<<endl;
    cout<<(xn = Radian(xn))<<endl;
    cout<<(xk = Radian(xk))<<endl;
    cout<<(dx = Radian(dx))<<endl;
    cout<<"Удаление целых периодов : "<<endl;
    cout<<(xn = Period(xn))<<endl;
    cout<<(xk = Period(xk))<<endl;
    cout<<(dx = Period(dx))<<endl;
 
    cout<<"|    x   |   N   |   Sum   |   cos   |   fabs(an)   |     e     |"<<endl;
    for(x = xn; x <= xk; x += dx)
    {
        y = Taylor(x, e, n, err);
        cout<<((x < 0) ? "|-" : "| ")
            <<setw(5)<<setprecision(3)<<fabs(x)<<" | "
            <<setw(5)<<n<<" | "
            <<((y < 0) ? "-" : " ")
            <<setw(6)<<setprecision(4)<<fabs(y)<<" | "
            <<((cos(x) < 0) ? "-" : " ")
            <<setw(6)<<setprecision(4)<<fabs(cos(x))<<" | "
            <<setw(8)<<setprecision(6)<<err<<" | "
            <<setw(8)<<setprecision(6)
            <<(cos(x) ? fabs((y - cos(x)) / cos(x)) : 0)<<" | "
            <<endl;
    }
    system("pause");//это нестандартная функция потому заменяйте на свою ф-цию стоп-точки
    return 0;
}
 
double Radian(double x)
{
    return x*asin(1.0) / 90;
}
 
double Period(double x)
{
    double period = 4*asin(1.0);
    while( x - period > 0)
        x -= period;
    return x;
}
 
double Taylor(double x, double e, long &n, double &err)
{
    double sum = 0;
    double an  = 1;
    n = 0;
    for(n = 0; e <= (err = fabs(an)); n++)
    {
        sum += an;
         an *= (-1)*x*x/((2*n + 1)*(2*n + 2));
    }
    return sum;
    
}

Не по теме:

valenti2, с тебя расписать как я вывел последний итератор, если сможешь то сам научишся их выводить

-=ЮрА=-
Заблокирован
Автор FAQ
25.03.2013, 17:59     Найти ошибку (Вычисление функции с помощью разложения в ряд) #13
Вот скрин на другой диапазон
Миниатюры
Найти ошибку (Вычисление функции с помощью разложения в ряд)  
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
31.03.2013, 20:33     Найти ошибку (Вычисление функции с помощью разложения в ряд)
Еще ссылки по теме:

Вычисление функции с использованием их разложения в ряд C++
Не могу найти ошибку в подсчете значения функции и разложения ее по ряду C++
C++ Вычисление функций с использованием их разложения в степенной ряд

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

Или воспользуйтесь поиском по форуму:
valenti2
0 / 0 / 0
Регистрация: 16.11.2012
Сообщений: 18
31.03.2013, 20:33  [ТС]     Найти ошибку (Вычисление функции с помощью разложения в ряд) #14
Ты вывел итератор почленно, циклом при an не меньше погрешности
Верно?

upd/ или вопрос был про функцию an?
Yandex
Объявления
31.03.2013, 20:33     Найти ошибку (Вычисление функции с помощью разложения в ряд)
Ответ Создать тему
Опции темы

Текущее время: 00:15. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru