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

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

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 13, средняя оценка - 4.85
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
#1

Цикл while, разложение ряда - C++

03.12.2011, 01:52. Просмотров 1789. Ответов 24
Метки нет (Все метки)

Составить программу для вычисления функции c помощью разложения в ряд:


Формула прикреплена


Вычислить и вывести на экран в виде таблицы значения функции на интервале от xn до xk с шагом dx с точностью eps.

Для вычисления последующего члена ряда воспользоваться рекуррентной формулой Cn+1=Cn*T, где T - некоторый множитель.

Предусмотреть максимальное количество итераций = 500.

Протестировать работу программы для следующих значений:

xn xk dx eps
3.5 5.8 0.2 1e-20
1.5 3.5 0.5 1e-15

Результат получить в виде:
-----------------------------------------------

| x | F | n | math F |

-----------------------------------------------

| -3.50 | 2.80254 | 11 | 2.81324 |

| -3.00 | 2.60142 | 13 | 2.60287 |
...
-----------------------------------------------
где: x - значение аргумента, F - значение функции, n - к-во просуммированных членов ряда, вычисление той же ф-ии спомощью библиотеки math
Для ввода/вывода используйте библиотеку stdio.


все из лабы сделал только это задание осталось вообще не представляю как делать(
Вот код программы похожей задачи, но там формула другая с арктангенсом

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
#include <iostream>
#include  <math.h>
#include <iomanip>
using namespace std;
 
doublexn, xk,  dx,  eps;
 
void main()
{
    cout<<"Vveditexn,xk,dx,eps"<<endl;
    cin>>xn>>xk>>dx>>eps;
cout<<"-------------------------------------"<<endl;
cout<<"|   x  |     f    |   n  |  math F  |"<<endl;
cout<<"-------------------------------------"<<endl;
double ck;
double S, T;
int n;
double x, mathF;
        x=xn;
while(x<=xk)
        {
                n=0;
            ck=-1/x;
                S=1.57;
while (fabs(ck)>eps&& n<=500)
                {
                        S=S+ck;
                        T=((-1.0)*(2.0*n+1.0)/(((2.0*n+3.0)*(x*x))));
ck=ck*T;
n++;
                }
mathF=atan(x);
            if (n<=500)
cout<<"|"<<setw(6)<<x<<"| "<<setw(8)<<S<<" |"<<setw(6)<<n<<"| "<<setw(8)<<mathF<<" |"<<endl;
                else
cout<<"|"<<setw(6)<<x<<"| "<<setw(8)<<S<<" |"<<setw(6)<<n-1<<"| "<<setw(8)<<mathF<<" |"<<endl; 
                x=x+dx;
        }
cout<<"-------------------------------------"<<endl;
system("PAUSE");
}
0
Миниатюры
Цикл while, разложение ряда  
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
03.12.2011, 01:52
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Цикл while, разложение ряда (C++):

Разложение экспоненты в ряд тейлора. Цикл в цикле - C++
Здравствуйте. Сама задача: составить программу для вычисления множества значений функции exp x на интервале −40,1 с шагом 0.5 не...

Цикл For, вычисление суммы ряда - C++
Всем привет, как правильно написать цикл для данного примера? S = sin(1) + sin(1+2) + sin(1+2+3) + ... + sin(1+2+...+n)

Вычислить сумму ряда (цикл while()) - C++
Помогите пожалуйста, напишите программу по примеру, прикреплённому ниже, через цикл While. по возможности, напишите пожалуйста так же,...

Используя цикл for вычислить сумму ряда - C++
8. Даны натуральное n и действительное х, Вычислить sinx+sin^2x+.....sin^nx . Использовать оператор for

Цикл с предусловием (Вычислить сумму ряда) - C++
Здравствуйте уважаемые форумчане. Вот задача. Не могу её реализовать, формулы больно пугают, теряюсь с чего начать. Помогите чем сможете....

Цикл: вычислить сумму знакопеременного ряда - C++
#include &quot;stdafx.h&quot; #include&lt;iostream&gt; #include &quot;math.h&quot; using namespace std; const double e=0.00001; int _tmain(int argc,...

24
-=ЮрА=-
Заблокирован
Автор FAQ
04.12.2011, 20:32 #16
Цитата Сообщение от mr_aleks Посмотреть сообщение
в формулу я подстовлял вашу (2.0*n - 1)/(2.0*n + 1)/(x*x) , посмотрите код на 1 стр, в самом низу!
- вот моя формула в коде без функций с ограниченим в 500 итераций выводящим всё в аккуратную табличку и ввёл mathF
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
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
int main()
{
    long n;
    double xn, xk, dx, e;
    cout<<"xn = ";cin>>xn;
    cout<<"xk = ";cin>>xk;
    cout<<"dx = ";cin>>dx;
    cout<<" e = ";cin>>e;
    double sum;
    double an, x = xn;
    cout<<"|  x   |  Sum   | mathF  |    N   |\n";
    while(x <= xk)
    {
        n   = 1;
        sum = 0;
        an  = 2/x;
        do
        {
            sum += an;
            n   += 1;
            an *= (2.0*n - 1)/((2.0*n + 1)*(x*x));
        }
        while(e <= fabs(an) && n <= 500);
        cout<<setw(6)<<setprecision(3)<<xn<<" | "
            <<setw(6)<<setprecision(3)<<sum<<" | "
            <<setw(6)<<setprecision(3)<<log((x + 1)/(x - 1))<<" | ";
        cout<<setw(4)<<n<<endl;
        x = x + dx;
    }
    system("pause");
    return 0;
}
1
Миниатюры
Цикл while, разложение ряда  
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
04.12.2011, 20:45  [ТС] #17
-только не пойму почему х всегда одинаковый выдает?
-все нашел, в выводе поставили xn надо x т.к. присваивали xn=x


И сверил с данными и расчетами вашими в MathCAD и не совсем сошлось, но да ладно может это и не принципиально
0
-=ЮрА=-
Заблокирован
Автор FAQ
04.12.2011, 21:09 #18
Цитата Сообщение от mr_aleks Посмотреть сообщение
присваивали xn=x
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
while(e <= fabs(an) && n <= 500);
cout<<setw(6)<<setprecision(3)<<xn<<" | "
- здесь нужно было <<x записать, я переделывал и не заметил

Добавлено через 2 минуты
Цитата Сообщение от mr_aleks Посмотреть сообщение
И сверил с данными и расчетами вашими в MathCAD и не совсем сошлось, но да ладно может это и не принципиально
- должно сходиться, почему есть расхождения не могу понять, выкладки проверял там придраться нечего. Самое непонятное что если записать ту сумму в MathCAD именно через суммирование то всё считается без погрешности
0
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
04.12.2011, 22:36  [ТС] #19
ну я уж точно не знаю, большое спасибо и на этом) думаю препод если есть погрешность скажет

Добавлено через 1 час 13 минут
Еще одна проблеммка)

посмотрел в след. лабе сказано
запрограммируйте с использованием следующих функций:
· функция T выполняет вычисление множителя из рекуррентной формулы;
· функция iter выполняет вычисление очередного члена ряда;


попробовал но n и сумму вообще не считает, суму считает но не правельно


поправте где надо, спасибо большое зарание

я скорее не правельно записал вычисление очередного члена ряда, можите сказать что за переменная "an" ?

Добавлено через 12 минут
Точнее вот так правельней будет, но все равно почемуто первый член ряда какойто заоблочный ((

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
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
 
 
double T(double n, double x)
{
    double a;
    a = (2.0*n - 1)/((2.0*n + 1)*(x*x));
    return a;
};
 
double inter (double an, double T)
{
    an=an*T;
    return an;
};
 
 
int main()
{
    long n;
    double xn, xk, dx, e,a;
    cout<<"xn = ";cin>>xn;
    cout<<"xk = ";cin>>xk;
    cout<<"dx = ";cin>>dx;
        cout<<" e = ";cin>>e;
        double sum;
        double an, x = xn;
    cout<<"-----------------------------------"<<endl;
    cout<<"|   x  |     f  |   n  |  math F  |"<<endl;
    cout<<"-----------------------------------"<<endl;
    while(x <= xk)
        {
                n   = 1;
                sum = 0;
                an  = 2/x;
                do
                {
                        sum=sum+an;
                            an=inter(an,a);
                            a=T(n,x);
                        n   += 1;
                        
                }
                while(e <= fabs(an) && n <= 500);
                cout<<"|"<<setw(5)<<setprecision(3)<<x<<" | "
            <<setw(6)<<setprecision(4)<<sum<<" | "
                        <<setw(4)<<n<<" | ";
                cout<<setw(8)<<setprecision(3)<<log((x + 1)/(x - 1))<<" | "<<endl;
        x = x + dx;
        }
    system("pause");
    return 0;
}
0
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
04.12.2011, 22:39  [ТС] #20
Вот проверил с вашими данными!
0
Миниатюры
Цикл while, разложение ряда  
-=ЮрА=-
Заблокирован
Автор FAQ
04.12.2011, 23:39 #21
mr_aleks, ничего ты не проверил, моих данных нет - есть формула для итератора, которую я вывел. Т - это и есть итератор!!! Далее - что значит 183 итерации при х = 1,5
Цитата Сообщение от mr_aleks Посмотреть сообщение
Cn+1=Cn*T,
- эта запись означает
a[n + 1] = a[n]*T
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
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
double T(double n, double x)
{
    double a = (2.0*n - 1)/((2.0*n + 1)*(x*x));
    return a;
};
 
int main()
{
    long n;
    double xn, xk, dx, e;
    cout<<"xn = ";cin>>xn;
    cout<<"xk = ";cin>>xk;
    cout<<"dx = ";cin>>dx;
    cout<<" e = ";cin>>e;
    double sum;
    double an, x = xn;
    cout<<"-----------------------------------\n";
    cout<<"|  x   |  Sum   | mathF  |    N   |\n";
    cout<<"-----------------------------------\n";
    while(x <= xk)
    {
        n   = 1;
        sum = 0;
        an  = 2/x;
        do
        {
            sum += an;
            n   += 1;
            an *= T(n, x);
        }
        while(e <= fabs(an) && n <= 500);
        cout<<setw(6)<<setprecision(3)<<xn<<" | "
            <<setw(6)<<setprecision(3)<<sum<<" | "
            <<setw(6)<<setprecision(3)<<log((x + 1)/(x - 1))<<" | ";
        cout<<setw(4)<<n<<endl;
        x = x + dx;
    }
    system("pause");
    return 0;
}
0
Миниатюры
Цикл while, разложение ряда  
-=ЮрА=-
Заблокирован
Автор FAQ
04.12.2011, 23:47 #22
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
a[n + 1] = a[n]*T
- следующий а равен предыдущему а умноженному на итератор, вот эта строка в коде
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
an *= T(n, x);
Добавлено через 1 минуту
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
an *= T(n, x);
- сокращённая запись аn = an*T будь массив было бы
a[n + 1] = a[n]*T.
0
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
04.12.2011, 23:51  [ТС] #23
просто в заднии сказано ( я про след. лабу) что надо отдельную функцию написать на вычисление множителя из рекуррентной формулы (T) и на вычесление очередного члена ряда (iter)!
на вычесление ряда есть а на очередного члена нет...
Спасибо за разьяснение)

Или просто
an *= T
эту формулу записать в функцию?
0
-=ЮрА=-
Заблокирован
Автор FAQ
04.12.2011, 23:54 #24
Цитата Сообщение от mr_aleks Посмотреть сообщение
эту формулу записать в функцию?
- ну если так

Цитата Сообщение от mr_aleks Посмотреть сообщение
что надо отдельную функцию написать на вычисление множителя из рекуррентной формулы (T) и на вычесление очередного члена ряда (iter)
- тогда угу

Добавлено через 13 секунд
Цитата Сообщение от mr_aleks Посмотреть сообщение
эту формулу записать в функцию?
- ну если так

Цитата Сообщение от mr_aleks Посмотреть сообщение
что надо отдельную функцию написать на вычисление множителя из рекуррентной формулы (T) и на вычесление очередного члена ряда (iter)
- тогда угу
0
mr_aleks
0 / 0 / 0
Регистрация: 13.10.2011
Сообщений: 131
05.12.2011, 00:12  [ТС] #25
Вот так будет ?

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
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
double T(double n, double x)
{
        double a = (2.0*n - 1)/((2.0*n + 1)*(x*x));
    return a;
};
 
double iter (double an, double T)
{
        an*=T;
        return an;
};
 
int main()
{
        long n;
    double xn, xk, dx, e,a;
    cout<<"xn = ";cin>>xn;
    cout<<"xk = ";cin>>xk;
    cout<<"dx = ";cin>>dx;
        cout<<" e = ";cin>>e;
        double sum;
        double an, x = xn;
        cout<<"-----------------------------------\n";
    cout<<"|  x   |  Sum   | mathF  |    N   |\n";
        cout<<"-----------------------------------\n";
    while(x <= xk)
        {
                n   = 1;
                sum = 0;
                an  = 2/x;
                do
                {
                        sum += an;
                        n   += 1;
                        a=T(n,x);
                        an=iter(an,a);
                        
                }
                while(e <= fabs(an) && n <= 500);
                cout<<setw(6)<<setprecision(3)<<x<<" | "
            <<setw(6)<<setprecision(3)<<sum<<" | "
                        <<setw(6)<<setprecision(3)<<log((x + 1)/(x - 1))<<" | ";
                cout<<setw(4)<<n<<endl;
        x = x + dx;
        }
    system("pause");
    return 0;
}
Добавлено через 3 минуты
считает правельно)
0
05.12.2011, 00:12
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.12.2011, 00:12
Привет! Вот еще темы с ответами:

Цикл с предусловием (while), вычислить сумму ряда - C++
Вычислить сумму ряда с точностью Е, общий член которого {a}_{n}=\frac{{(-1)}^{n-1}}{{n}^{n}}

Используя цикл for найти сумму ряда - C++
как сделать это с помощью for?))

Сумма ряда, разработать цикл для выполнения - C++
#include &lt;iostream&gt; #include &lt;conio.h&gt; using namespace std; int main() { float x, sum, n = 24,tx=0; for (x = 42; x &gt; 0;x...

Найти сумму ряда по заданной формуле, используя цикл for - C++
Помогите решить задачу с помощью цикла for С++ Дано вещественное число x (|x|&lt;1) и натуральное число n. Сначала выведите значение суммы...


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

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

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