[論文レビュー] Asymptotically Optimal Importance Sampling for Jackson Networks with a Tree Topology

Ali Devin Sezer|ArXiv.org|Aug 23, 2007

Probability and Risk Models参考文献 8被引用数 3

ひとこと要約

本稿は、木構造をとる安定なオープン・ジャクソンネットワークにおける希少なバッファオーバーフロー事象の漸近的に最適な重要度サンプリング（IS）アルゴリズムを、ハミルトニアン＝ヤコビ＝ベルマン（HJB）方程式の部分解法に基づき開発する。部分解法を用いたアプローチにより、効果的ノード効用とレートに基づく再帰的アルゴリズムを導入し、共有バッファ構造および個別バッファ構造の両方において漸近的最適性を達成する。

ABSTRACT

Importance sampling (IS) is a variance reduction method for simulating rare events. A recent paper by Dupuis, Wang and Sezer (Ann. App. Probab. 17(4):1306- 1346, 2007) exploits connections between IS and stochastic games and optimal control problems to show how to design and analyze simple and efficient IS algorithms for various overflow events for tandem Jackson networks. The present paper uses the same approach to build asymptotically optimal IS schemes for stable open Jackson networks with a tree topology. Customers arrive at the single root of the tree. The rare overflow event we consider is the following: given that initially the network is empty, the system experiences a buffer overflow before returning to the empty state. Two types of buffer structures are considered: 1) A single system-wide buffer of size $n$ shared by all nodes, 2) each node $i$ has its own buffer of size $β_i n$, $β_i \in (0,1)$.

研究の動機と目的

• 安定なオープン・ジャクソンネットワークに木構造をとる場合の希少バッファオーバーフロー事象に対する漸近的に最適なISスキームの設計。
• 直列構成を超えるネットワークに対し、HJB方程式および境界条件への部分解法アプローチの拡張。
• 1つのシステムワイドバッファ（サイズ $ n $）と、各ノードに個別バッファ（サイズ $ \beta_i n $）を設ける2種類のバッファ構造の取り扱い。
• 効果的ノード効用とレートに基づき、証明可能な最適性を有する滑らかで部分的な解を構築。
• 効用とサービスレートのヒューリスティック解釈を持つ部分解を生成する再帰的アルゴリズムの提供。

提案手法

• 先行研究の最適制御／部分解法フレームワークを用い、IS設計をHJB方程式の滑らかな部分解の構築問題に還元する。
• ノード状態から導かれる有効勾配をパrameterとする、アフィン関数の滑らかな最小値として部分解を計算する再帰的アルゴリズムを採用する。
• ノードバッファ状態（空または非空）に基づき、システムのダイナミクスを反映する有効勾配、効用、レートを定義する。
• 滑らかさとリプシッツ連続性を保証するため、重み付き指数平滑化（ソフトマックスに類似）変換をアフィン関数の組み合わせに適用する。
• 近似精度を制御し、部分解の性質を保証するパrameter化された関数族 $ W^{\epsilon,\delta}(x) $（$ \epsilon, \delta > 0 $）を用いる。
• 構築された部分解が[16]における一般最適性定理のすべての条件を満たすことを証明し、結果として得られるIS推定量の漸近的最適性を保証する。

実験結果

リサーチクエスチョン

• RQ1 木構造ジャクソンネットワークにおけるバッファオーバーフローに対して、部分解法を用いた漸近的最適ISを構築可能か？
• RQ2 バッファ状態がオーバーフロー確率に与える影響を反映するため、効果的ノード効用とレートをどのように定義できるか？
• RQ3 部分解の再帰的構築が、HJB方程式および境界条件を満たす滑らかで有効な部分解を生成するか？
• RQ4 得られたISアルゴリズムは、共有バッファおよび個別バッファ構成の両方において漸近的に最適か？
• RQ5 本手法は、DupuisとWang（2010）が提示したより一般的な枠組みとどのような関係にあるか？

主な発見

• 提案された再帰的アルゴリズムは、HJB方程式およびその境界条件の滑らかな部分解を構築する。これは、ISの漸近的最適性にとって必要かつ十分である。
• 部分解は $ W(x) = -\delta \log \sum_{l=1}^{L} \exp\left( -\frac{1}{\delta}(2\gamma - \alpha_l \epsilon + \langle q_l, x \rangle) \right) $ の形をとり、滑らかさと適切なスケーリングを保証する。
• 任意の $ x \in \mathbb{R}_+^d $ に対して、部分解は $ H_b(DW(x)) \geq -C_1 \exp(-\epsilon/\delta) $ を満たし、$ \epsilon, \delta \to 0 $ のとき近似的最適性が証明される。
• アルゴリズムにより、IS推定量は漸近的効率性を達成し、希少事象確率に対して指数的でない速度で分散が増加する。
• 第5および第6節の数値結果により、提案手法の実用的有効性と低分散性が確認される。
• 本手法は、直列ネットワークに関する先行結果を一般化し、任意の木構造ジャクソンネットワークに適用可能な体系的な枠組みを提供する。