Pn(x)=y0+k=1∑nΔky0⋅i=0∏k−1(x−xi), где:
Pn(x)— интерполяционный многочлен степени n;
y0 — значение функции в первом узле (x0);
Δky0 — k-я разделённая разность, вычисленная для узлов x0,…,xk;
∏i=0k−1(x−xi) — произведение линейных множителей для k-го члена.
Нулевая разность (значение функции)
\Delta^0 y[x_i] = y_i.
$
Первая разделённая разность
Δ1y[xi,xi+1]=xi+1−xiyi+1−yi. Общая формула ( k )-й разделённой разности
Δky[xi,xi+1,…,xi+k]=xi+k−xiΔk−1y[xi+1,…,xi+k]−Δk−1y[xi,…,xi+k−1].
Произведение линейных множителей
Для каждого члена интерполяционного многочлена:
i=0∏k−1(x−xi)=(x−x0)(x−x1)…(x−xk−1), где k — текущий порядок разности.
Итоговая структура
Общий вид многочлена:
Pn(x)=y0+Δ1y[x0,x1]⋅(x−x0)+Δ2y[x0,x1,x2]⋅(x−x0)(x−x1)+… +Δny[x0,x1,…,xn]⋅(x−x0)(x−x1)…(x−xn−1). Каждый член состоит из:
Коэффициента Deltaky[x0,…,xk];
Произведения линейных множителей ∏i=0k−1(x−xi).
Last updated