Вычисление суммы бесконечного ряда с заданной точностью.

1. Цель работы: ознакомиться с текстовым редактором С++; освоить простейшие язы­ковые конструкции ввода - вывода данных; приобрести практические навыки составле­ния и отладки программ с простыми переменными; знать основные предопределённые язы­ком типы данных и правила их описания; усвоить арифметические и логические операции; усвоить способы организации циклических вычислений и операторы управления.

2. Указания по подготовке к работе:

Необходимо ознакомиться с основными предопределёнными типами данных (целочисленные типы, вещественные типы, тип char и логический тип (bool)) и операторами управления (условными и цикла), а также операторами ввода-вывода. [см. лекцию 4, 5, 7, 8].

3. Теоретическая часть:

Часто в математических вычислениях возникает необходимость определить сумму бесконечного ряда с заданной точностью. На такие задачи налагается одно жесткое условие: ряд должен быть сходящимся, т.е. абсолютное значение текущего элемента ряда должно быть меньше абсолютного значения его предыдущего элемента.

В подобных случаях удобно пользоваться итерационными циклами, реализуемыми операторами цикла с предварительным и последующим условиями.

В противном случае условие выхода из цикла никогда не будет выполнено и цикл станет бесконечным.

Пример: написать программу вычисления суммы бесконечного ряда с точностью eps=0.001:

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

Этот ряд сходится при |x| < вычисление суммы бесконечного ряда с заданной точностью. - student2.ru , а так как он знакопеременный, то его предел равен 0 (значения его элементов стремятся к нулю).Для достижения заданной точности требуется суммировать члены ряда, абсолютная величина которых больше eps. Для сходящегося ряда модуль члена ряда C вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru при увеличении n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru стремится к нулю. При некотором n неравенство

| C вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru | вычисление суммы бесконечного ряда с заданной точностью. - student2.ru eps перестает выполняться и вычисления прекращаются.

Общий алгоритм решения этой задачи очевиден: требуется задать начальное значение суммы ряда, затем многократно вычислять очередной член ряда и добавлять его к ранее найденной сумме. Вычисления заканчиваются, когда величина очередного члена ряда станет меньше заданной точности. Алгоритм представим в виде блок-схемы (см. рис.1). Как видно из рис.1. каждый блок имеет один вход и один выход за исключением условного блока (ромб), который имеет два выхода. Выходы могут располагаться или по обе стороны блока, или любой из них по центру. Ошибки в этом нет. Выходы обязательно должны иметь названия: либо «истина» («да», «+»), либо «ложь» («нет», «-«). Это необходимо для однозначного понимания логики алгоритма.

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

Линии потока, указывающие последовательность выполнения действий, должны располагаться параллельно границам листа (только горизонтально или вертикально), не могут пересекаться друг с другом. Если линии указывают естественный ход выполнения действий (сверху вниз и слева направо), то их можно не обозначать стрелками, а если снизу вверх или справа налево, то стрелками обозначать обязательно.

Блоки, кроме блоков «начало» и «конец» («останов»), можно нумеровать, что также улучшает наглядность схемы.

До выполнения программы предсказать, сколько членов ряда потребуется просуммировать, невозможно. В цикле такого рода есть опасность, что он никогда не завершиться – как из-за возможных ошибок в вычислениях, так и из-за ограниченной области сходимости ряда. Поэтому для надежности программы нужно предусмотреть аварийный выход из цикла с печатью предупреждающего сообщения по достижении некоторого максимально допустимого количества итераций. В качестве признака сходимости ряда будем использовать булеву переменную flag. В случае, если ряд расходится, значение этой переменной обнулится (примет значение «ложь»). Напомним, что для расходящегося ряда значение суммы его элементов не может быть подсчитано, а, следовательно, результат не должен быть выведен. Однако, в ранних версиях компиляторов языка С++ нет такого типа данных, как булевый, поэтому вместо него можно использовать либо тип int или самый короткий тип char.

#include <stdio.h>

#include <conio.h>

#include <math.h>

void main(){

int k, const Max = 500; // максимально допустимое количество

// итераций

float sum = 0., eps = 1.0e-3, s, x;

bool flag = true; // char flag = 1;

printf (“Введите значение х<1 ”); // можно использовать

// cout<<…

scanf(“%f”, &x); // можно использовать cin>>…

k = 2; // согласно условию

s = 5*eps; // значение нужно для входа в цикл

while (fabs(s)>eps)

{s = pow(x,k)/(k-1.); //вычисление абсолютного значения

// текущего элемента ряда

if ((k-1)%2 = = 0) // определение знака элемента ряда

sum=sum+s;

else

sum=sum-s;

k++;

if (k>Max){

printf(“Ряд расходится!”);

flag = false; // flag = 0;

break;

}

}

if (flag) printf(“Сумма равна %f\n”, sum); // if (flag == 1)…

}

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

Для уменьшения количества выполняемых действий следует воспользоваться рекуррентной формулой получения последующего члена ряда через предыдущий

Rn = вычисление суммы бесконечного ряда с заданной точностью. - student2.ru , где R – некоторый множитель

тогда

Cn=R вычисление суммы бесконечного ряда с заданной точностью. - student2.ru Cn-1.

Подставив в эту формулу значения Cn и C вычисление суммы бесконечного ряда с заданной точностью. - student2.ru получим:

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

Множитель (-1)k-1 сокращается и в итоге получается просто «-», который мы, опустив эти вычисления, сразу указали в рекуррентном соотношении.

Вот эту формулу мы подставим в текст программы и получим:

. . .

k = 3; // увеличили номер члена ряда на 1

s = - x*x; // нашли по формуле значение 1-го члена ряда

while (fabs(s)>eps)

{sum+=s;

s*=-( x*(k-2.)/(k-1.)); //вычисление значения

// текущего элемента ряда по предыдущему

k++;

if (k>Max){

cout <<“\nРяд расходится!”;

flag = false; // или flag = 0;

break;

}

} . . .

Максимально допустимое количество итераций удобно задать с помощью именованной константы (Max). Для аварийного выхода из цикла применяется оператор break, который выполняет выход на первый после цикла оператор.

Поскольку выход как в случае аварийного, так и в случае нормального завершения цикла происходит в одну и ту же точку программы, вводится булева переменная flag, которая предотвращает печать значения функции после выхода из цикла, в случае, когда точность вычислений не достигнута. Создание подобных переменных-«флагов», принимающих значение «истина» в случае успешного окончания вычислений и «ложь» в противном случае, является распространенным приемом программирования.

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

Рис.1. Блок-схема нахождения суммы знакопеременного бесконечного ряда

Задание.

Для выданного преподавателем варианта задачи написать и отладить программу на языке С++ расчета суммы знакопеременного ряда с использованием рекуррентного соотношения.

Измените программу так, чтобы она печатала не только значения аргумента и функции, но и количество просуммированных членов ряда.

Реализуйте ту же логику без специальной булевой переменной, объединив проверку обоих вариантов выхода из цикла в заголовке.

5. Требования к отчету по лабораторной работе:

Отчет должен содержать:

1) расчеты рекуррентного множителя для Вашей задачи;

2) блок-схему алгоритма задачи

3) распечатки или тексты программ с комментариями.

6. Варианты индивидуальных заданий:

.

1. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-2.

2. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru = 10-3.

3. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

4. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

5. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-5.

6. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3 , X0=0.3 вычисление суммы бесконечного ряда с заданной точностью. - student2.ru .

7. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3 , X0= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru /9.

8. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n+1

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-8, X0=1.6 .

9. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n+1

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-8, X0=1.22.

10. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4, X0=0.6785. вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

11. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3.

12. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n+1 вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

13. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3.

14. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru n3n-1(1/4)n-1(-1)n

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3.

15. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-5.

16. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-5, в точке X0=2.44.

17. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-6.

18. Вычислить

S=4 вычисление суммы бесконечного ряда с заданной точностью. - student2.ru (-1)n вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-2.

19. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4 в точке X0=2.334.

20. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-3.

21. Вычислить сумму ряда

S=2 вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4 в точке X0=0.4357.

22. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-2.

23. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

24. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-6.

25. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-6.

26. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

27. Вычислить сумму ряда

S= вычисление суммы бесконечного ряда с заданной точностью. - student2.ru вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru =10-4.

28. Вычислить сумму ряда

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru , в точке вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

29. . Вычислить сумму ряда

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru , в точке вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

30. . Вычислить сумму ряда

вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

с точностью вычисление суммы бесконечного ряда с заданной точностью. - student2.ru

7. Контрольные вопросы.

1. Как ввести начальные значения переменных.

2. Какое место расположение операторов описания типа в программе.

3. Какие формы представления чисел существуют в С++.

4. Какая последовательность выполнения операций в выражениях.

5. Что понимается под термином « тело цикла ».

6. Какого типа могут быть метки в операторе swith.

7. Какие существуют формы записи условного оператора? В чём их различие?

Лабораторная работа № 2.

Наши рекомендации