如何在PARI/GP中计算任意mxn-矩阵的伪逆?有没有一种简单的方法,或者我必须完全对这个过程进行编程?
发布于 2015-04-22 02:58:14
Jörg Arndt在这里编写了(Moore-Penrose)伪逆的代码:
http://www.jjj.de/pari/
它出现在matsvd.gpi下的位置。
在矩阵具有完整列秩的简单(但常见)情况下,可以将其计算为
pseudoinverse(M) = my(ct=conj(M)~); (ct*M)^-1 * ct;发布于 2014-08-22 20:22:42
我想你是指摩尔-彭罗斯伪逆吧?
The tutorial and the manual for Pari/GP没有提到伪逆,因此您可能必须编写自己的解决方案。
Wikipedia entry可能help.You还可以在优秀的高级线性代数书籍中找到算法,例如Jonathan Golan的The Linear Algebra a Beginning Graduate Student Ought to Know。
https://stackoverflow.com/questions/24537661
复制相似问题