大矩阵的Matlab逆

Question

这是我要求解的方程式:

h = (X'*X)^-1*X'*y

其中 X 是矩阵，y 是向量((X'X)^-1 是 X 转置乘以 X 的倒数)。我在 Matlab 中将其编码为:

h = (X'*X)\X'*y

我认为这是正确的。问题是 X 大约是 10000x10000，并且试图计算逆函数会使 Matlab 在我能找到的最强大的计算机(16 核，24GB RAM)上崩溃。有没有什么方法可以将其拆分，或者有一个库是为进行如此大的反演而设计的吗？

谢谢。

Answer 1

这看起来像是伪逆。你也许只是在寻找

h = X \ y;