当前位置:   article > 正文

【生物信息学】奇异值分解(SVD)_np.linalg.svd

np.linalg.svd

目录

一、奇异值分解(SVD)

二、Python实现

1. 调包np.linalg.svd()

2. 自定义

三、SVD实现链路预测


一、奇异值分解(SVD)

        SVD分解核心思想是通过降低矩阵的秩来提取出最重要的信息,实现数据的降维和去噪。

ChatGPT:

        SVD(奇异值分解)是一种常用的矩阵分解方法,它可以将一个矩阵分解为三个矩阵的乘积:U、Σ和V^T。这里,U和V是正交矩阵,Σ是一个对角矩阵。SVD分解的详细过程如下:

        假设我们有一个m×n的矩阵A,其中m表示行数,n表示列数。SVD分解的目标是找到三个矩阵U、Σ和V^T,使得:A = UΣV^T

        其中,U是一个m×m的正交矩阵,Σ是一个m×n的对角矩阵,V^T是一个n×n的正交矩阵。正交矩阵满足U^TU = I和V^TV = I,其中I是单位矩阵。

对角矩阵Σ的对角线上的元素称为奇异值,通常用σ表示。奇异值按照从大到小的顺序排列,因此可以将它们写成σ₁ ≥ σ₂ ≥ ⋯ ≥ σᵣ,其中r是矩阵A的秩。

SVD分解的步骤如下:

  1. 计算矩阵A的转置A^T与A的乘积AA^T,得到一个m×m的对称矩阵。

  2. 对对称矩阵AA^T进行特征值分解,得到特征值和对应的特征向量。特征值按照从大到小排列。

  3. 根据特征值和特征向量,构造正交矩阵U。U的列向量是对应于AA^T的特征向量,按照特征值从大到小排列。

  4. 计算矩阵A的乘积A^TA,得到一个n×n的对称矩阵。

  5. 对对称矩阵A^TA进行特征值分解,得到特征值和对应的特征向量。特征值按照从大到小排列。

  6. 根据特征值和特征向量,构造正交矩阵V。V的列向量是对应于A^TA的特征向量,按照特征值从大到小排列。

  7. 从特征值中计算奇异值,奇异值的平方根即为特征值的平方根。

  8. 根据奇异值构造对角矩阵Σ。

最后,将矩阵A分解为A = UΣV^T。

二、Python实现

1. 调包np.linalg.svd()

在Python中,可以使用NumPy库来实现SVD分解。以下是一个示例代码:

  1. import numpy as np
  2. # 定义一个矩阵A
  3. A = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
  4. # 使用NumPy的svd函数进行SVD分解
  5. U, s, Vt = np.linalg.svd(A)
  6. # 打印分解结果
  7. print("矩阵U:")
  8. print(U)
  9. print("奇异值:")
  10. print(s)
  11. print("矩阵V的转置:")
  12. print(Vt)

        运行上述代码,将得到矩阵A的SVD分解结果。其中,U是一个正交矩阵,s是包含矩阵A的奇异值的一维数组,Vt是V的转置矩阵。

2. 自定义

  1. import numpy as np
  2. def svd_decomposition(A):
  3. # 计算 A 的转置与 A 的乘积
  4. AAT = np.dot(A, A.T)
  5. # 计算 A 的乘积与 A 的转置的乘积
  6. ATA = np.dot(A.T, A)
  7. # 计算 A 的转置与 A 的乘积的特征值和特征向量
  8. eigenvalues_U, eigenvectors_U = np.linalg.eig(AAT)
  9. # 计算 A 的乘积与 A 的转置的特征值和特征向量
  10. eigenvalues_V, eigenvectors_V = np.linalg.eig(ATA)
  11. # 对特征值进行排序,并获取排序索引
  12. sorted_indices_U = np.argsort(eigenvalues_U)[::-1]
  13. sorted_indices_V = np.argsort(eigenvalues_V)[::-1]
  14. # 获取奇异值
  15. singular_values = np.sqrt(np.sort(eigenvalues_U)[::-1])
  16. # 获取 U 矩阵
  17. U = eigenvectors_U[:, sorted_indices_U]
  18. # 获取 V 矩阵
  19. V = eigenvectors_V[:, sorted_indices_V]
  20. # 对 V 进行转置
  21. V = V.T
  22. return U, singular_values, V
  23. # 定义一个矩阵 A
  24. A = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])
  25. # 调用自定义的 SVD 分解函数
  26. U, s, Vt = svd_decomposition(A)
  27. # 打印分解结果
  28. print("矩阵 U:")
  29. print(U)
  30. print("奇异值:")
  31. print(s)
  32. print("矩阵 V 的转置:")
  33. print(Vt)

        自定义的 svd_decomposition 函数基于奇异值分解算法,通过计算矩阵 A 的转置与 A 的乘积以及 A 的乘积与 A 的转置的特征值和特征向量来实现 SVD 分解。

三、SVD实现链路预测

  1. import numpy as np
  2. # 定义邻接矩阵A,表示网络结构
  3. A = np.array([[0, 1, 0, 1],
  4. [1, 0, 1, 0],
  5. [0, 1, 0, 1],
  6. [1, 0, 1, 0]])
  7. # 进行奇异值分解
  8. U, S, V = np.linalg.svd(A)
  9. # 选择保留的前k个奇异值和对应的奇异向量
  10. k = 2 # 选择保留的奇异值个数
  11. U_k = U[:, :k]
  12. S_k = np.diag(S[:k])
  13. V_k = V[:k, :]
  14. # 进行链路预测
  15. A_pred = U_k @ S_k @ V_k
  16. # 保留两位小数
  17. A_pred = np.round(A_pred, decimals=2)
  18. # 输出链路预测结果
  19. print("链路预测结果:")
  20. print(A_pred)

        首先定义了一个邻接矩阵 A,然后使用 np.linalg.svd 函数进行奇异值分解。然后,我们选择保留前 k 个奇异值和对应的奇异向量(在示例中选择 k=2),并重新构造预测的邻接矩阵 A_pred。最后,输出链路预测的结果。

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/繁依Fanyi0/article/detail/713902
推荐阅读
相关标签
  

闽ICP备14008679号