numpy.random.noncentral_chisquare#
- random.noncentral_chisquare(df, nonc, size=None)#
从非中心卡方分布中抽取样本。
非中心 \(\chi^2\) 分布是 \(\chi^2\) 分布的推广。
注意
新代码应使用
Generator
实例的noncentral_chisquare
方法;请参阅快速入门。- 参数:
- df浮点数或浮点数数组
自由度,必须 > 0。
- nonc浮点数或浮点数数组
非中心参数,必须是非负数。
- sizeint 或 int 元组,可选
输出形状。如果给定形状为例如
(m, n, k)
,则抽取m * n * k
个样本。如果 size 为None
(默认),则当df
和nonc
均为标量时返回单个值。否则,抽取np.broadcast(df, nonc).size
个样本。
- 返回:
- outndarray 或标量
从参数化的非中心卡方分布中抽取的样本。
另请参阅
random.Generator.noncentral_chisquare
新代码应使用此方法。
注释
非中心卡方分布的概率密度函数为
\[P(x;df,nonc) = \sum^{\infty}_{i=0} \frac{e^{-nonc/2}(nonc/2)^{i}}{i!} P_{Y_{df+2i}}(x),\]其中 \(Y_{q}\) 是具有 q 个自由度的卡方分布。
参考
[1]维基百科,“非中心卡方分布” https://en.wikipedia.org/wiki/Noncentral_chi-squared_distribution
示例
从分布中抽取值并绘制直方图
>>> import matplotlib.pyplot as plt >>> values = plt.hist(np.random.noncentral_chisquare(3, 20, 100000), ... bins=200, density=True) >>> plt.show()
从非中心参数非常小的非中心卡方分布中抽取值,并与卡方分布进行比较。
>>> plt.figure() >>> values = plt.hist(np.random.noncentral_chisquare(3, .0000001, 100000), ... bins=np.arange(0., 25, .1), density=True) >>> values2 = plt.hist(np.random.chisquare(3, 100000), ... bins=np.arange(0., 25, .1), density=True) >>> plt.plot(values[1][0:-1], values[0]-values2[0], 'ob') >>> plt.show()
演示非中心参数的大值如何导致更对称的分布。
>>> plt.figure() >>> values = plt.hist(np.random.noncentral_chisquare(3, 20, 100000), ... bins=200, density=True) >>> plt.show()