本文介绍了用Python求解奇异值分解(SVD)的处理方法，对大家解决问题具有一定的参考价值，需要的朋友们下面随着小编来一起学习吧！
问题描述
我正在尝试将一个IDL程序翻译成Python。我必须通过以下方式解决SVD的结果
from scipy.linalg import svd
A = [[1,2,3],[4,5,6]]
b = [4,4,5]
u,w,v = svd(A)
这段代码运行得很好，并且可以很好地从IDL翻译过来。下一步是在IDL(！)
中x = svsol(u,w,v,b)
PYTHON和IDL中的u几乎相同(对于其他矩阵也是如此)。唯一的区别是维度，IDL的矩阵更大，但有很多零。从这个意义上看，Python的矩阵似乎压缩得更多。
有没有人知道类似的Python。
如果有人需要，这里是svsol的使用手册。
推荐答案
使用SVDC和SVSOL可以通过奇异值分解来解决线性最小二乘问题。这是由numpy.linalg.lstsq函数在numpy中完成的。(不需要先计算SVD分解，然后再反向求解。)
>>> import numpy as np
>>> A = np.array([[1,2,3],[4,5,6]])
>>> b = np.array([4,4])
>>> x, _, _, _ = np.linalg.lstsq(A,b)
>>> x
array([-2., 0., 2.])
>>> np.dot(A,x)
array([ 4., 4.])
b的长度必须与A的行数相同，因此您的示例是错误的。为了确保我正确解释IDL语义，下面是svsolreference manual：中的示例
>>> A = np.array(
... [[1.0, 2.0, -1.0, 2.5],
... [1.5, 3.3, -0.5, 2.0],
... [3.1, 0.7, 2.2, 0.0],
... [0.0, 0.3, -2.0, 5.3],
... [2.1, 1.0, 4.3, 2.2],
... [0.0, 5.5, 3.8, 0.2]])
>>> B = np.array([0.0, 1.0, 5.3, -2.0, 6.3, 3.8])
>>> x, _, _, _ = np.linalg.lstsq(A,B)
>>> print x
[ 1.00095058 0.00881193 0.98417587 -0.01009547]
这篇关于用Python求解奇异值分解(SVD)的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持编程学习网！
本站部分内容来源互联网,如果有图片或者内容侵犯您的权益请联系我们删除！