Python实现的各种常见分布算法示例
作者:yan456jie 发布时间:2021-12-06 18:26:15
标签:Python,分布算法
本文实例讲述了Python实现的各种常见分布算法。分享给大家供大家参考,具体如下:
#-*- encoding:utf-8 -*-
import numpy as np
from scipy import stats
import matplotlib.pyplot as plt
#####################
#二项分布
#####################
def test_binom_pmf():
'''
为离散分布
二项分布的例子:抛掷10次硬币,恰好两次正面朝上的概率是多少?
'''
n = 10#独立实验次数
p = 0.5#每次正面朝上概率
k = np.arange(0,11)#0-10次正面朝上概率
binomial = stats.binom.pmf(k,n,p)
print binomial#概率和为1
print sum(binomial)
print binomial[2]
plt.plot(k, binomial,'o-')
plt.title('Binomial: n=%i , p=%.2f' % (n,p),fontsize=15)
plt.xlabel('Number of successes')
plt.ylabel('Probability of success',fontsize=15)
plt.show()
def test_binom_rvs():
'''
为离散分布
使用.rvs函数模拟一个二项随机变量,其中参数size指定你要进行模拟的次数。我让Python返回10000个参数为n和p的二项式随机变量
进行10000次实验,每次抛10次硬币,统计有几次正面朝上,最后统计每次实验正面朝上的次数
'''
binom_sim = data = stats.binom.rvs(n=10,p=0.3,size=10000)
print len(binom_sim)
print "mean: %g" % np.mean(binom_sim)
print "SD: %g" % np.std(binom_sim,ddof=1)
plt.hist(binom_sim,bins=10,normed=True)
plt.xlabel('x')
plt.ylabel('density')
plt.show()
#####################
#泊松分布
#####################
def test_poisson_pmf():
'''
泊松分布的例子:已知某路口发生事故的比率是每天2次,那么在此处一天内发生4次事故的概率是多少?
泊松分布的输出是一个数列,包含了发生0次、1次、2次,直到10次事故的概率。
'''
rate = 2
n = np.arange(0,10)
y = stats.poisson.pmf(n,rate)
print y
plt.plot(n, y, 'o-')
plt.title('Poisson: rate=%i' % (rate), fontsize=15)
plt.xlabel('Number of accidents')
plt.ylabel('Probability of number accidents', fontsize=15)
plt.show()
def test_poisson_rvs():
'''
模拟1000个服从泊松分布的随机变量
'''
data = stats.poisson.rvs(mu=2, loc=0, size=1000)
print "mean: %g" % np.mean(data)
print "SD: %g" % np.std(data, ddof=1)
rate = 2
n = np.arange(0,10)
y = stats.poisson.rvs(n,rate)
print y
plt.plot(n, y, 'o-')
plt.title('Poisson: rate=%i' % (rate), fontsize=15)
plt.xlabel('Number of accidents')
plt.ylabel('Probability of number accidents', fontsize=15)
plt.show()
#####################
#正态分布
#####################
def test_norm_pmf():
'''
正态分布是一种连续分布,其函数可以在实线上的任何地方取值。
正态分布由两个参数描述:分布的平均值μ和方差σ2 。
'''
mu = 0#mean
sigma = 1#standard deviation
x = np.arange(-5,5,0.1)
y = stats.norm.pdf(x,0,1)
print y
plt.plot(x, y)
plt.title('Normal: $\mu$=%.1f, $\sigma^2$=%.1f' % (mu,sigma))
plt.xlabel('x')
plt.ylabel('Probability density', fontsize=15)
plt.show()
#####################
#beta分布
#####################
def test_beta_pmf():
'''
β分布是一个取值在 [0, 1] 之间的连续分布,它由两个形态参数α和β的取值所刻画。
β分布的形状取决于α和β的值。贝叶斯分析中大量使用了β分布。
'''
a = 0.5#
b = 0.5
x = np.arange(0.01,1,0.01)
y = stats.norm.pdf(x,a,b)
print y
plt.plot(x, y)
plt.title('Beta: a=%.1f, b=%.1f' % (a,b))
plt.xlabel('x')
plt.ylabel('Probability density', fontsize=15)
plt.show()
#####################
#指数分布(Exponential Distribution)
#####################
def test_exp():
'''
指数分布是一种连续概率分布,用于表示独立随机事件发生的时间间隔。
比如旅客进入机场的时间间隔、打进客服中心电话的时间间隔、中文 * 新条目出现的时间间隔等等。
'''
lambd = 0.5#
x = np.arange(0,15,0.1)
y =lambd * np.exp(-lambd *x)
print y
plt.plot(x, y)
plt.title('Exponential: $\lambda$=%.2f' % (lambd))
plt.xlabel('x')
plt.ylabel('Probability density', fontsize=15)
plt.show()
def test_expon_rvs():
'''
指数分布下模拟1000个随机变量。scale参数表示λ的倒数。函数np.std中,参数ddof等于标准偏差除以 $n-1$ 的值。
'''
data = stats.expon.rvs(scale=2, size=1000)
print "mean: %g" % np.mean(data)
print "SD: %g" % np.std(data, ddof=1)
plt.hist(data, bins=20, normed=True)
plt.xlim(0,15)
plt.title('Simulating Exponential Random Variables')
plt.show()
test_expon_rvs()
测试运行结果如下:
希望本文所述对大家Python程序设计有所帮助。
来源:https://blog.csdn.net/Yan456jie/article/details/52170481
0
投稿
猜你喜欢
- 以发布目录为例:<OBJECT ID="agobjOraSession" RUNAT=&quo
- 以前没见过这个效果,滚动纵向滚动条看看效果就明白了这样的效果,广告商应该比较喜欢。<!DOCTYPE html PUBLIC &quo
- php去掉数组的第一个值的两种删除方法:1、使用array_shift()函数删除数组的第一个值,语法“array_sh
- Tools -> Build System -> Build New System{ "cmd": [&qu
- 最近一个项目中遇到ASP对FoxPro库表(*.DBF)的操作问题。现实中确有许多应用软件使
- 最近 全栈数据工程师养成攻略 的微信群已经将近500人,开了二群之后为了打通不同微信群之间的消息,花了点时间做了个消息同步机器人,在任意群收
- slice 可以用来获取数组片段,它返回新数组,不会修改原数组。除了正常用法,slice 经常用来将 array-like 对象转换为 tr
- 函数的递归调用:是函数嵌套调用的一种特殊形式具体是指:在调用一个函数的过程中又直接或间接地调用到了本身# 直接调用本身def func():
- 如何在线压缩Access数据库?Access数据库可以在线压缩吗?可以的,代码和说明见下:compact.asp<%option&nb
- 前言对于刚刚下载好的pycharm,初学者使用会有一些问题,这里将介绍关于字体,背景,这些简单的设置将会提升编程的舒适度(下面以PyChar
- 大家都知道系统存储过程是无法用工具导出的(大家可以试试 >任务>生成SQL脚本) 因为系统存储过程一般是不让开发人员修改的。 需
- pre标签会原样保留HTML内容的格式,可是如果宽度过大会把页面撑坏,这时候需要自动换行来帮忙:Making preformate
- oracle数据库的快照是一个表,它包含有对一个本地或远程数据库上一个或多个表或视图的查询的结果。正因为快照是一个主表的查询子集,使用快照可
- 前言defer语句被用于预定对一个函数的调用。我们把这类被defer语句调用的函数称为延迟函数。而defer 延迟语句在其他编程语言里好像没
- 本文实例讲述了php基于协程实现异步的方法。分享给大家供大家参考,具体如下:github上php的协程大部分是根据这篇文章实现的:http:
- 互联网是一个飞速发展的行业,任何的止步不前都会导致被淘汰,只是时间早晚的问题,所以一个公司的学习与创新能力是非常重要的,特别是对于一个年轻的
- 最近,使用 golang 去管理本地应用的生命周期,期间有几个有趣的点,今天就一起看下。场景一我们来看看下面两个脚本会产
- 昨天我问过这个问题怎么用ADODB.Stream来读取或写入文件,而不是用fso,不过没人回答到点上,今天搞定了.贴出来给觉得有用的朋友,希
- 相信每个前端工程师都有自己喜爱的javascript框架,说情感也好,道信仰也罢,javascript框架带给人的不仅仅是便捷的开发,更有一
- PDO::beginTransactionPDO::beginTransaction 启动一个事务(PHP 5 >= 5.1.0, P