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

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

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

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

Формула

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

Схема Ромберг для x⁴: таблица уточнений

Слева трапеции, затем устранение h² и h⁴. Последняя оценка 0,2 совпадает с интегралом.

Обозначения

$Tk$: составная трапеция с k ячейками, единица интеграла
$Rk,j$: элемент таблицы экстраполяции, единица интеграла
$k$: номер уточнения сетки, начиная с нуля, 1
$j$: столбец экстраполяции, 1

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

Начальный столбец состоит из трапеций на одной и той же паре границ с шагами h,h/2,h/4. Для последовательного устранения членов ошибки нужна соответствующая гладкость функции и чётностепенное разложение остатка. Индекс j не превышает номер строки k.

Ограничения

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

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

Трапеции для гладкой функции имеют ошибку c₂h²+c₄h⁴+… . Два значения с шагами h и h/2 позволяют исключить c₂: (4T(h/2)−T(h))/3. Именно это выполняет первый столбец экстраполяции после начального.

Следующий столбец сравнивает уже уточнённые оценки с ведущей ошибкой четвёртого порядка. При делении шага пополам такой член меняется в 16 раз, поэтому знаменатель становится 16−1. Общий знаменатель 4ʲ−1 последовательно соответствует степеням 2j.

Новые трапеции можно вычислять экономно: половина предыдущей суммы сохраняется, добавляются только значения в серединах прежних ячеек. После строки k известны 2ᵏ+1 узлов; экстраполяционные операции сами новых обращений к функции не требуют.

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

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

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

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

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

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

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

Название Ромберга относится к табличной схеме интегрирования с последовательной экстраполяцией составных трапеций. Её механизм — устранение степенных членов ошибки. Запись в NIST DLMF отделяет начальные квадратуры от последующих алгебраических уточнений.

Пример

Интегрируем x⁴ на [0;1]. Трапеции дают T₁=0,5, T₂=0,28125 и T₄=0,220703125. Первый уточнённый элемент R₁,₁=(4·0,28125−0,5)/3=0,208333333. В третьей строке R₂,₁=(4·0,220703125−0,28125)/3=0,200520833. Затем R₂,₂=R₂,₁+(R₂,₁−R₁,₁)/15=0,2. Точный интеграл также 1/5=0,2. Проверка: на первой экстраполяции убран член h², на второй — h⁴; для этой полиномиальной задачи остаток исчезает. Число узлов исходного столбца равно пяти; новые значения для алгебраического заполнения двух столбцов не требовались. Все члены таблицы относятся к одним границам.

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

Начальный столбец должен состоять из трапеций с отношением шагов ровно два. Знаменатель 4ʲ−1 меняется по столбцам, а не одинаков для всей таблицы. Нельзя смешивать разные границы интегрирования. Добавление строк после потери значащих цифр не обязательно улучшает результат.

Практика

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

Первое уточнение

Условие. T(h)=1,4 и T(h/2)=1,1. Найдите первый экстраполированный результат.

Решение. (4·1,1−1,4)/3=3/3=1.

Ответ. 1.

Второй столбец

Условие. R₁,₁=0,208333333333 и R₂,₁=0,200520833333. Найдите R₂,₂.

Решение. Разность −0,0078125 делится на 15; поправка −0,000520833333.

Ответ. 0,2 с точностью округления входа.

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

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

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

Математика

Составная формула трапеций для вычисления интеграла

$T_n=h\left[\frac{f(a)+f(b)}2+\sum_{i=1}^{n-1}f(a+ih)\right],\quad h=\frac{b-a}{n}$

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

Математика

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

$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]$

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

Математика

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

$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}$

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