概率分布（python实现）

二项分布

代码：

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

p=0.5
n=5
X=np.arange(0,n+1,1)
plist=stats.binom.pmf(X,n,p)
plt.bar(X,plist,width=0.1)
plt.xticks(X)
plt.show()

图像：

伯努利分布

代码：

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

n = 10
p = 0.3
k = np.arange(0, 21)
binomial = binom.pmf(k, n, p)
plt.plot(k, binomial, 'o-')
plt.show()

图像：

泊松分布

代码：

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

mu=5
k=7
X=np.arange(0,k+1,1)
plist=stats.poisson.pmf(X,mu)
plt.bar(X,plist,width=0.1)
plt.xticks(X)
plt.show()

图像：

几何分布

代码：

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

k=5
p=0.6
X=np.arange(1,k+1,1)
plist=stats.geom.pmf(X,p)
plt.bar(X,plist,width=0.1)
plt.xticks(X)
plt.show()

图像：

正态分布

代码：

import numpy as np
from numpy.linalg import cholesky
import matplotlib.pyplot as plt

sampleNo = 1000
mu = 3
sigma = 0.1
np.random.seed(0)
s = np.random.normal(mu, sigma, sampleNo )
plt.subplot(141)
plt.hist(s, 30, normed=True)

图像：