numpy.random和scipy.stats中定义的幂律分布在数学意义上没有为负a定义,如

this question的答案中所解释的那样:由于奇点为零,它们不可规范化.所以,遗憾的是,数学说’不’.

您可以使用与x ^ {g-1}成比例的pdf定义分布,其中g

In [3]: def rndm(a, b, g, size=1):

"""Power-law gen for pdf(x)\propto x^{g-1} for a<=x<=b"""

...: r = np.random.random(size=size)

...: ag, bg = a**g, b**g

...: return (ag + (bg - ag)*r)**(1./g)

然后,你可以这样做,例如,

In [4]: xx = rndm(1, 2, g=-2, size=10000)

等等.

为完整起见,这里是pdf:

In [5]: def pdf(x, a, b, g):

ag, bg = a**g, b**g

....: return g * x**(g-1) / (bg - ag)

这一切都假设a< b和g!= 0.这些公式应该与numpy.power和scipy.stats.powerlaw一致,a = 0,b = 1和g> 0.

python幂函数无序分布_Python:从幂律分布中生成随机数相关推荐

  1. Python数据可视化系列之幂律分布

    1.幂律分布 首先要说的是中心极限定理--在复杂的多因素情况下,只要个体相互独立,集体效果就应该是正态分布.然而实际运用中,尤其是金融中,更多面对的是尖峰胖尾现象,比如下面这幅图描述的是标普500指数 ...

  2. 详解幂律分布,以及用于重尾分布的Python库powerlaw的使用

    详解幂律分布,以及用于重尾分布的Python包powerlaw 总述 幂律分布简要回顾 powerlaw库拟合效果 powerlaw库基本操作介绍 可视化 拟合范围 离散与连续数据 与其他分布比较 总 ...

  3. 幂律分布(python)

    幂律分布 1.幂律分布 首先要说的是中心极限定理--在复杂的多因素情况下,只要个体相互独立,集体效果就应该是正态分布.然而实际运用中,尤其是金融中,更多面对的是尖峰胖尾现象,比如下面这幅图描述的是标普 ...

  4. 【转载】关于幂律分布的一个笔记

    关于幂律分布的一个笔记 原文转自:http://blog.sina.com.cn/s/blog_55954cfb0100ps89.html 0:题外话或补记 最早知道二八法则,还是一本介绍犹太民族杰出 ...

  5. 关于幂律分布的一个笔记_哈克_新浪博客

    关于幂律分布的一个笔记_哈克_新浪博客 关于幂律分布的一个笔记     (2011-03-02 18:12:27)     转载▼     标签:     幂律     二八法则     杂谈     ...

  6. 对数正态分布_正态分布与幂律分布

    1.正态分布 英国生物统计学家法兰西斯·高尔顿做了一个实验.他在一块木板上画了一块等腰三角形,并在三角形区域内钉上n+1层钉子.第1层钉2个钉子,第2层钉3个钉子,下面每一层都比上一层增加一个钉子,上 ...

  7. Power law and Power law distribution(幂律和幂律分布)

    原文:<Power-law distribution in empirical data> 1. Introduction 有些分布可以很好的描述,比如成年男性的身高,某物体的重量等,它们 ...

  8. Python数据可视化:幂律分布

    1.公式推导  对幂律分布公式:    对公式两边同时取以10为底的对数:    令,且为常数,所以公式变为:  所以对于幂律公式,对X,Y取对数后,在坐标轴上为线性方程. 2.可视化  从图形上来说 ...

  9. [Python舆情分析] 二.时间间隔分布研究及幂律分布图绘制

    本文主要是作者学习舆情分析.情感分析.人类行为动力学分析的在线笔记,主要包括两方面内容,一是幂律特性,二是讲解时间间隔分布,三是Python绘制基于时间间隔分布的幂律特性图,四提供了另一种方法.基础性 ...

最新文章

  1. android sdk启动不了,windows server 2008下android sdk不能正常启动
  2. 基础的shell编程问题(二)
  3. unity基础知识笔记一(快捷方式、基础概念)
  4. 商场内自动扶梯的研究
  5. java中bitconverter_【Java】BitConverter(数字转字节数组工具类)
  6. 太实用了 Python 合成多张图片到PDF格式
  7. 【对抗攻击代码实战】对抗样本的生成——FGSM
  8. Java性能优化全攻略
  9. 电脑怎么设置计算机系统,细说电脑怎么设置wifi
  10. 怎么把ppt弄成链接的形式_ppt链接excle表格:如何在ppt中超级链接到指定的excel工作表...
  11. 聊聊 SAP 产品 UI 上的消息显示机制
  12. excel ctrl shift+键盘方向键的使用
  13. echarts 饼状图 (数据为0或很小的扇形显示问题、扇形间隔)
  14. 辗转相除法——求最大公约数(易懂详解)
  15. 【windows下直接把文件打包成war包】
  16. 【Unity】场景异步加载、卸载
  17. 【全网热点】打造全网最全爱心代码仓库【火速领取爱心】
  18. MVC依赖注入Unity应用学习
  19. 基于Python实现制作的塔防游戏
  20. 微软语音朗读引擎MSSPEECH打造自己的TTS应用

热门文章

  1. linux----------2--3----(无名)管道通信原理及管道编程实战
  2. 分销系统具有哪些大优势?
  3. X-NUCA'2019部分题目WP
  4. 时间与日期处理模块 题目解题报告
  5. Delphi Web前端开发教程(2):基于TMS WEB Core框架
  6. CSP-J1 CSP-S1 信奥 第1轮 初赛 数据分析 成绩及分数线汇总
  7. 洗地机性价比高的是哪款?性价比高适合入手的洗地机介绍
  8. 助力高效办公,TOM企业邮箱网页版全面升级
  9. NRF52832+DW1000通信系统架构设计
  10. 计算机的知识体系结构,计算机体系结构知识点