1. Вычисление определенных интегралов Справочная информация

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

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

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

П

Рис.1.

оэтому для приближённого вычисления определённого интеграла доста­точно подсчитать пло­щадь трапецииS. Можно указать множество способов её вычисления. Простейшие из них – формулы прямоугольников, трапеций и Симпсона.

Формула средних прямоугольников

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

Рис.2.

или в развёрнутом виде

,

где – абсцисса центральной точкиi-го участка разбиения [a, b], h = (b – a)/n – шаг равномерного разбиения отрезка, n – количество участков разбиения.

Формула трапеций

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

Рис.3.

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

.

Формула Симпсона (j.Gregory(Грегори)1668,Th.Simpson1743)

В

Рис.4.

формуле Симпсона площадь криволинейной трапеции рассчитывается как сумма площадей ряда криволинейных трапеций, у которых криволинейная сторона представляет собой участок параболы, проходящей через три соседние точки на графике подынтегральной функции. Это можно видеть на рис.4. Поэтому число участков разбиения отрезка [a, b] в отличие от предыдущих способов обязательно должно быть чётным. Исходя из этого, значение определённого интеграла приближённо вычисляется по формуле

.

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

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

.

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

.

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

• для формулы средних прямоугольников ;

• для формулы трапеций ;

• для формулы Симпсона .

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

,

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

.

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

,

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

Например, если требуется вычислить интеграл

с относительной погрешностью, не превышающей 0.5%, то для выполнения задания можно воспользоваться формулой средних прямоугольников, выбрав в качестве первоначального значения шага интегрирования h = 0.5, т.е. разбив отрезок [0, 1] на два участка с граничными точками 0, 0.5 и 1.0. Тогда приближённое значение интеграла будет

.

При уменьшении шага интегрирования вдвое (h = 0.25), т.е. при разбиении отрезка интегрирования на 4 участка, значение интеграла будет равно

.

Таким образом, абсолютная и относительная погрешности вычисления последнего значения интеграла могут быть вычислены по правилу Рунге

,

.

В связи с тем, что требуемая точность вычисления интеграла не достигнута, решение продолжается: шаг интегрирования уменьшается ещё в два раза (h = 0.125). В этом случае значение интеграла будет равно 0.6141, а его абсолютная и относительная погрешности составят 0.0046 и 0.75%, соответственно. Видно, что требуемая точность опять не достигнута, поэтому шаг интегрирования уменьшается ещё в два раза (h = 0.0625) и вычисления повторяются. С таким шагом значение интеграла получается равным 0.6152, его абсолютная погрешность – 0.0011, а относительная – 0.18%. Таким образом, при шаге интегрирования h = 0.0625 решение задачи, равное 0.6152, получено с заданной погрешностью.

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