Оценки ошибок формул трапеций и центральных прямоугольников

Будем предполагать, что функция f(x) имеет на отрезке интегрирования .вторую производную f//(x), и f//(x)непрерывна на., причёмпри всех x..

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

По формуле Тейлора, применённой к функции f(x) в точке , получаем для x.:

где .— некоторая точка, лежащая междуи.

Заметим, что

поскольку — середина отрезка интегрирования в этом интеграле. Получаем тогда, что

где h_i=x_i-x_i-1. Таким образом, суммируя по всем отрезкам разбиения, получаем оценку ошибки:

Если все отрезки разбиения имеют одинаковую длину h_i=, то получаем

или

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

Рассмотрим снова рис. 5. Прямая, соединяющая концы хорды графика, то есть точки (x_i-1;f(x_i-1)) и (x_i;f(x_i)), имеет уравнение

Действительно, это равенство задаёт линейную функцию, и легко проверить, что l_i(x_i-1)=f(x_i-1) и l_i(x_i)=f(x_i). Разность между площадью под графиком функции на отрезке ,и площадью трапецииS_i равняется тогда

Докажем, что стоящая под знаком последнего интеграла разность удовлетворяет оценке

при всех ,. Эта оценка получается как следствие такой теоремы.

Теорема (о погрешности линейной интерполяции) Пусть f(x) — функция, имеющая на отрезке непрерывную вторую производную f^//(x), а l(x) — линейная функция, такая что ;. Назовём функциюl(x) линейной интерполирующей функцией для f(x) на , а разность погрешностью линейной интерполяции.

Тогда найдётся такая точка , что

(1)

Доказательство. Очевидно, что при иравенство (1) выполняется, как бы ни была выбрана точка, поскольку и левая, и правая части равенства обращаются тогда в ноль. Пусть теперь точкане совпадает ни с, ни с. Рассмотрим вспомогательную функцию

зависящую от параметра , и выберем значениетак, чтобы было выполнено равенство. Легко видеть, что для этого нужно взятьТогда

Функция (x) обращается в 0 в трёх точках: и. Значит, по теореме Ролля, её производная^/(x) обращается в 0 в каких-то двух точках и. Применяя снова теорему Ролля, теперь уже к производной^/(x) на отрезке , получаем, что^//(x) обращается в 0 в некоторой точке . Однако функцию ^//(x) легко вычислить:

так как вторая производная линейной функции равна 0. Таким образом,,

Остаётся заметить, что точка выбиралась как произвольная точка интервала, и доказательство завершено.

Возвращаемся к изучению ошибки формулы трапеций и связанным с этим обозначениям.

Следствие 1. Приимеет место оценка

Доказательство. Применим формулу (1) к функции g=f и отрезку и получим:

Здесь мы заметили, что xi-1 и xi — корни квадратного трёхчлена (x-xi-1)(x-xi) так, что на отрезке между корнями квадратный трёхчлен не больше нуля, и поэтому

Ошибку на -м отрезке разбиения мы можем теперь оценить так:

(последний интеграл легко вычисляется).

Просуммируем теперь оценки ошибок на каждом отрезке разбиения и получим оценку ошибки всей квадратурной формулы трапеций:

Если длины всех отрезков разбиения взяты одинаковыми, равными , то полученная оценка даёт

Оценки ошибок и, как мы видим, отличаются ровно в два раза. Выше мы отмечали, что эти ошибки имеют противоположные знаки, если функция f//(x)сохраняет знак на отрезке интегрирования. Значит, на каждом отрезке знакопостоянства функции f//(x)ошибкиибудут примерно компенсировать друг друга. Подобно тому, как мы получили формулу трапеций из формул левых и правых прямоугольников, попробуем получить усреднённую квадратурную формулу, скомбинировав формулы центральных прямоугольников и трапеций; при этом нас подогревает надежда на то, что новая формула будет иметь существенно меньшую ошибку. Умножая IR на, для того чтобы уравновесить ошибки противоположных знаков, получаем:

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

Суммируя эти величины по всем отрезкам разбиения, получаем квадратурную формулу:

(2)

Эта формула в точности совпадает с формулой Симпсона.