Easy question with probability

3 次查看(过去 30 天)
Hi,
Let's say
P = [0.1 0.3 0.4 0.2]
X = [1 2 3 4]
where P(n) is the probability to select the X(n) element. I wish to make a function that select an "random" element of X according to its probability, like
f = myfun(P,X)
>> f = 2
thx a lot

采纳的回答

Andrei Bobrov
Andrei Bobrov 2011-12-7
function F = myfun(P,X)
x = cumsum([0 P(:).'/sum(P(:))]);
x(end) = 1e3*eps + x(end);
[a a] = histc(rand,x);
F = X(a);
  5 个评论
Oleg Komarov
Oleg Komarov 2011-12-7
Andrei's function should read:
function F = myfun(P,X)
p = cumsum([0; P(1:end-1).'; 1+1e3*eps]);
[a a] = histc(rand,p);
F = X(a);
Jyotshna Ramashire
Jyotshna Ramashire 2017-3-18
can you please explain the codes?

请先登录,再进行评论。

更多回答(1 个)

Roger Stafford
Roger Stafford 2017-3-18
Also a little late, this also works:
C = cumsum(P);
f = X(1+sum(C(end)*rand>C));

类别

Help CenterFile Exchange 中查找有关 Logical 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by