Вторая интерполяционная формула Ньютона

Первая интерполяционная формула Ньютона

Интерполяционный многочлен Ньютона для равноотстоящих узлов

Рассмотрим случай, когда h=x_i+1 – x_i=const (i=0, 1, …).

Рассмотрим конечные разности:

– конечные разности 1-го порядка – разности между значениями функции в соседних узлах.

– конечные разности 2-го порядка – разности между конечными разностями 1-го порядка.

– конечные разности 3-го порядка.

…

– конечные разности k-го порядка.

Конечные разности удобно вычислять в таблице:

Будем искать интерполяционный многочлен Ньютона в виде многочлена n-ой степени:

P_n(x) = a₀+ a₁(x-x₀) + a₂(x-x₀)(x-x₁) +…+ a_n(x-x₀)…(x-x_n-1) (3)

Коэффициенты a₀, a₁, …, a_n находятся из условия совпадения значения исходной функции f(x) и многочлена P_n(x) в узлах интерполяции: .

Пусть , тогда , соответственно

…

Подставив в формулу (3), получим:

– первая интерполяционная формула Ньютона.

Погрешность вычислений оценивается следующим образом:

Так при n=2

, где

Вторая интерполяционная формула Ньютона применяется, когда значение аргумента ближе к концу отрезка интерполяции (формула для интерполирования назад). Интерполяционный многочлен ищется в виде многочлена n-ой степени:

P_n(x) = a₀+ a₁(x-x_n) + a₂(x-x_n)(x-x_n-1) +…+ a_n(x-x_n)…(x-x₁) (4)

Коэффициенты a₀, a₁, …, a_n находятся из условия совпадения значения исходной функции f(x) и интерполяционного многочлена P_n(x) в узлах: .

Подставив a_k в (4) и перейдя к переменной , получим вторую интерполяционную формулу Ньютона:

.

Погрешность вычислений оценивается следующим образом:

.

Рассмотрим задачу субтабулирования (уплотнения таблицы) функции на отрезке. Введем следующие обозначения:

a, b – концы субтабулирования;

H0 – старый шаг таблицы;

H – новый шаг таблицы;

y1, y2, y3 – конечные разности 1-го, 2-го, 3-го порядка;

d – границы погрешности метода.

Для вычисления конечных разностей составляется таблица:

Блок-схема уплотнения таблиц функций:

Программа уплотнения таблиц функций (субтабулирования)

program subtab;

var a,b,d,h0,h,y,y0,y1,y2,y3,x,t : real;

begin write('Введите a, b, H0, H - '); readln(a,b,H0,H);

write('Введите Y0, конечные раз-ти Y1, Y2, y3 - ');

readln(Y0,Y1,Y2,y3);

writeln(' X Y D');

x:=a;

while x<=b do

begin t:=(x-a)/h0;

y:=y0+t*y1+t*(t-1)*y2/2;

d:=y3*t*(t-1)*(t-2)/6;

writeln(x:8:4, y:12:6, d:14:8);

x:=x+h

end;

readln;

end.