scipy.linalg.
cho_solve_banded#
- scipy.linalg.cho_solve_banded(cb_and_lower, b, overwrite_b=False, check_finite=True)[源代码][源代码]#
求解线性方程
A x = b
,已知带状厄米特矩阵A
的 Cholesky 分解。- 参数:
- (cb, lower)tuple, (ndarray, bool)
cb 是 A 的 Cholesky 分解,由 cholesky_banded 给出。lower 必须与传递给 cholesky_banded 的值相同。
- barray_like
右侧
- overwrite_bbool, 可选
如果为真,函数将覆盖 b 中的值。
- check_finitebool, 可选
是否检查输入矩阵是否仅包含有限数值。禁用可能会提高性能,但如果输入包含无穷大或NaN,可能会导致问题(崩溃、无法终止)。
- 返回:
- x数组
系统 A x = b 的解
参见
cholesky_banded
带状矩阵的Cholesky分解
注释
Added in version 0.8.0.
示例
>>> import numpy as np >>> from scipy.linalg import cholesky_banded, cho_solve_banded >>> Ab = np.array([[0, 0, 1j, 2, 3j], [0, -1, -2, 3, 4], [9, 8, 7, 6, 9]]) >>> A = np.diag(Ab[0,2:], k=2) + np.diag(Ab[1,1:], k=1) >>> A = A + A.conj().T + np.diag(Ab[2, :]) >>> c = cholesky_banded(Ab) >>> x = cho_solve_banded((c, False), np.ones(5)) >>> np.allclose(A @ x - np.ones(5), np.zeros(5)) True