2次のアーラン分布の系列の重ね合せ系列の2乗変動係数(2)

2次のアーラン分布の系列の重ね合せ系列の2乗変動係数(1)」の続きです。


今度は2次のアーラン分布を3つ重ねた時の2乗変動係数 c^2を求めます。
3つの系列の名前を系列1、2、3とします。系列1でイベントが発生した時刻のうち一つを選んで、その時刻をt=0と置きます。この次に系列1、2あるいは3のいずれかにおいて発生するイベントまでの間隔を考えます。t=0の時の系列2の状態は確率1/2で図1の状態1、確率1/2で状態2です。また、系列3の状態も確率1/2で図1の状態1、確率1/2で状態2です。


今、t=0で系列2と3の状態が状態1であると仮定します。そうすると、それぞれの系列において次のイベントが発生するまでの間隔は、3つとも同じ2次のアーラン分布になります。よって、系列1でも2でも3でも、次のイベントまでの間隔の確率密度関数

  • f(t;2,\lambda)=\lambda^2t\exp(-\lambda{t})・・・・・(2)

となります。分布関数F(t;2,\lambda)は、

  • F(t;2,\lambda)=-\lambda{t}\exp(-\lambda{t})-\exp(-\lambda{t})+1・・・・・(3)

になります。「min(X,Y)の確率密度関数」の式(1)を応用して考えると、系列1,2,3を区別せずにt=0から次のイベントまでの間隔の分布の分布関数G_{11}(t)]は

  • G_{11}=1-[1-F(t;2\lambda)]^3

となります。これをt微分して

  • g_{11}(t)=3f(t;2,\lambda)[1-F(t;2,\lambda)]^2=3f(t;2,\lambda)[\lambda{t}\exp(-\lambda(-\lambda{t})+\exp(-\lambda{t})\)^2
    • =3f(t;2,\lambda)\exp(-2\lambda{t})(\lambda^2t^2+2\lambda{t}+1)=3\lambda^2t(\lambda^2t^2+2\lambda{t}+1)\exp(-3\lambda{t})
    • =3(\lambda^4t^3+2\lambda^3t^2+\lambda^2t)\exp(-3\lambda{t})=3\left(\frac{6}{3^4}\frac{(3\lambda)^4t^3}{3!}+\frac{4}{3^3}\frac{(3\lambda)^3t^2}{2!}+\frac{1}{3^2}(3\lambda)^2t\right)\exp(-3\lambda{t})
    • =\left(\frac{2}{3^2}\frac{(3\lambda)^4t^3}{3!}+\frac{4}{3^2}\frac{(3\lambda)^3t^2}{2!}+\frac{3}{3^2}(3\lambda)^2t\right)\exp(-3\lambda{t})
    • =\frac{1}{9}[2f(t;4.3\lambda)+4f(t;3,3\lambda)+3f(t;2,3\lambda)]

よって

  • g_{11}(t)=\frac{1}{9}[2f(t;4,3\lambda)+4f(t;3,3\lambda)+3f(t;2,3\lambda)]・・・・・(13)


今度は、t=0で系列2が状態1、系列3が状態2であると仮定します。そうすると、系列1、2での次のイベントまでの間隔の確率密度関数f(t;2,\lambda)のままですが、系列3での次のイベントまでの間隔の確率密度関数

  • f(t;1,\lambda)=\lambda\exp(-\lambda{t})・・・・・(5)

になります。F(t;1,\lambda)は、

  • F(t;1,\lambda)=1-\exp(-\lambda{t})・・・・・(6)

よって、系列1,2,3を区別せずにt=0から次のイベントまでの間隔の分布の分布関数G_{12}(t)は「min(X,Y)の確率密度関数」の式(1)を応用して考えると、

  • G_{12}(t)=1-[1-F(t;2,\lambda)]^2[1-F(t;1,\lambda)]

となります。これをt微分して

  • g_{12}(t)=2[1-F(t;2,\lambda)]f(t;2,\lambda)[1-F(t;1,\lambda)]+[1-F(t;2,\lambda)]^2f(t;1,\lambda)
    • =2[\lambda{t}\exp(-\lambda{t})+\exp(-\lambda{t})]\lambda^2t\exp(-\lambda{t})\exp(-\lambda{t})+[\lambda{t}\exp(-\lambda{t})+\exp(-\lambda{t})]^2\lambda\exp(-\lambda{t})
    • =2(\lambda{t}+1)\lambda^2t\exp(-3\lambda{t})+(\lambda{t}+1)^2\lambda\exp(-3\lambda{t})
    • =[2(\lambda{t}+1)\lambda^2t+(\lambda{t}+1)^2\lambda]\exp(-3\lambda{t})
    • =[2\lambda^3t^2+2\lambda^2t+\lambda^3t^2+2\lambda^2t+\lambda]\exp(-3\lambda{t})=(3\lambda^3t^2+4\lambda^2t+\lambda)\exp(-3\lambda{t})
    • =\left[\frac{3\cdot{2}}{3^3}\frac{(3\lambda)^3t^2}{2}+\frac{4}{3^2}(3\lambda)^2t+\frac{1}{3}(3\lambda)\right]\exp(-3\lambda{t})=\frac{1}{9}\left[2\frac{(3\lambda)^3t^2}{2}+4(3\lambda)^2t+3(3\lambda)\right]\exp(-3\lambda{t})
    • =\frac{1}{9}[2f(t;3,3\lambda)+4f(t;2,3\lambda)+3f(t;1,3\lambda)]

よって

  • g_{12}(t)=\frac{1}{9}[2f(t;3,3\lambda)+4f(t;2,3\lambda)+3f(t;1,3\lambda)]・・・・・(14)

となります。
t=0で系列2が状態2、系列3が状態1である場合の、t=0から次のイベントまでの間隔の確率密度関数g_{21}(t)は、系列2と3が同等であることを考えれば

  • g_{21}(t)=g_{12}(t)・・・・・(15)

であることが分かります。

2次のアーラン分布の系列の重ね合せ系列の2乗変動係数(3)」に続きます。