numpy.linalg.solve возвращает решение для матрицы, которая не является полным рангом - PullRequest
Новая
       9

numpy.linalg.solve возвращает решение для матрицы, которая не является полным рангом

0 голосов
/ 19 марта 2019

Я изучал основы регрессии OLS, и мне сказали, что в модели не может быть идеальной мультиколинарности. У меня есть матрица, которая явно не имеет полного ранга: 

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

Но NumPy все еще может решить уравнение enter image description here

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

Это не должно быть аналитически разрешимо без матрицы полного ранга. Что здесь происходит?

...