Квадратурные формулы Ньютона-Котеса, Чебышева, Гаусса. 



Квадратурные формулы Ньютона-Котеса

 

 

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

Так для функции f(x),  вычисленной в каких-то точках x0 ,x1 ,...,xm промежутка [a,b], можно построить многочлен Лагранжа

                                

и в результате интегрирования получить квадратурную формулу

                              ,                  ( 7 )

где Rm[f] - ошибка квадратур и

                                                                        ( 8 )

Хотя идея квадратур связана с интерполяционным многочленом Лагранжа, можно обойтись без его интегрирования.

Заметим, что для многочленов степени не выше m квадратура (7) является точной, т.е. Rm[f] =0 .

Полагая   f(x)  = xL ( L=0, 1, 2, ... , m), подстановкой в (7) получаем систему линейных уравнений для определения коэффициентов квадратурной формулы на выбранной системе узлов

                    .           ( 9 )

Большинство известных квадратурных формул строится на задании системы точек a+kh (k=0,1,..., m), где h=(b-a)/ m. Если пределы интегрирования a и b входят в состав узлов квадратурной формулы (7) , то ее называют квадратурой открытого типа (в противном случае - замкнутого типа). Создаваемые таким образом формулы называются формулами Ньютона-Котеса.

Рассмотрим пример построения какой-то квадратуры.

Создадим формулу открытого типа с тремя внутренними узлами a+h, a+2h, a+3h , где h=(b-a)/4 . Построив и решив систему

                    ,

получаем формулу вида

          . (10 )         

Создадим формулу закрытого типа с тремя узлами a, (a+b)/2, b . Построив и решив систему

                    ,

получаем известную формулу парабол ( формулу Симпсона )

                 .               (11)       

Если взять два узла a, b , то получаем формулу трапеций

                      .                              (12 )

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

     

Если представить квадратуры Ньютона-Котеса в стандартных формах:

1) для закрытого типа

          ,            ( 13 )

2) для открытого типа

                 ( 14 )

то можно легко построить любую из них с помощью таблиц

Что касается остаточных членов квадратурных формул, то их оценки зависят от четности числа (n) узлов формулы и типа (s=0 - закрытого и s=1 - открытого)

                             - n  нечетно ;           ( 15 )

                               - n  четно                  ( 16 )

(преимущество формул с нечетным числом узлов очевидно). Значения P и Q приведены в табл..

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

 

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

     П.Л. Чебышевым была предложена группа квадратурных формул с неравноотстоящими узлами, но c одинаковым весом слагаемых

                                                              ( 17 )

Очевидно, что эта формула применима к любому интегралу, т.к.

.

  Требуя, чтобы эта формула была точной для многочленов 1, t, t2,.., tn , получаем систему нелинейных уравнений:

                                               ( 18 )

Очевидно, что A=2 / n . Доказано, что приведенная система имеет действительные решения при значениях n £ 7 и n=9 , причем имеет место симметрия tk = tn-k+1 ( k=1,2,...,n/2 ).

Остаточный член квадратурной формулы Чебышева

                                 ,                  ( 19 )

где p(n) / q(n) приведены в табл.

 

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

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

Как было показано выше, всякий интеграл по конечному промежутку можно свести к интегралу по промежутку от -1  до 1

. ( 20 )

Возьмем квадратурную формулу

                                                     ( 21 )

и потребуем, чтобы она была точной для многочленов до 2n-1-го порядка (например, 1, t, t2, ..., t2n-1). Выполнение этого требования порождает систему 2n нелинейных уравнений с 2n неизвестными, решение которой отнюдь нетривиально даже с учетом симметрии узлов и коэффициентов

                            tk = tn-k+1 , Сk = Сn-k+1 ( k=1,2,...,n/2 ).

( при n=1 и n=2 формула Гаусса совпадает с формулой Чебышева).

Остаточный член n-точечной формулы Гаусса имеет вид [1]:

         . ( 22 )

Сравнение этой оценки с оценками для других квадратур показывает преимущество формул Гаусса.

 


Дата добавления: 2018-02-18; просмотров: 581; ЗАКАЗАТЬ РАБОТУ