Математика / Начала анализа

Адаптивный Симпсон: локальная оценка ошибки

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

Опубликовано: 26 октября 2026 г.Обновлено: 26 октября 2026 г.

Формула

S_1=S(a,b),\ S_2=S(a,m)+S(m,b),\quad \varepsilon_{\mathrm{est}}=\frac{|S_2-S_1|}{15},\quad I\approx S_2+\frac{S_2-S_1}{15}

Схема Симпсон для x⁴: изменение площади

Кривые: функция x⁴ и парабола одного правила. Узлы 0;1/4;1/2;3/4;1 используются для уточнения.

Обозначения

$S(a,b)$: (b−a)[f(a)+4f((a+b)/2)+f(b)]/6, единица интеграла
$m$: середина [a,b], единица x
$S₁,S₂$: одна оценка и сумма двух оценок Симпсона, единица интеграла
$εest$: асимптотическая оценка ошибки уточнённой суммы, единица интеграла

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

Функция вычислима на концах, серединах и четвертях каждого участка. Оценка с делителем 15 предполагает гладкость и ведущую ошибку четвёртого порядка при делении шага. Для рекурсии заранее задают допуск, предельную глубину и минимально различимый шаг.

Ограничения

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

Подробное объяснение

Одно правило Симпсона использует три значения функции. Два правила на половинах используют пять, из которых три уже известны. Разность S₂−S₁ отражает, насколько локальная парабола перестала описывать интеграл после уточнения.

Если ведущая ошибка меняется как четвёртая степень шага, ошибка S₁ примерно в 16 раз больше ошибки S₂. Разность двух оценок содержит 15 таких малых ошибок, откуда берётся делитель 15. Добавление (S₂−S₁)/15 компенсирует ведущий член с его знаком.

Если индикатор меньше локального допуска, участок принимают. Иначе его делят и передают каждой половине часть допуска, например половину. Тогда суммарный бюджет принятых участков не превышает исходный бюджет при справедливости локальных оценок.

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

Как пользоваться формулой

На исходном участке вычислите одно правило Симпсона, затем сумму двух правил на половинах с сохранением общих узлов.
Найдите модуль разности и разделите его на 15; сравнивайте эту оценку с выделенным допуском именно участка.
При успехе примите уточнённую сумму с поправкой её знака, иначе разделите участок и бюджет допуска между половинами.
Контролируйте глубину рекурсии и различимость узлов; при исчерпании ресурса сохраните результат как недостигший допуска.

Историческая справка

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

В вычислительной математике XX–XXI веков этот подход стал способом экономить расчёты на гладких частях, оставляя больше точек у сложных участков. Fundamentals of Numerical Computation подробно выводит идею вложенных оценок и распределения допуска; NIST DLMF задаёт свойства исходной квадратуры. Современная рекурсивная реализация добавляет контроль глубины и переиспользование узлов. Адаптивный алгоритм нельзя целиком приписывать имени в названии локального правила: размещение узлов и управление погрешностью составляют отдельную часть метода.

Историческая линия формулы

Имя Симпсона относится к базовому параболическому правилу. Адаптивная схема добавляет сравнение вложенных оценок, распределение допуска и рекурсивное деление. Делитель 15 следует из четвёртого порядка ошибки и удвоения разрешения, а не из нового эмпирического закона.

Пример

Для f(x)=x⁴ на [0;1] одно правило даёт S₁=(0+4·(1/2)⁴+1)/6=5/24≈0,208333333. Сумма на двух половинах S₂=77/384≈0,200520833. Разность S₂−S₁=−1/128, её модуль после деления на 15 равен 1/1920≈0,000520833. При локальном допуске 0,001 участок принимается; при 0,0001 требуется дальнейшее деление. Скорректированная оценка S₂+(S₂−S₁)/15=0,2 совпадает с точным интегралом. Это удачное полиномиальное совпадение не превращает индикатор в строгую универсальную границу для других функций. На следующем шаге узлы переиспользуются.

Частая ошибка

Модуль нужен для сравнения с допуском, но поправка к интегралу сохраняет знак S₂−S₁. Нельзя дать каждой половине полный исходный допуск без учёта общей суммы. Одно правило Симпсона ещё не адаптивный алгоритм. Отсутствие новых делений из-за машинного шага не означает математической сходимости.

Практика

Задачи с решением

Проверка допуска

Условие. S₁=2,015, S₂=2,000, локальный допуск 0,002. Принимается ли участок?

Решение. |S₂−S₁|/15=0,015/15=0,001<0,002; скорректированный ответ 2−0,001=1,999.

Ответ. Да, по индикатору; оценка 1,999.

Разделение бюджета

Условие. Общий допуск 0,0008 разделён поровну между двумя половинами. Какой допуск каждой?

Решение. 0,0008/2=0,0004; сумма локальных бюджетов равна исходному.

Ответ. По 0,0004.

Дополнительные источники

https://fncbook.com/adaptive/
https://dlmf.nist.gov/3.5

Связанные формулы

Математика

Формула Симпсона: интеграл по параболическим дугам

$S_n=\frac h3\left[f_0+f_n+4\!\sum_{\substack{i=1\\i\text{ нечётное}}}^{n-1}f_i+2\!\sum_{\substack{i=2\\i\text{ чётное}}}^{n-2}f_i\right]$

Составная формула Симпсона объединяет соседние промежутки попарно и заменяет функцию параболой через три узла. Для равномерной сетки число промежутков должно быть чётным.

Математика

Метод Ромберга: таблица уточнений интеграла

$R_{k,0}=T_{2^k},\quad R_{k,j}=R_{k,j-1}+\frac{R_{k,j-1}-R_{k-1,j-1}}{4^j-1},\quad 1\le j\le k$

Метод Ромберга строит таблицу из составных трапеций с последовательным удвоением числа отрезков. Экстраполяция устраняет ведущие чётные степени шага и повышает точность без смены уже вычисленных узлов.

Математика

Составная формула средних прямоугольников

$M_n=h\sum_{i=0}^{n-1}f\!\left(a+\left(i+\tfrac12\right)h\right),\quad h=\frac{b-a}{n}$

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