如何用特征库计算零空间的基础？

Question

如何计算零基空间与特征库矩阵？

我试图找到明确的函数名计算零基础，也是，作为一种解决方法，找到方法矩阵的计算RREF（如我们能够从RREF获得零基础）。

但我找不到任何相关的函数名称。

我认为必须有解决方案，但我对Eigen库和Eigen的代码知之甚少，这也很难让我理解。

请建议我解决这个问题。

Answer 1

你可以使用Eigen::FullPivLU::kernel()方法零空间的基础：

FullPivLU<MatrixXd> lu(A); 
MatrixXd A_null_space = lu.kernel(); 

Answer 2

替代方案：用OpenCV的计算零空间：

` 
cv::Mat EpipolarConstraint::getNullSpace(cv::Mat p) 
{ 
    cv::SVD svd = cv::SVD(p, cv::SVD::FULL_UV); 
    cv::Mat vt_ = svd.vt; 
    int i; 
    for (i = 1; i <= 3; i++) 
    { 
     if (p.at<double>(i - 1, i - 1) == 0) 
     { 
      break; 
     } 
    } 
    cv::Mat result = vt_(cv::Rect(0, i-1, p.cols, vt_.rows-i+1)); 
    cv::Mat result_t; 
    cv::transpose(result, result_t); 
    return result_t; 
}`