


poisson.m simulates a homogeneous Poisson process


0001 % poisson.m simulates a homogeneous Poisson process 0002 0003 lambda=10; % arrival rate 0004 Tmax=3; % maximum time 0005 0006 clear T; 0007 T(1)=random('Exponential',1/lambda); 0008 i=1; 0009 0010 while T(i) < Tmax, 0011 T(i+1)=T(i)+random('Exponential',1/lambda); 0012 i=i+1; 0013 end 0014 0015 Y=zeros(1,i-1); 0016 plot(T(1:(i-1)),Y,'.');