3. Узлы и весовые коэффициенты квадратурных формул Гаусса.

Изучив свойства полиномов Лежандра, перейдем к решению основной задачи – определению узлов и весовых коэффициентов квадратурных формул Гаусса. Составим полином -ой степени

,

где - искомые узлы. Возьмем произвольный полиномстепени, помножим его на полиноми проинтегрируем произведение по отрезкус помощью квадратурной формулы . Поскольку это произведение представляет собой полином степени, формула Гаусса должна быть для него точной. В результате согласно получим:

.

Мы видим, что полином ортогонален к любому полиному степенив том числе и к полиномам Лежандра индекса. Это означает, что он с точностью до множителя совпадает с-ым полиномом Лежандра:. Отсюда следует вывод: узлы квадратурной формулы Гаусса являются корнями полинома Лежандра. Напомним, что корни полиномов Лежандра располагаются на интервалесимметрично относительно его средней точки.

Для того, чтобы подсчитать весовые коэффициенты , введем специальные полиномы

.

Каждый из них является полиномом степени . В числителе у него стоит полиномс опущенным множителем, в знаменателе - значение числителя в точке. В результате такой структуры полиномв точкахудовлетворяет соотношениям:

.

Для полинома квадратурная формула Гаусса должна быть точной. С учетом это дает

.

В результате получаем следующее интегральное выражение для весовых коэффициентов квадратурной формулы Гаусса:

.

4. Исследование квадратурной формулы.

Нам осталось решить последний вопрос – доказать, что квадратурная формула, у которой в качестве узлов берутся корни полинома Лежандра, а весовые коэффициентывычисляются по формулам , действительно решают задачу Гаусса, являясь точной для любого полинома степени.

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

.

Справедливость данного разложения вытекает из следующих соображений. Здесь левая и правая части равенства совпадают в точках. Но, если два полинома-ой степени совпадают вточках, то они тождественно равны.

Интегрируя равенство по отрезку , получим

.

Итак, для полиномов -ой степени утверждение доказано.

Теперь рассмотрим произвольный полином степени. Разделим его с остатком на полином Лежандраи представим в виде:

,

где иполиномы степени. Проинтегрировав равенство по отрезку, будем иметь:

Поясним выполненные преобразования. Интеграл опущен, поскольку полином Лежандраортогонален к любому полиному-ой степени. Оставшийся интеграл от полиномавычислен с помощью квадратурной формулы . Выше уже доказано, что для полиномов степениона является точной. Последний переход заключается в том, что в суммудобавлены слагаемые. Они не меняют значения суммы, поскольку все равны нулю: ведь узлами квадратурной формулы являются корни полинома Лежандра.

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

Задача 5.

Построить квадратурную формулу Гаусса с двумя и тремя узлами.

Выведем сначала квадратурную формулу с двумя узлами. Узлы определяются как корни второго полинома Лежандра, выражение для которого мы выписывали выше . В данном случае имеем:

,.

Узлы расположены симметрично относительно точки .

Весовые коэффициенты рассчитываются по формуле :

,

.

Они равны между собой, а их сумма, в соответствии с общим соотношением , равна двум. В результате искомая квадратурная формула принимает вид:

.

Она является точной для любого полинома третьей степени.

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

,,.

Остается подсчитать весовые коэффициенты:

В результате квадратурная формула Гаусса с тремя узлами запишется в виде:

.

Она является точной для любого полинома пятой степени.

Задача 6.

Вычислить по формулам Симпсона и Гаусса при интеграл:

.

Сравнить результаты численного интегрирования с точным значением интеграла и между собой.

Формулы Симпсона и Гаусса дают в данном случае следующие результаты:

,

,

,

.

Мы видим, что даже с двумя узлами формула Гаусса дает хороший ответ. Его точность выше точности ответа, полученного по формуле Симпсона.

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