scipy.linalg.

eigvals_banded#

scipy.linalg.eigvals_banded(a_band, lower=False, overwrite_a_band=False, select='a', select_range=None, check_finite=True)[源代码][源代码]#

求解实对称或复Hermitian带矩阵的特征值问题。

找到矩阵 a 的特征值 w:

a v[:,i] = w[i] v[:,i]
v.H v    = identity

矩阵 a 存储在 a_band 中，可以是下对角线或上对角线有序形式：

a_band[u + i - j, j] == a[i,j] (如果是上三角形式; i <= j) a_band[ i - j, j] == a[i,j] (如果是下三角形式; i >= j)

其中 u 是主对角线上方的带数。

a_band 的示例（a 的形状为 (6,6)，u=2）:

upper form:
*   *   a02 a13 a24 a35
*   a01 a12 a23 a34 a45
a00 a11 a22 a33 a44 a55

lower form:
a00 a11 a22 a33 a44 a55
a10 a21 a32 a43 a54 *
a20 a31 a42 a53 *   *

标记为 * 的单元格未被使用。

参数:

a_band(u+1, M) array_like

M 乘 M 矩阵 a 的带状部分。

下限bool, 可选

矩阵是否为下三角形式。（默认是上三角形式）

覆盖一个波段bool, 可选

丢弃 a_band 中的数据（可能会提高性能）

选择{‘a’, ‘v’, ‘i’}, 可选

计算哪些特征值

选择	计算的
‘a’	所有特征值
‘v’	区间 (min, max] 中的特征值
‘我’	特征值与索引 min <= i <= max

选择范围(最小值, 最大值), 可选

选定特征值的范围

check_finitebool, 可选

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

返回:

w(M,) ndarray: 特征值，按升序排列，每个特征值根据其重数重复。

Raises:

LinAlgError: 如果特征值计算不收敛。

参见

eig_banded: 对称/厄米带矩阵的特征值和右特征向量
eigvalsh_tridiagonal: 对称/Hermitian 三对角矩阵的特征值
eigvals: 一般数组的特征值
eigh: 对称/厄米数组的特征值和右特征向量
eig: 非对称数组的特征值和右特征向量

示例

>>> import numpy as np
>>> from scipy.linalg import eigvals_banded
>>> A = np.array([[1, 5, 2, 0], [5, 2, 5, 2], [2, 5, 3, 5], [0, 2, 5, 4]])
>>> Ab = np.array([[1, 2, 3, 4], [5, 5, 5, 0], [2, 2, 0, 0]])
>>> w = eigvals_banded(Ab, lower=True)
>>> w
array([-4.26200532, -2.22987175,  3.95222349, 12.53965359])