[논문 리뷰] On the first-visit-time problem for birth and death processes with catastrophes

Antonio Di Crescenzo, Virginia Giorno외 2명|ArXiv.org|2003. 07. 15.

Advanced Queuing Theory Analysis참고 문헌 16인용 수 3

한 줄 요약

이 논문은 재해가 발생하는 생애-죽음 과정에 대한 첫 방문 시간 분포를 조사하며, 임의의 상태 $k \in \{r, r+1, \ldots\}$로의 첫 방문 시간 밀도의 라플라스 변환을 유도하고, 평균과 분산을 계산한다. 재해가 없는 등가 과정과의 연관을 통해 저자들은 첫 방문 시간 분포와 첫 재해 발생 시간에 대한 정확한 표현을 도출하며, 다양한 매개변수 구간에서 기대값과 정상분포에 대한 명시적 공식을 제공한다.

ABSTRACT

For a birth-death process subject to catastrophes, defined on the state-space $S=\{r,r+1,r+2,...\}$, with $r$ a positive integer or zero, the first-visit time to a state $k\in S$ is considered and the Laplace transform of its probability density function is determined, use of which is then made to obtain mean and variance. The Laplace transform of the probability density function of the first effective catastrophe occurrence time and its expected value are also obtained. Some extensions to time-non-homogeneous processes are then provided. Finally, certain additional results concerning the determination of the steady-state distribution and the representation of the transition probabilities are worked out, while some applications to particular birth-death processes are shown in the Appendix.

연구 동기 및 목표

• 재해가 발생하는 생애-죽음 과정에서 사전에 지정된 영이 아닌 상태로의 첫 방문 시간 분포를 결정하는 것.
• 상태 $k$로의 첫 방문 시간의 밀도 함수에 대한 라플라스 변환을 유도하고, 그 평균과 분산을 계산하는 것.
• 첫 재해 발생 시간을 분석하여 그 라플라스 변환과 기대값을 도출하는 것.
• 시간 비균일 과정으로 결과를 확장하고, 과정의 정상분포를 유도하는 것.
• 전이 확률의 명시적 표현을 제공하고, 부록에서 특정 생애-죽음 과정을 통해 프레임워크를 설명하는 것.

제안 방법

• 저자들은 재해가 있는 생애-죽음 과정과 재해가 없는 해당 과정 사이의 기능적 관계를 수립하여, 분석적 취급의 용이성을 위해 후자를 기준으로 삼는다.
• 두 과정 간의 관계를 활용하여 유도된 적분 방정식을 해결함으로써 첫 방문 시간 밀도의 라플라스 변환을 도출한다.
• 첫 재해 발생 시간에 대해서는 유사한 축소 기법을 적용하며, 포아송 분포의 기억 없음 성질을 활용한다.
• 정상분포는 콘플루언트 하이퍼지오메트릭 함수와 베타 함수를 사용하여 도출되며, 생애, 죽음 및 재해율의 상대적 크기에 따라 달라진다.
• 기준 과정의 평균 공식에 매개변수를 대입하여 $N(t)$의 기대값에 대한 명시적 표현을 얻는다.
• 이전 연구에서 사용된 생성함수 접근법을 수정하여 시간 비균일 과정으로 방법을 확장한다.

실험 결과

연구 질문

• RQ1 재해가 발생하는 생애-죽음 과정이 주어진 비흡수 상태 $k \geq r$에 도달하는 첫 시간의 분포는 무엇인가?
• RQ2 재해가 존재할 경우, 목표 상태로의 첫 방문 예상 시간이 재해가 없는 과정과 비교해 어떻게 달라지는가?
• RQ3 첫 재해 발생 시간의 라플라스 변환은 무엇이며, 그 기대값은 얼마인가?
• RQ4 재해가 있는 과정에서 정상분포가 존재하는 조건는 무엇이며, 그 명시적 형태는 무엇인가?
• RQ5 과정이 시간 비균일적일 경우 전이 확률과 첫 방문 시간 분포는 어떻게 변화하는가?

주요 결과

• 상태 $k$로의 첫 방문 시간 밀도의 라플라스 변환을 닫힌 형태로 도출하여, 첫 방문 시간의 평균과 분산을 정확하게 계산할 수 있다.
• 일반적으로 $\alpha \neq \beta + \xi$ 인 경우, $N(0) = j$일 때 $N(t)$의 기대값은 $j e^{(\alpha - \beta - \xi)t} + \frac{\nu (e^{(\alpha - \beta - \xi)t} - 1)}{\alpha - \beta - \xi}$ 로 주어지며, $\xi = 0$일 때는 재해가 없는 경우로 축소된다.
• 정상분포는 $\alpha < \beta + \xi$ 일 때 존재하며, $q_n = \xi \widehat{\pi}_{0,n}(\xi)$ 로 주어지며, $\widehat{\pi}_{0,n}(\lambda)$ 는 매개변수 구간에 따라 하이퍼지오메트릭 함수와 베타 함수를 통해 표현된다.
• $\alpha = \beta$ 인 경우, 정상분포는 제2종 콘플루언트 하이퍼지오메트릭 함수를 포함하며, $q_0 = \left(\frac{\xi}{\alpha}\right)^{\nu/\alpha} e^{\xi/\alpha} \Gamma\left(1 - \frac{\nu}{\alpha}, \frac{\xi}{\alpha}\right)$ 로 주어진다.
• $\alpha > \beta$ 인 경우, 정상분포는 하이퍼지오메트릭 함수 $F$ 와 베타 함수를 사용하여 표현되며, 생애율이 사망율을 초월하는 것을 반영하도록 매개변수를 조정한다.
• 첫 재해 발생 시간의 기대값은 재해율 $\xi$ 의 역수로 도출되며, 포아송 과정의 기억 없음 성질과 일치한다.