linear-algebra - 用 PARI/GP 计算伪逆

Question

如何计算 PARI/GP 中任意 mxn- 矩阵的伪逆？有没有简单的方法，还是我必须完全编程这个过程？

Answer 1

Jörg Arndt 在这里为 (Moore-Penrose) 伪逆写了代码：

http://www.jjj.de/pari/

它出现在matsvd.gpi.

在矩阵具有完整列秩的简单（但常见）情况下，您可以将其计算为

pseudoinverse(M) = my(ct=conj(M)~); (ct*M)^-1 * ct;

Answer 2

我想你的意思是摩尔彭罗斯伪逆？

Pari/GP 的教程和手册没有提到伪逆，因此您可能必须编写自己的解决方案。

维基百科条目可能会有所帮助。您还可以在优秀的高级线性代数书籍中找到算法，例如 Jonathan Golan 的The Linear Algebra a Starting Graduate Student Aught to Know。