7. Метод гаусса

В предыдущих параграфах рассматривались методы чис­ленного интегрирования с произвольным разбиением интер­вала. Фактически разбиение производилось на равные отрезки.

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

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

Чтобы упростить выкладки, изменим пределы интегри­рования так, чтобы они стали равными (+1, -1). Для этого введем новую переменную

так что

Интеграл (1) после такой подстановки запишется в виде

(23)

где

Это означает, что соответствующей заменой переменных можно свести все интегралы к виду (23). (Под «всеми» подразумеваются интегралы с конечными пределами интег­рирования и с непрерывной подынтегральной функцией.)

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

для которой

(24)

Рис. 5. Геометрическое представление метода Гаусса с двумя

ординатами.

Интеграл, стоящий в левой части этого уравнения, пред­ставляет собой площадь трапеции на рис. 5. Пусть сумма площадей вертикально заштрихованных участков (между -1 и μ₀ и от μ₁ до +1) равна площади зачерненного участ­ка (между μ₀ и μ₁). Тогда площадь трапеции в точности равна площади под кривой y=φ(μ). Задача заключается в нахождении такой прямой линии, для которой достигается это равенство.

Положим

(25)

где необходимо определить А₀, A₁, μ₀ и μ₁. Так как в форму­ле имеются четыре параметра, то естественно предположить,

что формула даст точный результат при интегрировании кубической параболы

Перепишем эту функцию в виде

Если a₀ и a₁ должны удовлетворять уравнению (24), то μ₀ и μ₁ определяется из условия

Так как это равенство должно быть справедливо для любых β₀ и β₁ то необходимо потребовать выполнения двух следующих равенств:

Выполнив интегрирование, можно записать два послед­них равенства в виде

откуда следует, что

(26)

Теперь остается только найти A₀ и A₁ в формуле (25). Заметим, что

(27)

и из формул (25) и (26) следует, что

Так как значение выражения в правой части последней формулы должно быть равно значению интеграла в формуле (27) при всех a₀ и а₁, то для А₀ и A₁ получаем систему уравнений

из которой находим

А₀ = A₁ = 1. (28)

Уравнение (25) окончательно запишется в виде

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

где

(29)

Чтобы найти К, предположим, что При этом так что

(30)

Точное значение интеграла легко вычислить:

С другой стороны, из формул (29) и (30)

Поэтому K=1/135, и окончательная формула для ошибки ограничения запи­шется в следующем виде:

Можно вывести гауссовы формулы численного инте­грирования более высоких порядков, предусматривая большее количество ординат и вводя различные весовые коэффициенты А_i.

(31)

В общем случае, при (n + 1) ординате, получается точ­ная формула для нахождения интеграла от многочлена степени 2n + 1.

Оказывается, что значения μ_i в уравнении (31) являют­ся корнями полиномов Лежандра степени n. По этой причи­не вышеописанный метод численного интегрирования часто называют методом Лежандра — Гаусса. Полиномы Лежанд­ра P_n(μ) можно найти с помощью рекуррентных формул

(32)

Например, для m = 2

Заметим, что корни Р₂ (μ) равны , как мы уже опре­делили ранее.

Весовые коэффициенты в формуле (31) можно найти из следующего соотношения:

В качестве примера возьмем n = 2, так что μ ₀ = - , μ₁ = и Р’₂= 3μ. При этом

и совершенно аналогично A₁ = 1.

В общем случае ошибка ограничения определяется формулой

Таблица μ_i, и A_i для n = 2,...,6 приведена в приложе­нии. Заметим, что μ_i расположены симметрично относительно начала координат и что коэффициенты A_k одинаковы для μ_k и для — μ_k. Более подробная таблица вплоть до n = 48 приведена в приложении к книге В. И. Крылова «Прибли­женное вычисление интегралов».

Итак, метод Гаусса позволяет достичь большей точности, нежели формула Симпсона при том же количестве ординат, но зато точки, где следует определять эти ординаты, пол­ностью определены и совершенно не зависят от произвола программиста. Как и во многих других случаях, при чис­ленном интегрировании приходится делать выбор между простотой формулы Симпсона и возможной экономией машинного времени при использовании метода Гаусса. На практике чаще используется метод Симпсона.