首先,我会写这个答案,假设你将numpy导入为np,因为它清楚地区分了numpy函数与内置函数或python的数学和随机包的函数.

我认为没有必要回答您指定的问题,因为您的基本假设是错误的:

是的,泊松统计量具有等于方差的均值,但假设您使用常数林.但你没有.你输入高斯的y值,所以你不能指望它们是恒定的(它们是你的定义高斯!).

使用np.random.poisson(lam = 0.5)从泊松分布中获取一个随机值.但要小心,因为这个泊松分布甚至与高斯分布几乎不相同,因为你处于“低均值”区间,这两个区间显着不同,例如参见Wikipedia article about Poisson distribution.

你也在创建随机数,所以你不应该真正绘制它们,而是绘制它们的np.histogram.由于统计分布都与概率密度函数有关(见Probability density function).

之前,我已经提到你创建了一个带有常数lam的泊松分布,所以现在是时候讨论大小:你创建随机数,所以要近似真实的泊松分布,你需要绘制大量的随机数.其大小来自:np.random.poisson(lam = 0.5,size = 10000)例如创建一个10000个元素的数组,每个元素从泊松概率密度函数中绘制,平均值为0.5.

如果您在之前提到的维基百科文章中没有阅读它,则泊松分布根据定义仅给出无符号(> = 0)整数作为结果.

所以我想你想要做的是创建一个包含1000个值的高斯和泊松分布:

gaussian = np.random.normal(0.5, 2*np.sqrt(2*np.log(2)), 1000)

poisson = np.random.poisson(0.5, 1000)

然后绘制它,绘制直方图:

import matplotlib.pyplot as plt

plt.hist(gaussian)

plt.hist(poisson)

plt.show()

要从随机样本中获取统计数据,您仍然可以在高斯和泊松样本上使用np.var和np.mean.而这次(至少在我的样本运行中)他们给出了很好的结果:

print(np.mean(gaussian))

0.653517935138

print(np.var(gaussian))

5.4848398775

print(np.mean(poisson))

0.477

print(np.var(poisson))

0.463471

请注意高斯值几乎与我们定义的参数完全相同.另一方面,泊松均值和变量几乎相等.您可以通过增加上面的大小来提高均值和var的精度.

为什么泊松分布不接近原始信号

原始信号仅包含0到1之间的值,因此泊松分布仅允许正整数,标准差与平均值相关联.从高斯的平均值到目前为止,你的信号大约为0,因此泊松分布几乎总是为0.高斯有最大值时,值为1. 1的泊松分布看起来像这样(左边是信号泊松和右边的泊松分布值为1)

所以你会在那个地区得到很多0和1以及2.但是也有可能你将值绘制到7.这正是我提到的反对称性.如果你改变高斯的幅度(例如乘以1000),那么“拟合”要好得多,因为泊松分布在那里几乎是对称的:

python生成泊松分布_Python Numpy泊松分布相关推荐

  1. python生成等差数列_python numpy函数中的linspace创建等差数列详解

    前言 本文主要给大家介绍的是关于linspace创建等差数列的相关内容,分享出来供大家参考学习,下面话不多说了,来一起看看详细的介绍吧. numpy.linspace 是用于创建一个由等差数列构成的一 ...

  2. python生成时间戳_python生成13位或16位时间戳以及反向解析时间戳的实例

    我就废话不多说了,直接上代码吧! import datetime import time def get_float_time_stamp(): datetime_now = datetime.dat ...

  3. python生成数字_Python生成数字图片代码分享

    本文向大家分享了几段Python生成数字图片的代码,喜欢的朋友可以参考.具体如下: 最终版本 # -*- coding:utf-8 -*- from PIL import Image,ImageFon ...

  4. python npv 计算公式_Python numpy 中常用的数据运算

    Numpy 精通面向数组编程和思维方式是成为Python科学计算大牛的一大关键步骤.--<利用Python进行数据分析> Numpy(Numerical Python)是Python科学计 ...

  5. python生成圆形_python生成圆形图片的方法

    本文实例为大家分享了python生成圆形图片的具体代码,供大家参考,具体内容如下 # -*- coding: utf-8 -*- """ __author__= 'Du' ...

  6. python生成零矩阵_python 实现矩阵填充0的例子

    python 实现矩阵填充0的例子 需求: 原矩阵 [[1 2 3] [4 5 6] [7 8 9]] 在原矩阵元素之间填充元素 0,得到 [[1. 0. 2. 0. 3.] [0. 0. 0. 0. ...

  7. python如何生成等差数列_python numpy函数中的linspace创建等差数列详解

    前言 本文主要给大家介绍的是关于linspace创建等差数列的相关内容,分享出来供大家参考学习,下面话不多说了,来一起看看详细的介绍吧. numpy.linspace 是用于创建一个由等差数列构成的一 ...

  8. python维度变换_Python NumPy用法

    介绍 NumPy是Python数值计算最重要的基础包,大多数提供科学计算的包都是用NumPy的数组作为构建基础.NumPy本身并没有提供多么高级的数据分析功能,理解NumPy数组以及面向数组的计算,将 ...

  9. python数组初始化_Python Numpy 数组的初始化和基本操作

    Python 是一种高级的,动态的,多泛型的编程语言.Python代码很多时候看起来就像是伪代码一样,因此你可以使用很少的几行可读性很高的代码来实现一个非常强大的想法. 一.基础: Numpy的主要数 ...

  10. Python生成随机数和numpy生成随机数

    生成随机数 使用random函数 random.random() 产生一个包含0不包含1之间的小数 random.randint(1,5) 产生一个包含开始和结束的整数 random.uniform( ...

最新文章

  1. 张宏江:开源时代如何解决人的思维孤岛
  2. linux shell 检查端口开放
  3. python为什么从0开始_为什么大多数编程语言中,索引都是从0开始
  4. iOS App的几种多视图介绍
  5. 使用CuteFTP登陆FTP(serv-U)服务器后无法LIST目录和文件的解决方法
  6. java里冒泡排序编程案例_冒泡排序法-java案例详解
  7. 使用OData协议查询Windows日志
  8. linux之解决libipopt.so.1: Cannot open shared object file
  9. redis下key的过期时间详解
  10. 你有什么道理后悔没有早点知道?
  11. 基于IR2136的三相电机控制
  12. Python玩转图像格式转换
  13. 适用于高速公路的查询软件,数据稳定、免维护,可查询高速路况、路线规划、未来天气等信息
  14. 06-课堂笔记-包相关
  15. developerWorks 中国Java technology文档库Spring 系列
  16. 计算机研究生论文数学公式,研究生论文公式符号细则.doc
  17. 大学计算机网络技术基础--笔记大全
  18. 9-ICMP报文种类
  19. 【数据结构】代码重现
  20. 计算机网络应用押韵句,沙雕押韵很可爱的rap的句子 笑掉大牙的押韵句子

热门文章

  1. 移动播放器html,支持移动平台的Html5播放器
  2. 使用vue 开发chrome 插件
  3. 机器语言入门 w3c,编程语言
  4. java并发测试工具_Java并发测试工具类
  5. JMeter并发测试工具
  6. 解决百度地图生成器添加标注后图标不显示的问题
  7. 利用Excel爬取网页数据
  8. SILVACO 学习笔记第一章
  9. PlaySound的同步播放与异步播放
  10. 五金模具设计统赢外挂提升效率技巧、外挂模具设计流程、常见问题归纳