Eigen库-矩阵的伪逆（Matlab-pinv）

Question

我正在尝试使用Eigen库找到矩阵的伪逆。他们有一个实现它的类，但我不知道如何编写语法。

这是网站上显示的方式(https://eigen.tuxfamily.org/dox/classEigen_1_1CompleteOrthogonalDecomposition.html#ab2fd4c81aa1cd8bc917c7f135505cb7f):

const Inverse Eigen::CompleteOrthogonalDecomposition< MatrixType > :: pseudoInverse() const

- Hisham Mohammad

---

你的问题究竟是什么？请修改你的问题，提供一个最小可复现示例（MCVE）以及一个清晰的问题陈述。 - Baum mit Augen

Answer 1

这是CompleteOrthogonalDecomposition类的一种方法。因此，在使用它之前，您必须对矩阵执行该分解。例如：

#include <Eigen/QR>    

Eigen::MatrixXd A = ... // fill in A
Eigen::MatrixXd pinv = A.completeOrthogonalDecomposition().pseudoInverse();

- Sean

---

虽然在大多数情况下建议使用其他求解器。 - Adi Shavit

是的，绝对可以。 - Sean

@AdiShavit 为什么建议使用其他求解器？ - Dsel

@AdiShavit，你能给出参考资料吗？你推荐哪些求解器？ - ignacio

@ignacio：Eigen文档表示不建议使用它（不高效或数值稳定性差：请参见https://eigen.tuxfamily.org/dox/classEigen_1_1CompleteOrthogonalDecomposition.html#a3c89639299720ce089435d26d6822d6f）。请参见https://eigen.tuxfamily.org/dox/group__LeastSquares.html，了解在Eigen中解决线性系统的更好建议。 - Sandy

A.completeOrthogonalDecomposition().pseudoInverse() 的结果和 MATLAB 中的 pinv() 不同。请参考 https://gitlab.com/libeigen/eigen/-/issues/2462 获取更详细的信息。 - 劉大為