×

一种改进的ziggurat算法,用于生成指数分布和正态分布的伪随机数。 (英文) Zbl 1510.62066号

摘要:ziggurat算法是一种从统计分布生成伪随机数(PRN)的快速拒绝采样方法。在该算法中,矩形采样域层叠在一起(类似于锯齿形),以封装所需的概率密度函数。对这些层中的随机值进行采样,如果它们位于概率密度函数图的下方,则返回这些值。在这里,我们提出了一种实现,其中ziggurat层完全位于概率密度函数之下,从而消除了在ziggura层中进行任何拒绝测试的需要。在新算法中,概率密度的小悬垂段保留在每个锯齿层的右侧,可以使用三角形采样域对其进行有效采样。指数和正态变量的新算法的中位数运行时间分别减少到58%和53%(总体范围:41-93%)。提供了一个可访问的C库,以及对Python和MATLAB/Octave的扩展。

MSC公司:

62-08 统计问题的计算方法
65立方厘米 数值分析中的随机数生成
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] L'Ecuyer P.测试随机数生成器。In:冬季模拟会议。新泽西州皮斯卡塔韦:电气和电子发动机研究所;1992年,第305-313页。[谷歌学者]
[2] Marsaglia G,Tsang WW公司。一种快速、易于实现的方法,用于从递减或对称单峰密度函数中进行采样。SIAM科学统计与计算杂志。1984;5(2):349-359. doi:10.1137/0905026[Crossref],[Web of Science®],[Google学者]·Zbl 0573.65116号
[3] Marsaglia G,Tsang WW公司。生成随机变量的ziggurat方法。J统计软件。2000;5(8):1-7. [Crossref],[Google学者]
[4] 日本Doornik。一种生成正常随机样本的改进ziggurat方法。牛津:牛津大学;2005.[谷歌学者]
[5] Zhang G,Leong PHW,Lee DU,Villasenor JD,Cheung RC,Luk W.基于Ziggurat的硬件高斯随机数生成器。参加:现场可编程逻辑和应用国际会议。坦佩雷:IEEE;2005年,第275-280页。[谷歌学者]
[6] 鲁宾H,约翰逊BC。指数偏差和正态偏差的有效生成。J统计计算模拟。2006;76(6):509-518. doi:10.1080/10629360500108004[Taylor&Francis Online],[Web of Science®],[Google学者]·Zbl 1157.62329号
[7] Thomas DB、Luk W、Leong PH、Villasenor JD。高斯随机数生成器。ACM计算调查。2007;39(4):11个。可从以下位置购买:http://portal.acm.org/citation.cfm?doid=1287620.1287622。doi:10.1145/128762.1287622[Crossref],[Web of Science®],[Google学者]
[8] Oved I.计算超越函数。2003年[引用日期:2014年3月24日];可从以下位置获得:http://math.arizona.edu/aprl/teach/iriso/oversive.ps.[谷歌学者]
[9] de Schryver C,Schmidt D,Wehn N,Korn E,Marxen H,Korn R.一种新的基于硬件高效反演的非均匀分布随机数生成器。2010年可重构计算和FPGA国际会议(ReConFig)。墨西哥坎昆:IEEE;2010年,第190-195页。[谷歌学者]
[10] 史密斯WD。如何从概率分布中取样。2002年4月【引用2014年3月24日】;可从以下位置获得:http://scorevoting.net/WarrenSmithPages/homepage/sampling.ps。[谷歌学者]
[11] Saito M,Matsumoto M.Simd-oriented fast mersenne twister:一个128位伪随机数生成器。蒙特卡罗和准蒙特卡罗方法2006。柏林:斯普林格·弗拉格;2008年,第607-622页。[Crossref],[Google学者]·Zbl 1141.65319号
[12] Marsaglia G,Tsang WW公司。生成gamma变量的简单方法。ACM Trans数学软件。2000;26(3):363-372. doi:10.1145/358407.358414[Crossref],[Web of Science®],[Google学者]·Zbl 1365.65022号
[13] Story S,Tak P,Tang P。ia-64上改进超越函数的新算法。IEEE计算机算法研讨会。1999;4.【谷歌学者】
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。