解决Python中的奇异值分解(SVD)
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了解决Python中的奇异值分解(SVD),小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含1194字,纯文字阅读大概需要2分钟。
内容图文
![解决Python中的奇异值分解(SVD)](/upload/InfoBanner/zyjiaocheng/675/eeadcf131a4e49c39ce2fba1b9f87004.jpg)
我正在尝试将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很好地翻译.下一步是IN IDL(!)
x = svsol(u,w,v,b)
python和IDL中的u几乎相同(其他矩阵也是如此).唯一的区别是维度,其中IDL的矩阵较大,但有很多零.在这种情况下,似乎Python的矩阵更加压缩.
有谁知道类似的Python知识.
如果有人需要,这里是svsol的手册.
解决方法:
使用IDL中的SVDC和SVSOL,您可以通过SVD分解解决线性最小二乘问题.这是通过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语义,这是svsol reference 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)全部内容,希望文章能够帮你解决解决Python中的奇异值分解(SVD)所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。