写出求正规方程的c++算法

以下是用C++实现求解正规方程的简单算法：

#include <iostream>
#include <Eigen/Dense> // Eigen库用于矩阵运算

using namespace Eigen;

int main() 
{
    MatrixXd X(3,2); // 自变量矩阵X，大小为3*2
    VectorXd y(3); // 因变量向量y，大小为3*1
    X << 1, 2, 1, 3, 1, 4; // 填充自变量矩阵X
    y << 1, 2, 3; // 填充因变量向量y
    VectorXd w(2); // 系数向量w，大小为2*1
    w = (X.transpose() * X).inverse() * X.transpose() * y; // 求解正规方程
    std::cout << "系数向量w为：" << std::endl << w << std::endl; // 输出系数向量w
    return 0;
}

在上述代码中，我们使用了Eigen库来进行矩阵运算。其中，MatrixXd表示矩阵类型，VectorXd表示向量类型。我们首先定义了自变量矩阵X和因变量向量y，并分别对其进行填充。然后，我们定义了系数向量w，并使用正规方程的公式求解得到w的值。最后，我们输出系数向量w的值。