Postać Newtona wzoru interpolacyjnego

by Jerry Sky

2020-11-17



Postać Newtona

pn(x)=k=0nckqk(x)=k=0nf[x0,x1,,xk]j=0k1(xxj). 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(0in)pn+1(xi)=yi(0in+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(xx0)(xxn). 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).