[論文レビュー] Recursive equations based models of queueing systems

Nikolai Krivulin|arXiv (Cornell University)|Oct 28, 2012

Simulation Techniques and Applications参考文献 19被引用数 13

ひとこと要約

この論文では、到着時刻と退去時刻に基づく再帰的式を提示し、max、min、加法演算のみを用いてキューイングシステムをモデル化する。これにより、効率的なシミュレーション、性能評価、感度推定が可能になる。このアプローチは、単一サーバ、トランジット、クローズドシステムの間でモデル化を統一し、計算上の利点を提供するとともに、分散低減および無限小パラメータ変動解析（IPA）をサポートする。

ABSTRACT

An overview of the recursive equations based models and their applications in simulation based analysis and optimization of queueing systems is given. These models provide a variety of systems with a convenient and unified representation in terms of recursions for arrival and departure times of customers, which involves only the operations of maximum, minimum, and addition.

研究の動機と目的

• 到着時刻と退去時刻に基づく再帰的式を用いて、単一サーバ、トランジット、クローズド、有限バッファといった多様なキューイングシステムの表現を統一すること。
• max、min、加法演算のみを含む再帰的構造を活用して、キューイングシステムの効率的シミュレーションと性能評価を可能にすること。
• モデルの明示的な構造的性質を活用して、分散低減や無限小パラメータ変動解析（IPA）などの高度なシミュレーション技術を支援すること。
• 特に複数の顧客クラスや非同一のキューイング順序を扱えないという、このアプローチの限界を特定すること。

提案手法

• 到着時刻と退去時刻の再帰的式を用いてシステムのダイナミクスを定式化する：$ A_k = A_{k-1} + \alpha_k $, $ D_k = (A_k \vee D_{k-1}) + \tau_k $（$ G/G/1 $ キューに対して）。
• トランジットキューへ再帰的フレームワークを拡張する：$ D_k^n = (D_k^{n-1} \vee D_{k-1}^n) + \tau_k^n $、およびクローズドシステムへの応用では遅延インデックスを用いる（例：$ D_k^1 = (D_{k-K_1}^N \vee D_{k-1}^1) + \tau_k^1 $）。
• 同じ再帰的フレームワーク内に、ブロッキングルール（例：マニュファクチャリングブロッキング、通信ブロッキング）を用いて有限バッファシステムを表現する。
• 再帰的構造を活用して、システムパラメータに関する性能指標の単調性および凸性の性質を導出する。
• 特に並列プロセッサおよびベクトルプロセッサに適した、効率的なシミュレーションアルゴリズムの設計にモデルを応用する。
• サンプル性能指標の明示的形を活用して、分散低減（例：アンチシスティック変数、共通の乱数）およびIPAに基づく勾配推定を実現する。

実験結果

リサーチクエスチョン

• RQ1 到着時刻と退去時刻に基づく再帰的式は、どのように多様なキューイングシステムを統一的かつ効率的に表現できるか？
• RQ2 max、min、加法演算を用いることで、キューイングダイナミクスのモデリングにどのような計算的・アルゴリズム的利点が得られるか？
• RQ3 再帰的構造は、シミュレーションにおける分散低減および感度分析をどのように支援するか？
• RQ4 このアプローチには、複数の顧客クラスや複雑なキューイング順序を扱えないという、どのような本質的限界があるか？
• RQ5 再帰的式は、効率的なシミュレーションアルゴリズムおよびIPAに基づく勾配推定の開発をどのように促進するか？

主な発見

• 到着時刻と退去時刻に基づく再帰的式は、max、min、加法演算のみを用いて、$ G/G/1 $、トランジット、クローズド、有限バッファキューイングシステムを統一的かつ計算的に効率的なフレームワークでモデル化可能である。
• 再帰的構造のおかげで、特に並列およびベクトルプロセッサに適した、極めて効率的なシミュレーションアルゴリズムの設計が可能となり、従来のイベントスケジューリング手法と比較して時間的・メモリ的コストを削減できる。
• サンプル性能指標の明示的形は、アンチシスティック変数や共通の乱数といった分散低減技術を容易にする。これらの技術は単調性の性質を必要とし、再帰的式から解析可能である。
• モデルは、性能指標の勾配推定の不偏性および一貫性を保証する構造的根拠を提供するため、無限小パラメータ変動解析（IPA）アルゴリズムの開発を支援する。
• このアプローチにより、特にmax-plus演算を含むシステムにおいて、性能指標の単調性および凸性の解析的考察が可能になる。
• 限界として、複数の顧客クラスや異なるキューイング順序を扱えないことが挙げられる。これは、モデルが個々の顧客の識別や選択ルールを追跡せず、ノードを施設として扱うためである。