用 PARI/GP 计算伪逆

【问题标题】：calculating the pseudoinverse with PARI/GP用 PARI/GP 计算伪逆
【发布时间】：2014-08-23 15:07:59
【问题描述】：

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

【问题讨论】：

【解决方案1】：

Jörg Arndt 在此处编写了 (Moore-Penrose) 伪逆的代码：

它出现在matsvd.gpi 下的位置。

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

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

【讨论】：

【解决方案2】：

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

The tutorial and the manual for Pari/GP 没有提及伪逆，因此您可能必须编写自己的解决方案。

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

【讨论】：