scipy.linalg.cho_solve_banded¶
- scipy.linalg.cho_solve_banded(cb_and_lower, b, overwrite_b=False, check_finite=True)[源代码]¶
解线性方程组
A x = b
,给定带状厄米特的Cholesky因式分解A
。- 参数
- (Cb,更低)元组,(ndarray,bool)
cb 是A的Cholesky因子分解,如Cholesky_Bandded所给出的。 lower 必须与赋予cholesky_banded的值相同。
- barray_like
右手边
- overwrite_b布尔值,可选
如果为True,则函数将覆盖中的值 b 。
- check_finite布尔值,可选
是否检查输入矩阵是否仅包含有限个数字。禁用可能会带来性能提升,但如果输入确实包含无穷大或NAN,则可能会导致问题(崩溃、非终止)。
- 退货
- x阵列
系统Ax=b的解
参见
cholesky_banded
带状矩阵的Cholesky分解
注意事项
0.8.0 新版功能.
示例
>>> 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