//计算Ax=b;Eigen::Matrix3d A;A<<1,1,1,2,2,2,3,3,3;Eigen::Matrix<double, 3, 1> b;b<<1,2,3;cout<<"行列式为:"<<A.determinant()<<endl; //看一下行列式的值是否为0,确保可逆cout<<"x="<<A.inverse()*b<<endl; //直接求逆来解,速度会慢cout<<"x="<<A.colPivHouseholderQr().solve(b)<<endl; //QR分解,速度快,注意调用格式