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 (テキサス大学オースティン校)、John J. Hasenbein (テキサス大学オースティン校)、Guodong Pang (ライス大学)分類 : math.PR (確率論)発表日 : 2025年11月3日 (arXivプレプリント)論文リンク : https://arxiv.org/abs/2510.27226v1 本論文は、単一サーバーキューイングシステムにおいて到着間隔時間とサービス時間が顧客待機時間に線形かつ確率的に依存する場合を研究し、待機時間過程の標本路径中偏差原理(MDP)を確立する。待機時間は確率的係数を有する修正Lindley再帰として表現される。確率係数の自然なスケーリング変換の下で、著者らはワークロード過程の流体挙動を分析し、安定平衡点(ゼロまたは正の値)を導出する。中偏差スケール過程は安定平衡点の周りで中心化され、その後、2つの確率遊歩によって駆動される線形確率微分方程式として表現される。これに漸近的に無視可能な誤差項と、ゼロでの可能な反射が加わる。両方の場合のMDPのレート関数は明示的に特性化でき、ゼロ中心化の場合は線形一般化Skorokhod反射写像を含むが、正中心化の場合は含まない点が異なる。
現実のキューイングシステムでは、到着過程とサービス時間はしばしばシステムの混雑度または遅延状態に依存する:
医療システム : 救急外来が過度に混雑している場合、患者は受診を放棄する(balking); 集中治療室が過負荷の場合、医師は患者の流転を加速させる可能性があるその他の応用 : 生物システム、製造業、在庫管理、コンピュータネットワーク、保険など多くの分野で同様の負荷依存行動が存在する理論的意義 : 古典的キューイング理論を状態依存システムに拡張し、このクラスのモデルにおける中偏差理論の空白を埋める実用的価値 : 混雑システムの稀な事象を理解するための理論的ツールを提供し、リスク評価とシステム設計に有用である方法論的貢献 : 確率係数を有する反射自己回帰過程を分析するための新しい技術を開発する分布分析が主流 : 既存研究の多くは定常分布と変換方法に焦点を当てている(Boxma et al. 2007, 2016, 2021)標本路径レベルの結果が限定的 :
Whitt (1990)の関数型中心極限定理は極限拡散過程の明示的形式を与えていない 大偏差原理は特殊な場合にのみ確立されている(Vlasiou and Palmowski 2014) 中偏差原理が完全に欠けている : これが本論文が埋める重要な空白である待機時間依存キューイングシステムの最初の標本路径中偏差原理を確立する 大偏差と中心極限定理の間にある精密な確率近似を提供する より広い反射AR(1)過程に適用可能な分析技術を開発する 最初のMDP結果 : 待機時間依存キューイングシステムの標本路径中偏差原理を確立し、この分野の理論的空白を埋める完全な流体分析 :異なるパラメータ領域(過負荷/臨界負荷/低負荷、異なる状態依存強度)における流体極限の挙動を系統的に分析 すべての安定平衡点(ゼロまたは正の値)を特定(表1に要約) 明示的なレート関数 : 2つの中心化ケース(ゼロ中心化と正中心化)に対して明示的に計算可能なレート関数を導出(定理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 はψ ∈ A C 0 \psi \in AC_0 ψ ∈ A C 0 に対して、そうでなければ∞ \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.総合評価 : これは高品質な理論数学論文であり、待機時間依存キューイングシステムの標本路径中偏差原理を厳密に確立し、この分野の重要な理論的空白を埋める。方法は革新的で、結果は明示的で、証明は完全である。主な不足は数値検証と応用討論の欠落、および特定の技術的仮定の制限性である。キューイング理論、大偏差理論、確率過程理論の研究者にとって重要な参考価値があり、実際のシステムのリスク分析にも理論的ツールを提供する。今後の研究では数値研究を補足し、実際の応用を探索することを推奨する。