Polinomial Newton

Dalam analisis numerik, polinomial Newton adalah interpolasi polinomial untuk suatu himpunan titik data yang diketahui. Polinomial ini dinamai dari penemunya, Isaac Newton.Templat:R Terkadang, polinomial ini disebut interpolasi polinomial beda terbagi Newton karena koefisien dari polinomialnya dihitung menggunakan metode beda terbagi Newton.

Diberikan suatu himpunan $k+1$ titik data

$(x_0,y_0),\dots ,(x_i,y_i),\dots ,(x_k,y_k)$

dengan syarat $x_p\ne x_q$ apabila $p\ne q$, maka interpolasi polinomial Newton merupakan suatu kombinasi linear dari polinomial basis Newton. Secara matematis, hal tersebut dapat dituliskan sebagai

$N(x):=\sum _{i=0}^k{a}_i{n}_i(x)$

dengan polinomial basis Newton didefinisikan sebagai

$n_i(x):=\prod _{j=0}^{i-1}(x-x_j)$

untuk $i>0$ dan $n_0(x)=1$. Koefisien dari polinomial tersebut didefinisikan sebagai $a_i:=[y_0,\dots ,y_i]$ dengan $[y_0,\dots ,y_i]$ adalah notasi untuk beda terbagi. Dengan demikian, interpolasi polinomial Newton dapat ditulis sebagai $$N(x)=[y_0]+[y_0,y_1](x-x_0)+\dots +[y_0,\dots ,y_k](x-x_0)(x-x_1)\dots (x-x_{k-1})$$

Rumus interpolasi Newton-Gregory maju

Bentuk dari Polinomial Newton dapat disederhanakan apabila $x_0,x_1,\dots ,x_k$ diatur sehingga berjarak sama. Jika didefinisikan

• $h:=x_{i+1}-x_i$
• $x-x_0=sh$, untuk suatu bilangan riil $s$
• ${\mathrm{\Delta }}^0{y}_i:=y_i$
• ${\mathrm{\Delta }}^{j+1}{y}_i:=\left({\mathrm{\Delta }}^j{y}_{i+1}\right)-\left({\mathrm{\Delta }}^j{y}_i\right)$

dengan $i,j=0,1,\dots ,k-1$, perhatikan bahwa

$\begin{array}{cc}x-x_1& =x-x_0-x_1+x_0\\ & =x-x_0-(x_1-x_0)\\ & =sh-h\\ & =h(s-1)\\ \\ x-x_2& =x-x_0-x_2+x_1-x_1+x_0\\ & =x-x_0-(x_2-x_1)-(x_1-x_0)\\ & =sh-h-h\\ & =h(s-2)\\ & \dots \\ x-x_i& =h(s-i)\\ ---------& ---------------------\\ \left[y_0\right]& :=y_0\\ [y_0,y_1]& =\frac{y_1-y_0}{x_1-x_0}\\ & =\frac{\mathrm{\Delta }{y}_0}{h}\\ \\ [y_0,y_1,y_2]& =\frac{\frac{y_2-y_1}{x_2-x_1}-\frac{y_1-y_0}{x_1-x_0}}{x_2-x_0}\\ & =\frac{\frac{\mathrm{\Delta }{y}_1}{h}-\frac{\mathrm{\Delta }{y}_0}{h}}{2h}\\ & =\frac{{\mathrm{\Delta }}^2{y}_0}{2!h^2}\\ \\ [y_0,y_1,y_2,y_3]& =\frac{\frac{\frac{y_3-y_2}{x_3-x_2}-\frac{y_2-y_1}{x_2-x_1}}{x_3-x_1}-\frac{\frac{y_2-y_1}{x_2-x_1}-\frac{y_1-y_0}{x_1-x_0}}{x_2-x_0}}{x_3-x_0}\\ & =\frac{\frac{\frac{\mathrm{\Delta }{y}_2}{h}-\frac{\mathrm{\Delta }{y}_1}{h}}{2h}-\frac{\frac{\mathrm{\Delta }{y}_1}{h}-\frac{\mathrm{\Delta }{y}_0}{h}}{2h}}{3h}\\ & =\frac{\frac{{\mathrm{\Delta }}^2{y}_1}{2h^2}-\frac{{\mathrm{\Delta }}^2{y}_0}{2h^2}}{3h}\\ & =\frac{{\mathrm{\Delta }}^3{y}_0}{3!h^3}\\ & \dots \\ [y_0,y_1,\dots ,y_k]& =\frac{{\mathrm{\Delta }}^k{y}_0}{k!h^k}\end{array}$

Sehingga, interpolasi polinomial Newton menjadi

$\begin{array}{cc}N(x)& =[y_0]+[y_0,y_1](x-x_0)+[y_0,y_1,y_2](x-x_0)(x-x_1)+\dots +[y_0,\dots ,y_k](x-x_0)(x-x_1)\dots (x-x_{k-1})\\ & =y_0+\frac{\mathrm{\Delta }{y}_0}{h}(sh)+\frac{{\mathrm{\Delta }}^2{y}_0}{2!h^2}(sh)(h(s-1))+\dots +\frac{{\mathrm{\Delta }}^k}{k!h^k}(sh)(h(s-1))\dots (h(s-k+1))\\ & =y_0+s\left(\mathrm{\Delta }{y}_0\right)+\frac{s(s-1)}{2!}\left({\mathrm{\Delta }}^2{y}_0\right)+\dots +\frac{s(s-1)\dots (s-k+1)}{k!}\left({\mathrm{\Delta }}^k{y}_0\right)\\ & =(\genfrac{}{}{0ex}{}{s}{0})\left({\mathrm{\Delta }}^0{y}_0\right)+(\genfrac{}{}{0ex}{}{s}{1})\left({\mathrm{\Delta }}^1{y}_0\right)+(\genfrac{}{}{0ex}{}{s}{2})\left({\mathrm{\Delta }}^2{y}_0\right)+\dots +(\genfrac{}{}{0ex}{}{s}{k})\left({\mathrm{\Delta }}^k{y}_0\right)\\ & =\sum _{i=0}^k(\genfrac{}{}{0ex}{}{s}{i})\left({\mathrm{\Delta }}^i{y}_0\right)\end{array}$

Rumus ini disebut sebagai rumus beda terbagi Newton maju.Templat:Citation needed

Rumus interpolasi Newton-Gregory mundur

Jika titik data yang ada diubah urutannya menjadi $x_k,x_{k-1},\dots ,x_0$, bentuk polinomial Newtonnya menjadi

$N(x)=[y_k]+[y_k,y_{k-1}](x-x_k)+\dots +[y_k,\dots ,y_0](x-x_k)(x-x_{k-1})\dots (x-x_1)$

Dengan asumsi serupa, jika $x_k,x_{k-1},\dots ,x_0$ berjarak sama, serta didefinisikan

• $h:=x_{i+1}-x_i$
• $x-x_k=sh$, untuk suatu bilangan riil $s$
• ${\mathrm{\Delta }}^0{y}_i:=y_i$
• ${\mathrm{\Delta }}^{j+1}{y}_i:=\left({\mathrm{\Delta }}^j{y}_{i+1}\right)-\left({\mathrm{\Delta }}^j{y}_i\right)$

dengan $i,j=0,1,\dots ,k-1$, maka diperoleh

$\begin{array}{cc}x-x_{k-1}& =(x-x_k)+(x_k-x_{k-1})\\ & =sh+h\\ & =h(s+1)\\ \\ x-x_{k-2}& =(x-x_k)+(x_k-x_{k-1})+(x_{k-1}-x_{k-2})\\ & =sh+h+h\\ & =h(s+2)\\ \dots \\ x-x_{k-j}& =h(s+j)\\ \end{array}$

Sehingga, interpolasi polinomial Newton menjadi

$\begin{array}{cc}N(x)& =[y_k]+[y_k,y_{k-1}](x-x_k)+[y_k,y_{k-1},y_{k-2}](x-x_k)(x-x_{k-1})+\dots +[y_k,\dots ,y_0](x-x_k)(x-x_{k-1})\dots (x-x_1)\\ & =y_k+\frac{\mathrm{\Delta }{y}_{k-1}}{h}(sh)+\frac{{\mathrm{\Delta }}^2{y}_{k-2}}{2!h^2}(sh)(h(s+1))+\dots +\frac{{\mathrm{\Delta }}^k{y}_0}{k!h^k}(sh)(h(s+1))\dots (h(s+k-1))\\ & =y_k+s\left(\mathrm{\Delta }{y}_{k-1}\right)+\frac{s(s+1)}{2!}\left({\mathrm{\Delta }}^2{y}_{k-2}\right)+\dots +\frac{s(s+1)(s+2)\dots (s+k-1)}{k!}\left({\mathrm{\Delta }}^k{y}_0\right)\\ & =(\genfrac{}{}{0ex}{}{s-1}{0})\left({\mathrm{\Delta }}^0{y}_k\right)+(\genfrac{}{}{0ex}{}{s}{1})\left(\mathrm{\Delta }{y}_{k-1}\right)+(\genfrac{}{}{0ex}{}{s+1}{2})\left({\mathrm{\Delta }}^2{y}_{k-2}\right)+\dots +(\genfrac{}{}{0ex}{}{s+k-1}{k})\left({\mathrm{\Delta }}^k{y}_0\right)\\ & =\sum _{j=0}^k(\genfrac{}{}{0ex}{}{s+j-1}{j})\left({\mathrm{\Delta }}^j{y}_{k-j}\right).\end{array}$

Rumus ini disebut sebagai rumus beda terbagi Newton mundurTemplat:Citation needed

Bentuk Newton dan polinomial Taylor

Rumus Newton sangat menarik karena rumus Newton terlihat seperti versi beda hingga dari polinomial Taylor. Polinomial Taylor memberitahu nilai suatu fungsi berdasarkan nilai $y$ beserta turunan-turunannya (turunan pertama, turunan kedua, dst.) dari suatu titik $x$ tertentu. Rumus Newton adalah polinomial Taylor yang didasarkan pada beda hingga, dan bukan laju perubahan sesaat.

Limit dari polinomial Newton jika semua titik berhimpit menuju $x=a$ adalah polinomial Taylor di sekitar $x=a$, lantaran nilai dari koefisien beda terbagi akan menjadi turunan. Hal ini akan lebih mudah terlihat menggunakan rumus Newton-Gregory maju $$\begin{array}{cc}\underset{(x_0,\dots ,x_n)\to (a,\dots ,a)}{\lim }& [y_0]+[y_0,y_1]\cdot (x-x_0)+[y_0,y_1,y_2]\cdot (x-x_0)(x-x_1)+\dots +[y_0,\dots ,y_n]\cdot (x-x_0)(x-x_1)\dots (x-x_{n-1})\\ =\underset{(x_0,\dots ,x_n)\to (a,\dots ,a)}{\lim }& y_0+\frac{\mathrm{\Delta }{y}_0}{h}\cdot (x-x_0)+\frac{{\mathrm{\Delta }}^2{y}_0}{2!h^2}\cdot (x-x_0)(x-x_1)+\dots +\frac{{\mathrm{\Delta }}^n{y}_0}{n!h^n}\cdot (x-x_0)(x-x_1)\dots (x-x_{n-1})\\ & =f(a)+f^{\prime }(a)\cdot (x-a)+\frac{f^{″}(a)}{2!}\cdot (x-a{)}^2\dots +\frac{f^{(n)}(a)}{n!}\cdot (x-a{)}^n\end{array}$$

Ide utama

Penyelesaian masalah interpolasi akan mengarah ke permasalahan dalam aljabar linear mengenai penyelesaian sistem persamaan linear. Dengan menggunakan basis monomial, maka akan diperoleh matriks Vandermonde yang sangat rumit. Dengan memilih basis lain, basis Newton, maka akan diperoleh sistem persamaan linear dengan matriks segitiga bawah, yang lebih sederhana dan dapat diselesaikan dengan lebih cepat.

Jika diberikan $k+1$ titik data, maka dikonstruksikan basis Newton sebagai

$n_0(x):=1n_j(x):=\prod _{i=0}^{j-1}(x-x_i)j=1,2,\dots ,k$

Dengan menggunakan polinomial ini sebagai basis yang baru, maka penyelesaian dari

$\left[\begin{array}{ccccc}1& & \dots & & 0\\ 1& x_1-x_0& & & \\ 1& x_2-x_0& (x_2-x_0)(x_2-x_1)& & ⋮\\ ⋮& ⋮& & \ddots & \\ 1& x_k-x_0& \dots & \dots & \prod _{j=0}^{k-1}(x_k-x_j)\end{array}\right]\left[\begin{array}{c}{a}_0\\ \\ ⋮\\ \\ a_k\end{array}\right]=\left[\begin{array}{c}{y}_0\\ \\ ⋮\\ \\ y_k\end{array}\right]$

akan menjawab masalah interpolasi polinomial di awal.

Sistem persamaan ini dapat diselesaikan secara iteratif dengan menyelesaikan

$\sum _{i=0}^j{a}_i{n}_i(x_j)=y_{j}j=0,1,\dots ,k.$

Penambahan titik baru

Sama seperti rumus selisih lainnya, derajat dari interpolasi polinomial Newton dapat dinaikkan dengan menambahkan suku dan titik baru tanpa membuang suku yang sudah ada. Bentuk Newton memiliki kesederhanaan karena titik baru akan selalu ditambahkan di ujung: titik baru pada rumus Newton-Gregory maju dapat ditambahkan di kanan, sedangkan titik baru pada rumus Newton-Gregory mundur dapat ditambahkan di kiri.

Keunggulan dan kelemahan berbagai rumus

Untuk setiap himpunan titik data yang berhingga, hanya terdapat satu polinomial dengan derajat terendah yang melewati semua titik yang diberikan. Maka dari itu, kata "interpolasi polinomial Bentuk Newton" (atau bentuk Lagrange, dll.) tidak akan memiliki makna yang ambigu. Meskipun demikian, pencarian polinomial ini bisa saja memiliki efisiensi komputasi yang berbeda apabila digunakan metode yang berbeda.

Newton versus Lagrange

Templat:Further Dibandingkan interpolasi polinomial bentuk Newton, bentuk Lagrange memerlukan usaha yang lebih sedikit, dan terkadang direkomendasikan untuk masalah yang sudah diketahui (dari pengalaman sebelumnya) berapa banyak suku yang diperlukan agar tingkat akurasinya cukup.

Jika ditambahkan titik data yang baru, metode beda terbagi mempunyai keunggulan karena suku yang diperoleh dari data sebelumnya dapat digunakan kembali. Dengan rumus Lagrange yang biasa, pengerjaan soal dengan penambahan titik data akan membutuhkan perhitungan ulang seluruh soal.

Secara umum, rumus beda terbagi lebih serbaguna pada berbagai situasi permasalahan. Dengan interpolasi polinomial bentuk Newton, terdapat algoritma yang singkat dan efektif untuk mencari koefisien dari polinomialnya.

Penerapan

Seperti yang bisa dilihat dari bentuk umum dari interpolasi polinomial Newton, titik data yang baru dapat ditambahakan ke himpunan data yang sudah ada untuk membuat interpolasi polinomial yang baru tanpa perlu menghitung ulang koefisien yang lama. Kalaupun terdapat titik data yang berubah nilainya, maka seluruh koefisien tidak perlu dihitung ulang. Selain itu, jika $x_0,x_1,\dots ,x_n$ memiliki jarak yang sama antara satu dengan yang lain, perhitungan beda terbagi akan jauh lebih mudah. Maka dari itu, rumus beda terbagi lebih disukai dibandingkan bentuk Lagrange untuk diterapkan.

Contoh

Nilai beda terbagi dapat dituliskan dalam bentuk tabel. Sebagai contoh, untuk mencari suatu fungsi $f$ yang menginterpolasi titik $\{(x_0,y_0),(x_1,y_1),\dots ,(x_n,y_n)\}$, maka dapat ditulis sebagai berikut

$\begin{array}{cccccc}x& [y_i]& [y_i,y_{i+1}]& [y_i,y_{i+1},y_{i+2}]& [y_i,y_{i+1},y_{i+2},y_{i+3}]& \dots \\ ---& ---& ------& --------& -----------& ------\\ x_0& y_0& & & & \\ & & \frac{y_1-y_0}{x_1-x_0}& & & \\ x_1& y_1& & \frac{\frac{y_2-y_1}{x_2-x_1}-\frac{y_1-y_0}{x_1-x_0}}{x_2-x_0}& & \\ & & \frac{y_2-y_1}{x_2-x_1}& & \frac{\frac{\frac{y_3-y_2}{x_3-x_2}-\frac{y_2-y_1}{x_2-x_1}}{x_3-x_1}-\frac{\frac{y_2-y_1}{x_2-x_1}-\frac{y_1-y_0}{x_1-x_0}}{x_2-x_0}}{x_3-x_0}& \\ x_2& y_2& & \frac{\frac{y_3-y_2}{x_3-x_2}-\frac{y_2-y_1}{x_2-x_1}}{x_3-x_1}& & \\ & & \frac{y_3-y_2}{x_3-x_2}& & ⋮& \\ x_3& y_3& & ⋮& & \\ & & ⋮& & ⋮& \\ ⋮& & & ⋮& & \\ & & ⋮& & \\ x_n& y_n& & & & \end{array}$

Maka koefisien dari interpolasi polinomialnya diperoleh dari entri teratas pada setiap kolom.

Sebagai contoh, misalkan dikonstruksikan interpolasi polinomial dari $f(x)=\sin x$ dengan menggunakan metode beda terbagi pada titik berikut

$n$	$x_n$	$y_n$
$0$	$0$	$0$
$1$	$0.1$	$0.09983$
$2$	$0.3$	$0.29552$
$3$	$0.6$	$0.56464$
$4$	$1$	$0.84147$

Dengan menggunakan akurasi sebesar lima digit, maka diperoleh tabel beda terbagi sebagai berikut

$\begin{array}{cccccc}0& 0& & & & \\ & & 0.9983& & & \\ 0.1& 0.09983& & -0.06617& & \\ & & 0.97845& & -0.16102& \\ 0.3& 0.29552& & -0.16278& & 0.01651\\ & & 0.89706& & -0.14451& \\ 0.6& 0.56464& & -0.29284& & \\ & & 0.69208& & & \\ 1& 0.84147& & & & \end{array}$

Sehingga, interpolasi polinomialnya adalah

$\begin{array}{cc}N(x)& =0+0.9983(x-0)-0.06617(x-0)(x-0.1)-0.16102(x-0)(x-0.1)(x-0.3)+0.01651(x-0)(x-0.1)(x-0.3)(x-0.6)\\ & =0.999789x+0.0026957x^2-0.17753x^3+0.01651x^4\end{array}$

Lihat juga

• De numeris triangularibus et inde de progressionibus arithmeticis: Magisteria magna, karya dari Thomas Harriot yang menjelaskan metode interpolasi yang serupa, ditulis 50 tahun lebih awal dari karya Newton namun baru diterbitkan pada tahun 2009
• Deret Newton
• Algoritma Neville
• Interpolasi Polinomial
• Interpolasi polinomial Bentuk Lagrange
• Interpolasi polinomial Bentuk Bernstein
• Interpolasi Hermite
• Teorema Carlson

Pranala luar

• Modul tentang Polinomial Newton, karya John H. Mathews