polyfit C语言实现

在C语言中，可以使用最小二乘法来实现多项式拟合（Polyfit），可以使用线性代数库，比如GNU Scientific Library (GSL) 来实现。下面是一个简单的示例代码：

#include <stdio.h>
#include <gsl/gsl_fit.h>

#define N 5 // 数据点个数
#define ORDER 3 // 多项式拟合阶数

int main() {
    double x[N] = {1, 2, 3, 4, 5}; // 自变量 X 的数据
    double y[N] = {2, 3, 5, 6, 8}; // 因变量 Y 的数据
    double c[ORDER + 1]; // 存储拟合结果的系数

    double cov[ORDER * ORDER]; // 存储协方差矩阵
    double chisq; // 存储卡方值

    gsl_fit_linear(x, 1, y, 1, N, c, cov, &chisq); // 进行多项式拟合

    printf("拟合结果：\n");
    for (int i = 0; i <= ORDER; i++) {
        printf("c[%d] = %g\n", i, c[i]);
    }

    return 0;
}

这个示例代码使用了 GSL 库中的 `gsl_fit_linear` 函数进行多项式拟合。你可以根据需要修改自变量 x 和因变量 y 的数据，以及指定拟合的阶数 ORDER。拟合结果将存储在数组 c 中，每个元素对应一个拟合系数。

请注意，为了编译这个代码，你需要先安装 GSL 库，并将其链接到你的项目中。