About simluation of the poisson point process. Thank you

6 次查看(过去 30 天)
Dear friends,
I am new with poisson point process. I did one simluation as below. My intensity lambda = 50;
clear all;
lambda=50;
npoints = poissrnd(lambda);
pproc = rand(npoints, 2);
plot(pproc(:, 1), pproc(:, 2), '.');
Then I have plot,
However, the link
http://connor-johnson.com/2014/02/25/spatial-point-processes/
showed me that when intensity lamuda = 0.2, smaller than 1 , he got
The link also showed the code in Python.Please check it.
Here is my question, why intensity is smaller than 1, he still can plot something here? If I let my code's lamda = 0.2, there will be no value to plot. I think I miss something about poiion point process? or it's a programming problem?
Thank you so much for your help.

采纳的回答

sun
sun 2015-1-14
Guys, I think I got it. I was plotting the unit figure. but the 2nd figure is not plotting the unit but 0 - 20. Thank you.
la=0.4;
lala=0.4*100;
npoints = poissrnd(lala);
pproc = rand(npoints, 2);
plot(pproc(:, 1).*100, pproc(:, 2).*100, '.');

更多回答(1 个)

Syed
Syed 2016-10-10
Don't you think it should be something like below.
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')

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by