norm#
- scipy.linalg.norm(a, ord=None, axis=None, keepdims=False, check_finite=True)[源代码]#
矩阵或向量范数。
此函数能够返回八种不同的矩阵范数之一,或者无限多种向量范数之一(如下所述),具体取决于
ord
参数的值。 对于秩不同于 1 或 2 的张量,仅支持 ord=None。- 参数:
- a类数组
输入数组。如果 axis 为 None,则 a 必须为 1-D 或 2-D,除非 ord 为 None。如果 axis 和 ord 均为 None,则返回
a.ravel
的 2-范数。- ord{int, inf, -inf, ‘fro’, ‘nuc’, None}, 可选
范数的阶数(请参见
Notes
下的表格)。inf 表示 NumPy 的 inf 对象。- axis{int, 整数的 2 元组, None}, 可选
如果 axis 是一个整数,则指定 a 中计算向量范数的轴。如果 axis 是一个 2 元组,则指定包含 2-D 矩阵的轴,并计算这些矩阵的矩阵范数。如果 axis 为 None,则返回向量范数(当 a 为 1-D 时)或矩阵范数(当 a 为 2-D 时)。
- keepdimsbool, 可选
如果将其设置为 True,则在结果中保留进行范数计算的轴,其大小为 1。 使用此选项,结果将正确地广播到原始的 a。
- check_finitebool, 可选
是否检查输入矩阵是否仅包含有限数字。禁用可能会提高性能,但如果输入包含无穷大或 NaN,可能会导致问题(崩溃、无法终止)。
- 返回:
- nfloat 或 ndarray
矩阵或向量的范数。
备注
对于
ord <= 0
的值,严格来说,结果不是数学上的“范数”,但它可能仍然对各种数值目的有用。可以计算以下范数
ord
矩阵的范数
向量的范数
None
Frobenius 范数
2-范数
‘fro’
Frobenius 范数
–
‘nuc’
核范数
–
inf
max(sum(abs(a), axis=1))
max(abs(a))
-inf
min(sum(abs(a), axis=1))
min(abs(a))
0
–
sum(a != 0)
1
max(sum(abs(a), axis=0))
如下
-1
min(sum(abs(a), axis=0))
如下
2
2-范数(最大奇异值)
如下
-2
最小奇异值
如下
其他
–
sum(abs(a)**ord)**(1./ord)
Frobenius 范数由 [1] 给出
\(||A||_F = [\sum_{i,j} abs(a_{i,j})^2]^{1/2}\)
核范数是奇异值之和。
Frobenius 和核范数阶数仅为矩阵定义。
参考文献
[1]G. H. Golub 和 C. F. Van Loan,《矩阵计算》,巴尔的摩,马里兰州,约翰霍普金斯大学出版社,1985 年,第 15 页
示例
>>> import numpy as np >>> from scipy.linalg import norm >>> a = np.arange(9) - 4.0 >>> 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.]])
>>> norm(a) 7.745966692414834 >>> norm(b) 7.745966692414834 >>> norm(b, 'fro') 7.745966692414834 >>> norm(a, np.inf) 4.0 >>> norm(b, np.inf) 9.0 >>> norm(a, -np.inf) 0.0 >>> norm(b, -np.inf) 2.0
>>> norm(a, 1) 20.0 >>> norm(b, 1) 7.0 >>> norm(a, -1) -4.6566128774142013e-010 >>> norm(b, -1) 6.0 >>> norm(a, 2) 7.745966692414834 >>> norm(b, 2) 7.3484692283495345
>>> norm(a, -2) 0.0 >>> norm(b, -2) 1.8570331885190563e-016 >>> norm(a, 3) 5.8480354764257312 >>> norm(a, -3) 0.0