Провести
исследование внутренней сходимости численного интегрирования методами Симпсона
и трапеций различных функций, задаваемых с помощью функций языка C.
1.2 Подробное описание
задания
Постройте зависимости количества итераций от различный величин
критерия точности.
Постройте обратные зависимости критерия точности от количества
итераций.
Повторите все вышеуказанные исследования для случая, когда при
вычислении критерия точности разность значений интеграла на смежных итерациях
относится не к предыдущему значению, а к точному значению аналитически
вычисленного интеграла.
Исследуйте влияние увеличения верхнего предела интегрирования на
точность (при прочих неизменных величинах).
Предполагается,
что отрезок интегрирования [a,b] разбит на n равных частей системой точек (сеткой). Контроль внутренней
сходимости заключается в циклическом вычислении приближенных значений интеграла
для удваиваемого по сравнению со значением на предыдущем прохождении цикла
числа п. Отношение абсолютной величины разности этих значений к абсолютной
величине предыдущего приближенного значения принимается в качестве критерия достижения
точности вычисления интеграла.
Приведем
пример построения зависимости количества итераций от критерия точности для
подынтегральной функции 1/x, решенный при помощи метода трапеций, критерий точности
вычисляется как отношение разности значений интегралов, полученных на смежных
итерациях, к аналитическому значению.
Входные
данные: номер зависимости (в данном случае – 1), номер делителя (аналитическое
значение -2), номер метода, номер функции, верхний и нижний пределы,
коэффициент m,
равный 1.
Критерий
точности автоматически изменяется от 10-7 до 10-3, с
каждым шагом увеличиваясь в 101/24, таким образом, получается 25
точек, которых вполне хватает для исследования зависимости. На каждой итерации
кол-во отрезков, на которое делится отрезок интегрирования увеличивается в два
раза, т.е n
= 2iter, где itter – кол-во итераций. Циклически вычисляется до достижения заданного
критерия точности.
Вычисленные
значения сохраняются в файле “D:\Zavisimost1.txt”. После переноса значений в MS Excel создается таблица 1 и
строится график, изображенный на рисунке 1.
Таблица 1
критерий точности
1,00E-07
1,47E-07
2,15E-07
3,16E-07
4,64E-07
6,81E-07
1,00E-06
1,47E-06
2,15E-06
Кол-во итераций
22
21
20
20
19
19
18
18
17
критерий точности
3,16E-06
4,64E-06
6,81E-06
1,00E-05
1,47E-05
2,15E-05
3,16E-05
4,64E-05
6,81E-05
Кол-во итераций
17
16
16
15
15
14
14
13
13
критерий точности
1,00E-04
1,47E-04
2,15E-04
3,16E-04
4,64E-04
6,81E-04
1,00E-03
Кол-во итераций
12
12
11
11
11
10
10
Рисунок 1
Результаты
всех измерений приведены в таблицах 2, 3, 4, 5, 6.
Таблица 2 - Зависимости кол-ва итераций от критерия точности для
функций 1/x, 1/x^2, sin x