Python 进行 Cholesky分解

Posted 2021-07-31 yangbocsu

Python 进行 Cholesky分解

---

#yangbocsu 2021.05.24  科技园7栋
import numpy as np
from scipy import linalg
A = np.array([[9, 6, 3],
         [6, 13, 11],
         [3, 11, 35]])
L = linalg.cholesky(A, lower=True) # 默认计算 upper， 所以指定 lower = True
U = linalg.cholesky(A)

print("L = \\n",L)
print("\\nU = \\n",U)

# 验证  np.allclose()  print(np.dot(L, L.T))