在Matalb上,我希望生成跟随Poisson到达的事件以及事件之间的到达间隔时间呈指数分布?
我应该使用哪个命令?
答案 0 :(得分:2)
假设您有一个计数过程{ N(t), t > = 0},其中 N ( t )是 t 发生的事件数。
如果这些事件遵循泊松过程(PP),且价格为 lambda ,那么我们有许多免费的属性。这些帖子假设均质(或固定)的 PP 。有关时间相关情况,请参见here和here和here。
PP 中的事件时间以速率 lambda 指数分布。
方法1:
1.产生间隔时间
2.获取实际事件时间<---您想要的
Rate = 0.75; % lambda (events per unit time)
NumEvents = 20; % number of events to generate
% Approach 1a requires exprnd function
InterEventTimes = exprnd(1/Rate,NumEvents,1);
EventTimes = cumsum(InterEventTimes);
% Approach 1b using inverse transform of exponential distribution
fh =@(n) - log(rand(n,1) / Rate; % inverse transform of exponential distribution
InterEventTimes2 = fh(NumEvents);
EventTimes2 = cumsum(InterEventTimes2);
这两种方法的验证均为available here。
答案 1 :(得分:1)
我不确定是否有功能,或者我是否正确理解了这个问题。
泊松过程的事件时间公式为
tau=1/lambda*ln(1/r)
其中lambda
是流程的强度,r
是单位统一分布,例如由rand
生成的分布。