Формула Симпсона (формула парабол)
Для вычисления интеграла разобьем отрезок интегрирования на два равных отрезка:
[х0 , х1] и [х1 , х2] (х0 = а, х2 = b)
и заменим подынтегральную функцию по формуле квадратичного интерполирования
(12)
где t = .
Тогда:
.
Перейдем к новой переменной интегрирования, учитывая, что
х = х0 +ht, dx= hdt,
при х=х0 t=0
при х=х2 t=2
или
(13)
Формула (13) называется формулой Симпсона или формулой парабол.
Ее геометрический смысл состоит в следующем: на отрезке [х0, х2] кривая у = f(x) заменяется квадратной параболой — графиком интерполяционного многочлена. При вычислении по формуле (13) значение интеграла будет численно равно значению площади криволинейной трапеции, ограниченной сверху дугой параболы, проходящей через точки: [х0, f(х0)], [х1, f(х1)], [х2, f (х2)]
На рисунке сплошной линией изображен график функции f(x) пунктирной — график многочлена Р2(х).
Для получения более точного результата достаточно разбить отрезок интегрирования [а; b] на четное число (2n) частей и применить формулу (13) для каждой пары смежных отрезков разбиения:
(14)
Суммируя равенства (14), получим обобщенную формулу Симпсона (парабол):
(15)
Пример. Вычислим приближенное значение интеграла по формуле Симпсона. Разбив отрезок интегрирования на десять равных частей и используя данные, содержащиеся в таблице, получим:
Итак, .
Выше показали, что .
Абсолютная погрешность найденного значения не превосходит 0,000005.
Сравнение приближенных значений интеграла , вычисленных по разным формулам, показывает, что наиболее точное значение было получено по обобщенной формуле Симпсона и наименее точное — по формуле прямоугольников.
Погрешность r обобщенной формулы Симпсона можно вычислить по формуле
(16)
где а < ξ< b.
Для абсолютной погрешности обобщенной формулы Симпсона можно получить следующую оценку:
где (17)
Дата добавления: 2022-07-20; просмотров: 74;