我发现 scipy 的 weibull_min pdf 函数不直观且难以使用。我想生成一个规模= 30和形状= 2.5的weibull PDF。这应该是这样的:
这是我迄今为止尝试过的:
import matplotlib.pyplot as plt
from scipy.stats import weibull_max, uniform
import numpy as np
aoas = np.linspace(0, 8, 1000)
speeds = np.linspace(1, 80, 1000)
plt.fill_between(speeds, speeds * 0, weibull_max.pdf(speeds, 2.5, 30), facecolor='k')
plt.ylabel('Probability Density')
plt.savefig('speedDist.pdf')
plt.clf()
当然，我的代码是错误的。文档称 weibull_min 需要 pdf(x, c, loc=0, scale=1)，其中 x 是分位数， c 是形状因子，scale 是比例因子。但是，当我将代码更改为 weibull_max.pdf(speeds, 2.5, 0, 30) 时，输出 pdf 仅由零值组成。所以，我很困惑。如何生成所需的威 bool PDF？
最佳答案
您的代码中有两个问题。
weibull_min和weibull_max不是同一分布。您在问题的标题和描述中引用了weibull_min，但在代码中使用了weibull_max。要匹配您显示的绘图，请使用weibull_min。pdf方法的参数为(x, shape, loc, scale)。您编写了pdf(speeds, 2.5, 30)，它将 30 分配给loc参数。相反，您应该使用pdf(speeds, 2.5, scale=30)
这是您的脚本的变体:
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import weibull_min
speeds = np.linspace(0, 80, 1000)
p = weibull_min.pdf(speeds, 2.5, scale=30)
plt.plot(speeds, p, 'b', linewidth=1)
plt.fill_between(speeds, speeds * 0, p, facecolor='b', alpha=0.1)
plt.ylabel('Probability Density')
plt.show()
它生成这个图:
关于python - scipy:如何使用 weibull_min.pdf？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/53459599/