Главная | Обратная связь | Поможем написать вашу работу!
МегаЛекции

Лабораторная работа №7 Многоканальная СМО с ограниченной очередью




Теория.

Рассмотрим S-канальную СМО с ожиданием, на которую поступает поток заявок с интенсивностью . Интенсивность обслуживания µ (для одного канала), число мест в очереди m. Заявка, поступившая в момент, когда все каналы заняты и очередь заполнена, получает отказ на обслуживание.

Граф состояний такой системы выглядит следующим образом:

E0
λ
µ
E1
Es
Es+m
λ
λ
λ
λ
Ek
λ
λ
(k+1)µ
очереди нет

Рисунок 7.1 Граф состояний СМО < µ | µ | S | m >

 

Составив уравнения Колмогорова, прировняв их нулю, решив данную систему и преобразовав решения с помощью условия нормировки и формул геометрической прогрессии, получим формулы (7.1)-(7.2).

 

(7.1)

 

(7.2)

 

По формуле (7.1) можно высчитать вероятность отказа, т.е. вероятность того, что все каналы заняты и вся очередь наполнена:

(7.3)

 

 

Относительная пропускная способность системы q будет являться противоположным событием событию «отказа заявки на обслуживание», и высчитывается по формуле.

(7.4)

 

 

Абсолютная пропускная способность (А) - среднее число заявок, которое может обслужить СМО в единицу времени:

(7.5)

 

Среднее число занятых каналов будет находиться по формуле (7.6).

(7.6)

 

Среднее число находящихся в очереди заявок:

(7.7)

 

 

Среднее время пребывания заявки в очереди:

(7.8)

 

 

Среднее число приборов свободных от работы:

(7.9)

 

 

Цель. Решить оптимизационную задачу на примере СМО < µ | µ | S | m >.

Входные значения: параметр λ=10шт/сутки, µ=2шт/сутки.

 

Задача. В отделе нагрева металла в цехе крупной плавки часть работ происходит в режиме «копильника»: если слиток застает все печи занятыми, то он помещается в «копильник», где ему обеспечивается нужная температура. Если все печи и «копильники» заняты - он отправляется на склад. Но при этом для того, чтобы его заново разогреть, потребуются дополнительные растраты в размере 100$ на его разогрев. Поступающий поток простейший с λ=10шт/сутки. Интенсивность нагрева слитков перед ковкой (распределение показательное) µ=2шт/сутки. В цехе имеется 10 печей, из которых часть должна быть для плавки, а часть для «копильника». Цена одной печи -100000$. Текущие затраты 50$ в сутки на обслуживание одной печи. Затраты на простой печи 30$ сутки. Затраты на содержание слитков в «копильнике» 60$ сутки на один слиток. Фонд времени в году 6000 часов. Найти оптимальное количество «копильников»(m) и печей для плавки(s).

Алгоритм решения.

1. Формула для расчета затрат:

 

(7.10)

 

 

где – затраты на работу системы;

- коэффициент эффективности капиталовложения (0,15$/год);

- цена одной печи (100000$);

- текущие затраты на содержание одной печи (50$/сутки);

– текущие затраты на содержание одной печи в простое (30$/сутки);

– затраты на содержание слитка в «копильнике» (60$/сутки шт.);

– затраты на разогрев отосланных на склад слитков (100$/сутки шт.);

– средняя длина очереди ();

–среднее число приборов свободных от работы;

–годовой фонд рабочего времени (6000 часов).

1. Протабулировать функцию (7.10) при (m=10-S);

 

2. Построить график зависимости затрат от числа обслуживающих каналов, и найти оптимальное число обслуживающих каналов для данной задачи.

 

Не забывайте проверять условие существования стационарного режима для текущей итерации и единицы измерения.


ПРИЛОЖЕНИЕ 1
Таблица значений квантилей

n\ε 0,01 0,025 0,05 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0,95 0,975 0,99
  0,0002 0,0010 0,0039 0,0158 0,0642 0,1485 0,2750 0,4549 0,7083 1,0742 1,6424 2,7055 3,8415 5,0239 6,6349
  0,0201 0,0506 0,1026 0,2107 0,4463 0,7133 1,0217 1,3863 1,8326 2,4079 3,2189 4,6052 5,9915 7,3778 9,2103
  0,1148 0,2158 0,3518 0,5844 1,0052 1,4237 1,8692 2,3660 2,9462 3,6649 4,6416 6,2514 7,8147 9,3484 11,3449
  0,2971 0,4844 0,7107 1,0636 1,6488 2,1947 2,7528 3,3567 4,0446 4,8784 5,9886 7,7794 9,4877 11,1433 13,2767
  0,5543 0,8312 1,1455 1,6103 2,3425 2,9999 3,6555 4,3515 5,1319 6,0644 7,2893 9,2364 11,0705 12,8325 15,0863
  0,8721 1,2373 1,6354 2,2041 3,0701 3,8276 4,5702 5,3481 6,2108 7,2311 8,5581 10,6446 12,5916 14,4494 16,8119
  1,2390 1,6899 2,1673 2,8331 3,8223 4,6713 5,4932 6,3458 7,2832 8,3834 9,8032 12,0170 14,0671 16,0128 18,4753
  1,6465 2,1797 2,7326 3,4895 4,5936 5,5274 6,4226 7,3441 8,3505 9,5245 11,0301 13,3616 15,5073 17,5345 20,0902
  2,0879 2,7004 3,3251 4,1682 5,3801 6,3933 7,3570 8,3428 9,4136 10,6564 12,2421 14,6837 16,9190 19,0228 21,6660
  2,5582 3,2470 3,9403 4,8652 6,1791 7,2672 8,2955 9,3418 10,4732 11,7807 13,4420 15,9872 18,3070 20,4832 23,2093
  3,0535 3,8157 4,5748 5,5778 6,9887 8,1479 9,2373 10,3410 11,5298 12,8987 14,6314 17,2750 19,6751 21,9200 24,7250
  3,5706 4,4038 5,2260 6,3038 7,8073 9,0343 10,1820 11,3403 12,5838 14,0111 15,8120 18,5493 21,0261 23,3367 26,2170
  4,1069 5,0088 5,8919 7,0415 8,6339 9,9257 11,1291 12,3398 13,6356 15,1187 16,9848 19,8119 22,3620 24,7356 27,6882
  4,6604 5,6287 6,5706 7,7895 9,4673 10,8215 12,0785 13,3393 14,6853 16,2221 18,1508 21,0641 23,6848 26,1189 29,1412
  5,2293 6,2621 7,2609 8,5468 10,3070 11,7212 13,0297 14,3389 15,7332 17,3217 19,3107 22,3071 24,9958 27,4884 30,5779
  5,8122 6,9077 7,9616 9,3122 11,1521 12,6243 13,9827 15,3385 16,7795 18,4179 20,4651 23,5418 26,2962 28,8454 31,9999
  6,4078 7,5642 8,6718 10,0852 12,0023 13,5307 14,9373 16,3382 17,8244 19,5110 21,6146 24,7690 27,5871 30,1910 33,4087
  7,0149 8,2307 9,3905 10,8649 12,8570 14,4399 15,8932 17,3379 18,8679 20,6014 22,7595 25,9894 28,8693 31,5264 34,8053
  7,6327 8,9065 10,1170 11,6509 13,7158 15,3517 16,8504 18,3377 19,9102 21,6891 23,9004 27,2036 30,1435 32,8523 36,1909
  8,2604 9,5908 10,8508 12,4426 14,5784 16,2659 17,8088 19,3374 20,9514 22,7745 25,0375 28,4120 31,4104 34,1696 37,5662
  8,8972 10,2829 11,5913 13,2396 15,4446 17,1823 18,7683 20,3372 21,9915 23,8578 26,1711 29,6151 32,6706 35,4789 38,9322
  9,5425 10,9823 12,3380 14,0415 16,3140 18,1007 19,7288 21,3370 23,0307 24,9390 27,3015 30,8133 33,9244 36,7807 40,2894
  10,1957 11,6886 13,0905 14,8480 17,1865 19,0211 20,6902 22,3369 24,0689 26,0184 28,4288 32,0069 35,1725 38,0756 41,6384
  10,8564 12,4012 13,8484 15,6587 18,0618 19,9432 21,6525 23,3367 25,1063 27,0960 29,5533 33,1962 36,4150 39,3641 42,9798
  11,5240 13,1197 14,6114 16,4734 18,9398 20,8670 22,6156 24,3366 26,1430 28,1719 30,6752 34,3816 37,6525 40,6465 44,3141
  12,1981 13,8439 15,3792 17,2919 19,8202 21,7924 23,5794 25,3365 27,1789 29,2463 31,7946 35,5632 38,8851 41,9232 45,6417
  12,8785 14,5734 16,1514 18,1139 20,7030 22,7192 24,5440 26,3363 28,2141 30,3193 32,9117 36,7412 40,1133 43,1945 46,9629
  13,5647 15,3079 16,9279 18,9392 21,5880 23,6475 25,5093 27,3362 29,2486 31,3909 34,0266 37,9159 41,3371 44,4608 48,2782
  14,2565 16,0471 17,7084 19,7677 22,4751 24,5770 26,4751 28,3361 30,2825 32,4612 35,1394 39,0875 42,5570 45,7223 49,5879
  14,9535 16,7908 18,4927 20,5992 23,3641 25,5078 27,4416 29,3360 31,3159 33,5302 36,2502 40,2560 43,7730 46,9792 50,8922
  15,6555 17,5387 19,2806 21,4336 24,2551 26,4397 28,4087 30,3359 32,3486 34,5981 37,3591 41,4217 44,9853 48,2319 52,1914
  16,3622 18,2908 20,0719 22,2706 25,1478 27,3728 29,3763 31,3359 33,3809 35,6649 38,4663 42,5847 46,1943 49,4804 53,4858

 


 

Продолжение таблицы значений квантилей

  17,0735 19,0467 20,8665 23,1102 26,0422 28,3069 30,3444 32,3358 34,4126 36,7307 39,5718 43,7452 47,3999 50,7251 54,7755
  17,7891 19,8063 21,6643 23,9523 26,9383 29,2421 31,3130 33,3357 35,4438 37,7954 40,6756 44,9032 48,6024 51,9660 56,0609
  18,5089 20,5694 22,4650 24,7967 27,8359 30,1782 32,2821 34,3356 36,4746 38,8591 41,7780 46,0588 49,8018 53,2033 57,3421
  19,2327 21,3359 23,2686 25,6433 28,7350 31,1152 33,2517 35,3356 37,5049 39,9220 42,8788 47,2122 50,9985 54,4373 58,6192
  19,9602 22,1056 24,0749 26,4921 29,6355 32,0532 34,2216 36,3355 38,5348 40,9839 43,9782 48,3634 52,1923 55,6680 59,8925
  20,6914 22,8785 24,8839 27,3430 30,5373 32,9919 35,1920 37,3355 39,5643 42,0451 45,0763 49,5126 53,3835 56,8955 61,1621
  21,4262 23,6543 25,6954 28,1958 31,4405 33,9315 36,1628 38,3354 40,5935 43,1053 46,1730 50,6598 54,5722 58,1201 62,4281
  22,1643 24,4330 26,5093 29,0505 32,3450 34,8719 37,1340 39,3353 41,6222 44,1649 47,2685 51,8051 55,7585 59,3417 63,6907
  22,9056 25,2145 27,3256 29,9071 33,2506 35,8131 38,1055 40,3353 42,6506 45,2236 48,3628 52,9485 56,9424 60,5606 64,9501
  23,6501 25,9987 28,1440 30,7654 34,1574 36,7550 39,0774 41,3352 43,6786 46,2817 49,4560 54,0902 58,1240 61,7768 66,2062
  24,3976 26,7854 28,9647 31,6255 35,0653 37,6975 40,0496 42,3352 44,7063 47,3390 50,5480 55,2302 59,3035 62,9904 67,4593
  25,1480 27,5746 29,7875 32,4871 35,9743 38,6408 41,0222 43,3352 45,7336 48,3957 51,6389 56,3685 60,4809 64,2015 68,7095
  25,9013 28,3662 30,6123 33,3504 36,8844 39,5847 41,9950 44,3351 46,7607 49,4517 52,7288 57,5053 61,6562 65,4102 69,9568
  26,6572 29,1601 31,4390 34,2152 37,7955 40,5292 42,9682 45,3351 47,7874 50,5071 53,8177 58,6405 62,8296 66,6165 71,2014
  27,4158 29,9562 32,2676 35,0814 38,7075 41,4744 43,9417 46,3350 48,8139 51,5619 54,9056 59,7743 64,0011 67,8206 72,4433
  28,1770 30,7545 33,0981 35,9491 39,6205 42,4201 44,9154 47,3350 49,8401 52,6161 55,9926 60,9066 65,1708 69,0226 73,6826
  28,9406 31,5549 33,9303 36,8182 40,5344 43,3664 45,8895 48,3350 50,8660 53,6697 57,0786 62,0375 66,3386 70,2224 74,9195
  29,7067 32,3574 34,7643 37,6886 41,4492 44,3133 46,8638 49,3349 51,8916 54,7228 58,1638 63,1671 67,5048 71,4202 76,1539

 

 


ПРИЛОЖЕНИЕ 2
Таблица значений

n/ε 0.40 0.25 0.10 0.05 0.025 0.01 0.005 0.0005
  0.324920 1.000000 3.077684 6.313752 12.70620 31.82052 63.65674 636.6192
  0.288675 0.816497 1.885618 2.919986 4.30265 6.96456 9.92484 31.5991
  0.276671 0.764892 1.637744 2.353363 3.18245 4.54070 5.84091 12.9240
  0.270722 0.740697 1.533206 2.131847 2.77645 3.74695 4.60409 8.6103
  0.267181 0.726687 1.475884 2.015048 2.57058 3.36493 4.03214 6.8688
  0.264835 0.717558 1.439756 1.943180 2.44691 3.14267 3.70743 5.9588
  0.263167 0.711142 1.414924 1.894579 2.36462 2.99795 3.49948 5.4079
  0.261921 0.706387 1.396815 1.859548 2.30600 2.89646 3.35539 5.0413
  0.260955 0.702722 1.383029 1.833113 2.26216 2.82144 3.24984 4.7809
  0.260185 0.699812 1.372184 1.812461 2.22814 2.76377 3.16927 4.5869
  0.259556 0.697445 1.363430 1.795885 2.20099 2.71808 3.10581 4.4370
  0.259033 0.695483 1.356217 1.782288 2.17881 2.68100 3.05454 4.3178
  0.258591 0.693829 1.350171 1.770933 2.16037 2.65031 3.01228 4.2208
  0.258213 0.692417 1.345030 1.761310 2.14479 2.62449 2.97684 4.1405
  0.257885 0.691197 1.340606 1.753050 2.13145 2.60248 2.94671 4.0728
  0.257599 0.690132 1.336757 1.745884 2.11991 2.58349 2.92078 4.0150
  0.257347 0.689195 1.333379 1.739607 2.10982 2.56693 2.89823 3.9651
  0.257123 0.688364 1.330391 1.734064 2.10092 2.55238 2.87844 3.9216
  0.256923 0.687621 1.327728 1.729133 2.09302 2.53948 2.86093 3.8834
  0.256743 0.686954 1.325341 1.724718 2.08596 2.52798 2.84534 3.8495
  0.256580 0.686352 1.323188 1.720743 2.07961 2.51765 2.83136 3.8193
  0.256432 0.685805 1.321237 1.717144 2.07387 2.50832 2.81876 3.7921
  0.256297 0.685306 1.319460 1.713872 2.06866 2.49987 2.80734 3.7676
  0.256173 0.684850 1.317836 1.710882 2.06390 2.49216 2.79694 3.7454
  0.256060 0.684430 1.316345 1.708141 2.05954 2.48511 2.78744 3.7251
  0.255955 0.684043 1.314972 1.705618 2.05553 2.47863 2.77871 3.7066
  0.255858 0.683685 1.313703 1.703288 2.05183 2.47266 2.77068 3.6896
  0.255768 0.683353 1.312527 1.701131 2.04841 2.46714 2.76326 3.6739
  0.255684 0.683044 1.311434 1.699127 2.04523 2.46202 2.75639 3.6594
  0.255605 0.682756 1.310415 1.697261 2.04227 2.45726 2.75000 3.6460
0.253347 0.674490 1.281552 1.644854 1.95996 2.32635 2.57583 3.2905

ПРИЛОЖЕНИЕ 3
Алгоритм решения задачи моделирования работы СМО


/// <summary>

/// класс заявка

/// </summary>

class Item

{

//индекс

public int index;

//и ее значение

public double value;

 

public Item(int i, double v)

{

index = i;

value = v;

}

}

/// <summary>

/// функция, моделирующая выполнение и отклонение заявок

/// </summary>

/// <param name="T">поток заявок</param>

/// <param name="b">число мест в очереди</param>

/// <param name="tb">время обслуживания</param>

/// <param name="QTask">возвращаемый поток выполн.заявок</param>

/// <param name="queue">очередь</param>

/// <param name="WorkF">выполн.заявки</param>

/// <param name="leaveTask">отклоненные заявки</param>

void Work(List<double> T, int b, double tb, out List<Item> QTask,out Queue<Item> queue,out List<Item> WorkF, out List<Item> leaveTask)

{

QTask = null; queue = null; WorkF = null; leaveTask = null;

if (T == null) return;

if (T.Count == 0) return;

QTask = new List<Item>();//обработка заявки

queue = new Queue<Item>();// очередь

WorkF = new List<Item>(); //выполненные и вышедшие заявки

leaveTask = new List<Item>();//отклоненные заявки

//первая зааявка сразу уходит на выполнение

WorkF.Add(new Item(0, T[0] + tb));

QTask.Add(new Item(0, T[0]));

for (int i = 1; i < T.Count; i++)

{

//берем время вхождения текущей заявки в СМО

//если есть заявка, которая выполняется до ее прихода, выполняем ее

while (WorkF[WorkF.Count - 1].value <= T[i])

{

//проверяем, есть ли заявка в очереди

if (queue.Count!= 0)

{

//забираем ее

Item buff = queue.Dequeue();

//добавляем в массив на выполнение

//причем задаем время ее завершения

WorkF.Add(new Item(buff.index,

//если заявка давно висит в очереди

buff.value < WorkF[WorkF.Count - 1].value?

//она выполняется сразу после предыдущей заявки

WorkF[WorkF.Count - 1].value + tb:

//если она только пришла, а очередь пустовала,

//то она выполняется после своего прихода

buff.value + tb));

 

}

else break;

}

//если в очереди есть места, сразу добавляем в нее заявку

if (queue.Count < b)

{

queue.Enqueue(new Item(i, T[i]));

QTask.Add(new Item(i, T[i]));

}

else

{

//иначе заявка отклоняется

leaveTask.Add(new Item(i, T[i]));

}

}

//когда прошли по всем заявкам, а в очереди еще что-то осталось, завершаем и их

while (queue.Count!= 0)

{

Item buff = queue.Dequeue();

WorkF.Add(new Item(buff.index,

buff.value < WorkF[WorkF.Count - 1].value?

WorkF[WorkF.Count - 1].value + tb:

buff.value +tb));

 

}

}

 

Поделиться:





Воспользуйтесь поиском по сайту:



©2015 - 2024 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...