C++最小二乘法求多项式系数

void Multinomial(Mat_& F, double * x, double * y, int m, int n)
{
    Mat_ X = Mat(m, n, CV_64FC1, Scalar(0));
    Mat_ Y = Mat(m, 1, CV_64FC1, Scalar(0));
    double t = 1;
    for (int i = 0; i < m; i++)
    {
        t = 1;
        X(i, 0) = 1;
        for (int j = 1; j < n; j++)
        {
            t = t * x[i]; //多项式的j阶项
            X(i, j) = t;
        }
        Y(i, 0) = y[i];
    }
    F = (X.t() * X).inv() * (X.t() * Y);
    F = F.t();
}

void main()
{
    double x[1026];
    double y[1026];
    ReadFile("data.txt", 1026, x, y);
    Mat_ F;
    Multinomial(F, x, y, 1026, 10);
    cout << F << endl;

}

运行结果：输出10个系数值