Формулы Ньютона — Котса

Формула Ньютона — Котса для $n=2$

Формулы Ньютона — Котса (Котеса), называемые также правилами квадратуры Ньютона — Котса или просто правилами Ньютона — Котса, — это группа формул для численного интегрирования (называемых также квадратурами), основанных на вычислении интегрируемой функции в одинаково отстоящих друг от друга точках. Формулы названы именами Исаака Ньютона и Роджера Котса.

Формулы Ньютона — Котса полезны, когда заданы значения интегрируемой функции на точках, отстоящих друг от друга на одинаковом расстоянии. Если можно менять положение точек, могут оказаться более пригодными другие методы, такие как метод Гаусса и Шаблон:Нп5.

Описание

Предполагается, что значения функции Шаблон:Mvar определены на отрезке $[a,b]$ и известны в $n+1$ точке $a⩽x_0<x_1<\dots <x_n⩽b$, расположенных на одинаковых расстояниях друг от друга. Если $x_0=a$ и $x_n=b$, то есть используются значения функции на границах интервала, то функция называется квадратурой «замкнутого» типа, а если $x_0>a$ и $x_n<b$, то есть значения функции в крайних точках интервала не используются, то «открытого» типаШаблон:Sfn. Формулы Ньютона — Котса, использующие $n+1$ точек, могут быть определены (для обоих случаев) какШаблон:Sfn

$\underset{a}{\overset{b}{\int }}f(x)dx\approx \sum _{i=0}^n{A}_{i}f(x_i)$,

где

• для формулы замкнутого типа $x_i=a+ih$, где $h=\frac{b-a}{n}$,
• для формулы открытого типа $x_i=a+(i+1)h$, где $h=\frac{b-a}{n+2}$.

Число Шаблон:Mvar называется размером шага, а $A_i$ называется квадратурным коэффициентомШаблон:Sfn.

$A_i$ можно вычислить как интегралы от базисных многочленов Лагранжа, которые зависят только от $x_i$ и не зависят от функции Шаблон:Mvar. Пусть $L(x)$ — интерполяционный многочлен в форме Лагранжа для заданных точек $(x_0,f(x_0)),\dots ,(x_n,f(x_n))$, тогда

$\underset{a}{\overset{b}{\int }}f(x)dx\approx \underset{a}{\overset{b}{\int }}L(x)dx=\underset{a}{\overset{b}{\int }}(\sum _{i=0}^{n}f(x_i)l_i(x))dx=\sum _{i=0}^{n}f(x_i)\underset{A_i}{\underbrace{\underset{a}{\overset{b}{\int }}{l}_i(x)dx}}.$

Неустойчивость для высоких степеней

Можно построить формулы Ньютона — Котса любой степени Шаблон:Mvar. Однако для больших Шаблон:Mvar правило Ньютона — Котса может иногда страдать от феномена РунгеШаблон:Sfn, когда ошибка растёт экспоненциально для больших Шаблон:Mvar. Такие методы, как квадратура Гаусса или квадратура Кленшоу — Кёртиса — с неравными расстояниями между точками (имеющими бо́льшую плотность на концах интервала интегрирования) — устойчивы и более точны, а потому обычно более предпочтительны, чем квадратура Ньютона — Котса. Если эти методы нельзя использовать, то есть если значения интегрируемого выражения заданы только в фиксированной сетке с одинаковыми расстояниями, можно избежать феномена Рунге путём использования разбиения интервала, как разъяснено ниже.

Также устойчивые формулы Ньютона — Котса можно построить, если заменить интерполяцию на метод наименьших квадратов. Это позволяет записать численно устойчивые формулы даже для высоких степенейШаблон:RШаблон:R.

Формулы Ньютона — Котса замкнутого типа

Следующая таблица содержит перечисление некоторых формул Ньютона — Котса замкнутого типа. Для $0⩽i⩽n$ пусть $x_i=a+i\frac{b-a}{n}=a+ih$, а обозначение $f_i$ служит сокращением для $f(x_i)$.

Замкнутые формулы Ньютона — Котса

Шаблон:Mvar	Размер шага Шаблон:Mvar	Общее название	Формула	Ошибка
1	$b-a$	Метод трапеций	$\frac{h}{2}(f_0+f_1)$	$-\frac{1}{12}{h}^3{f}^{(2)}(\xi )$
2	$\frac{b-a}{2}$	Формула Симпсона	$\frac{h}{3}(f_0+4f_1+f_2)$	$-\frac{1}{90}{h}^5{f}^{(4)}(\xi )$
3	$\frac{b-a}{3}$	Формула Симпсона 3/8	$\frac{3h}{8}(f_0+3f_1+3f_2+f_3)$	$-\frac{3}{80}{h}^5{f}^{(4)}(\xi )$
4	$\frac{b-a}{4}$	Шаблон:Нп5	$\frac{2h}{45}(7f_0+32f_1+12f_2+32f_3+7f_4)$	$-\frac{8}{945}{h}^7{f}^{(6)}(\xi )$

Правило Буля иногда ошибочно называют правилом Боде, как результат типографической опечатки в книге Абрамовица и СтиганШаблон:SfnШаблон:R.

Степень размера сегмента h в ошибке показывает скорость, с которой убывает ошибка аппроксимации. Порядок производной функции f в ошибке даёт наименьшую степень многочлена, который не может быть вычислен точно (то есть с нулевой ошибкой) по этому правилу. Число $\xi$ должно быть взято из интервала (a, b).

Формулы Ньютона — Котса открытого типа

Таблица показывает некоторые формулы Ньютона — Котса открытого типа. Снова, $f_i$ служит сокращённой записью для $f(x_i)$, где $x_i=a+i\frac{b-a}{n+2}$.

Открытые формулы Ньютона — Котса

Шаблон:Mvar	Размер шага Шаблон:Mvar	Общее название	Формула	Ошибка
0	$\frac{b-a}{2}$	Сумма Римана или средняя сумма Римана	$2h{f}_1$	$\frac{1}{3}{h}^3{f}^{(2)}(\xi )$
1	$\frac{b-a}{3}$		$\frac{3}{2}h(f_1+f_2)$	$\frac{3}{4}{h}^3{f}^{(2)}(\xi )$
2	$\frac{b-a}{4}$	Формула Милна	$\frac{4}{3}h(2f_1-f_2+2f_3)$	$\frac{14}{45}{h}^5{f}^{(4)}(\xi )$
3	$\frac{b-a}{5}$		$\frac{5}{24}h(11f_1+f_2+f_3+11f_4)$	$\frac{95}{144}{h}^5{f}^{(4)}(\xi )$

Разбиение интервала

Чтобы формула Ньютона — Котса была более точной, нужно, чтобы длина Шаблон:Mvar была мала. Значит, интервал интегрирования $[a,b]$ сам должен быть маленьким, что в большинстве случаев не так. По этой причине обычно численное интегрирование осуществляется путём разбиения интервала $[a,b]$ на меньшие подынтервалы, на каждом из которых применяется формула Ньютона — Котса, после чего результаты складываются. См. статью Численное интегрирование.

См. также

• Квадратура (математика)
• Интерполяция
• Кубический сплайн

Литература

• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга
• Шаблон:Книга

Ссылки

• Шаблон:SpringerEOM
• Newton-Cotes formulas on www.math-linux.com
• Шаблон:MathWorld
• Newton-Cotes Integration, numericalmathematics.com