import numpy as np
def SolverEqCholesky(matrix_A, vector_b, n):
# 创建一个n*1的零向量作为结果
x = np.zeros((n, 1))
# Cholesky分解
L = np.linalg.cholesky(matrix_A)
# 前向替换 Ly = b
y = np.linalg.solve(L, vector_b)
# 后向替换 L^T x = y
x = np.linalg.solve(L.T, y)
return x
Cholesky分解法-python答案
于 2023-08-23 17:44:10 首次发布