[論文レビュー] A non-increasing Lindley-type equation

Maria Vlasiou|TU/e Research Portal (Eindhoven University of Technology)|Apr 22, 2014

Numerical methods for differential equations被引用数 3

ひとこと要約

この論文は、交互にサービスを行うシステムにおけるサーバー待機時間をモデル化する非増加型の Lindley 方程式 $ W_{n+1} = \max\{0, B_{n+1} - A_n - W_n\} $ を研究する。準備時間分布の有理ラプラス変換を含む十分条件のもとで、$ W $ の定常分布の閉形式表現を導出する。一般化された Wiener-Hopf 技法と収縮写像を用いて正確な解法が可能となり、安定性の解析も可能となる。

ABSTRACT

In this paper we study the Lindley-type equation $W=\max\{0, B - A - W\}$. Its main characteristic is that it is a non-increasing monotone function in its main argument $W$. Our main goal is to derive a closed-form expression of the steady-state distribution of $W$. In general this is not possible, so we shall state a sufficient condition that allows us to do so. We also examine stability issues, derive the tail behaviour of $W$, and briefly discuss how one can iteratively solve this equation by using a contraction mapping.

研究の動機と目的

• 非増加型 Lindley 方程式における待機時間 $ W $ の定常分布の閉形式表現を導出すること。
• 準備時間分布のラプラス変換が有理関数である場合に、その閉形式解が可能となる条件を確立すること。
• 交互にサービスを行うシステムにおける待機時間分布の安定性および尾部挙動を分析すること。
• 非単調性を示す再帰関係に対しても、収縮写像を用いた反復的解法の構成的手法を提供すること。
• キャロセルおよび二重キャロセル系に関する既存の結果を、指数分布サービス時間を有する一般の準備時間分布へと拡張すること。

提案手法

• システムを Lindley 型再帰 $ W_{n+1} = \max\{0, B_{n+1} - A_n - W_n\} $ としてモデル化し、$ B_n $ と $ A_n $ を非負の i.i.d. ランダム変数とする。
• 一般化された Wiener-Hopf 技法を適用し、核 $ \overline{F}_B(x+y) $ をラプラス変換が有理関数である成分に分解することで方程式を解く。
• 準備時間分布のラプラス変換が有理関数であるという仮定を用いて、未知定数および初期確率 $ \pi_0 $ の有限線形方程式系を導出する。
• 収縮写像を用いて、定常分布の存在および一意性を証明し、任意の初期状態からの収束を示す。
• 得られた方程式系の係数および $ \pi_0 $ を解くことにより、定常分布 $ F_W(x) $ を明示的に導出する。
• 具体的な例として、有理ラプラス変換を持つ非位相型分布 $ F_B $ を用い、$ F_W(x) $ に対する明示的閉形式表現が得られることを検証する。

実験結果

リサーチクエスチョン

• RQ1 非増加型 Lindley 方程式 $ W_{n+1} = \max\{0, B_{n+1} - A_n - W_n\} $ の定常分布に対して、どのような条件下で閉形式表現が導出可能か？
• RQ2 再帰関係の非単調性が、定常分布の存在性・一意性および収束性にどのように影響するか？
• RQ3 このシステムにおける待機時間分布の尾部挙動は何か？ また、古典的 Lindley 方程式と比較してどう異なるか？
• RQ4 準備時間のラプラス変換が有理関数である場合、一般化された Wiener-Hopf 法をこの非増加型再帰に適応可能か？
• RQ5 位相型分布でない場合でも、計算的に実行可能な方法で解を明示的に構成可能か？

主な発見

• 準備時間分布のラプラス変換が有理関数である場合、定常待機時間分布 $ F_W(x) $ の閉形式表現が得られ、位相型仮定なしに正確な解法が可能となる。
• 核 $ \overline{F}_B(x+y) $ の分解から得られる $ \pi_0 $ および係数 $ c_1, \dots, c_5 $ を含む線形方程式系を用いて解が構成される。
• 特定の例として、$ \overline{F}_B(x+y) $ に指数関数および三角関数項を含む場合、定常分布は明示的に $ F_W(x) = 1 - \frac{2\mu e^{-x}}{D(\mu)} \left( 5(720+744\mu+347\mu^2) + \dots \right) $ と与えられる。ここで $ D(\mu) = 10800+27000\mu+22353\mu^2+7940\mu^3 $ である。
• 非増加性を示す再帰に対しても、収縮写像により任意の初期状態から一意の平衡分布へ収束することが保証される。
• $ W $ の尾部挙動が解析され、無限大への指数的項および減衰する積分が一貫した方法で取り扱われる。
• 数値的に実装可能であり、近似を用いずに正確な結果が得られる。非位相型分布に対しても、$ \mu = 2 $ の例でその有効性が示された。