求伪逆矩阵c++代码（Eigen库）

最新推荐文章于 2023-08-05 17:38:33 发布

转载最新推荐文章于 2023-08-05 17:38:33 发布 · 1.7k 阅读

9 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/wxl845235800/p/8892681.html

文章标签：

#c/c++ #matlab #人工智能

本文介绍如何使用C++和Eigen库计算非方阵矩阵的伪逆，并提供了一个具体的实现示例。通过奇异值分解(SVD)的方法，能够有效地处理非方阵的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

非方阵的矩阵的逆矩阵 pseudoInverse

伪逆矩阵是逆矩阵的广义形式，广义逆矩阵

matlab中是pinv(A)--》inv(A)。

#include "stdafx.h"
#include<iostream> 
#include<Eigen/Core> 
#include<Eigen/SVD>   


template<typename _Matrix_Type_> 
_Matrix_Type_ pseudoInverse(const _Matrix_Type_ &a, double epsilon = 
    std::numeric_limits<double>::epsilon()) 
{  
    Eigen::JacobiSVD< _Matrix_Type_ > svd(a ,Eigen::ComputeThinU | Eigen::ComputeThinV);  
    double tolerance = epsilon * std::max(a.cols(), a.rows()) *svd.singularValues().array().abs()(0);  
    return svd.matrixV() *  (svd.singularValues().array().abs() > tolerance).select(svd.singularValues().array().inverse(), 0).matrix().asDiagonal() * svd.matrixU().adjoint(); 
}   
int main() 
{  
    Eigen::MatrixXd A(2,3);  
    A<< 1, 2, 3, 4, 5, 7;  
    std::cout<<A<<std::endl<<std::endl;  
    std::cout<<pseudoInverse(A)<<std::endl;   
    Eigen::MatrixXd B(3,2);  
    B<< 1, 2, 3, 4, 5, 7;  
    std::cout<<B<<std::endl<<std::endl;  
    std::cout<<pseudoInverse(B)<<std::endl;  
    getchar();
    return 0; 
}