在 Matlab 中使用离散帕累托分布进行采样

Question

我想在 matlab 中生成具有离散帕累托分布的 10000 数字。

我用 R 写了这段代码，但我想用 Matlab 写。

离散帕累托的 R 码：

N<-5 ##for instance,you need N large enough
i0<-(1:N)
p<-1/(i0^1.98); p<-p/sum(p)
X<-sample(i0,N,replace=TRUE,prob=p)
X

alpha=1.98 并且您在此链接中看到离散的 pPareto 分布 q=0， https://en.wikipedia.org/wiki/Zipf%E2%80%93Mandelbrot_law

我的问题是：如何在 Matlab 中编写这段代码？

我在matlab中使用randsample，并在matlab中翻译了这段代码。

例如：

N=5;
for i0=1:N
    p=1/(i0.^1.98);
    P=p/sum(p);`
    X=randsample(i0,N,true,P);
    X
end

但它有一个错误。

谁能给个建议？

Answer 1

以下是一些建议：

1. 您的代码中有错误： P=p/总和（p）；` 应该 P=p/sum(p);

（这可能是一个格式化的东西，只是需要编辑，但不清楚你的错误是什么：所以从这里开始）

1. 您可能需要 Statistics and Machine Learning Toolbox 才能使用 randsample，这需要额外的许可证。

2. 如果这些都不能解决问题，请尝试输入 帮助 randsample 在 MATLAB 终端上并验证您的设置。

这看起来是您感兴趣的部分（来自帮助）：

Y = randsample(N,K,true,W) 或 randsample(POPULATION,K,true,W) 返回一个 加权样本，使用正权重 W，带放回。 W是 通常是概率向量。此功能不支持加权 无需更换即可采样。

1. 我认为您的代码可能存在误解。看起来您认为 P 应该是一个分布，但这意味着 p 将是一个数组或向量。但是，从您的代码中

p=1/(i0.^1.98);

只是将评估为一个标量值（即p，对于 i0=1，评估为 1^(-1.98)，因此 P 将等于 1。

也许您正在寻找类似的东西？

N=5;
i0=1:N;
p=1./(i0.^1.98);
P=p/sum(p);
X=randsample(i0,N,true,P)  

（我省略了最后一个;，以便您可以将结果显示到终端）