我在MatLab中有一个支持(supp_epsilon)和一个概率质量函数(pr_mass_epsilon)，构造如下。

supp_epsilon=[0.005 0.01 0.015 0.02]; 

suppsize_epsilon=size(supp_epsilon,2);

pr_mass_epsilon=zeros(suppsize_epsilon,1);

alpha=1;
beta=4;

for j=1:suppsize_epsilon
    pr_mass_epsilon(j)=betacdf(supp_epsilon(j),alpha,beta)/sum(betacdf(supp_epsilon,alpha,beta));
end

请注意，pr_mass_epsilon的组件总和为1。现在，我想从pr_mass_epsilon中抽取n随机数。我怎么能这样做呢？我想要一个代码，任何suppsize_epsilon的工作。
换句话说：我想从supp_epsilon中随机抽取元素，每个元素的概率由pr_mass_epsilon给出。