[論文レビュー] The Maximum of an Asymmetric Simple Random Walk with Reflection

Steven R. Finch|arXiv (Cornell University)|Aug 6, 2018

Stochastic processes and statistical mechanics参考文献 4被引用数 6

ひとこと要約

本稿では、原点で反射を受ける非対称な単純ランダムウォークの最大値の平均および分散について、漸近的表現を導出している。強反射と弱反射のメカニズムを区別し、期待最大値が時間とともに対数的に増加すること、オイラー定数および周期的フラクチュエーションを含む補正項があることを示している。結果はシミュレーションおよび母関数解析により検証されており、ℓ=1 の場合の交通信号制御によるキューイング問題についても解を導出している。

ABSTRACT

Consider the extreme value of a Bernoulli random walk on the one-dimensional integer lattice, with reflection at 0, over a finite discrete time interval. Only the asymmetric (biased) case is discussed. Asymptotic mean/variance results are given as the time interval length approaches infinity. We similarly solve an elementary traffic light problem from queueing theory.

研究の動機と目的

• 0 における反射を伴う非対称な単純ランダムウォークの極値分布を、有限時間区間で分析すること。
• 時間の範囲 n → ∞ における期待最大値およびその分散の正確な漸近的表現を導出すること。
• 強反射と弱反射のメカニズムを区別し、それらが最大値に与える影響を定量化すること。
• 到着および出発制約を伴う周期的交通信号システムを含むキューイング理論問題に、この枠組みを適用すること。
• 大規模なモンテカルロシミュレーションを用いて理論的予測を検証し、実験的データと解析的予測を比較すること。

提案手法

• 強反射に対しては (k+1)×(k+1) の遷移行列 A_k を用いて定義される再帰関係で反射ランダムウォークをモデル化し、弱反射に対しても同様の行列を用いる。
• 行列幾何的手法を用いて、G_k(z) = ∑ₙ P(M_n ≤ k)zⁿ の母関数を構築し、有理関数 Q_k / R_k の形で表現する。
• 母関数の特異点解析を用いて係数の漸近的挙動を抽出し、特に支配的特異点 z_k に注目する。
• 離散ガブリー分布近似を用いて漸近的平均および分散を導出し、オイラー定数 γ および周期関数 φ(n), ψ(n), ω(n) を含む補正項を含める。
• 対称ケースからの結果を一般化し、ベース r = 2 を r = q/p > 1 に置き換えることで、非対称ウォークに拡張する。Louchard らの先行研究を拡張する。
• n = 10¹⁰ の大規模なシミュレーションを用い、両方のウォークタイプおよび交通信号モデルに対して 40,000 本のサンプルパスを用いて理論的予測を検証する。

実験結果

リサーチクエスチョン

• RQ1 非対称ランダムウォークが原点で反射を受ける場合、期待最大値 E[M_n] は時間 n → ∞ のときどのようにスケーリングするか？
• RQ2 強反射と弱反射のメカニズムにおける期待最大値の違いは何か？
• RQ3 偏りパラメータ p < q に依存する、最大値の漸近的平均および分散はどのように変化するか？
• RQ4 周期的キューイングシステム、例えば交互に緑と赤に切り替わる交通信号のような場合に、同じ解析枠組みを適用できるか？
• RQ5 交通信号モデルにおける ℓ=1 の理論的予測が ℓ>1 に一般化できる範囲はどの程度か？

主な発見

• 強反射の場合、E[M_n] ~ ln(n)/ln(q/p) + (γ + ln((1−2p)²/(2q²)))/ln(q/p) + 1/2 + φ(n) が成り立ち、φ(n) は log_{q/p}(n) を周期とする周期関数（周期 1）。
• 弱反射の場合、E[M_n] ~ ln(n)/ln(q/p) + (γ + ln(p(1−2p)²/q²))/ln(q/p) + 1/2 + ψ(n) が成り立ち、ψ(n) も同様に周期的。
• 両ケースの漸近的分散は V[M_n] ~ π²/(6 ln²(q/p)) + 1/12 + ω(n) であり、ω(n) は無視可能で周期的。
• ℓ=1 の交通信号問題では、E[M_n] ~ ln(n)/(2 ln(q/p)) + (γ + ln(p(1−2p)²/(2q³)))/(2 ln(q/p)) + 1/2 + ψ(n) であり、V[M_n] ~ π²/(24 ln²(q/p)) + 1/12 + ω(n) が成り立つ。
• n = 10¹⁰ および 40,000 パスのシミュレーションにより、両方のウォークタイプおよび ℓ=1 の交通信号モデルにおいて、実験的ヒストограмと理論的予測が良好に一致している。
• ℓ=1 の理論的予測は ℓ>1 に一般化できないことが示され、シミュレーションで顕著なずれが確認され、分布が ℓ に依存しないという仮説は否定された。