Разделы сайта

Вывод формулы Симпсона

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

Рассмотрим подынтегральную функцию на отрезке . Заменим эту подынтегральную функцию интерполяционным многочленом Лагранжа второй степени, совпадающим с в точках x0, x1 = x0 + h, x2 = x0 + 2h:

Проинтегрируем :

Формула:

называется формулой Симпсона.

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

Оценим теперь погрешность интегрирования по формуле Симпсона. Будем считать, что у на отрезке существуют непрерывные производные f ʹ, f ʹʹ, f ʹʹʹ, f ʹʹʹʹ.

Составим разность:

К каждому из этих двух интегралов уже можно применить теорему о среднем, поскольку непрерывна на и функция неотрицательна на первом интервале интегрирования и неположительна на втором ( то есть не меняет знака на каждом из этих интервалов). Поэтому:

(воспользовались теоремой о среднем, поскольку - непрерывная функция; ).

Дифференцируя дважды и применяя затем теорему о среднем, получим для другое выражение:

, где

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

, .

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

Формула Симпсона в общем виде:

Перейти на страницу: 1 2

Самое читаемое:

Автоматизация технологических процессов в условиях технологического комплекса КК-АДСК-МНЛЗ ПАО МК Азовсталь, г. Мариуполь
Автоматизация управления технологическими процессами и производством играет важную роль в непрерывном увеличении выпуска продукции в нашей стране. Интенсификация и усложнение технологических процессов, рост единичной мощности агрегатов и повышение требований к качеству готовой продукции в соответствии с международными стандартами де ...

www.techstages.ru : Все права защищены! 2025