python中的线性无关矩阵

Question

我有一个矩阵是15714x541，这个矩阵似乎不是线性独立的。如何删除非线性独立的列

我尝试使用此帖子中的解决方案： elimination the linear dependent columns of a non-square matrix in python

但是，它说： ValueError：项目长度541而不是15714错误

Q, R = np.linalg.qr(fixeff.T)
fixeff[np.abs(np.diag(R))>=1e-10]

fixeff是我描述的矩阵

Answer 1

你可以用scipy.linalg.qr来做：

from scipy.linalg import qr

Q, R, P = qr(A, mode="economic", pivoting=True)
inv = P.argsort() # reversed order (necessary here)
good_columns = (np.abs(np.diag(R)) > 1e-10)[inv]

A = A[:,good_columns]