In [1]:
import numpy as np
import matplotlib.pyplot as plt
ndata - počet hodnot, které budeme generovat nbis - počet binů histogramu v poli x budou hrancie binů histogramu v poli y budou počty hodnot, které padly do daného binu
In [3]:
ndata=1000
nbins=100
x=np.linspace(0,1,nbins)
y=np.zeros(nbins) #nejdřív sem dáme samé nuly
vygenerujeme ndata hodnot z rovnoměrného rozdělení (0,1) a vypočítámje bin histogramu ibin, do kterého daná hodnota patří
In [5]:
for i in range(ndata):
ibin=int(nbins*np.random.random())
y[ibin]=y[ibin]+1
normalizujeme histogram aby měl jednotkouvou plochu
In [7]:
y=y/(ndata)*nbins
nakreslíme graf
In [9]:
plt.step(x,y)
plt.plot([0,1],[1,1],c='red') #nakreslení červené čáry v y=1
plt.xlim(0,1) #nastavení limit na ose x
Out[9]:
(0.0, 1.0)
In [ ]: