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

hi

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

Best Answer

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)
Related Question