SciPy 的哪个线性函数用于求解 Hermitian 正定带状矩阵方程?
名为scipy.linalg.solveh_banded的线性函数用于求解Hermitian正定带状矩阵方程。此函数的形式如下-
scipy.linalg.solveh_banded(ab, b, overwrite_ab=False, overwrite_b=False, lower=False, check_finite=True)
这个线性函数将求解方程ax=bforx,其中a是Hermitian正定带状矩阵。
带状矩阵a以下或上对角线排序形式存储在ab中,如下所示-
ab[u + i - j, j] == a[i,j] (if upper form; i<=j) ab[ i - j, j] == a[i,j] (if lower form; i >= j)
上表中ab的示例如下-
* * a02 a13 a24 a35 * a01 a12 a23 a34 a45 a00 a11 a22 a33 a44 a55
这里a的形状是(6,6)并且u=2。
下面给出了ab的例子如下:
a00 a11 a22 a33 a44 a55 a10 a21 a32 a43 a54 * a20 a31 a42 a53 * *
这里a的形状是(6,6)并且u=2。
参数
下面给出函数scipy的参数。-linalg.solveh_banded()
ab-(u+1,M)array_like
该参数表示带状矩阵。
b-(M,)或(M,K)array_like
此参数代表右侧。
Overwrite_ab-bool,可选
该参数用于丢弃带状矩阵ab中的数据。它可以提高矩阵的性能。
Overwrite_b-布尔值,可选
该参数用于丢弃b中的数据。它可以提高矩阵的性能
下-布尔值,可选
此参数用于检查带状矩阵是否为下式。默认是带状矩阵的上形式。
check_finite-布尔值,可选
此参数用于检查输入矩阵是否仅包含有限数。禁用它后,我们可能会获得一些性能提升。如果输入不包含无穷大,则可能会导致问题。
退货
x-(M,)或(M,K)ndarray
它返回带状矩阵方程ax=b的解。输出的形状将取决于b的形状。