Метод Грегори

Метод Грегори — это алгоритм численного интегрирования, предназначенный для приближённого вычисления определённого интеграла. Относится к классу формул Ньютона — Котеса и является обобщением простейших квадратурных правил, таких как метод прямоугольников, трапеций и Симпсона. Метод основан на интерполяции подынтегральной функции многочленом, проходящим через заданное количество равноотстоящих узлов, и последующем интегрировании этого многочлена. Назван в честь шотландского математика Джеймса Грегори (1638—1675), который впервые описал подобный подход в своей работе «Vera Circuli et Hyperbolae Quadratura» (1667).

История

Разработка методов численного интегрирования была тесно связана с потребностями астрономии, навигации и физики XVII века, где требовалось вычислять площади сложных фигур и объёмы тел, не поддающиеся аналитическому интегрированию. Джеймс Грегори, работая над задачами квадратуры круга и гиперболы, предложил способ аппроксимации интеграла путём разложения подынтегральной функции в ряд и последующего почленного интегрирования. В своей книге «Vera Circuli et Hyperbolae Quadratura» он изложил принципы, которые позже легли в основу формул, известных как формулы Грегори.

Впоследствии метод был развит и систематизирован в работах других математиков, в частности Исаака Ньютона и Роджера Котеса, которые в начале XVIII века создали общую теорию квадратурных формул. Метод Грегори часто рассматривается как частный случай или как одна из первых реализаций формул Ньютона — Котеса. В XIX—XX веках, с развитием вычислительной техники, метод Грегори и его модификации стали широко применяться для инженерных расчётов, а затем и для компьютерных вычислений, где он используется как основа для построения адаптивных алгоритмов интегрирования.

Сущность метода

Метод Грегори основан на замене подынтегральной функции \(f(x)\) на отрезке \([a, b]\) интерполяционным многочленом, построенным по \(n+1\) равноотстоящим узлам: \(x_0 = a, x_1 = a + h, x_2 = a + 2h, \ldots, x_n = b\), где \(h = (b-a)/n\) — шаг интегрирования. Значение интеграла приближённо равно интегралу от этого многочлена.

Общая формула метода Грегори для \(n\) узлов (с порядком точности \(O(h^{n+1})\)) имеет вид:

\[ \int_a^b f(x) \, dx \approx h \sum_{k=0}^{n} C_k f(x_k), \]

где \(C_k\) — весовые коэффициенты, зависящие от \(n\) и определяемые из условия точного интегрирования многочленов степени \(n\). Эти коэффициенты могут быть найдены с помощью метода неопределённых коэффициентов или по формулам Котеса.

Частные случаи

В зависимости от количества узлов \(n\) метод Грегори порождает известные квадратурные правила:

• n = 1 (два узла, отрезок \([x_0, x_1]\)): формула трапеций.

\[ \int_{x_0}^{x_1} f(x) \, dx \approx \frac{h}{2} \left( f(x_0) + f(x_1) \right). \] Погрешность: \(O(h^3)\).

• n = 2 (три узла, отрезок \([x_0, x_2]\)): формула Симпсона (парабол).

\[ \int_{x_0}^{x_2} f(x) \, dx \approx \frac{h}{3} \left( f(x_0) + 4f(x_1) + f(x_2) \right). \] Погрешность: \(O(h^5)\).

• n = 3 (четыре узла, отрезок \([x_0, x_3]\)): правило трёх восьмых (правило Симпсона 3/8).

\[ \int_{x_0}^{x_3} f(x) \, dx \approx \frac{3h}{8} \left( f(x_0) + 3f(x_1) + 3f(x_2) + f(x_3) \right). \] Погрешность: \(O(h^5)\).

• n = 4 (пять узлов, отрезок \([x_0, x_4]\)): правило Милна — Боуза (формула Котеса пятого порядка).

\[ \int_{x_0}^{x_4} f(x) \, dx \approx \frac{2h}{45} \left( 7f(x_0) + 32f(x_1) + 12f(x_2) + 32f(x_3) + 7f(x_4) \right). \] Погрешность: \(O(h^7)\).

При \(n \geq 8\) весовые коэффициенты в формулах Ньютона — Котеса становятся отрицательными, что может приводить к численной неустойчивости для больших \(n\). Поэтому на практике метод Грегори редко используется с числом узлов более 7—8.

Алгоритм применения

На практике метод Грегори реализуется следующим образом:

1. Разбиение отрезка. Отрезок интегрирования \([a, b]\) разбивается на \(m\) подотрезков (панелей) одинаковой длины \(h = (b-a)/m\). Каждый подотрезок содержит \(n+1\) узел (обычно \(n\) — небольшое число, например 2 или 3).
2. Применение формулы. На каждом подотрезке вычисляется интеграл по соответствующей формуле Грегори (например, для \(n=2\) — по формуле Симпсона).
3. Суммирование. Полученные значения суммируются по всем подотрезкам.

Для повышения точности применяется метод двойного пересчёта (правило Рунге), при котором вычисления проводятся с шагом \(h\) и \(h/2\), а затем оценивается погрешность.

Свойства и особенности

• Точность. Метод Грегори имеет порядок точности \(O(h^{n+1})\) для гладких функций, что делает его эффективным при интегрировании многочленов степени не выше \(n\). Для функций с разрывами или особенностями точность может резко снижаться.
• Устойчивость. При малых \(n\) (до 7) метод устойчив к ошибкам округления. При \(n \geq 8\) отрицательные коэффициенты могут приводить к катастрофической потере точности.
• Простота реализации. Алгоритм легко программируется, так как требует только вычисления значений функции в равноотстоящих точках.
• Ограничения. Метод требует равномерной сетки, что не всегда удобно для функций с резкими изменениями на отдельных участках. В таких случаях предпочтительнее адаптивные методы (например, метод Гаусса — Кронрода).

Применение

Метод Грегори и его частные случаи (формула трапеций, Симпсона) широко используются в:

• Инженерных расчётах — для вычисления площадей, объёмов, моментов инерции, работы сил.
• Физике — при интегрировании уравнений движения, расчёте полей и потенциалов.
• Экономике и статистике — для численного нахождения определённых интегралов от эмпирических функций.
• Компьютерной графике — для расчёта площадей под кривыми, например, в алгоритмах рендеринга.
• Образовании — как классический пример численного метода в курсах вычислительной математики.

В современных вычислительных пакетах (MATLAB, Scilab, Python с библиотекой SciPy) метод Грегори не реализован как отдельная функция, но его частные случаи (например, scipy.integrate.simpson) доступны для использования.

Связь с другими методами

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

Критика и ограничения

Основным недостатком метода Грегори является его чувствительность к выбору числа узлов. При большом \(n\) (более 7) формула становится численно неустойчивой из-за осцилляций интерполяционного многочлена (феномен Рунге). Кроме того, метод не адаптируется к особенностям функции, что требует либо ручного выбора подходящего шага, либо использования составных (композитных) схем. Для функций с быстрыми изменениями или разрывами более эффективны методы с переменным шагом, такие как метод Ромберга или адаптивное интегрирование.

Источники

1. Грегори, Джеймс. «Vera Circuli et Hyperbolae Quadratura» (1667).
2. Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином. Лаборатория знаний, 2006.
3. Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
4. Формалев В. Ф., Ревизников Д. Л. Численные методы. — М.: Физматлит, 2004.
5. Abramowitz M., Stegun I. A. Handbook of Mathematical Functions. — Dover Publications, 1965.