matlab - 从 Matlab 中预先指定的概率质量函数中绘制随机数

Question

我在 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。现在，我想从中抽取n随机数pr_mass_epsilon。我怎样才能做到这一点？我想要一个适用于任何suppsize_epsilon.

标签： matlabrandomdistributionprobability-theory

Answer 1

使用统计工具箱

该randsample函数可以直接执行此操作：

result = randsample(supp_epsilon, n, true, pr_mass_epsilon);

不使用工具箱

手动方法：

1. 在区间 (0,1) 中生成n均匀随机变量的样本。
2. 将每个样本与分布函数（质量函数的累积和）进行比较。
3. 查看每个均匀样本位于分布函数的哪个区间。
4. 索引到可能值的数组

result = supp_epsilon(sum(rand(1,n)>cumsum(pr_mass_epsilon(:)), 1)+1);

---

对于您的示例，使用n=1e6这两种方法中的任何一种都会给出与此类似的直方图：

histogram(result, 'normalization', 'probability')