#抛10次硬币，求恰好两次正面朝上的概率
import numpy as np
from scipy import stats as sts
n=10
p=0.5
k=np.arange(0,11)   #总共有0-10次正面朝上的可能,arange其实是一个列表
binomial=sts.binom.pmf(k,n,p)
print('概率为：',binomial)  #输出的结果有11个，分别表示0-10次
print('概率的总和',sum(binomial))
print('两次正面朝上的概率是',binomial[2]) #竟然能用数组的形式

#泊松分布
rate=2 #λ的值
n=np.arange(0,11)
print('0-10次发生事故的概率分别是：',sts.poisson.pmf(n,rate))

#二项分布
k1=135
a=120
b=140
print('飞行135min的概率为：',sts.randint.pmf(k1,a,b))

#指数分布 sts.expon.pdf()
#正态分布 sts.norm.pdf(50,40,4) ,40是均值，4是方差，

##同样可求期望和方差