poisson.m

泊松随机过程

function [T,NT] = poisson(Na,lam)
tau = randexp(Na,1,lam);
T = cumsum(tau);
NT = [1:length(T)]';