请用C++实现matlab的spline函数

可以使用STL中的vector和map来实现matlab的spline函数。

首先，我们需要定义一个结构体表示样条插值的系数，包括a、b、c、d四个系数：

struct SplineCoef {
    double a, b, c, d;
};

然后，我们可以编写一个函数来计算样条插值的系数：

vector<SplineCoef> spline(vector<double>& x, vector<double>& y) {
    int n = x.size();
    vector<double> h(n - 1), alpha(n - 1), l(n), mu(n - 1), z(n);
    vector<SplineCoef> coef(n - 1);

    for (int i = 0; i < n - 1; i++) {
        h[i] = x[i + 1] - x[i];
        alpha[i] = 3.0 / h[i] * (y[i + 1] - y[i]) - 3.0 / h[i - 1] * (y[i] - y[i - 1]);
    }

    l[0] = 1.0;
    mu[0] = 0.0;
    z[0] = 0.0;

    for (int i = 1; i < n - 1; i++) {
        l[i] = 2.0 * (x[i + 1] - x[i - 1]) - h[i - 1] * mu[i - 1];
        mu[i] = h[i] / l[i];
        z[i] = (alpha[i - 1] - h[i - 1] * z[i - 1]) / l[i];
    }

    l[n - 1] = 1.0;
    z[n - 1] = 0.0;
    coef[n - 2].c = 0.0;

    for (int j = n - 3; j >= 0; j--) {
        coef[j].c = z[j] - mu[j] * coef[j + 1].c;
        coef[j].b = (y[j + 1] - y[j]) / h[j] - h[j] * (coef[j + 1].c + 2.0 * coef[j].c) / 3.0;
        coef[j].d = (coef[j + 1].c - coef[j].c) / (3.0 * h[j]);
        coef[j].a = y[j];
    }

    return coef;
}

最后，我们可以实现一个函数来使用样条插值计算函数值：

double spline_eval(vector<SplineCoef>& coef, double x) {
    int n = coef.size();
    if (x <= coef[0].a) return coef[0].a;
    if (x >= coef[n - 1].a) return coef[n - 1].a;

    int i = 0, j = n - 1;
    while (i < j) {
        int k = (i + j) / 2;
        if (x <= coef[k].a) j = k;
        else i = k + 1;
    }

    double dx = x - coef[i].a;
    return coef[i].a + coef[i].b * dx + coef[i].c * dx * dx + coef[i].d * dx * dx * dx;
}

这样，我们就可以使用spline和spline_eval函数来进行样条插值了。