Someone please help me, I have probability mass function as follwing, How can generate 100 random values of j from this probability mass function?

2 次查看（过去 30 天）
mohammed elmenshawy
编辑： Bruno Luong 2020-8-21
Note: λ>0 and 0≤α≤1/2 and 0≤P(X=j)≤1. (Take initial valus for λ, α and p(X=j) from thir domain as you like as example)
I want to generate randomly 100 numbers (values for j) with this probability mass function . But I really have no idea how and where to start.
Can somebody help me?
Thank you in advance

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

采纳的回答

Jeff Miller
Jeff Miller 2020-8-21
If you would like to generate 100 random j values, use the function KSSV gave you to generate the f values for the probability masses (they are not all 0's after the for loop, only before it). Then:
sump = cumsum(f);
r_unif = rand(100,1);
r_j = zeros(100,1);
for i=1:100
r_j(i) = find(sump>r_unif(i),1);
end
figure;
histogram(r_j)

更多回答（2 个）

KSSV
KSSV 2020-8-20
  1. DEfine alpha and lambda values.
  2. DEfine j value.
  3. WRite your formula and substitute those values.
  6 个评论
显示 4更早的评论

请先登录，再进行评论。

Bruno Luong
Bruno Luong 2020-8-21
编辑：Bruno Luong 2020-8-21
alpha=0.01; lambda=0.99;
f=@(j) ((1-alpha).^j-(-alpha).^j).*lambda.^j*exp(lambda*alpha)./(factorial(j)*(exp(lambda)-1));
p=f(1:22); % f is practically 0 beyond 10
c=cumsum(p); c=c/c(end);
n = 100; % 1e6;
[~,j] = histc(rand(1,n),[0 c]);

类别

在 Help Center 和 File Exchange 中查找有关 MATLAB 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by