雙數陣列的Akima插值

Question

假設我有一個雙精度數組，什麼是使用Akima interpolation來抽樣這個系列的好算法？我太愚蠢了，不能將數學描述翻譯成代碼。

// values is an array of doubles 
// idx is the index of the left-hand value for the current interpolation 
// t is the normalized parameter between values[idx] and values[idx+1] 
// Don't worry about array bounds, I'll handle that separately. 
public double InterpolateAkima(double[] values, int idx, double t) 
{ 
    ...? 
} 

Answer 1

轉發並擴充我對another SO question的回答，該回復作爲該問題的重複被關閉 - 正如對該問題的評論所暗示的那樣。

阿克瑪的原始紙： ``基於本地程序 '' 內插和平滑曲線擬合的新方法，雜誌ACM 17，4（1970）的，589-602

http://www.leg.ufpr.br/lib/exe/fetch.php/wiki:internas:biblioteca:akima.pdf

Ç實施

https://github.com/ampl/gsl/blob/master/interpolation/akima.c

C＃實現

https://gist.github.com/dreikanter/3526685

Delphi實現（見在Delphi/SRC/spline3.pas過程BuildAkimaSpline）

http://www.alglib.net/translator/re/alglib-2.6.0.delphi.zip

阿克瑪的的Fortran 66實施

http://cran.r-project.org/web/packages/akima/

Fortran 90的實施

http://miyoshi.googlecode.com/svn-history/r72/trunk/common/common.f90

Java實現

https://commons.apache.org/proper/commons-math/jacoco/org.apache.commons.math3.analysis.interpolation/AkimaSplineInterpolator.java.html

Lisp實現 「爲AutoCAD二維多段線」

http://autocad.xarch.at/code/candido/akima.lsp

matlab實現

http://www.mathworks.se/matlabcentral/fileexchange/1814-akima-interpolation

帕斯卡實現（program description）

http://jean-pierre.moreau.pagesperso-orange.fr/Pascal/akima_pas.txt

Python實現

http://www.lfd.uci.edu/~gohlke/code/akima.py.html

VB6實現（見在vb6/SRC/spline3子程序BuildAkimaSpline。BAS）

http://www.alglib.net/translator/re/alglib-2.6.0.vb6.zip

http://www.koders.com/cpp/fid1393B9D668316C1700966643DE0609660B9CB13A.aspx?s=%22Brian+Smith%22

Answer 2

在google代碼搜索上找到了一些hits，但這不是我知道的地方。第一個結果是Math.NET，這可能會引起一些興趣。