c++多样式拟合代码

多样式拟合是指通过数学方法，将一组离散的、具有多样式特征的数据点拟合成一个平滑的曲线或函数。下面以C语言为例，介绍一段多样式拟合的代码。

#include <stdio.h>
#include <math.h>

// 定义多项式拟合函数，根据输入的数据点和阶数，返回一个多项式函数的系数数组
void polyfit(double *x, double *y, int n, int order, double *coeff)
{
    int i, j;
    double temp;
    double mat[order + 1][order + 2] = {0};

    // 构造方程组的系数矩阵
    for (i = 0; i <= order; i++) {
        for (j = 0; j <= order; j++) {
            for (temp = 0, k = 0; k < n; k++) {
                temp += pow(x[k], i + j);
            }
            mat[i][j] = temp;
        }
    }

    // 构造方程组的常数项
    for (i = 0; i <= order; i++) {
        for (temp = 0, k = 0; k < n; k++) {
            temp += pow(x[k], i) * y[k];
        }
        mat[i][order + 1] = temp;
    }

    // 高斯消元解方程
    for (i = 0; i <= order; i++) {
        for (j = i + 1; j <= order; j++) {
            temp = mat[j][i] / mat[i][i];
            for (k = 0; k <= order + 1; k++) {
                mat[j][k] -= temp * mat[i][k];
            }
        }
    }

    // 回代求解多项式系数
    for (i = order; i >= 0; i--) {
        temp = mat[i][order + 1];
        for (j = i + 1; j <= order; j++) {
            temp -= mat[i][j] * coeff[j];
        }
        coeff[i] = temp / mat[i][i];
    }
}

int main()
{
    double x[] = {1, 2, 3, 4, 5}; // 数据点的x坐标
    double y[] = {3, 4, 2, 5, 6}; // 数据点的y坐标
    int n = sizeof(x) / sizeof(x[0]); // 数据点的个数
    int order = 2; // 拟合多项式的阶数
    double coeff[order + 1]; // 多项式系数

    polyfit(x, y, n, order, coeff);

    // 打印多项式拟合结果
    printf("拟合多项式为：");
    for (int i = order; i >= 0; i--) {
        printf("%lf ", coeff[i]);
    }
    printf("\n");

    return 0;
}

该段代码实现了多项式拟合，根据输入的数据点和阶数，通过最小二乘法求解多项式拟合的系数，并将结果打印出来。代码首先定义了一个polyfit函数，该函数用于求解多项式拟合的系数，然后再在主函数中调用该函数来完成拟合并输出结果。最后，通过GCC编译器将代码编译为可执行文件，运行后即可得到多项式拟合的结果。