length = np.linalg.matrix_rank(mul,tol=0.0000000000000001) ##求秩
谢谢大佬，我这样设置后，求得的秩是4.但是逆矩阵仍然报异常。
via1 = np.linalg.inv(mul) ##求逆
报： raise LinAlgError("Singular matrix")
numpy.linalg.LinAlgError: Singular matrix
【 在 snoopyzhao 的大作中提到: 】
:
: 我没有你说的这两个软件
: 你的原始数据精度是什么样子的？你 print 出来的是一个单精度的矩阵，而内部存贮的可能是双精度，这种精度损失可能会影响到你的结果
: ...................
--
FROM 120.242.253.*