4. Составные формулы прямоугольников, трапеций, Симпсона.

На практике, при вычислении интеграла (16) обычно требуется, чтобы применение квадратурной формулы давало приемлемую точность, (как правило, известную заранее). В оценках погрешности (22), (24), (26) присутствует как длина отрезка - так и точность, с которой функция может быть приближена многочленом степени . В частности, если длина отрезка интегрирования не очень мала, то соответствующие оценки становятся малосодержательны. Как повысить точность?

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

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

может расти с ростом степени многочлена.

68

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

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

. (27)

Ее погрешность, согласно (22) оценивается следующим образом:

. (28)

b. Составная формула трапеций.

, (29)

для погрешности, которой справедлива оценка:

. (30)

с. Для получения составной формулы Симпсона, отрезок разбивают на равных частей и применяют формулу Симпсона

на отрезке длины . В результате приходим к формуле

, (31)

где

. (32)

69

5. Квадратурные формулы Гаусса.

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

. (33)

Пусть формула (33) точна для . Тогда, подставляя в (33) поочередно указанные функции, получаем уравнение с неизвестными. Если , то такие формулы называются квадратурными формулами Гаусса. Для простоты выкладок обычно рассматривается базовый отрезок . Решать систему для нахождения узлов и коэффициентов при малых можно напрямую, а при увеличении - довольно сложно.

Пример. .

.

Тогда

Откуда , и получаем формулу центральных прямоугольников.

Для решения поставленной задачи в общем случае нам понадобятся полиномы Лежандра, которые определяются формулами

. (34)

В частности

. (35)

70

Для полиномов Лежандра справедливы следующие свойства:

1. Полином Лежандра является полиномом -й степени, обладающим той же четностью, что и :

. (36)

2. Полином Лежандра имеет на интервале (-1,1) простых корней, которые располагаются симметрично относительно точки .

3. Любой полином степени ортогонален к полиному Лежандра :

. (37)

На доказательстве этих свойств мы останавливаться не будем.

Перейдем к построению квадратурных формул Гаусса. Возьмем в качестве узлов - корни полинома Лежандра и построим по ним квадратурную формулу

, (38)

которая точна для любого полинома степени . Так как любой многочлен степени можно представить в следующем виде ,

где

,

то, применяя квадратурную формулу (38), имеем:

(39)

Следовательно, квадратурная формула Гаусса (38) точна для любого многочлена степени .

71

Можно получить формулы для нахождения весовых коэффициентов , однако на практике, при небольших , коэффициенты вычисляются решая систему уравнений :

. (36)

Пример.

.

Тогда

.

Откуда .

Заметим, что сумма коэффициентов всегда равна двум.

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

.

Можно показать, что для погрешности формулы Гаусса справедлива оценка

. (37)

Доказательство формулы (37) выходит за рамки нашего пособия.

Замечание.

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

72

Упражнения для самостоятельного решения

1. Методом неопределенных коэффициентов построить формулу численного дифференцирования вида

.

Для какой точки формула имеет наибольший порядок аппроксимации?

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

.

3. Вычислить приближенно , используя формулы интерполирования, если .

4. С какой погрешностью вычисляется по обобщенным формулам прямоугольников и трапеций с шагом .

5. С помощью формулы Гаусса с 3 узлами вычислить . Найти погрешность ( ) .

6. С каким шагом необходимо проводить вычисления по составной формуле Симпсона, чтобы погрешность при вычислении интеграла

не превышала ?

73