numpy.linalg.solve返回矩阵不完整的解决方案

Question

我一直在学习OLS回归的基础知识，并被告知模型中没有完美的多重共线性。我有一个矩阵显然不完全排名：

np.linalg.matrix_rank(Matrix) == Matrix.shape[1]
   False

但是，numpy仍然可以求解方程

`y = np.random.randn(len(Matrix))
np.linalg.solve(Matrix.T.dot(Matrix),Matrix.T.dot(y))
array([ 8.65809521e-04,  5.16272205e-03, ...

没有完整等级的矩阵，这应该是不可以分析解决的。这是怎么回事？