Лекция-5
ОСНОВЫЧИСЛЕННЫХ МЕТОДОВ
Методы численного интегрирования
Понятие о формулах численного интегрирования.
Пусть требуется вычислить определенный интеграл вида
. (2.5.1)
Для многих функций 
первообразные представляют собой достаточно сложные комбинации элементарных функций, либо вовсе не выражаются через них. В таких случаях использование формулы Ньютона-Лейбница на практике не представляется возможным. Во многих практических случаях достаточно получить значение интеграла с заданной точностью 
. Для вычисления приближенного значения интеграла существуют формулы численного интегрирования. Суть построения формул численного интегрирования состоит в следующем.
Разобьем отрезок 
на 
частей. Для простоты изложения положим эти части одинаковой длины 
:
Пронумеруем точки разбиения так, как показано на рис. 2.5.1. Имеем:
,
при этом
.
Рис. 2.5.1. К вопросу о численном интегрировании.
Исходный интеграл (2.5.1) может быть представлен в виде суммы интегралов по полученным в результате разбиения «малым» отрезкам:
. (2.5.2)
Интегралы
(2.5.3)
вычисляются по приближенным формулам.
Простейшие формулы для приближенного вычисления интегралов по отрезку называются квадратурными формулами. Рассмотрим некоторые из них ниже, а также изучим вопросы их точности. Порядок точности квадратурной формулы определяется степенью полинома (многочлена), для которой эта квадратурная формула точна.
Формула прямоугольников (формула «средних»).
Заменим на i -ом участке интегрируемую функцию постоянной величиной, например, равной ее значению в средней точке (рис. 2.5.2):
Рис. 2.5.2. К интегрированию по формуле прямоугольников.
, где 
. (2.5.4)
Тогда интеграл на отрезке заменяется площадью прямоугольника, т.е.
, (2.5.5)
и вычисление исходного интеграла сводится к вычислению суммы
. (2.5.6)
Кроме того, часто из практических соображений в качестве 
в формуле (2.5.6) берется 
, либо 
. В результате получаем:
(2.5.7)
– квадратурная формула «левых» прямоугольников;
(2.5.8)
– квадратурная формула «правых» прямоугольников.
Формулы (2.5.7) и (2.5.8) менее точные, чем (2.5.6), но иногда более удобные, например, при численном решении дифференциальных уравнений.
Точность вычисления. Как следует из построения квадратурные формулы прямоугольников дают точный результат интегрирования для функций, постоянных на i -ом участке (
). Квадратурная формула «средних» прямоугольников дает точный результат также и для линейных на i -ом отрезке функций. Это утверждение достаточно проверить для простейшей линейной функции 
.
При точном интегрировании получаем:
,
а при интегрировании по формуле «средних» прямоугольников
Как видно, результаты точного и численного интегрирования совпадают.
Формула трапеций.
Заменяем на i -ом участке интегрируемую функцию линейной функцией 
, принимающей в точках и 
значения (рис. 2.5.3):
, 
. (2.5.9)
Рис. 2.5.3. К интегрированию по формуле трапеций.
Тогда интеграл на отрезке заменятся площадью трапеции (
и 
– основания, – высота)
. (2.5.10)
Вычисление исходного интеграла сводится к вычислению суммы
, (2.5.11)
или
. (2.5.12)
Точность вычисления. Как следует из построения квадратурная формула трапеций дает точный результат интегрирования для функций, линейных на i -ом участке ().
Формула Симпсона.
Разобьем интервал интегрирования на четное число отрезков. Рассмотрим сдвоенный участок [ 
].
Построим параболу
,
принимающую в точках
значения (см. рис. 3.5.4)
, 
, 
.
Рис. 2.5.4. К интегрированию по формуле Симпсона.
Такая парабола может быть представлена формулой
, (2.5.13)
где
, при этом 
(2.5.14)
или 
. (2.5.15)
Делая замену переменных, вычислим приближенное значение интеграла
Параметры 
, 
определим из условий
В итоге получим квадратурную формулу Симпсона
. (2.5.16)
Общая формула для вычисления приближенного значения интеграла примет вид
. (2.5.17)
Суммирование ведется только по нечетным i. Если перегруппировать члены суммы, получим
т.е.
(3.5.18)
Точность вычисления. Как следует из построения квадратурные формулы Симпсона дают точный результат интегрирования для функций, имеющих вид квадратичной параболы на сдвоенном участке [ 
] (
– нечетное). При этом одинаковая длина сдвоенных участков вовсе не обязательна и использовалась здесь исключительно для упрощения промежуточных выкладок и вида результирующих формул (2.5.17)–(2.5.18).
В том случае, когда сдвоенные участки имеют одинаковую длину, т.е.
– средняя точка сдвоенного участка,
формула Симпсона точна для функций, имеющих вид кубической параболы на этих участках. Это утверждение достаточно проверить для 
.
При точном интегрировании имеем:
,
а при интегрирование по формуле Симпсона получаем:
Как видно, результаты точного и численного интегрирования совпадают.
Вычисление интеграла с заданной точностью.
Обозначим через 
вычисленное приближенное значение интеграла
при разбиении отрезка на частей. При необходимости вычислить результат с заданной точностью 
вычисления повторяют с последовательно уменьшающимся (вдвое) шагом до тех пор, пока не выполнится условие окончание счета
. (2.5.19)