numpy.linalg.norm#

linalg.norm(x, ord=None, axis=None, keepdims=False)[source]#

矩阵或向量的范数。

此函数能够返回八种不同的矩阵范数之一,或无限多种向量范数之一(如下所述),具体取决于 ord 参数的值。

参数:
xarray_like

输入数组。如果 axis 为 None,除非 ord 为 None,否则 x 必须是一维或二维的。如果 axisord 都为 None,则返回 x.ravel 的 2-范数。

ord{int, float, inf, -inf, ‘fro’, ‘nuc’}, 可选

范数的阶数(有关矩阵和向量分别支持的值,请参见“注释”下的表格)。inf 表示 NumPy 的 inf 对象。默认为 None。

axis{None, int, 2-元组 ints}, 可选。

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

keepdimsbool, 可选

如果将其设置为 True,则归一化的轴将保留在结果中,大小为 1。使用此选项,结果将与原始 x 正确广播。

返回值:
nfloat 或 ndarray

矩阵或向量(s) 的范数。

参见

scipy.linalg.norm

SciPy 中的类似函数。

注释

对于 ord < 1 的值,结果严格来说不是数学上的“范数”,但它对于各种数值用途仍然可能有用。

可以计算以下范数

ord

矩阵范数

向量范数

None

Frobenius 范数

2-范数

‘fro’

Frobenius 范数

‘nuc’

核范数

inf

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

max(abs(x))

-inf

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

min(abs(x))

0

sum(x != 0)

1

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

如下

-1

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

如下

2

2-范数(最大奇异值)

如下

-2

最小奇异值

如下

其他

sum(abs(x)**ord)**(1./ord)

Frobenius 范数由 [1] 给出

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

核范数是奇异值的总和。

Frobenius 范数和核范数阶数仅针对矩阵定义,当 x.ndim != 2 时会引发 ValueError。

参考文献

[1]

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

示例

>>> import numpy as np
>>> from numpy import linalg as LA
>>> a = np.arange(9) - 4
>>> a
array([-4, -3, -2, ...,  2,  3,  4])
>>> b = a.reshape((3, 3))
>>> b
array([[-4, -3, -2],
       [-1,  0,  1],
       [ 2,  3,  4]])
>>> LA.norm(a)
7.745966692414834
>>> LA.norm(b)
7.745966692414834
>>> LA.norm(b, 'fro')
7.745966692414834
>>> LA.norm(a, np.inf)
4.0
>>> LA.norm(b, np.inf)
9.0
>>> LA.norm(a, -np.inf)
0.0
>>> LA.norm(b, -np.inf)
2.0
>>> LA.norm(a, 1)
20.0
>>> LA.norm(b, 1)
7.0
>>> LA.norm(a, -1)
-4.6566128774142013e-010
>>> LA.norm(b, -1)
6.0
>>> LA.norm(a, 2)
7.745966692414834
>>> LA.norm(b, 2)
7.3484692283495345
>>> LA.norm(a, -2)
0.0
>>> LA.norm(b, -2)
1.8570331885190563e-016 # may vary
>>> LA.norm(a, 3)
5.8480354764257312 # may vary
>>> LA.norm(a, -3)
0.0

使用 axis 参数计算向量范数

>>> c = np.array([[ 1, 2, 3],
...               [-1, 1, 4]])
>>> LA.norm(c, axis=0)
array([ 1.41421356,  2.23606798,  5.        ])
>>> LA.norm(c, axis=1)
array([ 3.74165739,  4.24264069])
>>> LA.norm(c, ord=1, axis=1)
array([ 6.,  6.])

使用 axis 参数计算矩阵范数

>>> m = np.arange(8).reshape(2,2,2)
>>> LA.norm(m, axis=(1,2))
array([  3.74165739,  11.22497216])
>>> LA.norm(m[0, :, :]), LA.norm(m[1, :, :])
(3.7416573867739413, 11.224972160321824)