MATLAB: Erlang Random variable

ahmed raheem

Hi all, i need to generate a random variable with erlang distribution, so is there anyone can help me please?

thanks…

Best Answer

In theory Erlang is the sum of exponential distributed random variables, so you can just do the same. If you have Statistics Toolbox, you can use exprnd to draw exponentially distributed random number.

Related Solutions

MATLAB: Random

x   = [1     2     3     4     5
6     7     8     9    10
11    12    13    14    15
16    17    18    19    20];

EDIT after Mohammad's coment

[m,n] = size(x);
df = bsxfun(@plus,fullfact(m*ones(1,n)),(0:n-1)*m);
k = randperm(size(df,1));
out = x(df(k(1:min(100,numel(k))),:));

[m,n] = size(x);
xc = mat2cell(x,m,ones(n,1));
cc = cell(1,n);
[cc{:}] = ndgrid(xc{:});
k = randperm(numel(cc{1}));
out =  cell2mat(cellfun(@(x)x(k(1:min(100,numel(k)))'),cc,'un',0));

MATLAB: Erlang distribution n=40 μ=8; 9; k=1; 10;

Using the definition of erlang pdf from Wikipedia https://en.wikipedia.org/wiki/Erlang_distribution, the following program can be written

mu = [8 9];
k = [1 10];
erlang_fun = @(x, mu, k) x.^(k-1).*exp(-x/mu)/(mu.^k.*factorial(k-1));
x = linspace(0, 300);
figure;
for i=1:numel(mu)
    subplot(1,2,i)
    plot(x, erlang_fun(x, mu(i), k(i)));
    title(sprintf('\\mu=%d, k=%d', mu(i), k(i)));
end

Best Answer

Related Solutions

MATLAB: Random

MATLAB: Erlang distribution n=40 μ=8; 9; k=1; 10;

Related Question