matlab中的泊松点过程

时间:2015-01-14 05:03:37

标签: python matlab probability poisson

我是poisson点过程的新手。我做了一个simluation(matlab)如下。 我的强度lambda = 50 ;

clear all;
lambda=50;
 npoints = poissrnd(lambda);
  pproc = rand(npoints, 2);
  plot(pproc(:, 1), pproc(:, 2), '.');

然后我有情节, enter image description here

然而,链接

http://connor-johnson.com/2014/02/25/spatial-point-processes/

告诉我,当强度 lamuda = 0.2,小于1 时,他得到了

enter image description here

该链接还显示了Python中的代码。请检查它。

这是我的问题,为什么强度小于1 ,他还可以在这里绘制一些东西?如果我让我的代码的lamda = 0.2,那么绘制就没有价值了。我想我对Poisson点过程有所遗漏?或者这是一个编程问题?我想在matlab中模拟这个lambda = 0.2。

非常感谢你的帮助。

总lambda = 0.4。我想拥有100个单位。这意味着单位强度= 0.4 * 100 = 40.非常感谢你。

 la=0.4;
 lala=0.4*100;
  npoints = poissrnd(lala);

  pproc = rand(npoints, 2);

  plot(pproc(:, 1).*100, pproc(:, 2).*100, '.');

enter image description here

2 个答案:

答案 0 :(得分:1)

我认为你的解决方案应该遵循。 (注:绘制PPP空间分布的两种不同方法)

clear all;
clc;
close all;

lambda=50;

%Method 1
pproc  = poissrnd(lambda, 100, 2);
size(pproc)
plot(pproc(:, 1), pproc(:, 2), '.');
title('Poisson with poissrnd')

%Method 2
pproc2 = random('Poisson', lambda, 100, 2);
size(pproc2)
figure;
plot(pproc2(:, 1), pproc2(:, 2), '.');
title('Poisson with Random statement')

答案 1 :(得分:0)

总lambda = 0.4。我想拥有100个单位。这意味着单位强度= 0.4 * 100 = 40.非常感谢你。

 la=0.4;
 lala=0.4*100;
  npoints = poissrnd(lala);

  pproc = rand(npoints, 2);

  plot(pproc(:, 1).*100, pproc(:, 2).*100, '.');

请在最后检查我的问题中的数字。