三次样条插值java算法

有n个已知点，(x0,y0),(x1,y1),...,(xn-1,yn-1).

矩阵运算求解m矩阵，所以用到ujmp-core这个jar 包，直接在maven中加入依赖

        <dependency>
            <groupId>org.ujmp</groupId>
            <artifactId>ujmp-core</artifactId>
            <version>0.3.0</version>
        </dependency>

定义样条插值的类

public class SplineInp {
    private double[] a = new double[10];
    private double[] b = new double[10];
    private double[] c = new double[10];
    private double[] d = new double[10];
    //已知点的x，y
    private long[] x;
    private double[] y;
    //计算系数
    public SplineInp(long[] x, double[] y){
        this.x = x;
        this.y = y;
        int n = x.length;
        //求a，总共有n个a
        for (int i = 0;i < n; i++){
            this.a[i] = y[i];
        }
        //求h[]
        double[] h = SplineInp.diff(x);
        //求左边矩阵A