Формула Симпсона (формула парабол).

Численное интегрирование

Задача численного интегрирования состоит в том, чтобы найти численное значение определенного интеграла

, (1)

где - функция, непрерывная на отрезке интегрирования . Формулы для решения этой задачи называются квадратурными. Квадратурная формула позволяет вместо точного значения интеграла (1) найти некоторое его приближенное значение . Разность точного и приближенного значений интеграла называется абсолютной погрешностью квадратурной формулы (или численного метода),

.

Квадратурные формулы используют для вычисления интеграла (1) значения , , …, функции в точках отрезка . Квадратурная формула имеет вид

, (2)

где - некоторые коэффициенты, которые называют весовыми.

Напомним геометрический смысл определенного интеграла: выражает площадь соответствующей криволинейной трапеции (фигуры, ограниченной графиком функции , прямыми и осью ).

Рассмотрим два подхода к решению задачи численного интегрирования.

1) Разобьем отрезок интегрирования на частичных отрезков, вычислим интегралы на частичных отрезках. Интеграл на всем отрезке интегрирования равен сумме интегралов на частичных отрезках (свойство аддитивности определенного интеграла).

2) Вычислим , заменяя подынтегральную функцию на всем отрезке интегрирования интерполяционным полиномом Лагранжа , построенным на узлах .

Обозначим через максимальное по модулю значение производной - го порядка функции на отрезке :

Рассмотрим варианты решения данной задачи.

Формулы прямоугольников

Для простоты разобьем отрезок интегрирования на частей точками, равноудаленными друг от друга: так, что будет выполняться равенство: , где .

Аппроксимируем площадь под графиком функции суммой площадей прямоугольников с основанием и высотой , где .

Причем, если взять (левую крайнюю точку частичного отрезка), то получим формулу левых прямоугольников:

.

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

Абсолютная погрешность формулы левых прямоугольников определяется выражением:

,

где .

Мы видим, что погрешность метода левых прямоугольников имеет тот же порядок, что шаг интегрирования (первый порядок по ). Поскольку для функций вида , то для таких функций формула левых прямоугольников является точной.

А если взять (правую крайнюю точку частичного отрезка), то получим формулу правых прямоугольников:

.

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

В случае, когда мы берем среднюю точку , получаем формулу средних прямоугольников:

, .

Абсолютная погрешность формулы средних прямоугольников оценивается выражением:

,

Погрешность формулы средних прямоугольников имеет второй порядок по ().

Наиболее употребительной является формула средних прямоугольников.

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

Заменим площадь криволинейной трапеции суммой площадей прямолинейных трапеций, построенных на частичных отрезках , . Площадь элементарной прямолинейной трапеции равна:

,

а интеграл равен:

.

Оценка абсолютной погрешности на всем отрезке интегрирования определяется выражением:

Погрешность метода трапеций имеет тот же порядок, что и . Для функций вида (полиномов первой степени) формула трапеций является точной.

Формула Симпсона (формула парабол).

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

Абсолютная погрешность формулы Симпсона оценивается выражением:

Абсолютная погрешность формулы Симпсона имеет тот же порядок, что и (четвертый порядок точности). Формула Симпсона точна для полиномов степени .

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

Правило Рунге основано на соотношении:

, (3)

где , - приближенные значения определенного интеграла, вычисленные при разбиении отрезка интегрирования на n и 2n частей соответственно; - порядок метода; - заданная точность. При каждом последующем приближении число отрезков разбиения удваивается. Если условие выполнено, за приближенное значение интеграла принимается значение , т.е. .Так как оценка осуществляется после вычисления, то она является апостериорной.

Напомним порядки методов (по ):