Формула трапеций. В формуле трапеций площадь криволинейной трапеции S заменяется суммой площадей элементарных трапеций (см

В формуле трапеций площадь криволинейной трапеции S заменяется суммой площадей элементарных трапеций (см. рис.2). Такой подход приводит к вычислению определённого интеграла I по формуле

или, раскрывая сумму,

.

Реализация формулы трапеций

Основой реализации любого вычислительного метода в программе Excel является представление его алгоритма в виде таблицы. При этом часто используется следующий подход: столбцами таблицы являются последовательно получаемые промежуточные результаты расчёта при одном значении аргумента или на одной итерации, а строками – последовательность значений аргумента или итерации. Иногда из этой закономерности выпадает первая строка таблицы, так как её используют для ввода начальных значений.

Описанный выше подход к формированию расчётной таблицы не является догмой. Каждый пользователь вправе использовать свой подход, но обычно он в той или иной мере основывается на описанных выше принципах.

Пример. Методом трапеций вычислить в программе Excel интеграл

с шагом интегрирования h = 0.1.

Решение. Одна из возможных реализаций формулы трапеций представлена таблицей 1.

Таблица 1.

i	xi
	x 0
	x 1
…	….	………
n	xn

В этой таблице столбец «i» не несет никакой информационной нагрузки, кроме пояснения к используемым формулам, а потому при реализации в программе Excel может быть опущен. Таблица интегрирования будет выглядеть следующим образом (см. рис.3)

Рис.3.

Результат интегрирования располагается в ячейке C11, где он вычисляется с помощью встроенной функции СУММ. Выбор этой ячейки для записи результата может быть произвольным и делается только из соображений эстетики.

Оценка погрешности численного интегрирования

Точность вычисления любой величины определяется погрешностью, которая может быть представлена в абсолютной или относительной форме. Абсолютная погрешность величины есть модуль разницы между её точным и приближённым значениями. Например, между точным значением определённого интеграла и его значением, полученным выбранным численным способом при конкретном количестве участков разбиения отрезка [ a, b ]

.

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

.

Анализ формулы трапеций численного интегрирования непрерывно диф­­ференциру­емых на отрезке [ a, b ] подынтегральных функций f (x) позволяет получить следующую оценку абсолютной погрешностей вычисления интегралов

.

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

,

где m для метода трапеций равно 2-м, а коэффициент C_m включает в себя длину участка интегрирования, максимум модуля производной и соответствующий коэффициент. Использование этой зависимости при уменьшении шага интегрирования вдвое позволяет записать оценку погрешности вычисления интеграла в виде

.

Сравнение последних двух формул даёт основное соотношение правила Рунге, справедливое для всех способов приближённого вычисления интеграла

,

где S (h /2) и S (h) – приближённые значения интеграла, вычисленные при шагах разбиения отрезка [ a, b ], отличающихся друг от друга в два раза. Исходя из этого, для оценки погрешности вычисленного значения интеграла с выбранным шагом надо повторить вычисления, удвоив величину шага, и воспользоваться приведённым выше соотношением.

Например, если есть необходимость определения погрешности полученного с шагом h = 0.125 значения рассмотренного в примере выше интеграла S (h = 0.125) = 0.6188, то расчёты надо повторить с шагом h = 0.25. Их результатом будет S (h = 0.25) = 0.6280. Таким образом, абсолютная и относительная погрешности вычисления приближённого значения интеграла при шаге h = 0.125 могут быть вычислены по правилу Рунге

,

.

Если постановка задачи требует получение результатов с меньшей погрешностью, чем была получена, то необходимо уменьшение величины шага интегрирования. Однако этот процесс нельзя продолжать бесконечно. Он ограничивается точностью представления данных в ЭВМ: существует некоторое минимальное значение шага разбиения отрезка [ a, b ], дальнейшее уменьшение которого вызовет рост погрешности вычисления интеграла.