用cholesky求解Ax=b。用python自带的函数

可以使用numpy.linalg.cholesky函数来求解Ax=b中的x，其中A是一个正定对称矩阵。以下是一个示例代码：

import numpy as np

# 构造一个正定对称矩阵A和一个向量b
A = np.array([[4, 2, 1], [2, 5, 3], [1, 3, 6]])
b = np.array([1, 2, 3])

# 使用cholesky分解求解Ax=b中的x
L = np.linalg.cholesky(A)
y = np.linalg.solve(L, b)
x = np.linalg.solve(L.T, y)

print(x)  # 输出解向量x

请注意，这种方法只适用于正定对称矩阵A，否则会出现错误。