函數 > 數據分析 > 內插與預測 > 多項式內插
多項式內插
polyint(vx, vy, x) - 傳回使用多項式函數在 x 的內插值,以及預期誤差。
polyint 函數使用 Neville 演算法,在要求的點 x 上,對長度為 N 的數據集執行多項式內插。此函數會求解度數為 N – 1 且通過每個點的唯一多項式。
polycoeff(vx, vy) - 傳回內插多項式函數的係數。
polycoeff 函數會計算內插多項式的係數,以在後續計算中使用。
polyiter(vx, vy, x, N, e) - 使用多項式函數以及最大階數 N 與最大誤差 e,傳回 x 上的內插值。
polyiter 的輸出是向量,其中第一個元素是「收斂」標記 (1 = 已收斂,0 = 未收斂),第二個元素是符合指定公差所需的疊代數,而第三個元素是針對 x 的輸入值所取得的 y 疊代預估。
引數
vx, vy 是長度相同之數據值的實數向量。
x 是想要用以計算內插曲線的自變數值。為獲得最佳結果,x 應該位在 vx 值所包含的範圍中。
若輸入向量使用單位,則 x 必須與 vx 使用相同的單位。
N 是最大疊代數。N 也是多項式函數的最大階數,因為每次疊代計算之後,多項式的度數皆會增加一。
e 是輸入公差。
若輸入向量使用單位,則 e 必須與 vy 使用相同的單位。
其他資訊
polyiter 函數所使用的 Aitken-Neville 內插類似在 polyintpolycoeff 中執行的多項式內插。但由於內插疊代,polyiter 可讓您指定輸入公差 e,以及最大疊代數 N。若數據點的最後兩個疊代預估在公差 e 內取得一致,或若疊代數達到輸入引數 N,則此演算法會停止。例如,在定積分的 Romberg 積分中可使用疊代內插以更加突顯。數值積分是運算相當密集的流程。提早結束可節省處理時間,但代價是解答較不準確。Aitken-Neville 內插通常僅會用於求解一些具有某些特定公差的內插點。
polyintpolycoeff 常式以 Numerical Recipes in C, The Art of Scientific Computing (C 語言中的數值秘笈:科學計算的藝術) (Cambridge University Press,Copyright (C) 1987, 1988 Numerical Recipes Software) 一書的 polyint (第 109 頁) 與 polycoeff (第 121 頁) 為基礎,經授權使用。Thomas Richard McCalla (1967) 已對 polyiter 常式加以說明。書名為 Introduction to Numerical Methods and FORTRAN Programming (數值方法與 FORTRAN 程式設計的簡介) (John Wiley)。
這是否有幫助?