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 日間)
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 件のコメント
mohammed elmenshawy
mohammed elmenshawy 2020 年 8 月 21 日
Thank you very much for your cooperation.

サインインしてコメントする。


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 ExchangeSystems of Nonlinear Equations についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by