Квадратурная формула Чебышева

Формулы Ньютона-Котеса высших порядков

Производя соответствующие вычисления при n =3, получим из формулы (8.3) и из выражения значения и квадратурную формулу Ньютона:

(правило ).

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

Рассмотрим квадратурную формулу

, (8.10)

где - постоянные коэффициенты. Чебышев предположил выбрать абсциссы таким образом, чтобы:

1. коэффициенты были равны между собой;

2. квадратурная формула (8.10) являлась точной для всех полиномов до степени n включительно.

Покажем, как могут быть найдены в этом случае величины и . Полагаем . Учитывая, что при , будем иметь , получаем . Следовательно, квадратурная формула Чебышева имеет вид:

. (8.11)

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

, (8.12)

из которой могут быть определены неизвестные . Заметим, что система (8.12) при n =8 и n ³10 не имеет действительных решений.

Выведем формулу Чебышева с тремя ординатами (n =3).

Для определения абсцисс имеем систему уравнений:

(8.13)

Рассмотрим симметрические функции корней:

Из системы (8.13) имеем:

Отсюда заключаем по теореме Виета, что есть корни вспомогательного уравнения или . Следовательно, можно принять: .

Таким образом, соответствующая формула Чебышева имеет вид .

Чтобы применить квадратурную формулу Чебышева к интегралу вида , следует преобразовать его с помощью подстановки:

, переводящей отрезок в отрезок . Применяя к преобразованному интегралу формулу Чебышева, будем иметь

,

где и - корни системы (8.13).

В таблице приведены значения корней t_i системы (8.12) для n= 1,2…,7.

Таблица 8.1

Значения абсцисс t_i в формуле Чебышева

Пример 8.3. Вычислить интеграл из предыдущего примера по формуле Чебышева для четырех и для пяти точек в Mathcad.

Оценить точность вычислений.

Вычисление интеграла методом Чебышева для 5точек

Рис. 8.3. Решение примера 8.2 в Mathcad

Дата добавления: 2014-12-26; Просмотров: 1172; Нарушение авторских прав?; Мы поможем в написании вашей работы!