Интерполяционная формула Ньютона используется для нахождения многочлена, проходящего через заданные точки ( x 0 , y 0 ) , ( x 1 , y 1 ) , … , ( x n , y n ) (x_0, y_0), (x_1, y_1), \ldots, (x_n, y_n) ( x 0 , y 0 ) , ( x 1 , y 1 ) , … , ( x n , y n ) . Многочлен строится по форме, удобной для пошагового добавления новых узлов.
Общая формула:
P n ( x ) = y 0 + ∑ k = 1 n Δ k y 0 ⋅ ∏ i = 0 k − 1 ( x − x i ) , P_n(x) = y_0 + \sum_{k=1}^{n} \Delta^k y_0 \cdot \prod_{i=0}^{k-1} (x - x_i), P n ( x ) = y 0 + k = 1 ∑ n Δ k y 0 ⋅ i = 0 ∏ k − 1 ( x − x i ) , где Δ k y 0 \Delta^k y_0 Δ k y 0 — конечные разности (или разделённые разности), а x 0 , x 1 , … , x n x_0, x_1, \ldots, x_n x 0 , x 1 , … , x n — узлы.
Построение формулы Ньютона
Исходные данные : Таблица узлов:
x y x 0 y 0 x 1 y 1 ⋮ ⋮ x n y n \begin{array}{|c|c|}
x & y \\
\hline
x_0 & y_0 \\
x_1 & y_1 \\
\vdots & \vdots \\
x_n & y_n \\
\end{array} x x 0 x 1 ⋮ x n y y 0 y 1 ⋮ y n Вычисление разделённых разностей : Для k k k -го порядка разности:
Δ k y = Δ k − 1 y [ x 1 , … , x k ] − Δ k − 1 y [ x 0 , … , x k − 1 ] x k − x 0 . \Delta^k y = \frac{\Delta^{k-1} y[x_1, \ldots, x_k] - \Delta^{k-1} y[x_0, \ldots, x_{k-1}]}{x_k - x_0}. Δ k y = x k − x 0 Δ k − 1 y [ x 1 , … , x k ] − Δ k − 1 y [ x 0 , … , x k − 1 ] . Построение интерполяционного многочлена : Подставляем значения разностей в формулу Ньютона.
Пример
Условие: Найти интерполяционный многочлен для узлов:
x y 1 2 2 3 4 7 \begin{array}{|c|c|}
x & y \\
\hline
1 & 2 \\
2 & 3 \\
4 & 7 \\
\end{array} x 1 2 4 y 2 3 7 Шаг 1. Вычисляем разделённые разности:
Δ 0 y = y i (начальные значения) . \Delta^0 y = y_i \quad \text{(начальные значения)}. Δ 0 y = y i ( начальные значения ) . Δ 1 y [ x 0 , x 1 ] = y 1 − y 0 x 1 − x 0 = 3 − 2 2 − 1 = 1. \Delta^1 y[x_0, x_1] = \frac{y_1 - y_0}{x_1 - x_0} = \frac{3 - 2}{2 - 1} = 1. Δ 1 y [ x 0 , x 1 ] = x 1 − x 0 y 1 − y 0 = 2 − 1 3 − 2 = 1. Δ 1 y [ x 1 , x 2 ] = y 2 − y 1 x 2 − x 1 = 7 − 3 4 − 2 = 2. \Delta^1 y[x_1, x_2] = \frac{y_2 - y_1}{x_2 - x_1} = \frac{7 - 3}{4 - 2} = 2. Δ 1 y [ x 1 , x 2 ] = x 2 − x 1 y 2 − y 1 = 4 − 2 7 − 3 = 2. Δ 2 y [ x 0 , x 1 , x 2 ] = Δ 1 y [ x 1 , x 2 ] − Δ 1 y [ x 0 , x 1 ] x 2 − x 0 = 2 − 1 4 − 1 = 1 3 . \Delta^2 y[x_0, x_1, x_2] = \frac{\Delta^1 y[x_1, x_2] - \Delta^1 y[x_0, x_1]}{x_2 - x_0} = \frac{2 - 1}{4 - 1} = \frac{1}{3}. Δ 2 y [ x 0 , x 1 , x 2 ] = x 2 − x 0 Δ 1 y [ x 1 , x 2 ] − Δ 1 y [ x 0 , x 1 ] = 4 − 1 2 − 1 = 3 1 . Шаг 2. Формула Ньютона:
P 2 ( x ) = y 0 + Δ 1 y [ x 0 , x 1 ] ⋅ ( x − x 0 ) + Δ 2 y [ x 0 , x 1 , x 2 ] ⋅ ( x − x 0 ) ( x − x 1 ) . P_2(x) = y_0 + \Delta^1 y[x_0, x_1] \cdot (x - x_0) + \Delta^2 y[x_0, x_1, x_2] \cdot (x - x_0)(x - x_1). P 2 ( x ) = y 0 + Δ 1 y [ x 0 , x 1 ] ⋅ ( x − x 0 ) + Δ 2 y [ x 0 , x 1 , x 2 ] ⋅ ( x − x 0 ) ( x − x 1 ) . Подставляем значения:
P 2 ( x ) = 2 + 1 ⋅ ( x − 1 ) + 1 3 ⋅ ( x − 1 ) ( x − 2 ) . P_2(x) = 2 + 1 \cdot (x - 1) + \frac{1}{3} \cdot (x - 1)(x - 2). P 2 ( x ) = 2 + 1 ⋅ ( x − 1 ) + 3 1 ⋅ ( x − 1 ) ( x − 2 ) . Упрощаем:
P 2 ( x ) = 2 + ( x − 1 ) + 1 3 ⋅ ( x − 1 ) ( x − 2 ) . P_2(x) = 2 + (x - 1) + \frac{1}{3} \cdot (x - 1)(x - 2). P 2 ( x ) = 2 + ( x − 1 ) + 3 1 ⋅ ( x − 1 ) ( x − 2 ) . P 2 ( x ) = 2 + x − 1 + 1 3 ⋅ ( x 2 − 3 x + 2 ) . P_2(x) = 2 + x - 1 + \frac{1}{3} \cdot (x^2 - 3x + 2). P 2 ( x ) = 2 + x − 1 + 3 1 ⋅ ( x 2 − 3 x + 2 ) . P 2 ( x ) = 1 3 x 2 − x + 8 3 . P_2(x) = \frac{1}{3} x^2 - x + \frac{8}{3}. P 2 ( x ) = 3 1 x 2 − x + 3 8 . Итог
Многочлен Ньютона второго порядка:
P 2 ( x ) = 1 3 x 2 − x + 8 3 . P_2(x) = \frac{1}{3} x^2 - x + \frac{8}{3}. P 2 ( x ) = 3 1 x 2 − x + 3 8 . Last updated 5 months ago