scipy.linalg.lu_solve

scipy.linalg.lu_solve(lu_and_piv, b, trans=0, overwrite_b=False, check_finite=True)[源代码]

在给定a的LU因式分解的情况下,求解方程组a x=b

参数
(Lu,PIV)

系数矩阵a的因式分解,由lu_factor给出

b阵列

右手边

trans{0,1,2},可选

要解决的系统类型:

运输

系统

0

a x=b

1

a^T x=b

2

a^H x=b

overwrite_b布尔值,可选

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

check_finite布尔值,可选

是否检查输入矩阵是否仅包含有限个数字。禁用可能会带来性能提升,但如果输入确实包含无穷大或NAN,则可能会导致问题(崩溃、非终止)。

退货
x阵列

系统解决方案

参见

lu_factor

Lu分解矩阵

示例

>>> from scipy.linalg import lu_factor, lu_solve
>>> A = np.array([[2, 5, 8, 7], [5, 2, 2, 8], [7, 5, 6, 6], [5, 4, 4, 8]])
>>> b = np.array([1, 1, 1, 1])
>>> lu, piv = lu_factor(A)
>>> x = lu_solve((lu, piv), b)
>>> np.allclose(A @ x - b, np.zeros((4,)))
True