scipy.linalg.

cho_solve#

scipy.linalg.cho_solve(c_and_lower, b, overwrite_b=False, check_finite=True)[源代码][源代码]#

给定矩阵A的Cholesky分解,求解线性方程 A x = b。

参数:
(c, 小写)tuple, (array, bool)

Cholesky 分解 a,由 cho_factor 给出

b数组

右侧

overwrite_bbool, 可选

是否覆盖b中的数据(可能会提高性能)

check_finitebool, 可选

是否检查输入矩阵是否仅包含有限数值。禁用可能会提高性能,但如果输入包含无穷大或NaN,可能会导致问题(崩溃、无法终止)。

返回:
x数组

系统 A x = b 的解

参见

cho_factor

矩阵的Cholesky分解

示例

>>> import numpy as np
>>> from scipy.linalg import cho_factor, cho_solve
>>> A = np.array([[9, 3, 1, 5], [3, 7, 5, 1], [1, 5, 9, 2], [5, 1, 2, 6]])
>>> c, low = cho_factor(A)
>>> x = cho_solve((c, low), [1, 1, 1, 1])
>>> np.allclose(A @ x - [1, 1, 1, 1], np.zeros(4))
True