用PARI / GP计算伪逆

Question

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

Answer 1

JörgArndt在这里编写了（Moore-Penrose）伪逆的代码：

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

它出现在matsvd.gpi下。

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

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

Answer 2

我想你的意思是Moore-Penrose伪逆？

The tutorial and the manual for Pari/GP没有提到pseudoinverse，所以你可能需要编写自己的解决方案。

Wikipedia entry可能会有所帮助。您还可以在优秀的线性代数书籍中找到算法，例如Jonathan Golan的The Linear Algebra a Beginning Graduate Student Ought to Know。