Формула прямоугольников


Аппроксимируем подынтегральную функцию полиномом нулевой степени – константой f(x*), где х* – единственный узел на отрезке [xi–1, xi]. Тогда

.

В качестве х* обычно выбирают середину отрезка (известно, что симметрия повышает точность численных формул):

x* = (xi–1+xi)/2=xi–1/2 .

Таким образом, получаем формулу «центральных» прямоугольников для элементарного отрезка:

, где hi = xi–xi–1; yi–1/2 = f(xi–1/2).

Просуммировав, получаем составную (обобщенную) формулу «центральных» прямоугольников:

. (3.7)

Для равномерной сетки ( h=(b–a)/n ) формула примет вид:

.

В дальнейшем будем рассматривать равномерную сетку.

Теперь оценим погрешность формулы «центральных» прямоугольников.

Сначала попробуем это сделать с помощью погрешности интерполяции (n=0; w(x)=x – xi–1/2).

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

Оценим погрешность метода с помощью разложения в ряд Тейлора:

где .

Для составной формулы «центральных» прямоугольников

,

где .

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

 

Примечание 3. О формулах «левых» и «правых» прямоугольников.

Можно легко показать (например, с помощью погрешности интерполяции), что при ином выборе узловой точки х* квадратурные формулы «левых» прямоугольников и «правых» прямоугольников, соответственно

имеют первый порядок точности.

Конец примечания 3.

Формула трапеций

Аппроксимируем подынтегральную функцию полиномом первой степени, построенным по двум узлам: xi–1 и xi . Интерполяционный полином Лагранжа первой степени имеет вид

.

Тогда

.

Получили формулу трапеций для элементарного отрезка [xi–1, xi].

Суммируя, получаем составную формулу трапеций:

. (3.8)

Оценим погрешность формулы трапеции.

Тогда .

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

 

Примечание 4. О симметричных квадратурных формулах.

Отметим еще один интересный факт. Формула прямоугольников получена при использовании полинома нулевой степени, а формула трапеций – первой. Но порядок точности один и тот же, хотя при увеличении степени полинома следовало бы ожидать увеличения точности. Это проявление симметрии формулы прямоугольников (сравните для примера порядки точности несимметричной формулы левых прямоугольников и формулы трапеций – соответствие степени интерполяционного полинома и порядка точности налицо).

Признаки симметричности квадратурных формул:

1) n – четное;

2) узлы расположены симметрично относительно середины отрезка [a,b], то есть ;

3) коэффициенты симметричны: Ci=Cn– i .

Конец примечания 4.

Формула Симпсона

Аппроксимируем подынтегральную функцию полиномом второй степени, построенным по трем узлам: xi–1, xi–1/2, xi . Интерполяционный полином Лагранжа второй степени имеет вид

 

Тогда

Получили формулу Симпсона для элементарного отрезка [xi–1,xi].

Суммируя, получаем составную формулу Симпсона:

(3.9)

Оценим погрешность формулы Симпсона.

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

будет безрезультатным (проверьте!). Поэтому используем разложение в ряд Тейлора узловых значений функции в окрестности точки x=xi–1/2:

.

Для составной формулы Симпсона

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

 

Примечание 5. О формуле Симпсона без полуцелых точек.

Формулу Симпсона можно записать по-другому, если пронумеровать насквозь и целые, и полуцелые точки.

Тогда n=2m (то есть nчетное число), новый шаг интегрирования h1 = h/2=(b–a)/(2m), а формула имеет вид:

(3.10)

Погрешность

.

Конец примечания 5.

 



Дата добавления: 2020-10-25; просмотров: 470;


Поиск по сайту:

Воспользовавшись поиском можно найти нужную информацию на сайте.

Поделитесь с друзьями:

Считаете данную информацию полезной, тогда расскажите друзьям в соц. сетях.
Poznayka.org - Познайка.Орг - 2016-2024 год. Материал предоставляется для ознакомительных и учебных целей.
Генерация страницы за: 0.01 сек.