ケリーネットワークの積形式解の存在証明の試み(失敗)

ケリーネットワークへの助走」の続きです。
ケリーネットワークが積形式解を持つことの証明は、「ジャクソンネットワークの積形式解の存在(1)」「ジャクソンネットワークの積形式解の存在(2)」「ジャクソンネットワークの積形式解の存在(3)」で展開した証明をほぼなぞることで証明出来るのではないかと、「眼が覚めて眠れない」でひらめきました。
ところが、これを試してみたのですがうまくいきませんでした。以下は、その「失敗した」試みです。


待ち行列ネットワークのステーションの数をNとします。ステーションを区別するために1から順番に番号をつけていきます。i番目のステーションm_i台の装置から成り、装置の処理時間は指数分布であるとします。ジョブのクラスの数をQとします。クラスに1から順番に番号をつけていきます。外から各ステーションへの各クラスのジョブの到着の時刻の間隔も指数分布であるとします。
i番目のステーションを通るクラスkジョブスループット\theta_{ik}で表します。また、i番目のステーションに外から入ってくるクラスkジョブスループット\lambda_{ik}で表します。ステーションiを終えたクラスkジョブがクラスsになってステーションjに進む確率をr_{ikjs}で表します。そうすると、ステーションjにクラスsジョブの入ってくる量\theta_{js}

  • \theta_{js}=\lambda_{js}+\Bigsum_{i=1}^N\Bigsum_{k=1}^Q\theta_{ik}r_{ikjs}・・・・・・(1)

で表すことが出来ます。式(1)は\thetaについての連立一次方程式になっています。
これが\thetaについて解けると仮定します。これを、

  • \theta_{js}=f_{js}(\Lambda,R)・・・・・・(2)

と表わすことにします。ただし\Lambda\lambda_{ik}を簡略的に表わしたものであり、R(r_{ikjs})を簡略的に表わしたものとします。さてaを任意の定数として、式(1)で

  • \theta_{ik}a\theta_{ik}
  • \lambda_{ik}a\lambda_{ik}

で置き換えると、やはり式(1)が成り立つので式(2)について

  • f_{js}(a\Lambda,R)=af_{jm}(\Lambda,R)・・・・・・(3)

が成り立つことが分かります。一方、定義から

  • \Bigsum_{s=1}^Q\theta_{js}=\frac{m_ju_j}{t_{ej}}・・・・・・(4)

です。式(4)を変形して

  • \frac{1}{t_{ej}}=\frac{\Bigsum_{s=1}^Q\theta_{js}}{m_ju_j}・・・・・・(5)

となります。さらに

  • \Bigsum_{s=1}^Q\theta_{js}=\theta_j・・・・・・(6)

とおけば

  • \frac{1}{t_{ej}}=\frac{\theta_j}{m_ju_j}・・・・・・(7)

ネットワークの「状態」を各ステーションで待っている、あるいは処理中である各々のクラスのジョブの数の組[tex:*1]で表すことにします。そして、状態K(+(j,s))について以下のような局所平衡方程式が成立すると仮定します。

  • p(K(+(j,s)))\frac{min(k_j+1,m_j)}{t_{ej}}
    • =p(K)\lambda_{js}+\Bigsum_{i=1}^N\Bigsum_{k=1}^Qp(K(+(i,k)))\frac{min(k_i+1,m_i)}{t_{ei}}r_{ikjs}・・・・・・(8)


さあ、ここで詰まってしまいました。
状態K(+j,s)ではステーションjのクラスsジョブはいったいいくつ処理中でしょう?
ステーションjの装置台数はm_jです。処理中であるクラスsジョブ数は、可能性としては0からm_jまであります。0というのはつまり、m_j台の装置が全て他のクラスのジョブを処理している場合のことを指しています。
ですから式(8)の左辺、これは、ステーションjからクラスsジョブdtの間に終了する確率を表したつもりだったのですが、かならずしも

  • p(K(+(j,s)))\frac{min(k_j+1,m_j)}{t_{ej}}

に等しくないことが分かります。よって式(8)は正しくありません。式(8)が成り立たないならば、ここから先へ話が進みません。ここで立ち往生です。


複数クラスを持つM/M/1待ち行列(1)」に続きます。

*1:i,k),(j,s),.....)]で定義します。さらに状態((i,k),(j,s),.....)Kと略記することにします。また、ステーションjのクラスsジョブ数をk_{js}で表します。また、状態Kから、k_{js}だけを+1した状態を[tex:K(+(j,s