Generating integer random numbers from power-law distribution

6 次查看(过去 30 天)
Is there a code to generate integer random numbers through power-law distribution of the form p(x)=x^-(alpha+1)?
thanks in advance

采纳的回答

Image Analyst
Image Analyst 2013-12-13
Try this: RANDRAW

更多回答(3 个)

Miguel Ángel
Miguel Ángel 2013-12-13
编辑:Miguel Ángel 2013-12-13
I don't know if round,ceil or floor is correct, but in case this may be a way
n=1000;
p=rand(1,n);
alpha=0.1;
%inverse of distribution
%properties of logs %LOGa(n) = LOGb(N) / LOGb(a)
%~round/ceil/floor
x=round(log(p)/log(alpha+1));
  1 个评论
Sunstar
Sunstar 2013-12-13
Thank you for your reply. Unfortunately, I cannot simply round or floor or ceil. I need to generate real integers.

请先登录,再进行评论。


Roger Stafford
Roger Stafford 2013-12-13
The distribution you mention is equivalent to the zeta distribution with alpha+1 replaced by s. See the website
http://en.wikipedia.org/wiki/Zeta_distribution
As Image Analyst has mentioned, there is a File Exchange contribution called 'randraw' which is located at:
http://www.mathworks.com/matlabcentral/fileexchange/7309-randraw
which claims to generate random numbers in accordance with fifty various distributions including the zeta distribution. It might be what you need.

Weijia Chen
Weijia Chen 2019-9-18
I also meet the same problem. The minimum integer that zeta distribution can generate seems to be 1, but I hope 0 can be generated. Is there any other distributin can do this? Thank you

类别

Help CenterFile Exchange 中查找有关 Random Number Generation 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by