Unsplash前回の記事待ち行列理論を応用して10分カットが儲かる仕組みを解明してみた。収容人数に制限がある待ち行列を扱いました。床屋に入れる人数、つまり系の収容人数が有限の場合です。この場合、床屋を訪れようとする人数に制限はありません。今回は、これが有限の待ち行列を扱います。これを有限ソースの待ち行列と呼びます。床屋の例でいうと、近所に住んでいる男性30人の得意客しか、その床屋を訪れる可能性はないようなケースです。このようにいうと特殊なケースを想定しているように思えますが、そうでもありません。ある事務所に10人しか従業員がいないのであればコピー機に並ぶ人は最大10人ですし、従業員100人の会社の社員食堂であれば、どんなに混雑していても行列の人数は100人が最大です。今回はまず、このような有限ソースの待ち行列を解析する公式を紹介します。5台のフォークリフトがある物流センターにおいて、最低4台のフォークリフトが正常稼働する確率を推定します。解いているのですが、今回はその別解になります。今回の方が応用範囲が広い解法であることが、後ほど明らかになります。まとめ◆仕事や勉強の息抜きに。センターに5台のフォークリフトがあり、それぞれが1年間に平均λ回故障するとします。故障したフォークリフトを修理する修理工場には十分なキャパがあり、常に待たずに修理ができるものとします。修理は1か月にμ台できるとします。この状況は次のように図示できます。ここで、修理工場を系と考えます。すると、修理工場へのフォークリフトの到着率は5λとなるはずです。1年間で平均λ回故障する(物流センターを離れて修理工場に到着する)フォークリフトが、5台物流センターにあるためです。一方で、修理工場ではフォークリフトがありませんので、修理工場からフォークリフトの離脱率は0になります。これを上図に書き入れると次のようになります。次に、フォークリフトが1台故障して、修理工場で修理している場合を考えてみましょう。次のようになるはずです。同様にフォークリフトが2台故障した場合から5台故障した場合を図で表すと、次のようになります。以上を状態遷移図にまとめると、次のようになります。次に系内にj台ある確率を求める式を見てみましょう。修理工場にフォークリフトが0台ある確率Π0、1台ある確率Π1、、、Π5を求める式です。この公式の導出は少し難しいので、こんな式があることを知っておいて、必要な時に使えればOKです。kCjはk個からj個を選ぶ組み合わせの数を表し、kCj=k!/j!(k-j)!で計算できます。二項分布の公式の意味とExcelでグラフを描く2通りの方法を実演またΠjの定義から、であることは明らかです。次に先で求めた確率を用いて、系に滞在するフォークリフト台数の期待値を求めます。期待値の公式から、次のように計算できます。この結果を使うと、系にフォークリフトが滞在する時間Wはリトルの法則からW=L/λで求められます。先ほどの状態遷移図で見たように、到着率λは状態(系に滞在するフォークリフト台数)により異なります。よって、到着率λの期待値を求めておきます。先ほどと同様に期待値の公式を使うと、λの期待値λ*は次のように計算できます。このλの期待値を先ほどのリトルの法則式に代入することによって、系にフォークリフトが滞在する時間の期待値は次のように求められます。W=L/λ*=L/λ(K-L)フォークリフトが正常稼働する台数をシミュレーション物流センターAでは5台のフォークリフトを所有しています。平均車齢が7年と古いため、平均して1年に1度くらい故障します。故障すると修理が終わるまでに平均1か月かかります。修理期間中は少ないフォークリフト台数でしのぐことになります。修理工場を系と考えると、修理工場にフォークリフトが1台ある確率Π1を求める問題に帰着します。修理にかかる平均日数を求めることは、フォークリフトが系(修理工場)に滞在する日数の期待値Wを求めることと同じです。到着率λや離脱率μがフォークリフトの稼働台数で変わるのは、先の状態遷移図で見た通りです。Excelで計算Excelシートで計算すると次のようになります。クリックすると拡大しますΠ1が物流センターで稼働するフォークリフトが4台になる確率で、28%です。Π0が67%ということは、フォークリフトが5台すべて稼働する確率は2/3しかないことになります。一か月のうち10日間くらいは、フォークリフトが欠けた状態になることを覚悟しましょうということです。これは下記の記事で計算した結果と同じです。待ち行列理論を使ってフォークリフトが正常稼働する台数を推定してみるまた修理にかかる平均日数は30日で、これは修理工場で待ち時間がないため当然の結果ですね。修理工場のキャパに制限がある場合それでは修理工場で待ち行列ができる場合、つまり修理工場のキャパに制限がある場合はどうなるでしょうか?この方が現実的な設定ですね。ここでは、修理工場が同時に2台までしか修理できなくて、フォークリフトの故障頻度が1年に平均2回である場合でシミュレーションしてみましょう。各状態での到着率と離脱率がどうなるかを見てみると、次のように離脱率が2μで頭打ちになることがわかります。クリックすると拡大します従って、状態遷移図は次のようになります。このことに注意して先ほどと同じようにExcelで計算すると、次のようになります。このようにフォークリフトが4台稼働する(修理工場にフォークリフトが1台ある)確率は30%になります。フォークリフトが5台すべて稼働する確率は37%しかないことがわかり、これではとてもオペレーションが回せないことがわかります。1台の修理にかかる平均日数は54日もあることもわかります。これは、修理工場で1度に2台しか修理できるキャパがないために、待ち時間ができてしまうためです。まとめ以前の記事待ち行列理論を使ってフォークリフトが正常稼働する台数を推定してみるでは、物流センターを系として待ち行列理論を適用しました。今回は逆に、修理工場を系として待ち行列理論を適用しました。どちらでも解けますが、修理工場のキャパに制限がある場合には、今回のように修理工場を系として有限ソースの待ち行列の公式を適用する必要があります。ロジギークよろしければシェアお願いしますTwitterFacebookLinkedInB!HatenaLINELINE【直列型待ち行列モデル】センター出荷作業のサイクルタイムを求めてみるNext待ち行列理論を応用して10分カットが儲かる仕組みを解明してみた。...◆キーワードで探すイチローの安打数がポアソン分布にならず正規分布になる理由を考察してみた【マルコフ連鎖】標準体重の人が太りすぎになるまで平均何年かかる?待ち行列理論を応用して10分カットが儲かる仕組みを解明してみた。計算グラフモンテカルロ積分クラスタリングカーネルトリックPythonSimPy(離散事象シミュレーション)inventorize(在庫シミュレーション)SARIMAX(時系列解析)investpy(財務データ取得)PyPDF2(PDF操作)仕事や勉強の合間の息抜きに管理人ロジギークロジスティクス業界で6か国に勤務。