Интерполирование с кратными узлами
— задача о построении
многочлена
минимальной
степени
, принимающего в некоторых точках (узлах
интерполяции
) заданные значения, а также заданные значения
производных
до некоторого
порядка
.
Показывается, что существует единственный многочлен
P
n
(
x
)
{\displaystyle \ P_{n}(x)}
степени
n
{\displaystyle \ n}
, удовлетворяющий условиям:
P
n
(
k
)
(
x
i
)
=
f
i
,
k
,
i
=
1
,
⋯
,
m
;
k
=
0
,
⋯
,
n
i
−
1
{\displaystyle P_{n}^{(k)}(x_{i})=f_{i,k},i=1,\cdots ,m;k=0,\cdots ,n_{i}-1}
, где
n
1
+
n
2
+
⋯
+
n
m
=
n
+
1
{\displaystyle n_{1}+n_{2}+\cdots +n_{m}=n+1}
.
Этот многочлен называют многочленом с кратными узлами, или
многочленом Эрмита
. В общем виде:
P
n
(
x
)
=
∑
i
=
1
m
∑
k
=
0
n
i
−
1
l
i
,
k
(
x
)
f
i
,
k
{\displaystyle P_{n}(x)=\sum _{i=1}^{m}\sum _{k=0}^{n_{i}-1}l_{i,k}(x)f_{i,k}}
,
m
{\displaystyle \ m}
— количество узлов и
n
i
{\displaystyle \ n_{i}}
— кратность узла
x
i
{\displaystyle \ x_{i}}
.
Шарль Эрмит
показал, что
l
i
,
k
(
x
)
=
[
1
k
!
∏
j
=
1
m
(
x
−
x
j
)
n
j
(
x
−
x
i
)
n
i
]
∑
s
=
0
n
i
−
k
−
1
c
s
i
(
x
−
x
i
)
k
+
s
{\displaystyle l_{i,k}(x)=\left[{\frac {1}{k!}}{\frac {\prod _{j=1}^{m}(x-x_{j})^{n_{j}}}{(x-x_{i})^{n_{i}}}}\right]\sum _{s=0}^{n_{i}-k-1}c_{s}^{i}(x-x_{i})^{k+s}}
, где
c
s
i
{\displaystyle \ c_{s}^{i}}
— коэффициенты
ряда Тейлора
для функции
(
x
−
x
i
)
n
i
∏
j
=
1
m
(
x
−
x
j
)
n
j
=
∑
s
=
0
∞
c
s
i
(
x
−
x
i
)
s
{\displaystyle {\frac {(x-x_{i})^{n_{i}}}{\prod _{j=1}^{m}(x-x_{j})^{n_{j}}}}=\sum _{s=0}^{\infty }c_{s}^{i}(x-x_{i})^{s}}
.
Доказательство
Частные случаи
Если все
n
i
{\displaystyle \ n_{i}}
равны единице, то интерполяционный многочлен Эрмита совпадает с
интерполяционным многочленом Лагранжа
.
Если количество узлов интерполяции равно единице, то интерполяционный многочлен Эрмита совпадает с
многочленом Тейлора
.
Если количество узлов интерполяции равно двум и в каждом задано значение функции и значение её производной — имеем задачу о построении
кубического сплайна
.
Оценка остатка интерполяции
См. также
Литература
Бахвалов Н. С., Численные методы, М., 1973.