numpy.std#

numpy.std(a, axis=None, dtype=None, out=None, ddof=0, keepdims=<no value>, *, where=<no value>, mean=<no value>, correction=<no value>)[源代码]#

沿指定轴计算标准差。

返回数组元素的标准差，即分布的离散程度的度量。默认情况下，标准差是为展平的数组计算的，否则则沿指定的轴计算。

参数:
a类数组对象

计算这些值的标准差。

axisNone 或 int 或 int 元组，可选

计算标准差的轴。默认是计算展平数组的标准差。如果这是一个整数元组，则对多个轴执行标准差计算，而不是像以前那样对单个轴或所有轴进行计算。

dtypedtype, optional

用于计算标准差的类型。对于整数类型的数组，默认值为float64；对于浮点数类型的数组，它与数组类型相同。

outndarray，可选

用于放置结果的替代输出数组。它必须具有与预期输出相同的形状，但类型（计算出的值）将在必要时进行转换。有关更多详细信息，请参阅 输出类型确定。

ddof{int, float}，可选

Means Delta Degrees of Freedom。计算中使用的除数是 N - ddof，其中 N 表示元素的数量。默认情况下 ddof 为零。有关 ddof 用法的详细信息，请参阅说明。

keepdimsbool，可选

如果设置为 True，则被缩减的轴将保留在结果中，其大小为 1。使用此选项，结果将能够与输入数组正确广播。

如果传递默认值，则 keepdims 不会传递给 std 方法的子类 ndarray，但是任何非默认值都会传递。如果子类的方法未实现 keepdims，则将引发任何异常。

wherearray_like of bool，可选

包含在标准差中的元素。有关详细信息，请参阅 reduce。

版本 1.20.0 中新增。

meanarray_like, optional

提供平均值以防止其重新计算。平均值的形状应如同使用 keepdims=True 计算一样。计算平均值的轴应与调用此 std 函数时使用的轴相同。

版本 2.0.0 中新增。

correction{int, float}，可选

ddof 参数的 Array API 兼容名称。一次只能提供其中一个。

版本 2.0.0 中新增。

返回:
standard_deviationndarray, 请参阅上面的 dtype 参数。

如果 out 为 None，则返回一个包含标准差的新数组；否则，返回输出数组的引用。

另请参阅

var, mean, nanmean, nanstd, nanvar
输出类型确定

备注

数组标准差计算有几种常见的变体。假设输入 a 是一个一维 NumPy 数组，并且 mean 要么作为参数提供，要么计算为 a.mean()，NumPy 将数组的标准差计算为

N = len(a)
d2 = abs(a - mean)**2  # abs is for complex `a`
var = d2.sum() / (N - ddof)  # note use of `ddof`
std = var**0.5

参数 ddof 的不同值在不同情况下很有用。NumPy 的默认 ddof=0 对应于以下表达式：

\[\sqrt{\frac{\sum_i{|a_i - \bar{a}|^2 }}{N}}\]

这在统计学领域有时被称为“总体标准差”，因为它将标准差的定义应用于 a，就好像 a 是所有可能观测值的完整总体一样。

许多其他库以不同的方式定义数组的标准差，例如

\[\sqrt{\frac{\sum_i{|a_i - \bar{a}|^2 }}{N - 1}}\]

在统计学中，所得数量有时被称为“样本标准差”，因为如果 a 是来自更大总体的随机样本，则此计算提供了总体方差的无偏估计的平方根。分母中使用 \(N-1\) 通常称为“贝塞尔校正”，因为它校正了当 a 的样本均值用于替代总体的真实均值时引入的方差估计的偏差（趋向于较低值）。所得的标准差估计仍然有偏差，但偏差比没有校正时要小。对于此数量，请使用 ddof=1。

请注意，对于复数，std 在平方之前取绝对值，因此结果始终是实数且非负。

对于浮点输入，标准差的计算精度与输入相同。根据输入数据，这可能会导致结果不准确，尤其对于 float32（请参阅下面的示例）。使用 dtype 关键字指定更高的精度累加器可以缓解此问题。

示例

>>> import numpy as np
>>> a = np.array([[1, 2], [3, 4]])
>>> np.std(a)
1.1180339887498949 # may vary
>>> np.std(a, axis=0)
array([1.,  1.])
>>> np.std(a, axis=1)
array([0.5,  0.5])

单精度下的 std() 可能不准确

>>> a = np.zeros((2, 512*512), dtype=np.float32)
>>> a[0, :] = 1.0
>>> a[1, :] = 0.1
>>> np.std(a)
np.float32(0.45000005)

使用 float64 计算标准差更准确

>>> np.std(a, dtype=np.float64)
0.44999999925494177 # may vary

指定 where 参数

>>> a = np.array([[14, 8, 11, 10], [7, 9, 10, 11], [10, 15, 5, 10]])
>>> np.std(a)
2.614064523559687 # may vary
>>> np.std(a, where=[[True], [True], [False]])
2.0

使用 mean 关键字以节省计算时间

>>> import numpy as np
>>> from timeit import timeit
>>> a = np.array([[14, 8, 11, 10], [7, 9, 10, 11], [10, 15, 5, 10]])
>>> mean = np.mean(a, axis=1, keepdims=True)
>>>
>>> g = globals()
>>> n = 10000
>>> t1 = timeit("std = np.std(a, axis=1, mean=mean)", globals=g, number=n)
>>> t2 = timeit("std = np.std(a, axis=1)", globals=g, number=n)
>>> print(f'Percentage execution time saved {100*(t2-t1)/t2:.0f}%')

Percentage execution time saved 30%
在此页面