2025-11-13T20:07:10.997459

Sample Path Moderate Deviation Principle for Queues with Waiting-time Dependent Interarrival and Service Times

Feng, Hasenbein, Pang
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.
academic

Principio di Deviazione Moderata su Traiettorie Campionarie per Code con Tempi Interarrivo e Servizio Dipendenti dal Tempo di Attesa

Informazioni Fondamentali

  • ID Articolo: 2510.27226
  • Titolo: Sample Path Moderate Deviation Principle for Queues with Waiting-time Dependent Interarrival and Service Times
  • Autori: Chang Feng (University of Texas at Austin), John J. Hasenbein (University of Texas at Austin), Guodong Pang (Rice University)
  • Classificazione: math.PR (Teoria della Probabilità)
  • Data di Pubblicazione: 3 novembre 2025 (preprint arXiv)
  • Link dell'Articolo: https://arxiv.org/abs/2510.27226v1

Riassunto

Questo articolo studia sistemi di code a singolo servitore nei quali i tempi interarrivo e i tempi di servizio dipendono linearmente e stocasticamente dal tempo di attesa dei clienti. Gli autori stabiliscono il principio di deviazione moderata (MDP) su traiettorie campionarie per il processo dei tempi di attesa. Il tempo di attesa può essere rappresentato come una ricorsione di Lindley modificata con coefficienti di peso stocastici. Sotto una riscalatura naturale dei coefficienti stocastici, gli autori analizzano il comportamento fluido del processo di carico di lavoro, derivando punti di equilibrio stabile (che possono essere zero o positivi). I processi su scala di deviazione moderata sono centrati attorno ai punti di equilibrio stabile e quindi rappresentati come equazioni differenziali stocastiche lineari guidate da due passeggiate casuali, più termini di errore asintoticamente trascurabili e una possibile riflessione nello zero. Le funzioni di velocità dell'MDP in entrambi i casi possono essere caratterizzate esplicitamente, con la differenza che il caso centrato nello zero coinvolge la mappa di riflessione di Skorokhod generalizzata lineare, mentre il caso centrato positivamente no.

Contesto di Ricerca e Motivazione

Formulazione del Problema

Nei sistemi di code reali, i processi di arrivo e i tempi di servizio spesso dipendono dallo stato di congestione o ritardo del sistema:

  • Sistemi Medici: quando i pronto soccorso sono sovraccarichi, i pazienti rinunciano alle cure (balking); quando le unità di terapia intensiva sono al limite, i medici possono accelerare il trasferimento dei pazienti
  • Altre Applicazioni: sistemi biologici, manifatturiero, gestione dell'inventario, reti di computer e assicurazioni presentano tutti comportamenti simili dipendenti dal carico

Importanza

  1. Significato Teorico: estende la teoria classica delle code ai sistemi dipendenti dallo stato, colmando il vuoto nella teoria della deviazione moderata per questi modelli
  2. Valore Pratico: fornisce strumenti teorici per comprendere gli eventi rari nei sistemi congestioni, utile per la valutazione del rischio e la progettazione del sistema
  3. Contributo Metodologico: sviluppa nuove tecniche per l'analisi di processi autoregressivi riflessi con coefficienti stocastici

Limitazioni dei Metodi Esistenti

  • Predominanza dell'Analisi Distributiva: la ricerca esistente si concentra principalmente su distribuzioni stazionarie e metodi di trasformazione (Boxma et al. 2007, 2016, 2021)
  • Risultati Limitati a Livello di Traiettoria Campionaria:
    • Il teorema del limite centrale funzionale di Whitt (1990) non fornisce la forma esplicita del processo di diffusione limite
    • Il principio di grande deviazione è stabilito solo in casi speciali (Vlasiou e Palmowski 2014)
    • Assenza Completa del Principio di Deviazione Moderata: questo è il vuoto chiave che questo articolo colma

Motivazione della Ricerca

  1. Stabilire il primo principio di deviazione moderata su traiettoria campionaria per sistemi di code dipendenti dal tempo di attesa
  2. Fornire un'approssimazione probabilistica fine tra la grande deviazione e il teorema del limite centrale
  3. Sviluppare tecniche analitiche applicabili a processi autoregressivi riflessi più generali

Contributi Principali

  1. Primo Risultato MDP: stabilisce il principio di deviazione moderata su traiettoria campionaria per sistemi di code dipendenti dal tempo di attesa, colmando il vuoto teorico in questo campo
  2. Analisi Fluida Completa:
    • Analizza sistematicamente il comportamento del limite fluido in diverse regioni di parametri (sovraccarico/carico critico/sottocarico, diverse intensità di dipendenza dallo stato)
    • Identifica tutti i punti di equilibrio stabile (zero o positivi), riassunti nella Tabella 1
  3. Funzioni di Velocità Esplicite: per i due casi di centralizzazione (centrato nello zero e centrato positivamente) derivano funzioni di velocità esplicitamente calcolabili (Teorema 2.6):
    • Centralizzazione nello zero: coinvolge la mappa di riflessione di Skorokhod generalizzata lineare
    • Centralizzazione positiva: non coinvolge riflessione, forma della funzione di velocità più semplice
  4. Nuove Tecniche di Prova:
    • Sviluppa il metodo di analisi MDP per catene di Markov ricorsive lineari (Sezione 4)
    • Utilizza innovativamente tecniche di martingala per provare la scomparsa esponenziale dei termini di errore
    • Stabilisce un framework sistematico di argomentazione per la compattezza esponenziale e l'equivalenza esponenziale
  5. Approssimazione di Diffusione Complementare: nell'Appendice B prova il teorema del limite centrale funzionale con limite come processo OU o processo OU riflesso, completando il lavoro di Whitt (1990)

Dettagli Metodologici

Definizione del Compito

Considerare una sequenza di sistemi di code FIFO a singolo servitore (indicizzati da n):

  • Input: sequenza i.i.d. di vettori casuali {(Ain,Sin,Ain,Bin),iN0}\{(A^n_i, S^n_i, \mathcal{A}^n_i, \mathcal{B}^n_i), i \in \mathbb{N}_0\}
  • Meccanismo di Dipendenza dallo Stato:
    • Intervallo di arrivo effettivo: Ain=Ain+AinWinA'^{n}_i = A^n_i + \mathcal{A}^n_i W^n_i
    • Tempo di servizio effettivo: Sin=Sin+BinWinS'^{n}_i = S^n_i + \mathcal{B}^n_i W^n_i
  • Obiettivo: stabilire il principio di deviazione moderata su traiettoria campionaria per il processo dei tempi di attesa su scala di deviazione moderata W~n(t)=nbn(Wˉn(t)Wˉ)\tilde{W}^n(t) = \frac{\sqrt{n}}{b_n}(\bar{W}^n(t) - \bar{W}^*)

Dove:

  • bnb_n \to \infty, bn/n0b_n/\sqrt{n} \to 0 (scala di deviazione moderata)
  • Wˉn(t)=1nWntn\bar{W}^n(t) = \frac{1}{n}W^n_{\lfloor nt \rfloor} (scala fluida)
  • Wˉ\bar{W}^* è il punto di equilibrio stabile del limite fluido

Struttura Ricorsiva Centrale

Ricorsione di Lindley Modificata: Wi+1n=(CinWin+Xin)+W^n_{i+1} = (C^n_i W^n_i + X^n_i)^+

Dove:

  • Xin=SinAinX^n_i = S^n_i - A^n_i (incremento nominale)
  • Cin=1+BinAinC^n_i = 1 + \mathcal{B}^n_i - \mathcal{A}^n_i (coefficiente di peso stocastico)

Assunzione di Scala Chiave (Assunzione 2.1): Cin=11nΘiC^n_i = 1 - \frac{1}{n}\Theta_i

Dove {Θi}\{\Theta_i\} è una sequenza i.i.d. con media θ\theta e varianza σΘ2\sigma^2_\Theta. Questa scelta di scala assicura che:

  • Quando nn è grande, (Ci)1/neΘi/n1Θi/n(C_i)^{1/n} \approx e^{-\Theta_i/n} \approx 1 - \Theta_i/n
  • È compatibile con la scala di Whitt (1990) e Boxma et al. (2016)

Analisi Fluida (Sezione 3)

Attraverso somme telescopiche e introducendo termini di errore, si ottiene la rappresentazione su scala fluida: Wˉn(t)=Wˉ0n+1ni=0nt1Xin0tθWˉn(s)ds+ϵˉ1n(t)+ϵˉ2n(t)+1nLnt1n\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}

Limite Fluido (Teorema 3.2): Wˉ=Rθ(wˉ0+μe)\bar{W} = R_\theta(\bar{w}_0 + \mu e)

Dove RθR_\theta è la mappa di riflessione di Skorokhod generalizzata lineare, che soddisfa la forma differenziale: dWˉ(t)=μθWˉ(t)+dLˉ(t)d\bar{W}(t) = \mu - \theta\bar{W}(t) + d\bar{L}(t)

Analisi del Punto di Equilibrio Stabile (Tabella 1 riassume):

Condizione di CaricoDipendenza dallo StatoPunto di Equilibrio Stabile Wˉ\bar{W}^*
μ>0\mu > 0 (sovraccarico)θ>0\theta > 0μ/θ\mu/\theta
μ=0\mu = 0 (critico)θ0\theta \geq 00
μ<0\mu < 0 (sottocarico)qualsiasi θ\theta0

MDP per Catena di Markov Ricorsiva Lineare (Sezione 4)

Per gestire il sistema ausiliario senza riflessione: Vi+1n=CinVin+XinV^n_{i+1} = C^n_i V^n_i + X^n_i

Rappresentazione su Scala di Deviazione Moderata: V~n(t)=V~n(0)+R~Xn(t)VˉR~Θn(t)0tθV~n(s)ds+n(μnμ)t+ϵ~1,Vn(t)+ϵ~2,Vn(t)+ϵ~3,Vn(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)

Dove R~Xn\tilde{R}^n_X, R~Θn\tilde{R}^n_\Theta sono passeggiate casuali su scala di deviazione moderata.

Passaggi Tecnici Chiave:

  1. Limiti Esponenziali (Lemma 4.5): prova il limite casuale esponenziale del processo su scala fluida
  2. Analisi dei Termini di Errore (Lemma 4.6): utilizza tecniche di martingala per provare la compattezza esponenziale di ϵ~2,Vn\tilde{\epsilon}^{n}_{2,V}
  3. Compattezza Esponenziale (Teorema 4.7): combina tutti i termini per provare la compattezza esponenziale di {V~n}\{\tilde{V}^n\}
  4. Equivalenza Esponenziale (Lemma 4.9): prova che i termini di errore sono esponenzialmente equivalenti al processo zero

MDP per il Processo dei Tempi di Attesa (Sezione 5)

Rappresentazione su Scala di Deviazione Moderata: W~n(t)=W~n(0)+R~Xn(t)WˉR~Θn(t)0tθW~n(s)ds+n(μnμ)t+ϵ~1n(t)+ϵ~2n(t)+ϵ~3n(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)

Innovazioni Chiave:

  1. Limiti del Sistema Ausiliario (Lemma 5.1-5.2): costruisce sistemi di limite superiore e inferiore per controllare WnW^n
  2. Analisi del Termine di Riflessione:
    • Centralizzazione Positiva (Wˉ>0\bar{W}^* > 0): L~nP1/bn20\tilde{L}^n \xrightarrow{P^{1/b^2_n}} 0 (riflessione non attivata)
    • Centralizzazione Nello Zero (Wˉ=0\bar{W}^* = 0): L~n\tilde{L}^n come processo di aggiustamento della mappa di Skorokhod

Forma Esplicita della Funzione di Velocità (Teorema 2.6)

Caso di Centralizzazione Positiva (μ>0,θ>0,Wˉ=μ/θ\mu > 0, \theta > 0, \bar{W}^* = \mu/\theta): I(ϕ)=θ22(θ2σX2+μ2σΘ2)0T(ϕ˙(t)r+θϕ(t))2dtI(\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

Per ϕAC\phi \in AC e ϕ(0)=w0\phi(0) = w_0.

Caso di Centralizzazione Nello Zero (μ=0,θ0,Wˉ=0\mu = 0, \theta \geq 0, \bar{W}^* = 0): I(ϕ)=0T1{ϕ(t)>0}12σX2(ϕ˙(t)r+θϕ(t))2dt+12σX2r20T1{ϕ(t)=0}1{r>0}dtI(\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

Per ϕ\phi non negativo in ACAC e ϕ(0)=w0\phi(0) = w_0.

Punti di Innovazione Tecnica

1. Metodo Diretto su Traiettoria Campionaria

Diversamente dal metodo di convergenza debole di Dupuis e Johnson (2015), questo articolo:

  • Lavora direttamente nello spazio DTD_T
  • Stabilisce il risultato attraverso MDP di distribuzioni finite-dimensionali e compattezza esponenziale
  • Applica il principio di contrazione

2. Applicazione Innovativa di Tecniche di Martingala

Per il termine di errore ϵ~2n\tilde{\epsilon}^{n}_{2} (Lemma 4.6):

  • Costruisce una martingala: Zkn=i=0k(θΘi)(VˉinVˉ)Z^n_k = \sum^k_{i=0} (\theta - \Theta_i)(\bar{V}^n_i - \bar{V}^*)
  • Prova la proprietà di supermartingala: ζkn=exp{bnnZknbn2nKα2σΘ2k}\zeta^n_k = \exp\{\frac{b_n}{\sqrt{n}}Z^n_k - \frac{b^2_n}{n}K^2_\alpha \sigma^2_\Theta k\}
  • Applica la disuguaglianza di Doob per ottenere limiti esponenziali

3. Argomentazione Sistematica di Equivalenza Esponenziale

Stabilisce una catena di argomentazione progressiva (Lemma 4.5 → 4.6 → Teorema 4.7 → Corollario 4.8 → Lemma 4.9):

  • Ogni passo stabilisce un risultato più forte
  • Utilizza le conclusioni del passo precedente per semplificare le prove successive
  • Infine prova che tutti i termini di errore sono esponenzialmente equivalenti al processo zero

4. Tecnica di Limite del Sistema Ausiliario

Attraverso la costruzione del processo Υn\Upsilon^n (Lemma 5.1): Υin=max{0,Xi1n,Xi1n+Ci1nXi2n,,Xi1n++C1nX0n+Ci1nC0nW0n}\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\}

Prova che 0WinΥin0 \leq W^n_i \leq \Upsilon^n_i, quindi collega Υn\Upsilon^n a un sistema ricorsivo lineare costruito indipendentemente.

5. Caratterizzazione della Compattezza C-Esponenziale

Nell'Appendice C sistematizza:

  • Condizioni necessarie e sufficienti per la compattezza esponenziale (Teoremi C.6, C.7)
  • Preservazione sotto mappe continue (Lemma C.8)
  • Caratterizzazione della convergenza superesponenziale (Lemmi C.9, C.10)
  • Teorema di trasformazione di tempo casuale (Teorema C.11)

Configurazione Sperimentale

Nota: Questo articolo è un lavoro di matematica pura teorica e non contiene esperimenti numerici o simulazioni. Tutti i risultati sono teoremi matematici rigorosi e loro prove.

Framework di Verifica Teorica

L'articolo verifica la ragionevolezza della teoria nei seguenti modi:

  1. Coerenza con Risultati Noti:
    • Quando Cin=1α/nC^n_i = 1 - \alpha/n (deterministico), recupera i risultati di Boxma et al. (2016)
    • Il limite di diffusione (Appendice B) è coerente con l'approssimazione della distribuzione stazionaria di Whitt (1990) (Osservazione B.4)
  2. Verifica di Casi Estremi:
    • Quando μ<0\mu < 0 prova che W~nP1/bn20\tilde{W}^n \xrightarrow{P^{1/b^2_n}} 0 (stabilità del sistema molto forte)
    • Il comportamento del limite fluido in diverse regioni di parametri è coerente con l'intuizione (Figure 1-3)
  3. Coerenza Interna:
    • La funzione di velocità dell'MDP è coerente con la struttura di covarianza del limite di diffusione
    • La differenza di forma della funzione di velocità tra centralizzazione nello zero e positiva riflette l'effetto della mappa di riflessione

Risultati Sperimentali

Risultati Teorici Principali

Teorema 2.5 (Risultato MDP Principale): Sotto le Assunzioni 2.1 e 2.3, {W~n,nN}\{\tilde{W}^n, n \in \mathbb{N}\} soddisfa l'MDP in DTD_T con velocità bn2b^2_n e funzione di velocità:

  1. Centralizzazione Positiva (μ>0,θ>0,Wˉ=μ/θ\mu > 0, \theta > 0, \bar{W}^* = \mu/\theta): I(ϕ)=infψ1,ψ2DT,ϕ=Mθ(w0+ψ1μθψ2+re)[IX(ψ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)]
  2. Centralizzazione Nello Zero (μ=0,θ0,Wˉ=0\mu = 0, \theta \geq 0, \bar{W}^* = 0): I(ϕ)=infψ1DT,ϕ=Rθ(w0+ψ1+re)IX(ψ1)I(\phi) = \inf_{\psi_1 \in D_T, \phi = R_\theta(w_0 + \psi_1 + re)} I_X(\psi_1)

Dove IX(ψ)=12σX20Tψ˙(t)2dtI_X(\psi) = \frac{1}{2\sigma^2_X}\int^T_0 |\dot{\psi}(t)|^2 dt per ψAC0\psi \in AC_0, altrimenti \infty.

Teorema 2.6 (Funzione di Velocità Esplicita): I problemi di ottimizzazione possono essere risolti esplicitamente (vedi sopra "Forma Esplicita della Funzione di Velocità").

Comportamento del Limite Fluido (Sezione 3)

Teorema 3.2: WˉnWˉ\bar{W}^n \to \bar{W} u.o.c. in probabilità, dove Wˉ=Rθ(wˉ0+μe)\bar{W} = R_\theta(\bar{w}_0 + \mu e).

Stabilità in Diverse Regioni di Parametri (Tabella 1):

  • Sovraccarico + Forte Dipendenza dallo Stato (μ>0,θ>0\mu > 0, \theta > 0): stabile a valore positivo μ/θ\mu/\theta
  • Carico Critico + Dipendenza dallo Stato Non Negativa (μ=0,θ0\mu = 0, \theta \geq 0): stabile a 0
  • Sottocarico (μ<0\mu < 0): stabile a 0 (qualsiasi θ\theta)
  • Casi Instabili: μ>0,θ0\mu > 0, \theta \leq 0 o μ=0,θ<0\mu = 0, \theta < 0

Risultati Ausiliari

Teoremi 4.3-4.4 (MDP per Sistema Ricorsivo Lineare): Stabilisce MDP per il sistema senza riflessione VnV^n, con forma della funzione di velocità simile ma senza coinvolgere la mappa di riflessione.

Teorema B.3 (Teorema del Limite Centrale Funzionale): Provato nell'Appendice B:

  • Centralizzazione Positiva: W^nW^=Mθ(W^0+ηe+σX2+μ2θ2σΘ2B)\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) (processo OU)
  • Centralizzazione Nello Zero: W^nW^=Rθ(W^0+ηe+σXB)\hat{W}^n \Rightarrow \hat{W} = R_\theta(\hat{W}_0 + \eta e + \sigma_X B) (processo OU riflesso)

Dove W^n(t)=n(Wˉn(t)Wˉ)\hat{W}^n(t) = \sqrt{n}(\bar{W}^n(t) - \bar{W}^*) (scala di diffusione).

Lemmi Chiave

  1. Lemma 3.1: Termini di errore ϵˉ1n,ϵˉ2n0\bar{\epsilon}^{n}_{1}, \bar{\epsilon}^{n}_{2} \to 0 u.o.c. in probabilità
  2. Lemma 4.5: limKlim supn1bn2logP(VˉnT>K)=\lim_{K \to \infty} \limsup_{n \to \infty} \frac{1}{b^2_n} \log P(\|\bar{V}^n\|_T > K) = -\infty
  3. Lemma 4.6: {ϵ~2,Vn}\{\tilde{\epsilon}^{n}_{2,V}\} è esponenzialmente compatto in DTD_T
  4. Lemma 5.2: limKlim supn1bn2logP(WˉnT>K)=\lim_{K \to \infty} \limsup_{n \to \infty} \frac{1}{b^2_n} \log P(\|\bar{W}^n\|_T > K) = -\infty

Lavori Correlati

Sistemi di Code Dipendenti dal Tempo di Attesa

  • Lavori Iniziali: Harris (1967), Callahan (1973), Brill (1988) hanno studiato modelli fondamentali
  • Progressi Recenti:
    • Boxma e Vlasiou (2007): caso Bernoulli Ci{1,1}C_i \in \{-1, 1\}
    • Boxma et al. (2016): processo AR(1) riflesso con CiC_i deterministico, FCLT
    • Boxma et al. (2021), Huang (2023), Dimitriou e Fiems (2024): forme più generali di CiC_i
    • Vlasiou e Palmowski (2014): grande deviazione della coda della distribuzione stazionaria

Contributo di Questo Articolo: primo risultato MDP su traiettoria campionaria, colmando il vuoto nello spettro dei teoremi limite.

Teoria della Grande Deviazione e Deviazione Moderata

  • LDP su Traiettoria Campionaria:
    • Bazhba et al. (2025): LDP su scala sublineare per ricorsione di Lindley ordinaria (Ci=1C_i = 1)
    • Chen et al. (2024): LDP per ricorsione affine con coda pesante
  • MDP su Traiettoria Campionaria:
    • Puhalskii (1999): coda GI/GI/1
    • Chang et al. (1999): flusso fluido casuale con input a dipendenza lunga
    • Anugu e Pang (2024a): coda con infiniti servitori e tempo di servizio variabile nel tempo
    • Feng et al. (2025): coda GI/GI/1+GI

Posizionamento di Questo Articolo: estende la teoria MDP ai sistemi di code dipendenti dallo stato.

Perpetuità di Vervaat

La ricorsione senza riflessione Vi+1n=CinVin+XinV^n_{i+1} = C^n_i V^n_i + X^n_i è chiamata perpetuità di Vervaat nella letteratura:

  • Kesten (1973), Brandt (1986), Embrechts e Goldie (1994): teoria della distribuzione stazionaria
  • Vervaat (1979), Whitt (1990): teoremi del limite funzionale
  • Chen et al. (2024): LDP per caso con coda pesante

Contributo di Questo Articolo: stabilisce MDP per questo tipo di processo (Sezione 4), come passo intermedio per l'analisi del processo riflesso.

MDP per Sistemi Markov Ricorsivi

  • Dupuis e Johnson (2015): metodo di convergenza debole attraverso formula variazionale, consente alcune ricorsioni non lineari
    • Differenza: questo articolo sviluppa il metodo diretto, gestisce esplicitamente il termine di riflessione, lavora nello spazio DTD_T

Conclusioni e Discussione

Conclusioni Principali

  1. Completezza Teorica: stabilisce un framework completo di teoremi limite per sistemi di code dipendenti dal tempo di attesa (limite fluido, limite di diffusione, principio di deviazione moderata)
  2. Dicotomia della Funzione di Velocità:
    • Punto di equilibrio positivo: forma della funzione di velocità semplice, non coinvolge riflessione
    • Punto di equilibrio zero: funzione di velocità coinvolge la mappa di riflessione di Skorokhod, più complessa
  3. Contributo Metodologico: le tecniche sviluppate (metodo di martingala, argomentazione di compattezza esponenziale, limite del sistema ausiliario) possono essere applicate a processi casuali riflessi più generali
  4. Sensibilità ai Parametri: il comportamento del sistema è altamente sensibile al carico nominale μ\mu e all'intensità della dipendenza dallo stato θ\theta (riassunto nella Tabella 1)

Limitazioni

  1. Assunzioni del Modello:
    • Richiede relazione di dipendenza lineare (i sistemi reali possono essere non lineari)
    • Scelta di scala specifica Cin=1Θi/nC^n_i = 1 - \Theta_i/n (sebbene coerente con la letteratura esistente)
    • L'assunzione di indipendenza (Θi\Theta_i indipendente da XinX^n_i) può essere troppo forte
  2. Limitazione di Stabilità: analizza solo il caso di punto di equilibrio stabile, non copre il caso instabile (come μ>0,θ<0\mu > 0, \theta < 0)
  3. Limitazione a Singolo Servitore: il metodo non si estende a sistemi multi-servitore o di rete
  4. Condizioni Tecniche:
    • Richiede l'esistenza della funzione generatrice dei momenti (Assunzione 2.3(ii))
    • Richiede velocità di convergenza della condizione iniziale (Assunzione 2.3(i))
  5. Assenza di Verifica Numerica: come lavoro puramente teorico, non fornisce esempi numerici o verifiche di simulazione

Direzioni Future

Possibili Estensioni Proposte dall'Articolo:

  1. Dipendenza Non Lineare: estendere a dipendenza non lineare dallo stato Ai=Ai+f(Wi)A'_i = A_i + f(W_i)
  2. Sistemi Multi-Servitore: analizzare code multi-servitore dipendenti dallo stato nel regime di Halfin-Whitt
  3. Sistemi di Rete: estendere a reti di Jackson o reti di code più generali
  4. Distribuzioni con Coda Pesante: studiare il caso in cui XinX^n_i o Θi\Theta_i hanno distribuzioni con coda pesante
  5. Parametri Variabili nel Tempo: considerare μn(t)\mu_n(t), θn(t)\theta_n(t) che variano nel tempo
  6. Applicazioni di Ottimizzazione: utilizzare i risultati MDP per simulazione di eventi rari e gestione del rischio

Valutazione Approfondita

Punti di Forza

1. Rigore Teorico

  • Prove Complete: tutti i teoremi hanno prove dettagliate, dettagli tecnici sufficienti
  • Logica Chiara: dall'analisi fluida al sistema ausiliario al risultato principale, livelli ben definiti
  • Appendici Ricche: le Appendici A-C forniscono background necessario e lemmi tecnici

2. Innovazione Metodologica

  • Metodo Diretto su Traiettoria Campionaria: rispetto al metodo di convergenza debole è più intuitivo, più adatto a gestire la riflessione
  • Applicazione Astuta di Tecniche di Martingala: la costruzione di martingala nella prova della scomparsa esponenziale è creativa
  • Framework Sistematico: il framework di argomentazione di compattezza esponenziale stabilito ha applicabilità generale

3. Esplicitezza dei Risultati

  • Funzione di velocità completamente esplicita (Teorema 2.6), conveniente per ulteriore analisi e applicazione
  • Comportamento del limite fluido completamente caratterizzato (Tabella 1), copre tutte le regioni di parametri

4. Qualità della Scrittura

  • Struttura Chiara: introduzione fornisce motivazione sufficiente, corpo logico rigoroso, appendici supportano fortemente
  • Uso di Simboli Regolamentato: uso coerente di simboli (come xˉn\bar{x}^n per scala fluida, x~n\tilde{x}^n per scala di deviazione moderata)
  • Rassegna della Letteratura Completa: discussione sufficiente dei collegamenti e delle differenze con lavori correlati

Insufficienze

1. Limitazioni di Praticità

  • Assenza di Verifica Numerica: mancanza di esempi numerici che illustrino il significato pratico dei risultati teorici
  • Discussione di Applicazione Insufficiente: non discute come utilizzare i risultati MDP per l'analisi o la progettazione di sistemi reali
  • Problema di Stima dei Parametri: non affronta come stimare θ\theta, σΘ2\sigma^2_\Theta e altri parametri dai dati

2. Assunzioni Tecniche

  • Condizione della Funzione Generatrice dei Momenti (Assunzione 2.3(ii)) esclude distribuzioni con coda pesante, limitando l'applicabilità
  • Unicità della Scelta di Scala: non discute sufficientemente perché Cin=1Θi/nC^n_i = 1 - \Theta_i/n è la scelta "naturale"

3. Copertura di Portata

  • Assenza di Casi Instabili: μ>0,θ0\mu > 0, \theta \leq 0 e altri casi instabili non sono analizzati
  • Comportamento Transitorio: si concentra principalmente su comportamento asintotico, la qualità dell'approssimazione per nn finito non è discussa

4. Leggibilità

  • Alta Densità Tecnica: per lettori non specialisti, i dettagli delle prove possono essere difficili da seguire
  • Spiegazione Intuitiva Insufficiente: alcuni passaggi tecnici (come la costruzione del sistema ausiliario Υn\Upsilon^n) mancano di spiegazione intuitiva

Impatto

Contributo al Campo

  1. Colmamento di Vuoto Teorico: primo MDP per code dipendenti dallo stato, perfeziona il sistema di teoremi limite in questo campo
  2. Valore Metodologico: le tecniche sviluppate possono essere applicate all'analisi di altri processi casuali riflessi
  3. Ruolo di Ponte: collega la teoria delle code, la teoria della grande deviazione e la teoria dei processi casuali

Valore Pratico

  1. Valutazione del Rischio: l'MDP fornisce stima di probabilità di evento raro più fine rispetto alla grande deviazione
  2. Progettazione del Sistema: la funzione di velocità può essere utilizzata per ottimizzare i parametri del sistema al fine di controllare la probabilità di deviazione
  3. Accelerazione della Simulazione: i risultati MDP possono guidare tecniche di riduzione della varianza come il campionamento per importanza

Riproducibilità

  • Risultati Teorici: enunciati dei teoremi sono chiari, assunzioni esplicite, in linea di principio verificabili
  • Tecniche di Prova: passaggi dettagliati, riferimenti sufficienti, possono servire come esempio per imparare le tecniche MDP
  • Assenza di Codice: non fornisce codice di implementazione numerica (sebbene questo sia accettabile per un lavoro puramente teorico)

Scenari Applicabili

1. Applicazione Diretta

  • Sistemi Medici: pronto soccorso, unità di terapia intensiva e altri sistemi medici sensibili alla congestione
  • Centri di Contatto: scenari in cui il tempo di attesa influenza gli arrivi successivi e la velocità di servizio
  • Sistemi di Computer: pianificazione dei compiti dipendente dal carico e servizio

2. Estensione Teorica

  • Code Multi-Servitore: il metodo può estendersi a sistemi di tipo M/M/cM/M/c dipendenti dallo stato
  • Reti di Code: serve come blocco costitutivo per l'analisi di reti complesse
  • Matematica Finanziaria: strutture ricorsive simili appaiono in alcuni modelli finanziari

3. Prestito Metodologico

  • Processo di Diffusione Riflesso: le tecniche di martingala e l'argomentazione di compattezza esponenziale possono essere utilizzate per altri processi riflessi
  • Sistema Dipendente dallo Stato: il metodo del limite del sistema ausiliario ha applicabilità generale
  • Teoria della Deviazione Moderata: fornisce un esempio per l'analisi MDP di altri sistemi casuali

Raccomandazioni per Ricerca Successiva

  1. Studio Numerico: verificare le previsioni teoriche attraverso simulazione, esplorare la velocità di convergenza per nn finito
  2. Analisi di Sensibilità ai Parametri: studiare sistematicamente la dipendenza della funzione di velocità da θ\theta, σX2\sigma^2_X e altri parametri
  3. Estensione a Coda Pesante: ricercare se l'Assunzione 2.3(ii) può essere rilassata
  4. Caso di Studio Applicativo: selezionare uno scenario di applicazione concreta (come dati di pronto soccorso) per verificare il modello
  5. Metodo di Calcolo: sviluppare algoritmi numerici per il calcolo della funzione di velocità e del percorso ottimale

Riferimenti Bibliografici (Selezionati)

  1. Whitt, W. (1990). Queues with service times and interarrival times depending linearly and randomly upon waiting times. Queueing Systems, 6:335-351.
    • Lavoro classico esteso da questo articolo
  2. Boxma, O., Mandjes, M., and Reed, J. (2016). On a class of reflected AR(1) processes. Journal of Applied Probability, 53(3):818-832.
    • Risultato FCLT per dipendenza dallo stato deterministica
  3. Dupuis, P. and Johnson, D. (2015). Moderate Deviations for Recursive Stochastic Algorithms. Stochastic Systems, 5(1):87-119.
    • Metodo MDP correlato (approccio di convergenza debole)
  4. Puhalskii, A. A. (1999). Moderate deviations for queues in critical loading. Queueing Systems, 31(3):359-392.
    • Lavoro classico MDP per coda GI/GI/1
  5. 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.
    • LDP per sistema ricorsivo con coda pesante

Valutazione Complessiva: Questo è un articolo di matematica teorica di alta qualità che stabilisce rigorosamente il principio di deviazione moderata su traiettoria campionaria per sistemi di code dipendenti dal tempo di attesa, colmando un importante vuoto teorico in questo campo. Il metodo è innovativo, i risultati sono espliciti e le prove sono complete. Le principali insufficienze risiedono nella mancanza di verifica numerica e discussione di applicazione, nonché nelle limitazioni di alcune assunzioni tecniche. Ha importante valore di riferimento per i ricercatori in teoria delle code, teoria della grande deviazione e teoria dei processi casuali, e fornisce anche strumenti teorici per l'analisi del rischio di sistemi reali. Si raccomanda che i lavori futuri integrino la ricerca numerica ed esplorino applicazioni pratiche.