wolfram-mathematica - 随机泊松噪声

Question

我在 Mathematica 中寻找与以下 Matlab 函数等效的函数:

“R = poissrnd(lambda) 从具有平均参数 lambda 的泊松分布生成随机数。lambda 可以是向量、矩阵或多维数组。R 的大小是 lambda 的大小。”

下面的函数输出示例。

b = 95.7165   95.7165   95.7165   95.7165   95.7165   98.9772   98.9772   98.9772   98.9772    0.3876

poissrnd(b)

ans =100   115     81    90   109   106   104    87   104     2

我怎样才能在 Mathematica 8 中做类似的事情？

Answer 1

Poisson distribution仅为整数定义。所以你需要使用 RandomInteger与 PoissonDistribution像这样:

poissrnd[lambda_]:=RandomInteger[PoissonDistribution[lambda]]

用法:

b = {95.7165, 95.7165, 95.7165, 95.7165, 95.7165, 98.9772, 98.9772, 
  98.9772, 98.9772, 0.3876};

poissrnd /@ b

Out[1] = {104, 97, 67, 84, 96, 123, 93, 96, 100, 0}