Julia三次样条插值函数

# 本文将演示如何使用Julia语言实现三次样条插值
# 如果没有安装Dierckx包。
# 首先要Pkg.add("Dierckx"),然后Pkg.update()
# 当然也可以用Interpolations.jl,但对其语法不太熟悉
# 用下面代码做示例。(https://github.com/kbarbary/Dierckx.jl)
using Dierckx
x = linspace(0,pi)
v = sin.(x)
xq = linspace(pi/2,pi)

# 计算三次样条插值
spl = Spline1D(x, v; k = 1)  # k: order of interpolant; can be between 1-5
vq = spl(xq)  # 等价于vq = Spline1D(x, v;  k = 1)(xq)

# result as following:
################################################################

# 接下来使用Interpolations来作三次样条。(https://github.com/JuliaMath/Interpolations.jl)
# 如果没有安装，先安装此包。Pkg.add("Interpolations")
# 以下为示例

using Interpolations
f(x) = log(x)
xs = 1:0.2:10
A = [f(x) for x in xs]
# linear interpolation
interp_linear = LinearInterpolation(xs, A)
xt = 3:0.01:8
B = interp_linear( xt )

# cubic spline interpolation
interp_cubic = CubicSplineInterpolation(xs, A)
C = interp_cubic( xt )