Оценка погрешности полиномиальной интерполяции

Оценим остаточный член, т.е. разность Rn-(x)=f(x)-Ln(x)

Потребуем, чтобы в области [ a, b] изменения x функция f(x) была дифференцируема (n+1) раз, т.е. f’(x),…, f(n+1) (x).

Ясно, что в узлах интерполяции эта погрешность Rn(x)=0.

Оценим ее в любой точке x Î [ a, b].

Рассмотрим вспомогательную функцию

u(x)=f(x)-Ln(x)-kПn+1(x), где

- K - постоянный коэффициент, который выберем ниже.

U(x) имеет (n+1) нуль (корень) в т. x0, x1, …, xn.

Подберем k так, чтобы u(x) имела n+2 корень в любой фиксированной точке x’ отрезка [ a, b ], не совпадающий с узлами интерполирования. Т.е. u(x’)=0.

Положим f(x’)=Ln(x’)-kПn+1(x’).

Тогда

При этом значении k функция u(x) будет иметь (n+2) корня в
узлах x’,x0, x1, …, xn.

Применяя теорему Ролля на каждом из отрезков [x0, x1],[x1, x2],…,
[xi, x’], [x’
, xi],…,[xn-1, xn] убеждаемся, что производная u’(x) имеет не менее (n+1) корня на [ a, b].

U’’(x) – n раз обращается в нуль на [ a, b].

………………………………………………

U (n+1)(x) - имеет 1 корень

Пусть U (n+1)(x)=0, xÎ[ a, b].

Так как Ln(n+1)(x)=0 и Пn+1(n+1)(x)=(n+1)!, имеем

При x=x, получаем

®

Т.о. Rn(x)=f(x)-Ln(x)=?

В силу того, что U(x’)=0, имеем

Т.к. x’ - произвольно, то

, где xÎ[ a, b].

Обозначим через

Получим формулу для погрешности интерполяционной формулы Лагранжа

Если f(x) - многочлен n -ой степени, то f(x)=Ln(x) и Rn(x)=0 в силу единственности.

Пример:

С какой точностью можно вычислить Ö115 с помощью Ln(x) для y=Öx, если x0=100, x1=121, x2=144.

Y’=0.5x-1/2, y’’=-1/4x-3/2, y’’’=3/8x-5/2. M3=max|y’’|=3/8 1/Ö1005=3/2 10-5.

|R2|£3/8 10-5 1/3! |(115-100)(1150121)(115-144)|=1.6 10-3.


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: