双数阵列的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，这可能会引起一些兴趣。