Cvičení 7 Systémy hromadné obsluhy s omezenou kapacitou 1. Systém M/M/n/m/FIFO Vstupní proud zákazníků je Poissonův proces s parametrem λ, doba obsluhy se řídí rozložením ( )µEx , v systému je n linek obsluhy, kapacita systému je omezená (je rovna m) a frontový režim je „první vstupuje, první je obsloužen“. Označme µ λ =β , n β =ρ . Systém se může stabilizovat vždy. Stacionární rozložení:        +=ρ = β = m,1,njproa n! n n,1,2,jproa !j a 0 j n 0 j j K K , kde ∑ ∑ − = = ρ+ β = 1n 0j m nj j nj0 !n n !j 1 a . Charakteristiky stabilizovaného systému: Pravděpodobnost, že přicházející zákazník bude odmítnut: 0 m n Z a !n n P ρ= . Pravděpodobnost, že přicházející zákazník bude čekat ve frontě: ( )     =ρ− ≠ρ ρ− ρ− = − 1pronma 1pro 1 1 a P n nm n Q . Střední hodnota počtu přijatých zákazníků za jednotku času: ( )ZP P1−λ=λ . Střední hodnota počtu odmítnutých zákazníků za jednotku času: ZZ Pλ=λ . Střední hodnota počtu zákazníků ve frontě: ( ) ( )∑+= −= m 1nj jQ anjNE . Střední hodnota počtu obsluhovaných zákazníků: ( ) ( )ZS P1NE −β= . Využití systému: ( )ZP1−ρ=κ . Ostatní charakteristiky dostaneme pomocí Littleova vzorce. Charakteristiky tohoto systému počítá funkce odmitani.m Příklad 1.: V autoservisu jsou 3 mycí rampy a jeden pracovník, jemuž mytí auta trvá v průměru 12 min. Za 1 h přijedou průměrně 3 auta. Jsou-li však v okamžiku příjezdu auta všechny rampy obsazeny, auto nečeká a vrací se později. a) Jaká je pravděpodobnost, že v autoservisu budou 0, 1, 2, 3 auta? b) Vypočtěte střední hodnotu počtu zákazníků v autoservisu a ve frontě. c) Vypočtěte střední hodnotu doby čekání ve frontě. d) Jaká je pravděpodobnost, že bude volná aspoň jedna rampa? e) Vypočtěte využití systému. Výsledky: Jde o systém M/M/n/m/FIFO, kde n = 1, m = 3, 5 3 , 5 3 ,5,3 =ρ=β=µ=λ ad a) 272 27 a, 272 45 a, 272 75 a, 272 125 a 3210 ==== ad b) ( ) ( ) 272 99 NE, 272 246 NE Q == ad c) ( ) s18min7 272 33 WE Q == ad d) 272 245 a1 3 =− ad e) 54,0 272 147 ==κ Návod na řešení pomocí MATLABu: lambda=3;mi=5;n=1;m=3; [a,PZ,PQ,lambdaP,lambdaZ,kappa,ENS,ENQ,EN,EWS,EWQ,EW]=odmitani(lambda,mi,n,m) 2. Uzavřený systém M/M/n/m/FIFO V systému je m zákazníků, přičemž mohou čekat v omezené frontě délky m – n ≥ 0. Zákazníci po ukončení obsluhy opouštějí systém, ale později se do něj vracejí s novým požadavkem. Doba pobytu každého zákazníka mimo systém se řídí rozložením ( )λEx , doba obsluhy každé linky se řídí rozložením ( )µEx . Označme µ λ =β , n β =ρ . Stacionární rozložení: ( )       ++=ρ − =β      = m,,2n,1njproa !jm!n !mn n,,2,1jproa j m a 0 j n 0 j j K K , kde ∑= −= m 1j j0 a1a Pravděpodobnost, že přicházející zákazník bude čekat ve frontě: ( )ρ− β = 1!n aP n 0Q Charakteristiky stabilizovaného systému: Střední hodnota počtu zákazníků v systému: ( ) ∑= = m 0j jjaNE . Střední hodnota počtu obsluhovaných zákazníků: ( ) ∑∑ = − = += m nj j 1n 0j jS anjaNE . Střední hodnota počtu zákazníků mimo systém: ( ) ( )NEmNE R −= . Střední hodnota počtu zákazníků přicházejících za jednotku času (intenzita vstupního proudu): ( )RR NEλ=λ . Využití systému: ( )RNEρ=κ . Ostatní charakteristiky dostaneme pomocí Littleova vzorce. Charakteristiky tohoto systému počítá funkce uzavreny.m Příklad 2.: V uzavřeném systému M/M/n/m/FIFO, kde n = 1, m = 2, 2,3 =µ=λ vypočtěte stacionární rozložení a základní charakteristiky systému. Úlohy řešte pomocí funkce uzavreny.m. Výsledek: Stacionární rozložení: ( )5294,0;3529,0;1176,0 17 9 , 17 6 , 17 2 =      =a Střední hodnota počtu zákazníků v systému je 1,41. Střední hodnota počtu obsluhovaných zákazníků je 0,88. Střední hodnota počtu zákazníků ve frontě je 0,53. Střední hodnota počtu zákazníků mimo systém je 0,59. Intenzita vstupního proudu zákazníků je 1,76. Systém je využit na 88 %. Návod na řešení pomocí MATLABu: lambda=3;mi=2;n=1;m=2; function[a,ENS,ENR,EN,lambdaR,kappa]=uzavreny(lambda,mi,n,m) Příklad 3.: Skupinu pěti stejných strojů má na starosti jeden údržbář. Doba bezporuchového provozu stroje má exponenciální rozložení se střední hodnotou 1/2 směny a doba opravy má rovněž exponenciální rozložení se střední hodnotou 1/20 směny. a) Jaká je pravděpodobnost, že všechny stroje pracují? b) Jaká je pravděpodobnost, že budou současně vyřazeny aspoň dva stroje? Úlohy řešte pomocí funkce uzavreny.m. Výsledek: Jde o uzavřený systém M/M/n/m/FIFO, kde n = 1, m = 5, 10 1 , 10 1 ,20,2 =ρ=β=µ=λ Ad a) 0,564 Ad b) 0,154 Návod na řešení pomocí MATLABu: lambda=2;mi=20;n=1;m=5; function[a,ENS,ENR,EN,lambdaR,kappa]=uzavreny(lambda,mi,n,m) Příklad k samostatnému řešení: Do univerzitního bufetu s kapacitou osm stolků po čtyřech místech přichází v průměru 25 studentů za hodinu. V průměru se student zdrží 30 minut. Můžeme předpokládat, že vstupní proud studentů je Poissonův proces a doba pobytu má exponenciální rozložení. Jaké jsou pravděpodobnosti, že bufet bude prázdný neb o naopak plně obsazen? Na kolik procent je bufet vyžíván? Jaký je průměrný počet volných míst k sezení? Jaký je průměrný počet studentů, kteří si nemají kam sednout za hodinu provozu bufetu? Výsledky: Pravděpodobnost, že bufet bude prázdný, je 3,2767*10-6. Pravděpodobnost, že bufet bude plně obsazen, je 1,787*10-6 . Bufet je využíván na 39,1 %. Počet volných míst k sezení je asi 19. Průměrný počet studentů, kteří si nemají kam sednout za hodinu provozu bufetu, je asi 44,675*10-6 .