Postać Newtona wzoru interpolacyjnego

by Jerry Sky

2020-11-17



Postać Newtona

pn(x)=∑k=0nckqk(x)=∑k=0nf[x0,x1,…,xk]∏j=0k−1(x−xj). p_n(x) = \sum_{k=0}^n c_k q_k(x) = \sum_{k=0}^n f[x_0, x_1, \dots, x_k] \prod_{j=0}^{k-1} (x - x_j).

Współczynniki c0=f[x0], c1=f[x0,x1], …, cn=f[x0,x1,…,xn]c_0 = f[x_0],\, c_1 = f[x_0, x_1],\, \dots,\, c_n = f[x_0, x_1, \dots, x_n] wyznaczamy, konstruując tablicę trójkątną ilorazów różnicowych. Konstrukcja wielomianu interpolacyjnego pn+1(x)p_{n+1}(x): pn(xi)=yi(0≤i≤n)pn+1(xi)=yi(0≤i≤n+1) \begin{aligned} p_n(x_i) &= y_i \qquad (0 \le i \le n)\\ p_{n+1}(x_i) &= y_i \qquad (0 \le i \le n+1) \end{aligned}

Wielomian pn+1(x)p_{n+1}(x) można przedstawić w postaci pn+1(x)=pn(x)+cn+1(x−x0)⋯(x−xn). p_{n+1}(x) = p_n(x) + c_{n+1}(x - x_0) \dotsb (x - x_n).

Zatem wyznaczenie pn+1(x)p_{n+1}(x) na podstawie pn(x)p_n(x) wymaga wyznaczenia cn+1=f[x0,…,xn+1]c_{n+1} = f[x_0, \dots, x_{n+1}] (lokalna modyfikacja).