範例:Thiele 內插
在數據集上執行 Thiele 連續分數內插。
Thielecoeff 與 Thiele
使用 thielecoeff 與 thiele 函數內插來自有理函數的數據集。
1. 定義內插點數。
2. 定義有理函數,例如 Lorentzian 函數。
3. 取樣函數並將其 x 與 y 值記錄到向量中。
4. 使用 Thielecoeff 函數求解擬合樣本點之連續分數的係數。
5. 呼叫 Thiele 函數執行內插。
Thiele 函數成功,但 rationalint 函數在此內插範圍中的某幾點上失敗。這是因為數據對稱,因而導致在 x 的某些中間值上發生除以零的錯誤:
6. 繪製原始函數、樣本點及內插曲線。
7. 繪製表示 Thiele 內插之殘差的圖形。
點的順序
若兩個鄰接點具有相同的 y 值,則此演算法無法傳回正確的值。您必須重新排序數據,才可取得正確的結果。
1. 定義一組點。
2. 求解 Thiele 函數的係數。
係數太大。
3. 重新排序數據以取得適當的解。
4. 重複步驟 2。
5. 定義內插函數。
6. 繪製內插函數與原始數據點。