Mathematical discipline
待ち行列理論(数学的 確率論の 一分野) において 、 ジャクソン ネットワーク ( ジャクソンニアン ネットワークとも呼ばれる [1] )は、ネットワークが 積形式解を持つため 均衡分布の 計算が特に簡単な待ち行列ネットワークの一種である。これは 待ち行列ネットワーク 理論における最初の重要な発展であり 、この定理のアイデアを一般化して適用し、他のネットワークで同様の積形式解を探す研究は、 インターネットの開発に使用されたアイデアを含め、多くの研究の対象となってきた [2] 。 [3] このネットワークは、 James R. Jacksonによって初めて特定され [4] [5] 、彼の論文は Management Science 誌の 「Ten Most Influential Titles of Management Sciences First Fifty Years」 に再掲載された。 [6]
ジャクソンはバーク とライヒ の研究に触発された が [7] 、 ジャン・ウォーランド は「積形式の結果は…ジャクソン自身が基礎論文で信じていたほど出力定理の直接的な結果ではない」と指摘している [8] 。
以前の積形式の解法は、RRPジャクソンによって タンデムキュー (各顧客が順番に各キューを訪問しなければならない有限のキューの連鎖)と循環ネットワーク(各顧客が順番に各キューを訪問しなければならないキューのループ)に対して発見されました。 [9]
ジャクソンネットワークは複数のノードで構成され、各ノードはキューを表します。キューのサービスレートはノード依存(ノードごとにサービスレートが異なる)と状態依存(キューの長さに応じてサービスレートが変化する)の両方になります。ジョブは固定のルーティングマトリクスに従ってノード間を移動します。各ノードのすべてのジョブは単一の「クラス」に属し、同じサービス時間配分と同じルーティングメカニズムに従います。したがって、ジョブの処理に優先順位はなく、各ノードのすべてのジョブは 先着 順で処理されます。
有限個のジョブが閉じたネットワークを巡回するジャクソンネットワークも、ゴードン・ニューウェル定理 によって記述される積形式解を持つ 。 [10]
ジャクソンネットワークの必要条件 m個 の相互接続されたキューのネットワークは、 次の条件を満たす場合、 ジャクソンネットワーク [11] または ジャクソンネットワーク [12] と呼ばれます。
ネットワークが開いている場合、ノード iへの外部からの到着は ポアソン過程 を形成し 、 すべてのサービス時間は指数分布しており、すべてのキューでのサービス規律は 先着順 です。 キューi でサービスを完了した顧客は、 確率 で 新しいキュー j に移動する、または確率 でシステムを離れる。 オープンネットワークの場合、この確率はキューのいくつかのサブセットではゼロではない。 P i j {\displaystyle P_{ij}} 1 − ∑ j = 1 m P i j {\displaystyle 1-\sum _{j=1}^{m}P_{ij}} すべてのキューの使用率は 1 未満 です 。
定理 m個の M/M/1キュー のオープンジャクソンネットワークにおいて、各キューの 利用率 が1未満の 場合、均衡状態 確率分布が 存在し、状態は 個々のキューの均衡分布の積で与えられる。 ρ i {\displaystyle \rho _{i}} ( k 1 , k 2 , … , k m ) {\displaystyle \scriptstyle {(k_{1},k_{2},\ldots ,k_{m})}}
π ( k 1 , k 2 , … , k m ) = ∏ i = 1 m π i ( k i ) = ∏ i = 1 m [ ρ i k i ( 1 − ρ i ) ] . {\displaystyle \pi (k_{1},k_{2},\ldots ,k_{m})=\prod _{i=1}^{m}\pi _{i}(k_{i})=\prod _{i=1}^{m}[\rho _{i}^{k_{i}}(1-\rho _{i})].} この結果は、ステーションに c i サーバーがあり 、使用率要件が である M/M/c モデル ステーション にも当てはまります 。 π ( k 1 , k 2 , … , k m ) = ∏ i = 1 m π i ( k i ) {\displaystyle \pi (k_{1},k_{2},\ldots ,k_{m})=\prod _{i=1}^{m}\pi _{i}(k_{i})} i th {\displaystyle i^{\text{th}}} ρ i < c i {\displaystyle \rho _{i}<c_{i}}
意味 オープンネットワークでは、ジョブは外部から ポアソン過程 に従って到着し、その速度は です 。各到着ジョブは独立して、 確率 および で ノード jにルーティングされます。ノード i でのサービスが完了すると 、ジョブは 確率 で別のノード j に移動するか、確率 でネットワークから離脱します 。 α > 0 {\displaystyle \alpha >0} p 0 j ≥ 0 {\displaystyle p_{0j}\geq 0} ∑ j = 1 J p 0 j = 1 {\displaystyle \sum _{j=1}^{J}p_{0j}=1} p i j {\displaystyle p_{ij}} p i 0 = 1 − ∑ j = 1 J p i j {\displaystyle p_{i0}=1-\sum _{j=1}^{J}p_{ij}}
したがって、 外部到着と内部遷移の両方を含む
ノード i への全体的な到着率は次のようになります。 λ i {\displaystyle \lambda _{i}}
λ i = α p 0 i + ∑ j = 1 J λ j p j i , i = 1 , … , J . ( 1 ) {\displaystyle \lambda _{i}=\alpha p_{0i}+\sum _{j=1}^{J}\lambda _{j}p_{ji},i=1,\ldots ,J.\qquad (1)} (各ノードでの利用率は 1 未満であり、均衡分布、つまり長期平均動作を見ているため、 jから i に移行するジョブのレートは、 jへの 到着 レート の一部によって制限され、上記 ではサービス レートは無視されます 。) μ j {\displaystyle \mu _{j}}
を定義する と、 を解くことができます 。 a = ( α p 0 i ) i = 1 J {\displaystyle a=(\alpha p_{0i})_{i=1}^{J}} λ = ( I − P T ) − 1 a {\displaystyle \lambda =(I-P^{T})^{-1}a}
すべてのジョブもポアソン過程に従って各ノードから出発し、ノード i にジョブが ある場合のノード i のサービス率として定義されます 。 μ i ( x i ) {\displaystyle \mu _{i}(x_{i})} x i {\displaystyle x_{i}}
時刻 tにおけるノード i のジョブ数を 、 とすると 、 の 均衡 分布は 次 のバランス方程式系によって決定されます。 X i ( t ) {\displaystyle X_{i}(t)} X = ( X i ) i = 1 J {\displaystyle \mathbf {X} =(X_{i})_{i=1}^{J}} X {\displaystyle \mathbf {X} } π ( x ) = P ( X = x ) {\displaystyle \pi (\mathbf {x} )=P(\mathbf {X} =\mathbf {x} )}
π ( x ) ∑ i = 1 J [ α p 0 i + μ i ( x i ) ( 1 − p i i ) ] = ∑ i = 1 J [ π ( x − e i ) α p 0 i + π ( x + e i ) μ i ( x i + 1 ) p i 0 ] + ∑ i = 1 J ∑ j ≠ i π ( x + e i − e j ) μ i ( x i + 1 ) p i j . ( 2 ) {\displaystyle {\begin{aligned}&\pi (\mathbf {x} )\sum _{i=1}^{J}[\alpha p_{0i}+\mu _{i}(x_{i})(1-p_{ii})]\\={}&\sum _{i=1}^{J}[\pi (\mathbf {x} -\mathbf {e} _{i})\alpha p_{0i}+\pi (\mathbf {x} +\mathbf {e} _{i})\mu _{i}(x_{i}+1)p_{i0}]+\sum _{i=1}^{J}\sum _{j\neq i}\pi (\mathbf {x} +\mathbf {e} _{i}-\mathbf {e} _{j})\mu _{i}(x_{i}+1)p_{ij}.\qquad (2)\end{aligned}}} ここで は 単位ベクトル を表します 。 e i {\displaystyle \mathbf {e} _{i}} i th {\displaystyle i^{\text{th}}}
定理 それぞれが 確率質量関数 を持つ 独立確率変数のベクトルを仮定する 。 ( Y 1 , … , Y J ) {\displaystyle (Y_{1},\ldots ,Y_{J})} Y i {\displaystyle Y_{i}}
P ( Y i = n ) = p ( Y i = 0 ) ⋅ λ i n M i ( n ) , ( 3 ) {\displaystyle P(Y_{i}=n)=p(Y_{i}=0)\cdot {\frac {\lambda _{i}^{n}}{M_{i}(n)}},\quad (3)} ここで、 ieが適切に定義されている 場合 、オープンジャクソンネットワークの均衡分布は次の積の形を持ちます。 M i ( n ) = ∏ j = 1 n μ i ( j ) {\displaystyle M_{i}(n)=\prod _{j=1}^{n}\mu _{i}(j)} ∑ n = 1 ∞ λ i n M i ( n ) < ∞ {\displaystyle \sum _{n=1}^{\infty }{\frac {\lambda _{i}^{n}}{M_{i}(n)}}<\infty } P ( Y i = 0 ) = ( 1 + ∑ n = 1 ∞ λ i n M i ( n ) ) − 1 {\displaystyle P(Y_{i}=0)=\left(1+\sum _{n=1}^{\infty }{\frac {\lambda _{i}^{n}}{M_{i}(n)}}\right)^{-1}}
π ( x ) = ∏ i = 1 J P ( Y i = x i ) . {\displaystyle \pi (\mathbf {x} )=\prod _{i=1}^{J}P(Y_{i}=x_{i}).} すべての人のために 。⟩ x ∈ Z + J {\displaystyle \mathbf {x} \in {\mathcal {Z}}_{+}^{J}}
この定理は、各ノードの状態依存のサービス率を許容することで、上記の定理を拡張する。これは、 の分布を 独立 変数 のベクトルで関連付けるものである 。 X {\displaystyle \mathbf {X} } Y {\displaystyle \mathbf {Y} }
例 3ノードのオープンジャクソンネットワーク グラフに示す 3 つのノードを持つジャクソン ネットワークがあるとします。係数は次のようになります。
α = 5 , p 01 = p 02 = 0.5 , p 03 = 0 , {\displaystyle \alpha =5,\quad p_{01}=p_{02}=0.5,\quad p_{03}=0,\quad } P = [ 0 0.5 0.5 0 0 0 0 0 0 ] , μ = [ μ 1 ( x 1 ) μ 2 ( x 2 ) μ 3 ( x 3 ) ] = [ 15 12 10 ] for all x i > 0 {\displaystyle P={\begin{bmatrix}0&0.5&0.5\\0&0&0\\0&0&0\end{bmatrix}},\quad \mu ={\begin{bmatrix}\mu _{1}(x_{1})\\\mu _{2}(x_{2})\\\mu _{3}(x_{3})\end{bmatrix}}={\begin{bmatrix}15\\12\\10\end{bmatrix}}{\text{ for all }}x_{i}>0} すると定理により次のように計算できます。
λ = ( I − P T ) − 1 a = [ 1 0 0 − 0.5 1 0 − 0.5 0 1 ] − 1 [ 0.5 × 5 0.5 × 5 0 ] = [ 1 0 0 0.5 1 0 0.5 0 1 ] [ 2.5 2.5 0 ] = [ 2.5 3.75 1.25 ] {\displaystyle \lambda =(I-P^{T})^{-1}a={\begin{bmatrix}1&0&0\\-0.5&1&0\\-0.5&0&1\end{bmatrix}}^{-1}{\begin{bmatrix}0.5\times 5\\0.5\times 5\\0\end{bmatrix}}={\begin{bmatrix}1&0&0\\0.5&1&0\\0.5&0&1\end{bmatrix}}{\begin{bmatrix}2.5\\2.5\\0\end{bmatrix}}={\begin{bmatrix}2.5\\3.75\\1.25\end{bmatrix}}} の定義によれば 、次のようになります。 Y {\displaystyle \mathbf {Y} }
P ( Y 1 = 0 ) = ( ∑ n = 0 ∞ ( 2.5 15 ) n ) − 1 = 5 6 {\displaystyle P(Y_{1}=0)=\left(\sum _{n=0}^{\infty }\left({\frac {2.5}{15}}\right)^{n}\right)^{-1}={\frac {5}{6}}} P ( Y 2 = 0 ) = ( ∑ n = 0 ∞ ( 3.75 12 ) n ) − 1 = 11 16 {\displaystyle P(Y_{2}=0)=\left(\sum _{n=0}^{\infty }\left({\frac {3.75}{12}}\right)^{n}\right)^{-1}={\frac {11}{16}}} P ( Y 3 = 0 ) = ( ∑ n = 0 ∞ ( 1.25 10 ) n ) − 1 = 7 8 {\displaystyle P(Y_{3}=0)=\left(\sum _{n=0}^{\infty }\left({\frac {1.25}{10}}\right)^{n}\right)^{-1}={\frac {7}{8}}} したがって、各ノードに 1 つのジョブがある確率は次のようになります。
π ( 1 , 1 , 1 ) = 5 6 ⋅ 2.5 15 ⋅ 11 16 ⋅ 3.75 12 ⋅ 7 8 ⋅ 1.25 10 ≈ 0.00326 {\displaystyle \pi (1,1,1)={\frac {5}{6}}\cdot {\frac {2.5}{15}}\cdot {\frac {11}{16}}\cdot {\frac {3.75}{12}}\cdot {\frac {7}{8}}\cdot {\frac {1.25}{10}}\approx 0.00326} ここでのサービス率は状態に依存しないため、は 単純に 幾何分布 に従います。 Y i {\displaystyle Y_{i}}
一般化ジャクソンネットワーク 一般 化ジャクソンネットワークは、ポアソン過程である必要のない 更新到着過程 と、独立かつ同一分布に従う非指数分布のサービス時間を許容する。一般に、このネットワークは 積形定常分布 を持たない ため、近似値が求められる。 [13]
ブラウン運動近似 ある穏やかな条件下では、開放型一般化ジャクソンネットワークのキュー長過程 [ 要説明 ] は、と定義される 反射ブラウン運動 で近似できる。 ここで、 は過程のドリフト、 は 共分散行列 、 は反射行列である。これは、均質流体ネットワークを持つ一般化ジャクソンネットワークと反射ブラウン運動の関係から得られる2次の近似である。 Q ( t ) {\displaystyle Q(t)} RBM Q ( 0 ) ( θ , Γ ; R ) . {\displaystyle \operatorname {RBM} _{Q(0)}(\theta ,\Gamma ;R).} θ {\displaystyle \theta } Γ {\displaystyle \Gamma } R {\displaystyle R}
反射ブラウン運動過程のパラメータは次のように指定されます。
θ = α − ( I − P T ) μ {\displaystyle \theta =\alpha -(I-P^{T})\mu } Γ = ( Γ k ℓ ) with Γ k ℓ = ∑ j = 1 J ( λ j ∧ μ j ) [ p j k ( δ k ℓ − p j ℓ ) + c j 2 ( p j k − δ j k ) ( p j ℓ − δ j ℓ ) ] + α k c 0 , k 2 δ k ℓ {\displaystyle \Gamma =(\Gamma _{k\ell }){\text{ with }}\Gamma _{k\ell }=\sum _{j=1}^{J}(\lambda _{j}\wedge \mu _{j})[p_{jk}(\delta _{k\ell }-p_{j\ell })+c_{j}^{2}(p_{jk}-\delta _{jk})(p_{j\ell }-\delta _{j\ell })]+\alpha _{k}c_{0,k}^{2}\delta _{k\ell }} R = I − P T {\displaystyle R=I-P^{T}} ここで、記号は次のように定義されます。
近似式における記号の定義 シンボル 意味 α = ( α j ) j = 1 J {\displaystyle \alpha =(\alpha _{j})_{j=1}^{J}} 各ノードへの到着率を指定する J ベクトル。 μ = ( μ ) j = 1 J {\displaystyle \mu =(\mu )_{j=1}^{J}} 各ノードのサービスレートを指定する J ベクトル。 P {\displaystyle P} ルーティング マトリックス。 λ j {\displaystyle \lambda _{j}} ノードの有効な到着 。 j th {\displaystyle j^{\text{th}}} c j {\displaystyle c_{j}} ノードでのサービス時間の変動 。 j th {\displaystyle j^{\text{th}}} c 0 , j {\displaystyle c_{0,j}} ノードにおける到着間隔の変動 。 j th {\displaystyle j^{\text{th}}} δ i j {\displaystyle \delta _{ij}} ノード間の相関を指定するための係数。
参照
参考文献 ^ Walrand, J. ; Varaiya, P. (1980). 「ジャクソンネットワークにおける滞在時間と追い越し条件」. 応用確率論の進歩 . 12 (4): 1000–1018 . doi :10.2307/1426753. JSTOR 1426753. ^ Kelly, FP (1976年6月). 「キューのネットワーク」. 応用確率論の進歩 . 8 (2): 416– 432. doi :10.2307/1425912. JSTOR 1425912. ^ ジャクソン、ジェームズ・R.(2004年12月)「『ジョブショップ型待ち行列システム』に関するコメント:背景」 マネジメントサイエンス 誌 50 (12): 1796–1802 . doi :10.1287/mnsc.1040.0268. JSTOR 30046150. ^ ジャクソン、ジェームズ・R. (1963年10月). 「ジョブショップ型待ち行列システム」. 経営科学 . 10 (1): 131– 142. doi :10.1287/mnsc.1040.0268. JSTOR 2627213. 1963年1月のバージョンは、http://www.dtic.mil/dtic/tr/fulltext/u2/296776.pdf で入手できます。2018年4月12日に Wayback Machineでアーカイブされました。 ^ Jackson, JR (1957). 「待ち行列のネットワーク」. オペレーションズ・リサーチ . 5 (4): 518– 521. doi :10.1287/opre.5.4.518. JSTOR 167249. ^ ジャクソン、ジェームズ・R. (2004年12月). 「ジョブショップ型待ち行列システム」. マネジメントサイエンス . 50 (12): 1796– 1802. doi :10.1287/mnsc.1040.0268. JSTOR 30046149. ^ ライヒ, エドガー (1957年9月). 「行列がタンデムになっている場合の待ち時間」. Annals of Mathematical Statistics . 28 (3): 768– 773. doi : 10.1214/aoms/1177706889 . JSTOR 2237237. ^ ウォーランド、ジーン(1983年11月)「準可逆キューネットワークの確率的考察」 IEEE Transactions on Information Theory . 29 (6): 825– 831. doi :10.1109/TIT.1983.1056762. ^ Jackson, RRP (1995). 「書評:待ち行列ネットワークと積形式:システムアプローチ」. IMA Journal of Management Mathematics . 6 (4): 382– 384. doi :10.1093/imaman/6.4.382. ^ Gordon, WJ; Newell, GF (1967). 「指数関数型サーバーを備えた閉キューイングシステム」. オペレーションズ・リサーチ . 15 (2): 254. doi :10.1287/opre.15.2.254. JSTOR 168557. ^ Goodman, Jonathan B.; Massey, William A. (1984年12月). 「非エルゴード的ジャクソンネットワーク」. Journal of Applied Probability . 21 (4): 860– 869. doi :10.2307/3213702. JSTOR 3213702. ^ Walrand, J.; Varaiya, P. (1980年12月). 「ジャクソンネットワークにおける滞在時間と追い越し条件」. 応用確率論の進歩 . 12 (4): 1000–1018 . doi :10.2307/1426753. JSTOR 1426753. ^ チェン・ホン、ヤオ・デイビッド・D. (2001). 『待ち行列ネットワークの基礎:性能、漸近論、最適化 』 シュプリンガー. ISBN 0-387-95166-0 。
単一キューイングノード 到着プロセス キューイングネットワーク サービスポリシー 重要な概念 極限定理 拡張機能 情報システム