Метод Симпсона

Метод Симпсона (также известный как формула Симпсона или параболический метод) — это численный метод вычисления определённого интеграла, основанный на аппроксимации подынтегральной функции параболой на каждом отрезке разбиения. Относится к классу квадратурных формул Ньютона — Котеса и обеспечивает более высокую точность по сравнению с методами прямоугольников и трапеций при одинаковом количестве точек разбиения, особенно для гладких функций.

История

Метод был назван в честь британского математика Томаса Симпсона (1710—1761), который опубликовал его в 1743 году в своей работе «Mathematical Dissertations». Однако исторически формула была известна ещё раньше: её частный случай (для трёх точек) использовал немецкий математик и астроном Иоганн Кеплер в 1615 году для вычисления объёмов бочек. В честь Кеплера этот частный случай иногда называют «формулой Кеплера» или «правилом бочек». В русскоязычной литературе метод Симпсона часто называют «формулой парабол», подчёркивая геометрическую интерпретацию.

Математическая формулировка

Пусть требуется вычислить определённый интеграл от функции \( f(x) \) на отрезке \([a, b]\):

\[ I = \int_a^b f(x) \, dx. \]

Отрезок \([a, b]\) разбивается на чётное число \( n \) равных частей (шагов) длиной \( h = \frac{b-a}{n} \). Точки разбиения: \( x_i = a + i h \), где \( i = 0, 1, \dots, n \). Значения функции в этих точках: \( y_i = f(x_i) \).

Основная идея метода заключается в том, что на каждом двойном отрезке \([x_{2k}, x_{2k+2}]\) (где \( k = 0, 1, \dots, \frac{n}{2} - 1 \)) подынтегральная функция заменяется интерполяционным многочленом второй степени (параболой), проходящим через три точки: \((x_{2k}, y_{2k})\), \((x_{2k+1}, y_{2k+1})\) и \((x_{2k+2}, y_{2k+2})\).

Формула Симпсона

Общая формула для приближённого значения интеграла имеет вид:

\[ I \approx \frac{h}{3} \left[ y_0 + y_n + 4 \sum_{i=1,3,5,\dots}^{n-1} y_i + 2 \sum_{i=2,4,6,\dots}^{n-2} y_i \right]. \]

Или, в развёрнутом виде:

\[ I \approx \frac{h}{3} \left( y_0 + 4y_1 + 2y_2 + 4y_3 + 2y_4 + \dots + 4y_{n-1} + y_n \right). \]

Здесь:

• \( y_0 \) и \( y_n \) — значения функции на концах отрезка;
• коэффициенты 4 — для точек с нечётными индексами (внутренние точки каждого двойного отрезка);
• коэффициенты 2 — для точек с чётными индексами, кроме концов.

Погрешность метода

Остаточный член (погрешность) формулы Симпсона для одного двойного отрезка длины \( 2h \) выражается как:

\[ R = -\frac{h^5}{90} f^{(4)}(\xi), \]

где \( \xi \) — некоторая точка внутри отрезка \([x_{2k}, x_{2k+2}]\), а \( f^{(4)} \) — четвёртая производная функции.

Для всего отрезка \([a, b]\) глобальная погрешность имеет порядок \( O(h^4) \). Это означает, что при уменьшении шага \( h \) в два раза погрешность уменьшается примерно в 16 раз. Метод Симпсона является точным для многочленов степени не выше третьей (включая кубические параболы), так как четвёртая производная таких многочленов равна нулю.

Условия применимости

Метод Симпсона требует, чтобы подынтегральная функция была достаточно гладкой — имела непрерывную четвёртую производную на всём отрезке интегрирования. Для функций с разрывами, резкими осцилляциями или особенностями (например, вблизи полюсов) точность метода может резко снижаться, и требуется применение адаптивных алгоритмов или других методов (например, квадратур Гаусса).

Геометрическая интерпретация

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

Алгоритм реализации

На практике вычисление по методу Симпсона выполняется в несколько шагов:

1. Выбрать чётное число разбиений \( n \) (обычно \( n \geq 2 \)).
2. Вычислить шаг \( h = \frac{b-a}{n} \).
3. Определить значения функции \( y_i \) во всех точках разбиения.
4. Вычислить сумму:

• \( S_0 = y_0 + y_n \) — концы отрезка;
• \( S_1 = \sum_{i=1,3,5,\dots}^{n-1} y_i \) — нечётные индексы;
• \( S_2 = \sum_{i=2,4,6,\dots}^{n-2} y_i \) — чётные индексы (кроме концов).

1. Итоговое значение: \( I \approx \frac{h}{3} (S_0 + 4S_1 + 2S_2) \).

Пример вычисления

Вычислим интеграл \( \int_0^2 x^2 \, dx \) с \( n = 4 \) (шаг \( h = 0.5 \)):

• Точки: \( x_0 = 0, y_0 = 0 \); \( x_1 = 0.5, y_1 = 0.25 \); \( x_2 = 1, y_2 = 1 \); \( x_3 = 1.5, y_3 = 2.25 \); \( x_4 = 2, y_4 = 4 \).
• \( S_0 = 0 + 4 = 4 \); \( S_1 = 0.25 + 2.25 = 2.5 \); \( S_2 = 1 \).
• \( I \approx \frac{0.5}{3} (4 + 4 \cdot 2.5 + 2 \cdot 1) = \frac{0.5}{3} (4 + 10 + 2) = \frac{0.5}{3} \cdot 16 = \frac{8}{3} \approx 2.6667 \).

Точное значение: \( \int_0^2 x^2 \, dx = \frac{8}{3} \approx 2.6667 \). Погрешность равна нулю, так как подынтегральная функция — многочлен второй степени, а метод Симпсона точен для многочленов до третьей степени.

Разновидности и обобщения

Составная формула Симпсона

Стандартная реализация, описанная выше, называется составной (или композитной) формулой Симпсона. Она применяется, когда число разбиений \( n \) чётно.

Правило 3/8 Симпсона

Существует вариант метода, называемый «правилом 3/8», который использует аппроксимацию кубическим многочленом на отрезке, разбитом на три равные части. Формула имеет вид:

\[ I \approx \frac{3h}{8} (y_0 + 3y_1 + 3y_2 + y_3). \]

Этот вариант точен для многочленов степени не выше третьей, но требует, чтобы число разбиений \( n \) было кратно трём. На практике правило 3/8 часто комбинируют с обычной формулой Симпсона для обработки отрезков с нечётным числом интервалов.

Адаптивный метод Симпсона

Для повышения эффективности используются адаптивные алгоритмы, которые автоматически изменяют шаг разбиения в зависимости от локальной гладкости функции. Например, на участках с быстрым изменением функции шаг уменьшается, а на гладких — увеличивается. Это позволяет достичь заданной точности при минимальном количестве вычислений.

Применение

Метод Симпсона широко применяется в вычислительной математике, физике, инженерии и экономике для приближённого вычисления интегралов, когда аналитическое интегрирование затруднено или невозможно. Он используется:

• в расчётах площадей и объёмов (например, при определении водоизмещения судов по ватерлинии);
• в обработке сигналов (спектральный анализ);
• в статистике (вычисление моментов распределений);
• в численных методах решения дифференциальных уравнений (например, метод Рунге — Кутты использует квадратурные формулы).

Сравнение с другими методами

Метод Симпсона при одинаковом числе точек разбиения даёт существенно меньшую погрешность для гладких функций, но требует большего объёма вычислений на каждом шаге. Для функций с высокой осцилляцией или разрывами более эффективными могут оказаться методы Гаусса или адаптивные квадратуры.

Интересные факты

• Формула Симпсона является частным случаем квадратурных формул Ньютона — Котеса с \( n = 2 \) (три точки).
• В некоторых источниках метод Симпсона называют «правилом парабол», а его вариант с четырьмя точками — «правилом 3/8».
• В XIX веке русский математик Пафнутий Чебышёв разработал квадратурные формулы, оптимизированные по точности для заданного числа узлов, которые в ряде случаев превосходят метод Симпсона.

Источники

• Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином. Лаборатория знаний, 2006.
• Калиткин Н. Н. Численные методы. — М.: Наука, 1978.
• Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
• Формальский А. М. Численные методы. — М.: Изд-во МГУ, 2004.