[논문 리뷰] Inter-arrival time distribution for the non-homogeneous Poisson process

G. Yakovlev, John B. Rundle외 2명|arXiv (Cornell University)|2005. 07. 27.

Advanced Frequency and Time Standards인용 수 11

한 줄 요약

이 논문은 유한 또는 무한 시간 간격에서 비균일한 포아송 과정(NHPP)의 도착 간격 분포에 대한 정확한 해석적 표현을 유도한다. 평균값 함수 Λ(t)를 활용하고 도착 시간 밀도 함수를 통합하여, 관측 창 효과를 고려한 닫힌 형태의 확률밀도함수(pdf)를 도출한다. 이는 과정의 유한 시간 메모리로 인해 지수 분포에서 벗어나는 현상을 보여준다.

ABSTRACT

We derive an analytical expression of the inter-arrival time distribution for a non-homogeneous Poisson process (NHPP). This expression is exact and is applicable to any time interval, finite or infinite. As an illustration, we present simulation results for three different intensity functions.

연구 동기 및 목표

• 유한 또는 무한 시간 간격에서 비균일 포아송 과정(NHPP)의 도착 간격 분포에 대한 정확한 해석적 표현을 유도하기 위해.
• 특히 유한 관측 창 하에서 도착 간격 분포에 대한 정확한 해가 부족한 문제를 해결하기 위해.
• 유한 관측 간격이 도입하는 기억 효과가 균일 포아송 과정에서 관찰되는 지수 분포와 어떻게 다를지 정량화하기 위해.
• 세 가지 다른 강도 함수(일정, 감소, 증가)에 대해 수치 시뮬레이션을 통해 유도된 공식을 검증하기 위해.

제안 방법

• 불연속 간격에서의 증분의 독립성을 이용하여, 가능한 모든 사건 수 n과 도착 시간 t_k에 대해 합을 구함으로써 생존함수 G(x) = P(S > x)를 유도한다.
• 수치 실험에서 도착 간격을 시뮬레이션하기 위해 역변환 방법을 사용하여 NHPP 강도 함수 λ(t)와 일관성을 확보한다.
• 시간을 누적 사건 빈도로 매핑하기 위해 평균값 함수 Λ(t) = ∫₀ᵗ λ(u)du를 적용하여 시간 간격 내 사건 확률의 정확한 계산을 가능하게 한다.
• 생존함수 G(x)를 미분하여 도착 간격 pdf를 유도하며, 이는 λ(t)와 Λ(t)를 포함한 이重적분 표현식이 된다.
• 유한 간격의 경계 효과를 고려하기 위해 첫 번째 도착 간격에 대해 별도의 항을 포함하여 분포의 정규화 및 완전성을 확보한다.
• 10⁸개의 실현값에 대해 세 가지 강도 함수 λ(t) = 1/c, c/(1+t), 및 ct에 대해 이론적 결과와 시뮬레이션 결과를 비교하여 분석적 해를 검증한다.

실험 결과

연구 질문

• RQ1 유한 시간 간격 T에서 관측된 비균일 포아송 과정에 대해 도착 간격의 정확한 분포는 무엇인가요?
• RQ2 유한 관측 창 T가 무한 시간 경우와 비교해 도착 간격 분포에 어떤 영향을 미치나요?
• RQ3 유한 간격으로 인해 도입되는 기억 효과가 균일 포아송 과정에서 관찰되는 지수 분포와 얼마나 다를까요?
• RQ4 시간에 따라 변하는 강도 함수와 관측 창의 경계 제약 조건을 모두 고려한 분석적 표현을 도출할 수 있을까요?

주요 결과

• 유한 간격 T에 대한 정확한 도착 간격 확률밀도함수(pdf)는 pdf(x) = (1/Λ(T)) ∫₀ᵀ⁻ˣ λ(y)λ(y+x)exp{−[Λ(y+x)−Λ(y)]}dy + [λ(x)exp{−Λ(x)}]/Λ(T)로 주어지며, 이는 경계 효과를 고려한다.
• 균일한 경우 λ(t) = 1/c에 대해 유도된 pdf는 pdf(x) = (T + c − x)/(cT) exp(−x/c)로 간소화되며, 이는 T가 유한할 경우 지수 분포에서 벗어나는 현상을 보여준다.
• λ(t) = 1/c, c/(1+t), 및 ct에 대해 시뮬레이션이 분석적 pdf와 우수한 일치를 보이며, 테스트된 모든 강도 함수에서 이론과 수치 결과 간의 일치가 뛰어나다.
• 유한 관측 창은 과정에 기억을 도입한다. 시간 y 이후에는 T−y보다 긴 도착 간격은 관측될 수 없기 때문이다. 이는 일정한 λ(t)일지라도 지수 분포가 아닐 수 있음을 의미한다.
• 유도된 표현은 정규화되어 있으며, G(0) = 1 및 G(T) = 0를 만족하여 생존함수 정의와 일관성을 확보한다.
• 이 방법은 균일 포아송 과정의 지수 도착 간격 분포를 임의의 결정적 강도 함수를 가진 비균일 경우로 일반화하는 데 성공했다.