[
Імітаційне моделювання на основі попередньо встановлених залежностей
]
<div style="font-size:16px;text-align:justify;"> <div> КАФЕДРА <br /> Економічної кібернетики <br /> <a href="Імітаційне_моделювання" title="Імітаційне моделювання">імітаційне моделювання</a> на основі попередньо встановлених залежностей. <br /> м. Жовті Води 2004 <br /><br clear=all style=page-break-before:always> <b style=mso-bidi-font-weight:normal>Зміст</b> <br /> 1. Постановка <a href="Завдання_та_функції_діяльності_прокуратури" title="Завдання та функції діяльності прокуратури">завдання _________________________________</a> 3 <br /> 2. Метод рішення задачі <a href="Гарантії_місцевого_самоврядування_5" title="Гарантії місцевого самоврядування 5">_________________________________ 5</a> <br /> 3. Програмне <a href="Рішення_довільніх_систем_лінійніх_рівнянь" title="Рішення довільніх систем лінійніх рівнянь">рішення _________________________________</a> 8 <br /> 4. Керівництво пользователя_______________________________ 11 <br /> 5. Лістинг програми <a href="_Входные_барьеры_на_рынок_отрасли" title=" Входные барьеры на рынок отрасли">_________________________________</a> 12 <br /><br clear=all style=page-break-before:always> <b style=mso-bidi-font-weight:normal>Постановка завдання</b> <b style=mso-bidi-font-weight:normal> </b><br /> У сучасному світі гарантією ефективної <a href="Роботи" title="Роботи">роботи</a> будь-якого підприємства є раціональне використання грошових коштів і трудового фактора. Так для розрахунку економічного ефекту роботи касового залу необхідно провести імітаційне <a href="Моделювання" title="Моделювання">моделювання</a> на основі попередньо встановлених залежностей. <br /> Припустимо, що клієнти в зал прибувають з інтервалом, що обчислюються у хвилинах (див. рис. 1). <br /><o:lock v:ext=edit aspectratio=t /><o:lock v:ext=edit text=t /><img width=396 height=12 src=dopb174739.zip v:shapes=_x0000_s1026><o:lock v:ext=edit shapetype=t /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> t <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> <sup><u>Серпень</u></sup> 2000 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1026 DrawAspect=Content ObjectID=_1336184532></o:OLEObject><br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1028 DrawAspect=Content ObjectID=_1336184533></o:OLEObject><br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1030 DrawAspect=Content ObjectID=_1336184534></o:OLEObject><br /></div></td></tr></table><img width=300 height=167 src=dopb174740.zip v:shapes="_x0000_s1027 _x0000_s1028 _x0000_s1029 _x0000_s1030 _x0000_s1031 _x0000_s1032 _x0000_s1033 _x0000_s1034 _x0000_s1035 _x0000_s1036 _x0000_s1037 _x0000_s1038 _x0000_s1039"><table cellpadding=0 cellspacing=0><tr><td width=42 height=36 style=vertical-align:top><table cellpadding=0 cellspacing=0 width=100%><tr><td><div v:shape=_x0000_s1040 style="padding:3.6pt 7.2pt 3.6pt 7.2pt"><img width=19 height=23 src=dopb174741.zip v:shapes=_x0000_i1032><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1032 DrawAspect=Content ObjectID=_1336184538></o:OLEObject></xml><br /></div></td></tr></table></td></tr></table> SHAPE \ * MERGEFORMAT <img width=394 height=222 src=dopb174742.zip v:shapes=_x0000_i1033><br /> Рис. 1 - «Прихід клієнтів в зал» <br /> Прихід клієнтів в зал описується пуассоновским потоком з інтенсивністю r, який визначається наступним чином: <br /><img width=12 height=23 src=dopb174743.zip v:shapes=_x0000_i1034><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1034 DrawAspect=Content ObjectID=_1336184539></o:OLEObject></xml><img width=124 height=44 src=dopb174744.zip v:shapes=_x0000_i1035><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1035 DrawAspect=Content ObjectID=_1336184540></o:OLEObject></xml> (1.1) <br /> де: r - інтенсивність потоку; <br /> k - час між парафіями клієнтів. <br /> Параметр k може приймати дискретні значення від нуля до нескінченності. Причому k = 0 означає прихід відразу двох клієнтів. <br /> Припустимо, в залі є N кас. <a href="Математика" title="Математика"> Математичне</a> сподівання обслуговування клієнтів у банку позначимо <img width=16 height=17 src=dopb174745.zip v:shapes=_x0000_i1036><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1036 DrawAspect=Content ObjectID=_1336184541></o:OLEObject></xml> . Обслуговування клієнтів біля кас відбувається за експоненціальним законом розподілу випадкової величини <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1037><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1037 DrawAspect=Content ObjectID=_1336184542></o:OLEObject></xml> ( <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1038><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1038 DrawAspect=Content ObjectID=_1336184543></o:OLEObject></xml> - Час обслуговування клієнтів) з щільністю розподілу <img width=33 height=23 src=dopb174747.zip v:shapes=_x0000_i1039><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1039 DrawAspect=Content ObjectID=_1336184546></o:OLEObject></xml> : <br /><img width=47 height=23 src=dopb174748.zip v:shapes=_x0000_i1040><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1040 DrawAspect=Content ObjectID=_1336184547></o:OLEObject></xml><img width=85 height=51 src=dopb174749.zip v:shapes=_x0000_i1041><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1041 DrawAspect=Content ObjectID=_1336184548></o:OLEObject></xml> (1.2) <br /> Примітка: <br /> Якщо в залі є вільні каси, то клієнт стає на обслуговування до найближчої з них (тобто до каси з мінімальним номером). Якщо всі каси зайняті - клієнт стає в чергу до тієї каси, де черга мінімальна. Якщо черги однакові, то клієнт стає в будь-яку з них. <br /> Для вирішення поставленого завдання необхідно розробити алгоритм імітаційного <a href="Моделювання" title="Моделювання">моделювання</a> роботи касового залу за 24-годинний робочий день. А також визначити час простою кас і кількість клієнтів у черзі не обслужених на момент закриття банку. <br /><br clear=all style=page-break-before:always><h1 style="margin:0cm;margin-bottom:.0001pt;text-indent:35.45pt;line-height: 150%;mso-list:l2 level1 lfo2;tab-stops:list 0cm"> 1. Метод рішення задачі </h1> <a href="Імітаційне_моделювання" title="Імітаційне моделювання">Імітаційне моделювання</a> на ЕОМ <a href="Процес" title="Процес">процесу</a> функціонування автоматизованої <a href="Системи_Управління" title="Системи Управління">системи управління</a> роботою касового залу дозволяє отримати чисельний розв'язок поставленої задачі. Суть даного методу наближеного рішення полягає в проведенні низки випадкових випробувань ймовірнісної моделі досліджуваної системи та отриманні сукупності реалізацій випадкових <a href="Процес" title="Процес">процесів</a> зміни стану. <br /> У результаті багаторазового реалізації випадкових процесів визначаються оцінки ймовірності тих чи інших подій і <a href="Середні_Значення" title="Середні Значення">середні значення</a> випадкових величин. Імітаційне моделювання пов'язане з необхідністю відтворення випадкових подій і величин, розподілених за довільним законом. Існує кілька способів генерації випадкових величин і формування їх розподілів. Модель системи <a href="Управління" title="Управління">управління</a> роботою касового залу включає в себе: <br /> · Прихід клієнтів в зал <img width=23 height=23 src=dopb174750.zip v:shapes=_x0000_i1042><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1042 DrawAspect=Content ObjectID=_1336184549></o:OLEObject></xml> ; <br /> · Час обслуговування клієнтів біля кас <img width=23 height=23 src=dopb174751.zip v:shapes=_x0000_i1043><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1043 DrawAspect=Content ObjectID=_1336184550></o:OLEObject></xml> . <br /> За умовою поставленого завдання прихід клієнтів у зал описується пуассоновским потоком з інтенсивністю r. Для кращого <a href="Розуміння" title="Розуміння">розуміння</a> суті розподілу Пуассона необхідно знати основні визначення: <br /> Інтенсивність потоку - середнє число подій, що з'являється в одиницю часу. <br /> Потік - послідовність подій, які настають у випадкові моменти часу. <br /> Закон розподілу Пуассона виражається формулою (1.1). <br /> Будемо моделювати інтервал часу між двома послідовно зайшли в <a href="Банк" title="Банк">банк</a> клієнтами методом Монте-Карло з <a href="Датчики" title="Датчики">датчиком</a> випадкових чисел на інтервалі [0 - 1]. <br /> Сукупність <img width=100 height=24 src=dopb174752.zip v:shapes=_x0000_i1044><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1044 DrawAspect=Content ObjectID=_1336184551></o:OLEObject></xml> незалежних випадкових подій, що утворюють повну групу, характеризується ймовірностями появи кожного з подій <img width=40 height=24 src=dopb174753.zip v:shapes=_x0000_i1045><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1045 DrawAspect=Content ObjectID=_1336184553></o:OLEObject></xml> , Причому <img width=80 height=45 src=dopb174754.zip v:shapes=_x0000_i1046><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1046 DrawAspect=Content ObjectID=_1336184554></o:OLEObject></xml> . Для моделювання цієї сукупності випадкових подій використовується <a href="Генератор_випадкових_чисел" title="Генератор випадкових чисел">генератор випадкових чисел</a>, рівномірно розподілених в інтервалі [0 - 1]. При розподілі відрізка [0 - 1] на n частин, чисельно рівних <img width=40 height=24 src=dopb174753.zip v:shapes=_x0000_i1047><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1047 DrawAspect=Content ObjectID=_1336184555></o:OLEObject></xml> , Виникнення події <img width=16 height=24 src=dopb174755.zip v:shapes=_x0000_i1048><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1048 DrawAspect=Content ObjectID=_1336184556></o:OLEObject></xml> встановлюється шляхом визначення знаходження випадкового числа Х в межах інтервалу при перевірці умови <img width=157 height=45 src=dopb174756.zip v:shapes=_x0000_i1049><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1049 DrawAspect=Content ObjectID=_1336184557></o:OLEObject></xml> , Де <img width=13 height=19 src=dopb174757.zip v:shapes=_x0000_i1050><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1050 DrawAspect=Content ObjectID=_1336184558></o:OLEObject></xml> змінюється від нуля до n. При <img width=36 height=19 src=dopb174758.zip v:shapes=_x0000_i1051><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1051 DrawAspect=Content ObjectID=_1336184559></o:OLEObject></xml> маємо <img width=93 height=23 src=dopb174759.zip v:shapes=_x0000_i1052><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1052 DrawAspect=Content ObjectID=_1336184560></o:OLEObject></xml> ; При <img width=39 height=19 src=dopb174760.zip v:shapes=_x0000_i1053><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1053 DrawAspect=Content ObjectID=_1336184561></o:OLEObject></xml> маємо <img width=173 height=23 src=dopb174761.zip v:shapes=_x0000_i1054><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1054 DrawAspect=Content ObjectID=_1336184562></o:OLEObject></xml> і так далі. При підстановці <img width=16 height=24 src=dopb174755.zip v:shapes=_x0000_i1055><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1055 DrawAspect=Content ObjectID=_1336184563></o:OLEObject></xml> у формулу (1.1) отримаємо: <br /><img width=39 height=19 src=dopb174762.zip v:shapes=_x0000_i1056><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1056 DrawAspect=Content ObjectID=_1336184565></o:OLEObject></xml><img width=56 height=25 src=dopb174763.zip v:shapes=_x0000_i1057><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1057 DrawAspect=Content ObjectID=_1336184566></o:OLEObject></xml> ; <br /><img width=36 height=19 src=dopb174758.zip alt="*" v:shapes=_x0000_i1058><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1058 DrawAspect=Content ObjectID=_1336184567></o:OLEObject></xml><img width=124 height=41 src=dopb174764.zip v:shapes=_x0000_i1059><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1059 DrawAspect=Content ObjectID=_1336184568></o:OLEObject></xml> ; <br /><img width=39 height=19 src=dopb174760.zip alt="*" v:shapes=_x0000_i1060><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1060 DrawAspect=Content ObjectID=_1336184569></o:OLEObject></xml><img width=208 height=44 src=dopb174765.zip v:shapes=_x0000_i1061><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1061 DrawAspect=Content ObjectID=_1336184570></o:OLEObject></xml> і так далі. <br /> Причому <img width=136 height=24 src=dopb174766.zip v:shapes=_x0000_i1062><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1062 DrawAspect=Content ObjectID=_1336184571></o:OLEObject></xml> (Мін.) - максимальна кількість очікування клієнтів. <br /> Так як досвід проводиться багаторазово, то, очевидно, що частота потрапляння випадкових чисел на кожний з відрізків, що визначаються їх довжиною, і <a href="Відповідь" title="Відповідь">відповідає</a> отриманим ймовірностями <br /> <a href="Моделювання" title="Моделювання"> Моделювання</a> часу обслуговування клієнтів біля кас відбувається за експоненціальним законом розподілу, формула якого представлена вище (формула (1.2)). <br /> Час обслуговування клієнтів <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1063><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1063 DrawAspect=Content ObjectID=_1336184572></o:OLEObject></xml> , Як і будь-яка інша випадкова величина, описується <a href="Функції" title="Функції">функцією</a> розподілу <img width=36 height=21 src=dopb174767.zip v:shapes=_x0000_i1064><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1064 DrawAspect=Content ObjectID=_1336184573></o:OLEObject></xml> , Обумовлена як ймовірність <img width=16 height=17 src=dopb174768.zip v:shapes=_x0000_i1065><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1065 DrawAspect=Content ObjectID=_1336184574></o:OLEObject></xml> випадкової події, яка полягає в тому, що час обслуговування клієнтів <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1066><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1066 DrawAspect=Content ObjectID=_1336184575></o:OLEObject></xml> менше деякого заданого часу <img width=9 height=16 src=dopb174769.zip v:shapes=_x0000_i1067><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1067 DrawAspect=Content ObjectID=_1336184576></o:OLEObject></xml> : <br /><img width=101 height=23 src=dopb174770.zip v:shapes=_x0000_i1068><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1068 DrawAspect=Content ObjectID=_1336184577></o:OLEObject></xml><br /> Ця ймовірність розглядається як <a href="функція" title="функція">функція</a> <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1069><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1069 DrawAspect=Content ObjectID=_1336184578></o:OLEObject></xml> у всьому діапазоні можливих значень величини <img width=9 height=16 src=dopb174769.zip v:shapes=_x0000_i1070><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1070 DrawAspect=Content ObjectID=_1336184579></o:OLEObject></xml> .<a href="Функції" title="Функції"> Функція</a> розподілу будь-якої випадкової величини є неубутною функцією часу <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1071><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1071 DrawAspect=Content ObjectID=_1336184580></o:OLEObject></xml> . Зразковий вид <a href="Функції" title="Функції">функції</a> <img width=36 height=21 src=dopb174767.zip v:shapes=_x0000_i1072><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1072 DrawAspect=Content ObjectID=_1336184581></o:OLEObject></xml> дан на малюнку 3. <br /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1074 DrawAspect=Content ObjectID=_1336184582></o:OLEObject><br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1076 DrawAspect=Content ObjectID=_1336184583></o:OLEObject><br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 align=left><tr><td width=144 height=6></td></tr><tr><td></td><td><img width=360 height=185 src=dopb174771.zip v:shapes="_x0000_s1041 _x0000_s1042 _x0000_s1043 _x0000_s1044 _x0000_s1045 _x0000_s1046 _x0000_s1047 _x0000_s1048 _x0000_s1049"></td></tr></table><br /><br style=mso-ignore:vglayout clear=ALL> Рис. 3 - «Функція розподілу експоненціального закону» <br /> Так як значення <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1077><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1077 DrawAspect=Content ObjectID=_1336184585></o:OLEObject></xml> не можуть бути негативними, то <img width=61 height=21 src=dopb174772.zip v:shapes=_x0000_i1078><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1078 DrawAspect=Content ObjectID=_1336184586></o:OLEObject></xml> . При <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1079><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1079 DrawAspect=Content ObjectID=_1336184587></o:OLEObject></xml><img width=36 height=15 src=dopb174773.zip v:shapes=_x0000_i1080><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1080 DrawAspect=Content ObjectID=_1336184588></o:OLEObject></xml> величина <img width=36 height=21 src=dopb174767.zip v:shapes=_x0000_i1081><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1081 DrawAspect=Content ObjectID=_1336184589></o:OLEObject></xml> прагне до одиниці. Таким чином,<a href="Функції" title="Функції"> функція</a> розподілу часу обслуговування клієнтів: <br /><img width=97 height=24 src=dopb174774.zip v:shapes=_x0000_i1082><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1082 DrawAspect=Content ObjectID=_1336184590></o:OLEObject></xml> (1.3) <br /> де <img width=15 height=19 src=dopb174775.zip v:shapes=_x0000_i1083><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1083 DrawAspect=Content ObjectID=_1336184591></o:OLEObject></xml> - Параметр розподілу (середній час обслуговування клієнтів біля каси). <br /> <a href="Відповідь" title="Відповідь"> Відповідно</a> <a href="Щільність" title="Щільність">щільність</a> розподілу: <br /><img width=84 height=24 src=dopb174776.zip v:shapes=_x0000_i1084><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1084 DrawAspect=Content ObjectID=_1336184592></o:OLEObject></xml> (1.4) <br /> Для моделювання часу обслуговування клієнта біля каси проінтегруємо функцію розподілу <img width=36 height=21 src=dopb174767.zip v:shapes=_x0000_i1085><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1085 DrawAspect=Content ObjectID=_1336184593></o:OLEObject></xml> : <br /><img width=411 height=51 src=dopb174777.zip v:shapes=_x0000_i1086><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1086 DrawAspect=Content ObjectID=_1336184594></o:OLEObject></xml><br /><img width=97 height=24 src=dopb174774.zip v:shapes=_x0000_i1087><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1087 DrawAspect=Content ObjectID=_1336184596></o:OLEObject></xml> (1.5) <br /> Від датчика випадкових чисел рівномірно розподілених на інтервалі [0 - 1] отримуємо чергове число Х, яке підставляємо у формулу (1.5) і обчислюємо <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1088><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1088 DrawAspect=Content ObjectID=_1336184597></o:OLEObject></xml> : <br /><img width=83 height=48 src=dopb174778.zip v:shapes=_x0000_i1089><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1089 DrawAspect=Content ObjectID=_1336184598></o:OLEObject></xml><br /><img width=99 height=41 src=dopb174779.zip v:shapes=_x0000_i1090><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1090 DrawAspect=Content ObjectID=_1336184599></o:OLEObject></xml> (1.6) <br /> Зі співвідношення (1.6) знайдемо <img width=13 height=15 src=dopb174746.zip v:shapes=_x0000_i1091><xml><o:OLEObject Type=Embed ProgID=Equation.3 ShapeID=_x0000_i1091 DrawAspect=Content ObjectID=_1336184600></o:OLEObject></xml> <a href="Відповідь" title="Відповідь"> відповідне</a> Х, яке будемо приймати за випадкове число, що означає час обслуговування даної касою. <br /><br clear=all style=page-break-before:always><h1 style="margin:0cm;margin-bottom:.0001pt;text-indent:35.45pt;line-height: 150%;mso-list:l2 level1 lfo2;tab-stops:list 36.0pt"> 2. ПРОГРАМНЕ РІШЕННЯ </h1> Програма імітаційного моделювання роботи касового залу написана на мові C за допомогою середовища розробки Borland C + + 3.1. <br /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> kassa <br /></div></td></tr></table><img width=102 height=42 src=dopb174780.zip alt="Блок-схема: знак завершення: KASSA" v:shapes=_x0000_s1050> Блок-схема імітаційного моделювання роботи касового залу <br /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Введення кол-ва кас m <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Введення lambda <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Введення tau_min і tau_max <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> хв. і мак. час обслуго-вування клієнтів біля кас <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Параметр розподілу Пуассона <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 align=left><tr><td width=179 height=4></td></tr><tr><td></td><td><img width=461 height=155 src=dopb174781.zip v:shapes="_x0000_s1052 _x0000_s1058 _x0000_s1064 _x0000_s1051 _x0000_s1057 _x0000_s1063 _x0000_s1059 _x0000_s1060 _x0000_s1061 _x0000_s1062 _x0000_s1054 _x0000_s1055 _x0000_s1056 _x0000_s1053"></td></tr></table><br /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i:: = 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Введення hours <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Введення MO і SKO <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Час моделювання роботи касового залу в ч. <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> <a href="Математика" title="Математика">Математичне</a> очіку-ня і середн. кв. від-ня <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Створення і обнулення масиву і змінних <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> current_time <= minutes <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> current_time <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i <= m <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> tau [i] <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> tau [i] - <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> t_free + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 2 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 2 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> enter_time_ <= current_time <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> enter + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> enter_time_ = time_to () <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i = 1; INDPAS = 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> tau [i] <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> tau [i] = _tau () <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> INDPAS = 0 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i <= m & & INDPAS = = 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> INDPAS = = 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> refuse + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> INDPAS = = 0 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> serve + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> current_time + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 3 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 3 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 align=left><tr><td width=0 height=0></td><td width=192></td><td width=60></td><td width=102></td><td width=204></td><td width=119></td></tr><tr><td height=207></td><td></td><td colspan=4 align=left valign=top><img width=485 height=207 src=dopb174782.zip v:shapes="_x0000_s1070 _x0000_s1076 _x0000_s1069 _x0000_s1072 _x0000_s1077 _x0000_s1066 _x0000_s1067 _x0000_s1068 _x0000_s1065 _x0000_s1073 _x0000_s1074 _x0000_s1075 _x0000_s1071 _x0000_s1078"></td></tr><tr><td height=2></td></tr><tr><td height=52></td><td colspan=2></td><td align=left valign=top><img width=102 height=52 src=dopb174783.zip v:shapes="_x0000_s1080 _x0000_s1079"></td></tr><tr><td height=1></td></tr><tr><td height=1268></td><td colspan=4 align=left valign=top><img width=558 height=1268 src=dopb174784.zip v:shapes="_x0000_s1083 _x0000_s1087 _x0000_s1089 _x0000_s1088 _x0000_s1086 _x0000_s1091 _x0000_s1099 _x0000_s1098 _x0000_s1097 _x0000_s1103 _x0000_s1102 _x0000_s1106 _x0000_s1112 _x0000_s1111 _x0000_s1105 _x0000_s1104 _x0000_s1110 _x0000_s1109 _x0000_s1108 _x0000_s1107 _x0000_s1117 _x0000_s1115 _x0000_s1114 _x0000_s1113 _x0000_s1116 _x0000_s1119 _x0000_s1118 _x0000_s1096 _x0000_s1095 _x0000_s1101 _x0000_s1100 _x0000_s1125 _x0000_s1124 _x0000_s1093 _x0000_s1092 _x0000_s1123 _x0000_s1094 _x0000_s1090 _x0000_s1084 _x0000_s1085 _x0000_s1122 _x0000_s1120 _x0000_s1134 _x0000_s1126 _x0000_s1133 _x0000_s1132 _x0000_s1135 _x0000_s1136 _x0000_s1131 _x0000_s1130 _x0000_s1137 _x0000_s1140 _x0000_s1138 _x0000_s1141 _x0000_s1139 _x0000_s1150 _x0000_s1148 _x0000_s1149 _x0000_s1147 _x0000_s1145 _x0000_s1144 _x0000_s1146 _x0000_s1151 _x0000_s1152 _x0000_s1153 _x0000_s1161 _x0000_s1156 _x0000_s1155 _x0000_s1154 _x0000_s1157 _x0000_s1159 _x0000_s1160 _x0000_s1162 _x0000_s1142 _x0000_s1143 _x0000_s1158 _x0000_s1164 _x0000_s1163 _x0000_s1167 _x0000_s1172 _x0000_s1166 _x0000_s1171 _x0000_s1165 _x0000_s1170 _x0000_s1169 _x0000_s1168 _x0000_s1173 _x0000_s1185 _x0000_s1184 _x0000_s1176 _x0000_s1175 _x0000_s1174 _x0000_s1183 _x0000_s1182 _x0000_s1181 _x0000_s1180 _x0000_s1179 _x0000_s1178 _x0000_s1177 _x0000_s1186 _x0000_s1189 _x0000_s1188 _x0000_s1187 _x0000_s1129 _x0000_s1128 _x0000_s1121 _x0000_s1081 _x0000_s1082 _x0000_s1127"></td></tr></table><br /><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> tau [i]> 0 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> stand + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Так <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Ні <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i + + <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> i <m <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 4 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> 4 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> serve-= stand; <br /> t_free_av = (float) t_free / m; serve_av = (float) serve / enter; <br /> <a href="Refusenik" title="Refusenik"> refuse_av</a> = (float) refuse / enter; <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Висновок результатів <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> КІНЕЦЬ <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> Коррек. к-ть обслужених <br /> Обчислимо СР час простою <br /> Обчислимо СР % Обслугованих <br /> Обчислити. СР % Необслуженной <br /><br /></div></td></tr></table><table cellpadding=0 cellspacing=0 width=100%><tr><td><div> stand = 0; i = 1 <br /></div></td></tr></table><table cellpadding=0 cellspacing=0 align=left><tr><td width=0 height=0></td><td width=108></td><td width=118></td><td width=64></td><td width=240></td><td width=133></td></tr><tr><td height=14></td><td colspan=3 align=left valign=top><img width=290 height=14 src=dopb174785.zip v:shapes="_x0000_s1199 _x0000_s1192"></td></tr><tr><td height=5></td></tr><tr><td height=55></td><td colspan=2></td><td colspan=2 align=left valign=top><img width=304 height=55 src=dopb174786.zip v:shapes="_x0000_s1190 _x0000_s1191 _x0000_s1193"></td></tr><tr><td height=1></td></tr><tr><td height=402></td><td></td><td colspan=4 align=left valign=top><img width=555 height=402 src=dopb174787.zip v:shapes="_x0000_s1205 _x0000_s1198 _x0000_s1203 _x0000_s1204 _x0000_s1202 _x0000_s1201 _x0000_s1200 _x0000_s1206 _x0000_s1209 _x0000_s1197 _x0000_s1196 _x0000_s1208 _x0000_s1207 _x0000_s1210 _x0000_s1211 _x0000_s1214 _x0000_s1194 _x0000_s1195 _x0000_s1213 _x0000_s1212 _x0000_s1219 _x0000_s1218 _x0000_s1217 _x0000_s1221 _x0000_s1222 _x0000_s1224 _x0000_s1223 _x0000_s1220 _x0000_s1215 _x0000_s1216"></td></tr></table><br /><br style=mso-ignore:vglayout clear=ALL> Оригінальний <a href="Текст" title="Текст">текст</a> програми складається з одного файлу Kas1.c який містить реалізацію таких функцій програми: <br /> - Float RND_DIG (void) - <a href="функція" title="функція">Функція</a> повертає СЧ в діапазоні [0, 1]; <br /> - Void massive (void) - Функція виділяє пам'ять під <a href="Масив" title="Масив">масив</a> часу обслуговування в кожної з кас; <br /> - Float _tau (void) - Функція повертає час обслуговування біля каси; <br /> - Float time_to (void) - Функція визначає час входу наступного клієнта. Проміжки між входами розподілені за законом Пуассона з параметром lambda; <br /> - Void inf (void) - Функція виведення інформації користувачеві. <br /><br clear=all style=page-break-before:always> <b style=mso-bidi-font-weight:normal>КЕРІВНИЦТВО користувач я</b> <br /> 1. <a href="Програма_імітаційного_моделювання_роботи_банку" title="Програма імітаційного моделювання роботи банку">Програма імітаційного моделювання роботи банку</a> розташована за адресою: <br /> A: \ Kas1.exe <br /> 2. На запит програми: <br /> «Введіть кількість кас в касовому залі ж \ \ д вокзалу:» <br /> Вводимо передбачуване (дане) кількість кас обслуговують клієнтів. <br /> 3. На запит програми: <br /> «Введіть параметр розподілу Пуассона для <a href="Визначення_часу" title="Визначення часу">визначення часу</a> між входами пасажирів у зал:» <br /> Вводимо lambda - Параметр розподілу Пуассона <br /> 4. На запит програми: <br /> «Введіть мінімальну тривалість обслуговування біля кас» <br /> Вводимо tau_min - мінімальну тривалість обслуговування. <br /> 4. На запит програми: <br /> «Введіть максимальну тривалість обслуговування біля кас» <br /> Вводимо tau_max - максимальну тривалість обслуговування. <br /> 5. На запит програми: <br /> «Введіть час моделювання роботи залу (в ч.):» <br /> Вводимо hours - час моделювання роботи залу. <br /> 6. На запит програми: <br /> «Введіть <a href="Математика" title="Математика">математичне</a> очікування часу обслуговування біля кас:» <br /> Вводимо MO - математичне очікування. <br /> 7. На запит програми: <br /> «Введіть середньоквадратичне відхилення часу обслуговування біля кас:» <br /> Вводимо SKO - середньоквадратичне відхилення. <br /> 8. Після перегляду результатів, натисніть будь-яку клавішу для завершення роботи. <br /><br clear=all style=page-break-before:always> <b style=mso-bidi-font-weight:normal>5.</b> <b style=mso-bidi-font-weight:normal>Лістинг програми</b> <br /> / * Включає модулі * / <br /> # Include <stdio.h> / / Функції потокового вводу - виводу <br /> # Include <conio.h> / / Формування екракна <br /> # Include <math.h> / / <a href="Математика" title="Математика">Математичні</a> функції <br /> # Include <alloc.h> / / Функції по роботі з динам. пам'яттю <br /> / / Прототипи функцій <br /> float RND_DIG (void); / / Ф-я генерує случ. числа <br /> void inf (void); / / Ф-я повідомляє про приз-ии програми <br /> void massive (void); / / Ф-я виділяє пам'ять під масив <br /> float _tau (void); / / Розрахунок часу обсл-я клієнта біля каси <br /></div><div> float time_to (void); / / Розрахунок часу входу пасажира в зал <br /> / / Оголошення змінних <br /> float <br /> tau_min, / / Мін. час обсл-я біля кас <br /> tau_max, / / Макс. час обсл-я біля кас <br /> * Tau, / / Тек. час обсл-я біля кас <br /> current_time, / / Тек. час <br /> minutes, / / Заг. час роботи у хвилинах <br /> lambda, / / Пар-р закону распр-я Пуассона <br /> enter_time_, / / Час входу чергового пасажира <br /> MO, / / Мат. ожид-е часу обсл-я <br /> SKO, / / Среднекв. вимк-е часу обсл-я біля кас <br /> t_free, / / Заг. час простою кас <br /> t_free_av, / / СР час простою <br /> serve_av, / / СР % Обслуг-я пасажирів <br /> refuse_av, / / СР <a href="Методи_пошуку_відмов" title="Методи пошуку відмов">% Відмов</a> у обсл-ии <br /> N0 = 7836, / / 1-й множник для розрахунку случ. числа <br /> NI = 4792, / / 2-й множник для розрахунку случ. числа <br /> N, / / Прозведеніе N0 і NI <br /> ND; / / Середні 4 <a href="Цифри" title="Цифри">цифри</a> з числа N <br /> int <br /> m, / / Кількість кас <br /> i, / / Тек. номер каси <br /> enter, / / Кількість увійшли <br /> hours, / / Час роботи залу в годинах <br /> serve, / / Кількість обслуг-х пасажирів <br /> refuse, / / Кількість необслуж-х пасажирів <br /> stand, / / Кількість стоять у черзі на момент обслуговування <br /> INDPAS; / / Прапор обслуговування пасажирів <br /> / / Головна функція програми <br /> void main (void) <br /> { <br /> inf (); / / Повідомлення про призначення програми <br /> textcolor (14); / / установка кольору і фону <br /> textbackground (1); <br /> clrscr (); / / <a href="Очищення" title="Очищення">очищення</a> екрана <br /> / * Отримання даних від користувача * / <br /> printf ("\ n █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> for (i = 1; i <= 21; i + +) <br /> printf ("█ █"); <br /> printf ("█ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> gotoxy (10,7); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте кількість кас в касовому" <br /> "Залі ж \ \ д вокзалу:"); <br /> scanf ("% d", & m); <br /> } <br /> while ((m <= 0) | | (m> = 100)); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте параметр розподілу Пуассона для" <br /> "Визначення"); <br /> printf ("\ n █ \ t \ tвремені між входами пасажирів у зал:"); <br /> scanf ("% f", & lambda); <br /> } While (lambda <= 0); <br /> clrscr (); <br /> / / Висновок рамки <br /> printf ("\ n █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> for (i = 1; i <= 21; i + +) <br /> printf ("█ █"); <br /> printf ("█ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> gotoxy (5,3); <br /> do <br /> { <br /> printf ("\ n █ \ n █ \ tВведіте мінімальну тривалість обслуговування"); <br /> printf ("\ n █ \ t \ t \ tу касами:"); <br /> scanf ("% f", & tau_min); <br /> } While (tau_min <= 0); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте максимальну тривалість обслуговування"); <br /> printf ("\ n █ \ t \ t \ tу кас:"); <br /> scanf ("% f", & tau_max); <br /> } While (tau_max <= tau_min); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте час моделювання роботи залу (в ч.):"); <br /> scanf ("% d", & hours); <br /> } While (hours <= 0); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте математичне сподівання часу обслуговування"); <br /> printf ("\ n █ \ t \ tу кас:"); <br /> scanf ("% f", & MO); <br /> } While (MO <= 0); <br /> do <br /> { <br /> printf ("\ n █ \ tВведіте середньоквадратичне відхилення часу"); <br /> printf ("\ n █ \ tобслужіванія біля кас:"); <br /> scanf ("% f", & SKO); <br /> } While (SKO <= 0); <br /> massive (); / / створення і обнулення масиву <br /> minutes = floor (hours * 60); / / переведемо час роботи в хвилини <br /> enter = 0; / / ніхто не входив <br /> serve = 0; / / нікого не обслужили <br /> refuse = 0; / / нікому не відмовили <br /> t_free = 0; / / каси ще не простоювали <br /> current_time = 0; / / початок відліку часу з нуля <br /> enter_time_ = time_to (); / / коли зайде наступний <br /> / * Цикл моделювання роботи залу <br /> for (current_time = 0; current_time <= minutes; current_time + +) <br /> { <br /> / * Коригу-ка часу обсл-я у кожної каси * / <br /> if (current_time) <br /> for (i = 1; i <= m; i + +) <br /> if (tau [i]) / / Якщо каса зайнята, то зменшимо час обсл-я <br /> tau [i] -; / / у ній на 1 хв. <br /> else / / Якщо каса вільна, нарощуємо загальний час <br /> t_free + +; / / простою кас. <br /> if (enter_time_ <= current_time) <br /> { <br /> enter + +; / / Якщо входить пасажир, відзначаємо це <br /> enter_time_ = time_to (); / / Дізнаємося, коли зайде наступний? <br /> i = 1; / / Увійшовши підходить до 1-ї касі <br /> INDPAS = 1; / / Його ще не обслуговують <br /> do <br /> { <br /> if (tau [i] <= 0) / / Якщо каса вільна, ставимо пасажира на <br /> {/ / Обсл-е і дізнаємося, скільки йому стояти біля каси <br /> tau [i] = _tau (); <br /> INDPAS = 0; / / Зазначимо, що пасажира обслуговують <br /> } <br /> i + +; / / Підходимо до наступної касі <br /> } While ((i <= m) & & (INDPAS == 1)); <br /> if (INDPAS == 1) / / Якщо клієнт не обслужений, відзначаємо це <br /> refuse + +; <br /> if (INDPAS == 0) / / Якщо клієнт обслужений, відзначаємо це <br /> serve + +; <br /> } <br /> } <br /> stand = 0; / / Поки в черзі ніхто не варто <br /> i = 1; / / Починаємо з першої каси <br /> do <br /> { <br /> if (tau [i]> 0) / / Якщо час обсл-я не дорівнює нулю, значить <br /> stand + +; / / у каси хтось стоїть <br /> i + +; / / Дивимося у наступної каси <br /> } While (i <= m); <br /> serve-= stand; / / Коректуємо кол-во обслужених <br /> t_free_av = (float) t_free / m; / / Обчислимо СР час простою <br /> serve_av = (float) serve / enter; / / Обчислимо СР % Обслугованих <br /> refuse_av = (float) refuse / enter; / / Обчислимо ср % Необслуженной <br /> textcolor (1); / / Установка параметрів тексту <br /> textbackground (2); <br /> clrscr (); <br /> _setcursortype (_NOCURSOR); <br /><br /> / * Висновок отриманих результатів <br /> / / Висновок рамки <br /> printf ("\ n █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> for (i = 1; i <= 21; i + +) <br /> printf ("█ █"); <br /> printf ("█ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ █ "); <br /> gotoxy (3,3); <br /> printf ("\ n █ \ n █ \ t ╔ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ╗ "); <br /> printf ("\ n █ \ t ║% d залізничних кас працювало протягом:% d годин.", m, hours); <br /> gotoxy (70,6); <br /> printf ("║"); <br /> printf ("\ n █ \ t ║ У зал зайшло% d відвідувачів.", enter); <br /> gotoxy (70,7); <br /> printf ("║"); <br /> printf ("\ n █ \ t ║"); <br /> gotoxy (70,8); <br /> printf ("║"); <br /> printf ("\ n █ \ t ║ \ n █ \ t ║ Було обслуговано% d чоловік (а).", serve); <br /> gotoxy (70,9); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ Пішли необслуженной% d чоловік (а).", refuse); <br /> gotoxy (70,10); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ У черзі залишилося стояти% d чоловік (а).", stand); <br /> gotoxy (70,11); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ Середній відсоток <a href="Обслуговування_пасажирів,_що_вилітають_в_агентстві_та_аеропорту" title="Обслуговування пасажирів, що вилітають в агентстві та аеропорту">обслуговування пасажирів%</a> .2 f%%", serve_av * 100); <br /> gotoxy (70,12); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ Середній відсоток відмов в обслуговуванні% .2 f%%", refuse_av * 100); <br /> gotoxy (70,13); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ Всього каси простоювали:% .0 FЧ.% .0 f хв.", floor (t_free/60), <br /> t_free-floor (t_free/60) * 60); <br /> gotoxy (70,14); <br /> printf ("║"); <br /> printf ("\ n █ \ n █ \ t ║ Середній час простою 1 каси:% .0 FЧ.% .0 f хв.", <br /> floor (t_free_av/60), t_free_av-floor (t_free_av/60) * 60); <br /> gotoxy (70,15); <br /> printf ("║"); <br /> gotoxy (70,16); <br /> printf ("║"); <br /> printf ("\ n █ \ t ╚ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ═ ╝ "); <br /> getch (); <br /> } <br /> / / Функція повертає СЧ в діапазоні від 0 до 1 <br /> float RND_DIG (void) <br /> { <br /> float <br /> x; / / СЧ в діапазоні від 0 до 1 <br /> N = (NI * N0); <br /> ND = floor (N/100); <br /> ND = (ND/10000-floor (ND/10000)) * 10000; / / Відкинемо перші 2 цифри <br /> N0 = NI; <br /> NI = ND; <br /> x = ND/9999; <br /> return x; <br /> } <br /> / / Ф-я виділяє пам'ять під масив часу обсл-я у кожної з кас <br /> void massive (void) <br /> { <br /> int i; <br /> tau = (float *) malloc ((m +1) * sizeof (float)); <br /> / * Обнулення елементів масиву * / <br /> for (i = 1; i <= m; i + +) <br /> tau [i] = 0; <br /> } <br /> / / Ф-я возвр. час обсл-я біля каси <br /> float _tau (void) <br /> { <br /> float <br /> x, / / Випадкове число x <br /> y, / / Випадкове число y <br /> z, / / Випадкове число z <br /> tau_; / / Час обслуговування <br /> / / Перевірка на парність <br /> //--------------------- <br /> if (floor ((float) i / 2) == ceil ((float) i / 2)) <br /> { <br /> x = RND_DIG (); / / Запит случ. числа від ДСЧ <br /> y = RND_DIG (); / / Запит случ. числа від ДСЧ <br /> if (x> =. 0001) / / Якщо х не дуже маленьке, вважаємо z <br /> z = sqrt (-2 * log (x)) * sin (2 * M_PI * y); <br /> else <br /> { <br /> if (sin (2 * M_PI * y) <= 0) <br /> z = 0; / / Якщо обидва множника занадто <a href="Малі" title="Малі">малі</a>, то z = 0 <br /> else <br /> z = (float) minutes / 2; / / Якщо 2-й множ-ль нормальний, то <br /> } / / Вів-на z дорівнює половині часу роботи залу. <br /> tau_ = MO + z * SKO; / / Обчислимо нормально розподілене час обсл-я <br /> if (tau_ <0) / / Якщо час негативно, воно дорівнює 1 хв. <br /> tau_ = 1; <br /> } <br /> else <br /> { <br /> x = RND_DIG (); / / Запит случ. величини від ДСЧ <br /> / * Обчислимо рівномірно розподілене час обсл-я * / <br /> //------------------------------------------------ - <br /> tau_ = floor (tau_min + x * (tau_max-tau_min) + .5); <br /> } <br /> return floor (tau_ + .5); / / Виведемо в місце запиту <br /> } <br /> / / Ф-я визначає час входу слід-го клієнта. Проміжки <br /> / / Між входами розподілені за законом Пуассона з пар-м lambda <br /> / / ================================================ ============ <br /> float time_to (void) <br /> { <br /> int <br /> j; / / Пар-р закону распр-я Пуассона <br /> float <br /> x, / / Случ. вів-на <br /> tau, / / Час до входу <br /> a, / / Пар-р закону распр-я <br /> n, / / Половина загального часу роботи <br /> S; / / Пар-р закону распр-я <br /> n = floor (minutes / 2 +0.5); <br /> x = RND_DIG (); / / Приймемо случ. число <br /> tau = 0; <br /> j = 0; <br /> a = exp ((lambda * (-1))); / / Розрахунок поч. знач-я пар-ра а <br /> S = 0; <br /> do <br /> { <br /> if (S> = x & & x <S + a) / / Перевірка х на пренадл-ть проміжку [S; S + a) <br /> break; <br /> tau + +; <br /> j + +; <br /> S + = a; <br /> a *= lambda / j; <br /> } <br /> while (tau <= n); <br /> return current_time + tau; / / Повернемо час до входу слід-го клієнта <br /> } <br /></div> </div>
Будь ласка, не зберігайте тестовий текст.
Ваш ip: 3.129.13.201 буде збережений.
категорії
за типом
за алфавітом
завантажені
© Усі права захищені
написати до нас
