We consider a single-server queue where interarrival and service times depend linearly and randomly on customer waiting times, and establish a sample-path moderate deviation principle (MDP) for the waiting time process. The waiting times for the queue can be written as a modified Lindley recursion with a random weight coefficient. Under a natural scaling of the random coefficients, we analyze the fluid behavior of the workload process and derive the stable equilibrium point, which can be zero or a positive value. The moderate-deviation-scaled process is centered around the stable equilibrium point and then represented as a linear stochastic differential equation driven by two random walks together with additional asymptotically negligible error terms and possibly a reflection at zero. The rate functions of MDPs in the two scenarios can be characterized explicitly, and they differ in that the case with zero centering term involves the linearly generalized Skorokhod reflection mapping while the case with positive centering term does not (similar to the corresponding diffusion limits). Our analysis involves the MDP for the associated linearly recursive Markov chains, invoking a perturbation of two independent random walks, and employing martingale techniques to prove the asymptotically exponentially vanishing error terms.
논문ID : 2510.27226제목 : Sample Path Moderate Deviation Principle for Queues with Waiting-time Dependent Interarrival and Service Times저자 : Chang Feng (University of Texas at Austin), John J. Hasenbein (University of Texas at Austin), Guodong Pang (Rice University)분류 : math.PR (확률론)발표일시 : 2025년 11월 3일 (arXiv 사전인쇄본)논문링크 : https://arxiv.org/abs/2510.27226v1 본 논문은 단일 서버 대기열 시스템에서 도착간격시간과 서비스시간이 고객 대기시간에 선형적이고 확률적으로 의존하는 경우를 연구하며, 대기시간 과정의 표본경로 중편차원리(MDP)를 확립한다. 대기시간은 확률적 가중계수를 갖는 수정된 Lindley 재귀로 표현될 수 있다. 확률계수의 자연 스케일 변환 하에서, 저자들은 작업부하 과정의 유체 거동을 분석하여 안정 평형점(영 또는 양수일 수 있음)을 도출한다. 중편차 스케일 과정은 안정 평형점 주위로 중심화되며, 두 개의 확률 보행으로 구동되는 선형 확률미분방정식으로 표현되고, 점근적으로 무시할 수 있는 오차항과 영점에서의 가능한 반사를 더한다. 두 경우 모두에서 MDP의 속도함수는 명시적으로 특성화될 수 있으며, 차이점은 영 중심화 경우가 선형 일반화된 Skorokhod 반사 사상을 포함하는 반면, 양 중심화 경우는 포함하지 않는다는 것이다.
현실의 대기열 시스템에서 도착 과정과 서비스시간은 종종 시스템 혼잡도 또는 지연 상태에 의존한다:
의료 시스템 : 응급실이 과도하게 혼잡할 때 환자들이 진료를 포기(balking)하고, 중환자실이 과부하일 때 의사들이 환자 이동 속도를 높일 수 있다기타 응용 : 생물 시스템, 제조업, 재고 관리, 컴퓨터 네트워크 및 보험 등의 분야에서도 유사한 부하 의존 행동이 존재한다이론적 의의 : 고전적 대기열 이론을 상태 의존 시스템으로 확장하고, 이러한 모델에서 중편차 이론의 공백을 메운다실용적 가치 : 혼잡 시스템의 희귀 사건을 이해하기 위한 이론적 도구를 제공하며, 위험 평가 및 시스템 설계에 도움이 된다방법론적 기여 : 확률적 계수를 갖는 반사 자기회귀 과정을 분석하기 위한 새로운 기법을 개발한다분포 분석 주도 : 기존 연구는 주로 정상 분포 및 변환 방법에 초점을 맞춤(Boxma et al. 2007, 2016, 2021)표본경로 수준 결과 제한 :
Whitt (1990)의 함수형 중심극한정리는 극한 확산 과정의 명시적 형태를 제공하지 않음 대편차원리는 특수한 경우에만 확립됨(Vlasiou and Palmowski 2014) 중편차원리 완전히 부재 : 이것이 본 논문이 메우는 핵심 공백대기시간 의존 대기열 시스템의 첫 표본경로 중편차원리 확립 대편차와 중심극한정리 사이의 정교한 확률 근사 제공 더 광범위한 반사 AR(1) 과정에 적용 가능한 분석 기법 개발 첫 MDP 결과 : 대기시간 의존 대기열 시스템의 표본경로 중편차원리를 확립하여 해당 분야의 이론적 공백을 메운다완전한 유체 분석 :다양한 매개변수 영역(과부하/임계 부하/저부하, 다양한 상태 의존 강도)에서 유체 극한의 거동을 체계적으로 분석 모든 안정 평형점(영 또는 양수) 식별(표1 참조) 명시적 속도함수 : 두 가지 중심화 경우(영 중심 및 양 중심)에 대해 명시적으로 계산 가능한 속도함수 도출(정리 2.6):영 중심화: 선형 일반화된 Skorokhod 반사 사상 포함 양 중심화: 반사 미포함, 속도함수 형태 더 단순 새로운 증명 기법 :선형 재귀 Markov 연쇄의 MDP 분석 방법 개발(제4절) 오차항의 지수 소멸성 증명을 위한 마팅게일 기법의 혁신적 사용 지수 타이트성 및 지수 동치성의 체계적 논증 프레임워크 확립 보완적 확산 근사 : 부록 B에서 함수형 중심극한정리를 증명하며, 극한은 OU 과정 또는 반사 OU 과정(Whitt (1990)의 작업 보완)단일 서버 FIFO 대기열 시스템 수열(인덱스 n)을 고려:
입력 : i.i.d. 확률 벡터 수열 { ( A i n , S i n , A i n , B i n ) , i ∈ N 0 } \{(A^n_i, S^n_i, \mathcal{A}^n_i, \mathcal{B}^n_i), i \in \mathbb{N}_0\} {( A i n , S i n , A i n , B i n ) , i ∈ N 0 } 상태 의존 메커니즘 :
실제 도착간격: A i ′ n = A i n + A i n W i n A'^{n}_i = A^n_i + \mathcal{A}^n_i W^n_i A i ′ n = A i n + A i n W i n 실제 서비스시간: S i ′ n = S i n + B i n W i n S'^{n}_i = S^n_i + \mathcal{B}^n_i W^n_i S i ′ n = S i n + B i n W i n 목표 : 중편차 스케일 대기시간 과정 W ~ n ( t ) = n b n ( W ˉ n ( t ) − W ˉ ∗ ) \tilde{W}^n(t) = \frac{\sqrt{n}}{b_n}(\bar{W}^n(t) - \bar{W}^*) W ~ n ( t ) = b n n ( W ˉ n ( t ) − W ˉ ∗ ) 의 표본경로 MDP 확립여기서:
b n → ∞ b_n \to \infty b n → ∞ , b n / n → 0 b_n/\sqrt{n} \to 0 b n / n → 0 (중편차 스케일)W ˉ n ( t ) = 1 n W ⌊ n t ⌋ n \bar{W}^n(t) = \frac{1}{n}W^n_{\lfloor nt \rfloor} W ˉ n ( t ) = n 1 W ⌊ n t ⌋ n (유체 스케일)W ˉ ∗ \bar{W}^* W ˉ ∗ 는 유체 극한의 안정 평형점수정된 Lindley 재귀 :
W i + 1 n = ( C i n W i n + X i n ) + W^n_{i+1} = (C^n_i W^n_i + X^n_i)^+ W i + 1 n = ( C i n W i n + X i n ) +
여기서:
X i n = S i n − A i n X^n_i = S^n_i - A^n_i X i n = S i n − A i n (명목 증분)C i n = 1 + B i n − A i n C^n_i = 1 + \mathcal{B}^n_i - \mathcal{A}^n_i C i n = 1 + B i n − A i n (확률적 가중계수)핵심 스케일 가정 (가정 2.1):
C i n = 1 − 1 n Θ i C^n_i = 1 - \frac{1}{n}\Theta_i C i n = 1 − n 1 Θ i
여기서 { Θ i } \{\Theta_i\} { Θ i } 는 i.i.d. 수열, 평균 θ \theta θ , 분산 σ Θ 2 \sigma^2_\Theta σ Θ 2 . 이러한 스케일 선택은:
n n n 이 클 때, ( C i ) 1 / n ≈ e − Θ i / n ≈ 1 − Θ i / n (C_i)^{1/n} \approx e^{-\Theta_i/n} \approx 1 - \Theta_i/n ( C i ) 1/ n ≈ e − Θ i / n ≈ 1 − Θ i / n Whitt (1990) 및 Boxma et al. (2016)의 스케일과 호환 망원급 합과 오차항 도입을 통해 유체 스케일 표현을 얻음:
W ˉ n ( t ) = W ˉ 0 n + 1 n ∑ i = 0 ⌊ n t ⌋ − 1 X i n − ∫ 0 t θ W ˉ n ( s ) d s + ϵ ˉ 1 n ( t ) + ϵ ˉ 2 n ( t ) + 1 n L ⌊ n t ⌋ − 1 n \bar{W}^n(t) = \bar{W}^n_0 + \frac{1}{n}\sum^{\lfloor nt \rfloor -1}_{i=0} X^n_i - \int^t_0 \theta \bar{W}^n(s)ds + \bar{\epsilon}^{n}_{1}(t) + \bar{\epsilon}^{n}_{2}(t) + \frac{1}{n}L^n_{\lfloor nt \rfloor -1} W ˉ n ( t ) = W ˉ 0 n + n 1 ∑ i = 0 ⌊ n t ⌋ − 1 X i n − ∫ 0 t θ W ˉ n ( s ) d s + ϵ ˉ 1 n ( t ) + ϵ ˉ 2 n ( t ) + n 1 L ⌊ n t ⌋ − 1 n
유체 극한 (정리 3.2):
W ˉ = R θ ( w ˉ 0 + μ e ) \bar{W} = R_\theta(\bar{w}_0 + \mu e) W ˉ = R θ ( w ˉ 0 + μ e )
여기서 R θ R_\theta R θ 는 선형 일반화된 Skorokhod 반사 사상이며, 미분 형태를 만족:
d W ˉ ( t ) = μ − θ W ˉ ( t ) + d L ˉ ( t ) d\bar{W}(t) = \mu - \theta\bar{W}(t) + d\bar{L}(t) d W ˉ ( t ) = μ − θ W ˉ ( t ) + d L ˉ ( t )
안정 평형점 분석 (표1 요약):
부하 조건 상태 의존 안정 평형점 W ˉ ∗ \bar{W}^* W ˉ ∗ μ > 0 \mu > 0 μ > 0 (과부하)θ > 0 \theta > 0 θ > 0 μ / θ \mu/\theta μ / θ μ = 0 \mu = 0 μ = 0 (임계)θ ≥ 0 \theta \geq 0 θ ≥ 0 0 μ < 0 \mu < 0 μ < 0 (저부하)임의의 θ \theta θ 0
반사 없는 보조 시스템을 처리하기 위해:
V i + 1 n = C i n V i n + X i n V^n_{i+1} = C^n_i V^n_i + X^n_i V i + 1 n = C i n V i n + X i n
중편차 스케일 표현 :
V ~ n ( t ) = V ~ n ( 0 ) + R ~ X n ( t ) − V ˉ ∗ R ~ Θ n ( t ) − ∫ 0 t θ V ~ n ( s ) d s + n ( μ n − μ ) t + ϵ ~ 1 , V n ( t ) + ϵ ~ 2 , V n ( t ) + ϵ ~ 3 , V n ( t ) \tilde{V}^n(t) = \tilde{V}^n(0) + \tilde{R}^n_X(t) - \bar{V}^* \tilde{R}^n_\Theta(t) - \int^t_0 \theta \tilde{V}^n(s)ds + \sqrt{n}(\mu_n - \mu)t + \tilde{\epsilon}^{n}_{1,V}(t) + \tilde{\epsilon}^{n}_{2,V}(t) + \tilde{\epsilon}^{n}_{3,V}(t) V ~ n ( t ) = V ~ n ( 0 ) + R ~ X n ( t ) − V ˉ ∗ R ~ Θ n ( t ) − ∫ 0 t θ V ~ n ( s ) d s + n ( μ n − μ ) t + ϵ ~ 1 , V n ( t ) + ϵ ~ 2 , V n ( t ) + ϵ ~ 3 , V n ( t )
여기서 R ~ X n \tilde{R}^n_X R ~ X n , R ~ Θ n \tilde{R}^n_\Theta R ~ Θ n 는 중편차 스케일의 확률 보행.
핵심 기술 단계 :
지수 경계 (보조정리 4.5): 유체 스케일 과정의 지수 확률 유계성 증명오차항 분석 (보조정리 4.6): 마팅게일 기법을 사용하여 ϵ ~ 2 , V n \tilde{\epsilon}^{n}_{2,V} ϵ ~ 2 , V n 의 지수 타이트성 증명지수 타이트성 (정리 4.7): 모든 항을 조합하여 { V ~ n } \{\tilde{V}^n\} { V ~ n } 의 지수 타이트성 증명지수 동치성 (보조정리 4.9): 오차항이 영 과정과 지수 동치임을 증명중편차 스케일 표현 :
W ~ n ( t ) = W ~ n ( 0 ) + R ~ X n ( t ) − W ˉ ∗ R ~ Θ n ( t ) − ∫ 0 t θ W ~ n ( s ) d s + n ( μ n − μ ) t + ϵ ~ 1 n ( t ) + ϵ ~ 2 n ( t ) + ϵ ~ 3 n ( t ) + L ~ n ( t ) \tilde{W}^n(t) = \tilde{W}^n(0) + \tilde{R}^n_X(t) - \bar{W}^* \tilde{R}^n_\Theta(t) - \int^t_0 \theta \tilde{W}^n(s)ds + \sqrt{n}(\mu_n - \mu)t + \tilde{\epsilon}^{n}_{1}(t) + \tilde{\epsilon}^{n}_{2}(t) + \tilde{\epsilon}^{n}_{3}(t) + \tilde{L}^n(t) W ~ n ( t ) = W ~ n ( 0 ) + R ~ X n ( t ) − W ˉ ∗ R ~ Θ n ( t ) − ∫ 0 t θ W ~ n ( s ) d s + n ( μ n − μ ) t + ϵ ~ 1 n ( t ) + ϵ ~ 2 n ( t ) + ϵ ~ 3 n ( t ) + L ~ n ( t )
핵심 혁신 :
보조 시스템 경계 (보조정리 5.1-5.2): 상하 경계 시스템 구성을 통해 W n W^n W n 제어반사항 분석 :
양 중심화 (W ˉ ∗ > 0 \bar{W}^* > 0 W ˉ ∗ > 0 ): L ~ n → P 1 / b n 2 0 \tilde{L}^n \xrightarrow{P^{1/b^2_n}} 0 L ~ n P 1/ b n 2 0 (반사 미활성화)영 중심화 (W ˉ ∗ = 0 \bar{W}^* = 0 W ˉ ∗ = 0 ): L ~ n \tilde{L}^n L ~ n 은 Skorokhod 사상의 조정 과정양 중심화 경우 (μ > 0 , θ > 0 , W ˉ ∗ = μ / θ \mu > 0, \theta > 0, \bar{W}^* = \mu/\theta μ > 0 , θ > 0 , W ˉ ∗ = μ / θ ):
I ( ϕ ) = θ 2 2 ( θ 2 σ X 2 + μ 2 σ Θ 2 ) ∫ 0 T ( ϕ ˙ ( t ) − r + θ ϕ ( t ) ) 2 d t I(\phi) = \frac{\theta^2}{2(\theta^2\sigma^2_X + \mu^2\sigma^2_\Theta)} \int^T_0 (\dot{\phi}(t) - r + \theta\phi(t))^2 dt I ( ϕ ) = 2 ( θ 2 σ X 2 + μ 2 σ Θ 2 ) θ 2 ∫ 0 T ( ϕ ˙ ( t ) − r + θϕ ( t ) ) 2 d t
ϕ ∈ A C \phi \in AC ϕ ∈ A C 이고 ϕ ( 0 ) = w 0 \phi(0) = w_0 ϕ ( 0 ) = w 0 인 경우.
영 중심화 경우 (μ = 0 , θ ≥ 0 , W ˉ ∗ = 0 \mu = 0, \theta \geq 0, \bar{W}^* = 0 μ = 0 , θ ≥ 0 , W ˉ ∗ = 0 ):
I ( ϕ ) = ∫ 0 T 1 { ϕ ( t ) > 0 } 1 2 σ X 2 ( ϕ ˙ ( t ) − r + θ ϕ ( t ) ) 2 d t + 1 2 σ X 2 r 2 ∫ 0 T 1 { ϕ ( t ) = 0 } 1 { r > 0 } d t I(\phi) = \int^T_0 \mathbb{1}_{\{\phi(t) > 0\}} \frac{1}{2\sigma^2_X}(\dot{\phi}(t) - r + \theta\phi(t))^2 dt + \frac{1}{2\sigma^2_X}r^2 \int^T_0 \mathbb{1}_{\{\phi(t)=0\}} \mathbb{1}_{\{r > 0\}} dt I ( ϕ ) = ∫ 0 T 1 { ϕ ( t ) > 0 } 2 σ X 2 1 ( ϕ ˙ ( t ) − r + θϕ ( t ) ) 2 d t + 2 σ X 2 1 r 2 ∫ 0 T 1 { ϕ ( t ) = 0 } 1 { r > 0 } d t
비음 ϕ ∈ A C \phi \in AC ϕ ∈ A C 이고 ϕ ( 0 ) = w 0 \phi(0) = w_0 ϕ ( 0 ) = w 0 인 경우.
Dupuis and Johnson (2015)의 약 수렴 방법과 달리, 본 논문은:
D T D_T D T 공간에서 직접 작업유한 차원 분포의 MDP와 지수 타이트성을 통해 결과 확립 축약 원리(contraction principle) 적용 오차항 ϵ ~ 2 n \tilde{\epsilon}^{n}_{2} ϵ ~ 2 n 에 대해(보조정리 4.6):
마팅게일 구성: Z k n = ∑ i = 0 k ( θ − Θ i ) ( V ˉ i n − V ˉ ∗ ) Z^n_k = \sum^k_{i=0} (\theta - \Theta_i)(\bar{V}^n_i - \bar{V}^*) Z k n = ∑ i = 0 k ( θ − Θ i ) ( V ˉ i n − V ˉ ∗ ) 상마팅게일 성질 증명: ζ k n = exp { b n n Z k n − b n 2 n K α 2 σ Θ 2 k } \zeta^n_k = \exp\{\frac{b_n}{\sqrt{n}}Z^n_k - \frac{b^2_n}{n}K^2_\alpha \sigma^2_\Theta k\} ζ k n = exp { n b n Z k n − n b n 2 K α 2 σ Θ 2 k } Doob 부등식 적용으로 지수 경계 획득 재귀적 논증 체인 확립(보조정리 4.5 → 4.6 → 정리 4.7 → 추론 4.8 → 보조정리 4.9):
각 단계에서 더 강한 결과 확립 이전 단계의 결론을 활용하여 후속 증명 단순화 최종적으로 모든 오차항이 영 과정과 지수 동치임을 증명 과정 Υ n \Upsilon^n Υ n 구성을 통해(보조정리 5.1):
Υ i n = max { 0 , X i − 1 n , X i − 1 n + C i − 1 n X i − 2 n , … , X i − 1 n + ⋯ + C 1 n X 0 n + C i − 1 n ⋯ C 0 n W 0 n } \Upsilon^n_i = \max\{0, X^n_{i-1}, X^n_{i-1} + C^n_{i-1}X^n_{i-2}, \ldots, X^n_{i-1} + \cdots + C^n_1 X^n_0 + C^n_{i-1}\cdots C^n_0 W^n_0\} Υ i n = max { 0 , X i − 1 n , X i − 1 n + C i − 1 n X i − 2 n , … , X i − 1 n + ⋯ + C 1 n X 0 n + C i − 1 n ⋯ C 0 n W 0 n }
0 ≤ W i n ≤ Υ i n 0 \leq W^n_i \leq \Upsilon^n_i 0 ≤ W i n ≤ Υ i n 를 증명한 후, Υ n \Upsilon^n Υ n 을 독립적으로 구성된 선형 재귀 시스템과 연결.
부록 C에서 체계적으로 정리:
지수 타이트성의 필요충분조건(정리 C.6, C.7) 연속 사상 하에서의 보존성(보조정리 C.8) 초지수 수렴의 특성화(보조정리 C.9, C.10) 확률적 시간 변환 정리(정리 C.11) 주 : 본 논문은 순수 이론 수학 논문으로 수치 실험이나 시뮬레이션을 포함하지 않는다. 모든 결과는 엄격한 수학 정리 및 그 증명이다.
논문은 다음 방식으로 이론의 합리성을 검증:
알려진 결과와의 일치성 :C i n = 1 − α / n C^n_i = 1 - \alpha/n C i n = 1 − α / n (결정론적)일 때, Boxma et al. (2016)의 결과 회복확산 극한(부록 B)은 Whitt (1990)의 정상 분포 근사와 일치(주석 B.4) 극단 경우 검증 :μ < 0 \mu < 0 μ < 0 일 때 W ~ n → P 1 / b n 2 0 \tilde{W}^n \xrightarrow{P^{1/b^2_n}} 0 W ~ n P 1/ b n 2 0 증명(시스템 안정성 매우 강함)다양한 매개변수 영역의 유체 극한 거동이 직관과 부합(그림 1-3) 내재적 일관성 :MDP 속도함수와 확산 극한의 공분산 구조 일치 영 중심화와 양 중심화의 속도함수 형태 차이는 반사 사상의 작용을 반영 정리 2.5(주요 MDP 결과) :
가정 2.1과 2.3 하에서, { W ~ n , n ∈ N } \{\tilde{W}^n, n \in \mathbb{N}\} { W ~ n , n ∈ N } 는 D T D_T D T 에서 속도 b n 2 b^2_n b n 2 의 MDP를 만족하며, 속도함수는:
양 중심화 (μ > 0 , θ > 0 , W ˉ ∗ = μ / θ \mu > 0, \theta > 0, \bar{W}^* = \mu/\theta μ > 0 , θ > 0 , W ˉ ∗ = μ / θ ):
I ( ϕ ) = inf ψ 1 , ψ 2 ∈ D T , ϕ = M θ ( w 0 + ψ 1 − μ θ ψ 2 + r e ) [ I X ( ψ 1 ) + I Θ ( ψ 2 ) ] I(\phi) = \inf_{\psi_1, \psi_2 \in D_T, \phi = M_\theta(w_0 + \psi_1 - \frac{\mu}{\theta}\psi_2 + re)} [I_X(\psi_1) + I_\Theta(\psi_2)] I ( ϕ ) = inf ψ 1 , ψ 2 ∈ D T , ϕ = M θ ( w 0 + ψ 1 − θ μ ψ 2 + re ) [ I X ( ψ 1 ) + I Θ ( ψ 2 )] 영 중심화 (μ = 0 , θ ≥ 0 , W ˉ ∗ = 0 \mu = 0, \theta \geq 0, \bar{W}^* = 0 μ = 0 , θ ≥ 0 , W ˉ ∗ = 0 ):
I ( ϕ ) = inf ψ 1 ∈ D T , ϕ = R θ ( w 0 + ψ 1 + r e ) I X ( ψ 1 ) I(\phi) = \inf_{\psi_1 \in D_T, \phi = R_\theta(w_0 + \psi_1 + re)} I_X(\psi_1) I ( ϕ ) = inf ψ 1 ∈ D T , ϕ = R θ ( w 0 + ψ 1 + re ) I X ( ψ 1 ) 여기서 I X ( ψ ) = 1 2 σ X 2 ∫ 0 T ∣ ψ ˙ ( t ) ∣ 2 d t I_X(\psi) = \frac{1}{2\sigma^2_X}\int^T_0 |\dot{\psi}(t)|^2 dt I X ( ψ ) = 2 σ X 2 1 ∫ 0 T ∣ ψ ˙ ( t ) ∣ 2 d t (AC0 _0 0 에서 ψ \psi ψ 에 대해), 그 외는 ∞ \infty ∞ .
정리 2.6(명시적 속도함수) :
최적화 문제는 명시적으로 해결 가능(위의 "속도함수의 명시적 형태" 참조).
정리 3.2 : W ˉ n → W ˉ \bar{W}^n \to \bar{W} W ˉ n → W ˉ u.o.c. in probability, 여기서 W ˉ = R θ ( w ˉ 0 + μ e ) \bar{W} = R_\theta(\bar{w}_0 + \mu e) W ˉ = R θ ( w ˉ 0 + μ e ) .
다양한 매개변수 영역의 안정성 (표1):
과부하+강 상태 의존 (μ > 0 , θ > 0 \mu > 0, \theta > 0 μ > 0 , θ > 0 ): 양수 μ / θ \mu/\theta μ / θ 에서 안정임계 부하+비음 상태 의존 (μ = 0 , θ ≥ 0 \mu = 0, \theta \geq 0 μ = 0 , θ ≥ 0 ): 0에서 안정저부하 (μ < 0 \mu < 0 μ < 0 ): 0에서 안정(임의의 θ \theta θ )불안정 경우 : μ > 0 , θ ≤ 0 \mu > 0, \theta \leq 0 μ > 0 , θ ≤ 0 또는 μ = 0 , θ < 0 \mu = 0, \theta < 0 μ = 0 , θ < 0 정리 4.3-4.4(선형 재귀 시스템의 MDP) :
반사 없는 시스템 V n V^n V n 에 대해 MDP를 확립하며, 속도함수 형태는 유사하지만 반사 사상을 포함하지 않음.
정리 B.3(함수형 중심극한정리) :
부록 B에서 증명:
양 중심화 : W ^ n ⇒ W ^ = M θ ( W ^ 0 + η e + σ X 2 + μ 2 θ 2 σ Θ 2 B ) \hat{W}^n \Rightarrow \hat{W} = M_\theta(\hat{W}_0 + \eta e + \sqrt{\sigma^2_X + \frac{\mu^2}{\theta^2}\sigma^2_\Theta}B) W ^ n ⇒ W ^ = M θ ( W ^ 0 + ηe + σ X 2 + θ 2 μ 2 σ Θ 2 B ) (OU 과정)영 중심화 : W ^ n ⇒ W ^ = R θ ( W ^ 0 + η e + σ X B ) \hat{W}^n \Rightarrow \hat{W} = R_\theta(\hat{W}_0 + \eta e + \sigma_X B) W ^ n ⇒ W ^ = R θ ( W ^ 0 + ηe + σ X B ) (반사 OU 과정)여기서 W ^ n ( t ) = n ( W ˉ n ( t ) − W ˉ ∗ ) \hat{W}^n(t) = \sqrt{n}(\bar{W}^n(t) - \bar{W}^*) W ^ n ( t ) = n ( W ˉ n ( t ) − W ˉ ∗ ) (확산 스케일).
보조정리 3.1 : 오차항 ϵ ˉ 1 n , ϵ ˉ 2 n → 0 \bar{\epsilon}^{n}_{1}, \bar{\epsilon}^{n}_{2} \to 0 ϵ ˉ 1 n , ϵ ˉ 2 n → 0 u.o.c. in probability보조정리 4.5 : lim K → ∞ lim sup n → ∞ 1 b n 2 log P ( ∥ V ˉ n ∥ T > K ) = − ∞ \lim_{K \to \infty} \limsup_{n \to \infty} \frac{1}{b^2_n} \log P(\|\bar{V}^n\|_T > K) = -\infty lim K → ∞ lim sup n → ∞ b n 2 1 log P ( ∥ V ˉ n ∥ T > K ) = − ∞ 보조정리 4.6 : { ϵ ~ 2 , V n } \{\tilde{\epsilon}^{n}_{2,V}\} { ϵ ~ 2 , V n } 는 D T D_T D T 에서 지수 타이트보조정리 5.2 : lim K → ∞ lim sup n → ∞ 1 b n 2 log P ( ∥ W ˉ n ∥ T > K ) = − ∞ \lim_{K \to \infty} \limsup_{n \to \infty} \frac{1}{b^2_n} \log P(\|\bar{W}^n\|_T > K) = -\infty lim K → ∞ lim sup n → ∞ b n 2 1 log P ( ∥ W ˉ n ∥ T > K ) = − ∞ 초기 연구 : Harris (1967), Callahan (1973), Brill (1988)이 기본 모델 연구최근 진전 :
Boxma and Vlasiou (2007): C i ∈ { − 1 , 1 } C_i \in \{-1, 1\} C i ∈ { − 1 , 1 } Bernoulli 경우 Boxma et al. (2016): 결정론적 C i C_i C i 의 반사 AR(1) 과정, FCLT Boxma et al. (2021), Huang (2023), Dimitriou and Fiems (2024): 더 일반적인 C i C_i C i 형태 Vlasiou and Palmowski (2014): 정상 분포 꼬리 확률의 대편차 본 논문의 기여 : 첫 표본경로 MDP 결과로 극한 정리 스펙트럼의 공백을 메운다.
표본경로 LDP :Bazhba et al. (2025): 표준 Lindley 재귀(C i = 1 C_i = 1 C i = 1 )의 부선형 속도 LDP Chen et al. (2024): 중꼬리 아핀 재귀의 LDP 표본경로 MDP :Puhalskii (1999): GI/GI/1 대기열 Chang et al. (1999): 장거리 의존 입력의 확률적 유체 대기열 Anugu and Pang (2024a): 시변 서비스시간의 무한 서버 대기열 Feng et al. (2025): GI/GI/1+GI 대기열 본 논문의 위치 : MDP 이론을 상태 의존 대기열 시스템으로 확장.
반사 없는 재귀 V i + 1 n = C i n V i n + X i n V^n_{i+1} = C^n_i V^n_i + X^n_i V i + 1 n = C i n V i n + X i n 는 문헌에서 Vervaat 영속성이라 불림:
Kesten (1973), Brandt (1986), Embrechts and Goldie (1994): 정상 분포 이론 Vervaat (1979), Whitt (1990): 함수형 극한 정리 Chen et al. (2024): 중꼬리 경우의 LDP 본 논문의 기여 : 이러한 과정에 대해 MDP 확립(제4절), 반사 과정 분석의 중간 단계로 작용.
Dupuis and Johnson (2015) : 변분 공식을 통한 약 수렴 방법, 특정 비선형 재귀 허용
차이점 : 본 논문은 직접 방법을 개발하여 반사항을 명시적으로 처리하고 D T D_T D T 공간에서 작업이론적 완전성 : 대기시간 의존 대기열 시스템에 대한 완전한 극한 정리 프레임워크 확립(유체 극한, 확산 극한, 중편차원리)속도함수의 이분성 :양 평형점: 속도함수 형태 간결, 반사 미포함 영 평형점: 속도함수가 Skorokhod 반사 사상 포함, 더 복잡 방법론적 기여 : 개발된 기법(마팅게일 방법, 지수 타이트성 논증, 보조 시스템 경계)을 더 광범위한 반사 확률 과정에 적용 가능매개변수 민감성 : 시스템 거동이 명목 부하 μ \mu μ 와 상태 의존 강도 θ \theta θ 에 매우 민감(표1 요약)모델 가정 :선형 의존 관계 필요(실제 시스템은 비선형일 수 있음) 특정 스케일 선택 C i n = 1 − Θ i / n C^n_i = 1 - \Theta_i/n C i n = 1 − Θ i / n (기존 문헌과 일치하지만) 독립성 가정(Θ i \Theta_i Θ i 와 X i n X^n_i X i n 독립)이 과할 수 있음 안정성 제한 : 안정 평형점 경우만 분석, 불안정 경우(예: μ > 0 , θ < 0 \mu > 0, \theta < 0 μ > 0 , θ < 0 ) 미포함단일 서버 제한 : 방법이 다중 서버 또는 네트워크 시스템으로 확장되지 않음기술적 조건 :적률 생성함수 존재 필요(가정 2.3(ii)) 초기 조건의 수렴 속도 요구(가정 2.3(i)) 수치 검증 부재 : 순수 이론 작업으로 수치 예제 또는 시뮬레이션 검증 미제공논문에서 제시된 잠재적 확장 :
비선형 의존 : A i ′ = A i + f ( W i ) A'_i = A_i + f(W_i) A i ′ = A i + f ( W i ) 형태의 비선형 상태 의존으로 확장다중 서버 시스템 : Halfin-Whitt 체제에서 상태 의존 다중 서버 대기열 분석네트워크 시스템 : Jackson 네트워크 또는 더 일반적인 대기열 네트워크로 확장중꼬리 분포 : X i n X^n_i X i n 또는 Θ i \Theta_i Θ i 가 중꼬리 분포를 갖는 경우 연구시변 매개변수 : μ n ( t ) \mu_n(t) μ n ( t ) , θ n ( t ) \theta_n(t) θ n ( t ) 가 시간에 따라 변하는 경우 고려최적화 응용 : MDP 결과를 활용한 희귀 사건 시뮬레이션 및 위험 관리증명 완전성 : 모든 정리가 상세한 증명을 가지며 기술적 세부사항이 충분논리 명확성 : 유체 분석에서 보조 시스템을 거쳐 주요 결과로 진행되는 층차가 분명부록 풍부 : 부록 A-C가 필요한 배경 및 기술 보조정리 제공직접 표본경로 방법 : 약 수렴 방법 대비 더 직관적이며 반사 처리에 적합마팅게일 기법의 교묘한 응용 : 지수 소멸성 증명 시 마팅게일 구성이 창의적체계적 프레임워크 : 확립된 지수 타이트성 논증 프레임워크가 보편적 적용 가능속도함수 완전 명시(정리 2.6)로 추가 분석 및 응용 용이 유체 극한 거동 완전 특성화(표1)로 모든 매개변수 영역 포함 구조 명확 : 서론이 충분한 동기 제공, 본문이 논리 엄밀, 부록이 강력한 지원기호 규범 : 기호 사용 일관(예: x ˉ n \bar{x}^n x ˉ n 은 유체 스케일, x ~ n \tilde{x}^n x ~ n 은 중편차 스케일)문헌 검토 포괄 : 관련 연구와의 연결 및 차이점을 충분히 논의수치 검증 부재 : 이론 결과의 실제 의미를 설명하는 수치 예제 부족응용 논의 부족 : MDP 결과를 실제 시스템 분석 또는 설계에 활용하는 방법 미논의매개변수 추정 문제 : 데이터로부터 θ \theta θ , σ Θ 2 \sigma^2_\Theta σ Θ 2 등 매개변수 추정 방법 미포함적률 생성함수 조건 (가정 2.3(ii))이 중꼬리 분포를 배제하여 응용 범위 제한스케일 선택의 유일성 : C i n = 1 − Θ i / n C^n_i = 1 - \Theta_i/n C i n = 1 − Θ i / n 이 "자연스러운" 선택인 이유 충분히 논의 부족불안정 경우 부재 : μ > 0 , θ ≤ 0 \mu > 0, \theta \leq 0 μ > 0 , θ ≤ 0 등 불안정 경우 미분석과도 거동 : 주로 점근 거동에 초점, 유한 n n n 의 근사 품질 미논의기술 밀도 높음 : 비전문가 독자에게 증명 세부사항이 따라가기 어려울 수 있음직관 설명 부족 : 보조 시스템 Υ n \Upsilon^n Υ n 구성 등 일부 기술 단계의 직관적 설명 부족이론적 공백 해소 : 상태 의존 대기열의 첫 MDP로 해당 분야 극한 정리 체계 완성방법론적 가치 : 개발된 기법을 다른 반사 확률 과정 분석에 적용 가능교량 역할 : 대기열 이론, 대편차 이론, 확률 과정 이론 연결위험 평가 : MDP가 대편차보다 정교한 희귀 사건 확률 추정 제공시스템 설계 : 속도함수를 활용한 편차 확률 제어 매개변수 최적화시뮬레이션 가속 : MDP 결과가 중요도 샘플링 등 분산 축소 기법 지도이론 결과 : 정리 진술 명확, 가정 명시, 원칙적으로 검증 가능증명 기법 : 단계 상세, 인용 충분, MDP 기법 학습 범례로 활용 가능코드 부재 : 순수 이론 작업으로 수치 구현 코드 미제공(이는 수용 가능)의료 시스템 : 응급실, ICU 등 혼잡 민감 의료 서비스 시스템콜센터 : 대기시간이 후속 도착 및 서비스 속도에 영향하는 시나리오컴퓨터 시스템 : 부하 의존 작업 스케줄링 및 서비스다중 서버 대기열 : 방법이 M / M / c M/M/c M / M / c 형태 상태 의존 시스템으로 확장 가능대기열 네트워크 : 복잡 네트워크 분석의 구성 요소로 활용금융 수학 : 유사 재귀 구조가 특정 금융 모델에 나타남반사 확산 과정 : 마팅게일 기법 및 지수 타이트성 논증을 다른 반사 과정에 적용상태 의존 시스템 : 보조 시스템 경계 방법의 보편적 적용성중편차 이론 : 다른 확률 시스템의 MDP 분석 범례 제공수치 연구 : 시뮬레이션으로 이론 예측 검증, 유한 n n n 의 수렴 속도 탐색매개변수 민감성 분석 : 속도함수의 θ \theta θ , σ X 2 \sigma^2_X σ X 2 등 매개변수 의존성 체계적 연구중꼬리 확장 : 가정 2.3(ii) 완화 가능성 연구응용 사례 : 구체적 응용 분야(예: 응급실 데이터) 선택하여 모델 검증계산 방법 : 속도함수 및 최적 경로 계산 수치 알고리즘 개발Whitt, W. (1990) . Queues with service times and interarrival times depending linearly and randomly upon waiting times. Queueing Systems , 6:335-351.Boxma, O., Mandjes, M., and Reed, J. (2016) . On a class of reflected AR(1) processes. Journal of Applied Probability , 53(3):818-832.Dupuis, P. and Johnson, D. (2015) . Moderate Deviations for Recursive Stochastic Algorithms. Stochastic Systems , 5(1):87-119.Puhalskii, A. A. (1999) . Moderate deviations for queues in critical loading. Queueing Systems , 31(3):359-392.Chen, B., Rhee, C.-H., and Zwart, B. (2024) . Sample-path large deviations for a class of heavy-tailed Markov additive processes. Electron. J. Probab. , 29(1):1-44.종합 평가 : 이는 고품질의 이론 수학 논문으로, 대기시간 의존 대기열 시스템의 표본경로 중편차원리를 엄밀하게 확립하여 해당 분야의 중요한 이론적 공백을 메운다. 방법론이 혁신적이고 결과가 명시적이며 증명이 완전하다. 주요 부족점은 수치 검증 및 응용 논의 부재와 일부 기술 가정의 제한성이다. 대기열 이론, 대편차 이론, 확률 과정 이론 연구자에게 중요한 참고 가치가 있으며, 실제 시스템의 위험 분석에도 이론적 도구를 제공한다. 향후 연구에서 수치 연구를 보충하고 실제 응용을 탐색할 것을 권장한다.