matlab构造矩阵代码-PySPQR:用于SuiteSparseQR中的稀疏QR分解的Python包装器

matlab构造矩阵代码适用于SuiteSparseQR的Python包装器 该模块包装了分解函数以与一起使用。 这是Matlab的稀疏[Q,R,E] = qr() 。 由于某种原因，没有人包装过适用于Python的SuiteSparseQR的功能。 对于A x = b且稀疏A和稠密或稀疏b的情况，还包装了SuiteSparseQR解算器。 这对于解决最小二乘意义上的稀疏超定线性系统特别有用。 在此， A的大小为m×n， b为m×k（存储k不同的右侧向量，每个向量分别考虑）。 用法 import numpy import scipy . sparse . linalg import sparseqr # QR decompose a sparse matrix M such that Q R = M E # M = scipy . sparse . rand ( 10 , 10 , density = 0.1 ) Q , R , E , rank = sparseqr . qr ( M ) print ( abs ( Q * R - M * sparseqr . permutation