monte carlo random number generation

17 次查看（过去 30 天）
stavros boutsikos
stavros boutsikos 2021-3-29
评论： stavros boutsikos 2021-3-29
i want to use a function that generates a random number from [1,2,3,4,5,6,7,8,9,10] but the probabilities of each number are different. For exaample the pof drawing one is 0.2 but p of drawing 3 is 0.06? is there a function or a code that i can use in order to generate a random number by taking into consideration that each number has different weighting?

请先登录，再回答此问题。

采纳的回答

Andreas Apostolatos
Hello,
It seems that what you are seeking for can be utilized using function 'randsample' when using the variant that takes as argument variable 'population', which in your case is vector '[1,2,3,4,5,6,7,8,9,10]', and an additional vector of weights 'w' that has the same size as variable 'population', see the following documentation link for more information,
https://www.mathworks.com/help/stats/randsample.html#d123e694964
Note that function 'randsample' is using a uniform probability distribution and it is provided as part of the Statistics and Machine Learning Toolbox. The topic is also discussed in the following MATLAB Answers thread where additional workflows are discussed,
https://www.mathworks.com/matlabcentral/answers/502828-generating-a-weighted-random-number
I hope that this information helps.
Kind Regards,
Andreas
  1 个评论
stavros boutsikos
stavros boutsikos 2021-3-29
thank you very much Andrea, it was really helpfull and it helped me realise that my code was working and i was just overthinking it.

请先登录，再进行评论。

更多回答（1 个）

Image Analyst
Image Analyst 2021-3-29
You can simply get the cumulative sum and compare your random number to that. For example
r = rand
c = cumsum(p)
output = find(r < c, 1, 'last')
If that doesn't work for you, then give us all the probability values.
  1 个评论
stavros boutsikos
stavros boutsikos 2021-3-29
thank you very much, it actually helped me understand that the code i had was actually working!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Uniform Distribution (Continuous) 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by