Shortcuts

torch.linalg.cond

torch.linalg.cond(A, p=None, *, out=None) 张量

计算矩阵相对于矩阵范数的条件数。

K\mathbb{K}R\mathbb{R}C\mathbb{C}条件数 κ\kappa 的矩阵 AKn×nA \in \mathbb{K}^{n \times n} 定义为

κ(A)=ApA1p\kappa(A) = \|A\|_p\|A^{-1}\|_p

矩阵 A 的条件数衡量了线性系统 AX = B 相对于矩阵范数的数值稳定性。

支持输入 float、double、cfloat 和 cdouble 数据类型。 还支持矩阵的批处理,如果 A 是矩阵的批处理,则输出具有相同的批处理维度。

p 定义了所计算的矩阵范数。支持以下范数:

p

矩阵范数

2-范数(最大奇异值)

‘fro’

Frobenius 范数

‘nuc’

核范数

无穷大

max(sum(abs(x), dim=1))

-无穷大

min(sum(abs(x), dim=1))

1

max(sum(abs(x), dim=0))

-1

min(sum(abs(x), dim=0))

2

最大奇异值

-2

最小奇异值

其中 inf 指的是 float(‘inf’)、NumPy 的 inf 对象,或任何等效的对象。

对于 p(‘fro’, ‘nuc’, inf, -inf, 1, -1) 之一,此函数使用 torch.linalg.norm()torch.linalg.inv()。 因此,在这种情况下,矩阵(或批次中的每个矩阵)A 必须是方阵且可逆。

对于 p(2, -2) 中,此函数可以根据奇异值 σ1σn\sigma_1 \geq \ldots \geq \sigma_n

κ2(A)=σ1σnκ2(A)=σnσ1\kappa_2(A) = \frac{\sigma_1}{\sigma_n}\qquad \kappa_{-2}(A) = \frac{\sigma_n}{\sigma_1}

在这些情况下,它是使用 torch.linalg.svdvals() 计算的。对于这些范数,矩阵(或批次中的每个矩阵)A 可以具有任意形状。

注意

当输入在CUDA设备上时,如果p(‘fro’, ‘nuc’, inf, -inf, 1, -1)之一,此函数会同步该设备与CPU。

另请参阅

torch.linalg.solve() 用于解决方阵线性系统的函数。

torch.linalg.lstsq() 用于解决一般矩阵的线性系统函数。

Parameters
  • A (张量) – 形状为 (*, m, n) 的张量,其中 * 表示零个或多个批次维度, 对于 p(2, -2) 中,并且形状为 (*, n, n) 的张量,其中每个矩阵 对于 p(‘fro’, ‘nuc’, inf, -inf, 1, -1) 中都是可逆的。

  • p (int, inf, -inf, 'fro', 'nuc', 可选) – 计算中使用的矩阵范数类型(见上文)。默认值: None

Keyword Arguments

输出 (张量, 可选) – 输出张量。如果为,则忽略。默认值:

Returns

一个实值张量,即使 A 是复数。

Raises

RuntimeError – 如果 p(‘fro’, ‘nuc’, inf, -inf, 1, -1) 之一, 并且 A 矩阵或批次中的任何矩阵 A 不是方阵 或不可逆。

示例:

>>> A = torch.randn(3, 4, 4, dtype=torch.complex64)
>>> torch.linalg.cond(A)
>>> A = torch.tensor([[1., 0, -1], [0, 1, 0], [1, 0, 1]])
>>> torch.linalg.cond(A)
tensor([1.4142])
>>> torch.linalg.cond(A, 'fro')
tensor(3.1623)
>>> torch.linalg.cond(A, 'nuc')
tensor(9.2426)
>>> torch.linalg.cond(A, float('inf'))
tensor(2.)
>>> torch.linalg.cond(A, float('-inf'))
tensor(1.)
>>> torch.linalg.cond(A, 1)
tensor(2.)
>>> torch.linalg.cond(A, -1)
tensor(1.)
>>> torch.linalg.cond(A, 2)
tensor([1.4142])
>>> torch.linalg.cond(A, -2)
tensor([0.7071])

>>> A = torch.randn(2, 3, 3)
>>> torch.linalg.cond(A)
tensor([[9.5917],
        [3.2538]])
>>> A = torch.randn(2, 3, 3, dtype=torch.complex64)
>>> torch.linalg.cond(A)
tensor([[4.6245],
        [4.5671]])