每次我 运行 代码时，情节看起来都不一样

Question

我的代码有问题。每次我 运行 它看起来都不一样。有任何想法吗？我看不出有什么问题。我从 2 小时开始查看这段代码，但找不到问题...

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats


x = np.arange(0,24, 1)
y = stats.poisson.pmf(x, mu=13)
a =stats.poisson.rvs(mu=13, size=5000)

#plt.stem(a,x)
plt.hist(a,bins=x,density=True,edgecolor="red")
plt.title("Poisson Verteilung mit Erwartungswert 13")
plt.xlabel("Anzahl M.")
plt.ylabel("Wahrscheinlichkeit")
"""
i=10
o=30
while i != o:
  if y[i]*100<0.5:
    #print(i)
    break
  i+=1
"""

#plot to specific x value
plt.xlim(0, 23)

plt.plot()
plt.plot(13, y[13], marker='x', markersize=5, color="black",label="Varianz")
plt.plot(13, y[13], marker='x', markersize=5, color="black",label="Median")
plt.legend(loc="upper left")

Answer 1

您正在使用 random statistical distribution，每次您 运行 您的代码都会创建一个新的随机分布。
为了获得可重复性（当你 运行 它时，你的代码总是选择相同的随机分布）你必须在代码的开头 set a seed:

np.random.seed(42)

你可以选择你喜欢的种子，42是常见的。

完整代码

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats


np.random.seed(42)
x = np.arange(0,24, 1)
y = stats.poisson.pmf(x, mu=13)
a =stats.poisson.rvs(mu=13, size=5000)

#plt.stem(a,x)
plt.hist(a,bins=x,density=True,edgecolor="red")
plt.title("Poisson Verteilung mit Erwartungswert 13")
plt.xlabel("Anzahl M.")
plt.ylabel("Wahrscheinlichkeit")
"""
i=10
o=30
while i != o:
  if y[i]*100<0.5:
    #print(i)
    break
  i+=1
"""

#plot to specific x value
plt.xlim(0, 23)

plt.plot()
plt.plot(13, y[13], marker='x', markersize=5, color="black",label="Varianz")
plt.plot(13, y[13], marker='x', markersize=5, color="black",label="Median")
plt.legend(loc="upper left")

plt.show()

情节