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

Замена подынтегральной функции интерполяционным полиномом приводит к погрешности вычисления его значения R=|I1 – I|, где

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

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

· при использованииформул левых или правых прямоугольников

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

· при использованииформулы средних прямоугольников

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

· при использованииформулы трапеций

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

· при численном интегрировании по формуле Симпсона:

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

В приведенных выше формулах: a, bграницы интервала интегрирования; h=(b-a)/nшаг интегрирования; Оценка погрешности численного интегрирования - №6 - открытая онлайн библиотека . Оценка погрешности численного интегрирования - №7 - открытая онлайн библиотекаи Оценка погрешности численного интегрирования - №8 - открытая онлайн библиотека – среднее арифметическое, соответственно, первых, вторых и четвертых конечных разностей.

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

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

Формулы для оценки погрешности могут быть также использованы для выбора числа разбиений n или шага интегрирования h, необходимых для обеспечения заданной точности. Однако практическое использование этих формул ограничено в связи с трудоемкостью их вычислений, поэтому при реализации численных методов на ПК используется прием, позволяющий получить оценку погрешности в неявном виде. Этот прием основан на двукратном вычислении значения интеграла вначале с шагом h(где h=(b-a)/n), а затем с шагом h/2. Полученные значения интегралов IhиIh/2 могут быть применены для оценки погрешности интегрирования по формуле:

Оценка погрешности численного интегрирования - №9 - открытая онлайн библиотека (1.4.5-1)

где: k=1–для формул левых и правых прямоугольников;

k=2–для формул трапеции и средних прямоугольников;

k=4–для формулы Симпсона.

Если полученная погрешность не удовлетворяет требуемой точности, то вычисляется значение интеграла при h=h/4 и снова оценивается погрешность, и т.д. до тех пор, пока не окажется, что погрешность стала меньше заданной точности. Это правило называется правилом Рунге (или правилом двойного просчета).

Пример 1.4.5-1.Вычислить значение определенного интеграла Оценка погрешности численного интегрирования - №10 - открытая онлайн библиотека

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

x 0.0 0.1 0.2 0.3 0.4 0.5
f(x) 1.0 0.99005 0.960789 0.913913 0.852144 0.778801
0.6 0.7 0.8 0.9 1.0
0.697676 0.612626 0.527292 0.44858 0.367879

Используем формулы правыхилевых прямоугольников, считая, что h = 0.2, а n=(b-a)/h=5, имеем:

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

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

Используем формулу трапеций:

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

Используем формулы средних прямоугольников:

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

Используем формулу Симпсонаприm=2∙n=10(2∙5)и шаге h=0.1:

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

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

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

Следовательно, Оценка погрешности численного интегрирования - №18 - открытая онлайн библиотека

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