y = randsample(population,k)

Question

riad didou 2020-8-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/578751-y-randsample-population-k

编辑： Bruno Luong 2020-8-12

Please, I need to know the name of algorithm used by MATLAB 'randsample' function for the random selection without replacement

y = randsample(population,k)

3 个评论
显示 1更早的评论隐藏 1更早的评论

riad didou 2020-8-12

编辑：riad didou 2020-8-12

But i can't find the name of the algorithm for the randsample (i didn't mean the RNG function)!

Adam Danz 2020-8-12

编辑：Adam Danz 2020-8-12

But my answer tells you how to find that.

The answer explains that the syntax you're using in randsample is using randperm which uses a uniform pseudorandom number generator of which there are several types. See the "To determine which one you're using..." section to learn how to find which type you're using.

请先登录，再进行评论。

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

Answer 1

Bruno Luong 2020-8-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/578751-y-randsample-population-k#answer_478836

编辑：Bruno Luong 2020-8-12

TMW never mentions it officially, but it's Fisher-Yates shuffle algorithm

(source https://uk.mathworks.com/matlabcentral/answers/12248-what-does-randperm-miss, P. Perkin's post)

4 个评论
显示 2更早的评论隐藏 2更早的评论

Adam Danz 2020-8-12

编辑：Adam Danz 2020-8-12

在 MATLAB Online 中打开

Thanks for pointing to that discussion, Bruno Luong. Note that this algorithm still depends on the random number generator state settings as I tried to explain in this answer but knowing about the Fisher-Yates-Durstenfeld shuffle (FYD) is useful.

For OP, the FYD can be understood as randomly drawing k values from a hat containing n values where k<n. The "random" part uses the rng settings.

For example,

% Use the default Mersenne Twister(mt19937ar) with seed 9999
rng(9999, 'twister')
randsample(100,8)'
% ans =
%     43    61    83    55     4    16    22    19
% You can repeate those two lines and will get the same values. 
% Now change the random number generator to Combined multiple recursive (mrg32k3a) with the same seed
rng(9999, 'combRecursive')
randsample(100,8)'
% ans =
%     80     9    63    42    78    13    89    71
% You can repeate those two lines and will get the same values but note how they differ from the Mersenne Twister.
% Reset generator back to default
rng('default')

Maybe the OP didn't intend to dig that deeply.

Bruno Luong 2020-8-12

Read the pseudo code on the wikipedia link I posted above.

Adam Danz 2020-8-12

编辑：Adam Danz 2020-8-12

@riad didou, if you plan on using Matlab's randsample function, you will be relying on a random number generator.

I don't know how else to explain it other than this:

If you call randsample(n,k) or randsample(population, k), that function uses the randperm function whether you want it to or not. According to Matlab's document which I mentioned in the original thread, "The sequence of numbers produced by randperm is determined by the internal settings of the uniform pseudorandom number generator".

"Even if we don't use a random number generator the algorithm still working"

That's not true. The link Bruno provided vaguely describes that values from your array are "randomly" chosen. That requires you to use a random process which is what the random number generator is. Even if you write your own version that doesn't use rng() you'll still need to replace it with some kind of random processes.

Maybe it would help to know your end goal. If the goal it to write about it in a methods section of a paper, I wonder what bits of information you're missing. Many people just write that the function randsample was used in Matlab and they they cite Mathworks but you'd need to document the generator type and seed for reproducibility.

请先登录，再进行评论。