poisson.m

来自「泊松随机过程」· M 代码 · 共 4 行

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