利用三階多項式插值法,產生 x 範圍為 0.0 起之至 20.0 的 200 個均勻插值點,以 pgplot 繪出一曲線通過下表列點:

  xi       yi
-------------------
0.000000 0.000000
2.000000 0.9092974
4.000000 -0.7568025
6.000000 -0.2794155
8.000000 0.9893582
10.00000 -0.5440211
12.00000 -0.5365729
14.00000 0.9906074
16.00000 -0.2879033
18.00000 -0.7509872
20.00000 0.9129453

表列點與數值的檔案 tab_gen.dat

提示:

hunt 副程式來找尋每個插值點左側之表列點。

用公式 k = min( max(j-(m-1)/2,1) , N+1-m ) 來決定每次插值時之第一個表列點的指標,其中 N 是總表列值數,m 是多項式階數,j 是最接近插值點左側之表列點的指標。

 

三階多項式

 

 

補充參考:若採用其他階數時所畫出的圖形效果

四階多項式

 

五階多項式