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))