scipy.sparse.linalg.

norm#

scipy.sparse.linalg.norm(x, ord=None, axis=None)[源代码]#

稀疏矩阵的范数

此函数能够返回七种不同的矩阵范数之一,具体取决于 ord 参数的值。

参数:
x一个稀疏数组

输入的稀疏数组。

ord{非零整数, inf, -inf, ‘fro’}, 可选

范数的阶数(请参阅 Notes 下的表格)。inf 表示 numpy 的 inf 对象。

axis{整数, 整数的 2 元组, None}, 可选

如果 axis 是一个整数,则它指定计算向量范数的 x 的轴。如果 axis 是一个 2 元组,则它指定包含 2-D 矩阵的轴,并计算这些矩阵的矩阵范数。如果 axis 为 None,则返回向量范数(当 x 为 1-D 时)或矩阵范数(当 x 为 2-D 时)。

返回:
n浮点数或 ndarray

说明

一些 ord 没有实现,因为某些相关函数(例如 _multi_svd_norm)对于稀疏数组尚不可用。

此文档字符串基于 numpy.linalg.norm 修改。numpy/numpy

可以计算以下范数

ord

稀疏数组的范数

None

Frobenius 范数

‘fro’

Frobenius 范数

inf

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

-inf

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

0

abs(x).sum(axis=axis)

1

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

-1

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

2

谱范数(最大的奇异值)

-2

未实现

其他

未实现

Frobenius 范数由 [1] 给出

\(||A||_F = [\sum_{i,j} abs(a_{i,j})^2]^{1/2}\)

参考文献

[1]

G. H. Golub 和 C. F. Van Loan,《矩阵计算》,巴尔的摩,马里兰州,约翰·霍普金斯大学出版社,1985 年,第 15 页

示例

>>> from scipy.sparse import csr_array, diags_array
>>> import numpy as np
>>> from scipy.sparse.linalg import norm
>>> a = np.arange(9) - 4
>>> a
array([-4, -3, -2, -1, 0, 1, 2, 3, 4])
>>> b = a.reshape((3, 3))
>>> b
array([[-4, -3, -2],
       [-1, 0, 1],
       [ 2, 3, 4]])
>>> b = csr_array(b)
>>> norm(b)
7.745966692414834
>>> norm(b, 'fro')
7.745966692414834
>>> norm(b, np.inf)
9
>>> norm(b, -np.inf)
2
>>> norm(b, 1)
7
>>> norm(b, -1)
6

矩阵 2-范数或谱范数是最大的奇异值,它是近似计算的,并且存在局限性。

>>> b = diags_array([-1, 1], [0, 1], shape=(9, 10))
>>> norm(b, 2)
1.9753...