ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 02.04.2024

Просмотров: 229

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

68

потоками требований. Они определяются моментами поступлений τi и

количеством требований γi , поступающих в момент τi . При этом γi и τi в

общем случае случайны.

Рассмотрим класс потоков, называемых рекуррентными. У рекуррентного потока требований γi =1 для всех i=1,2,... и промежутки времени ξi = τi+1 − τi , i=1,2,..., между последовательными требованиями являются статистически независимыми положительными и одинаково распределенными случайными величинами.

Особенно важен частный случай, когда все ξi одинаковое экспоненциальное распределение с параметром λ. Такой поток называется пуассоновским потоком требований с интенсивностью λ, так как случайное количество требований, поступающих в промежутке времени длительности t, подчиняется пуассоновскому распределению с параметром λt.

Пуассоновский поток требований имеет следующие свойства:

1)стационарность: вероятность появления в r непересекающихся промежутках времени соответственно ν1 ,ν2 ,...,νr требований зависит только от числа и длин промежутков, но не зависит от их расположения на временной оси;

2)отсутствие последействия: количество требований, поступивших в промежутке времени [s, s+t), стохастически не зависит от поведения потока до момента времени s;

3)ординарность: вероятность P >1 ( t) того, что в промежутке времени

длиной t поступит более одного требования, имеет следующее свойство:

lim

P>1 ( t)

= 0, то есть P>1 ( t) = ο( t) .

t

t0

 

Совместное проявление этих трех свойств – стационарности, отсутствия последействия и ординарности – характеризует пуассоновский поток.

Последовательность длительностей обслуживания на любом приборе,


69

как правило, предполагается последовательностью независимых одинаково распределенных положительных случайных величин. Последовательность интервалов обслуживания, естественно, не содержит промежутков времени, когда прибор обслуживания свободен от требований. Особенно важен случай, когда длительности обслуживания имеют одинаковое экспоненциальное распределение с параметром μ.

При рассмотрении систем массового обслуживания будем использовать следующие обозначения:

λ – интенсивность входящего потока требований; μ – интенсивность обслуживания требований одним прибором; æ – число обслуживающих приборов в системе;

ψ – коэффициент использования обслуживающих приборов системы; b - математическое ожидание (м.о.) числа требований в очереди;

g - м.о. числа свободных приборов в СМО;

Pn - стационарная вероятность пребывания в СМО точно n требований.

4.2. Постановка задачи и метод решения

Рассмотрим задачу. На автомобильном заводе сложилась такая ситуация, что нельзя предоставить рабочим одновременно все необходимые для работы инструменты. Разнообразные инструменты имеются на складе. Поэтому при получении инструментов наблюдается образование очередей рабочих у склада. Очевидно, следует уменьшить время ожидания рабочих в очереди, так как оно потеряно для производства. Один кладовщик может обслуживать одновременно только одного рабочего. Если кладовщиков слишком много, то очереди рабочих не будет, но невыгодно платить «простаивающим» кладовщикам. Если кладовщиков недостаточно, то будут образовываться длинные очереди. Таким образом, возникает следующая задача: определить оптимальное количество кладовщиков в том смысле, чтобы время, потерянное рабочими, с одной стороны, и кладовщиками, с

70

другой, приводило бы к минимальным затратам, если себестоимость часа рабочего равна 6 руб., а себестоимость часа кладовщика – 3 руб.

Исследование начинается с определения статистических характеристик процесса поступления рабочих (требований) в кладовую и времени, затрачиваемого кладовщиками (приборами) на их обслуживание. Изучение входящего потока производится следующим образом. Каждые 10 минут в течение 100 последовательных десятиминутных интервалов отмечается число рабочих, пришедших в инструментальную кладовую за получением инструментов. Вычисляются частоты, соответствующие наблюдаемым числам. Результаты приведены в табл.1.

 

 

 

 

 

Таблица 1.

 

 

 

 

 

 

 

Число поступив ших требован ий за 10 мин. интервал времени

Наблюда емая частота

Частота по закону Пуассона

Число поступив ших требован ий за 10 мин. интервал времени

Наблюда емая частота

Частота по закону Пуассона

 

 

 

 

 

 

 

 

5

0,01

0,001

16

0,12

0,099

 

6

0

0,002

17

0,08

0,093

 

7

0,01

0,006

18

0,09

0,083

 

8

0,02

0,012

19

0,07

0,069

 

9

0,01

0,021

20

0,05

0,055

 

10

0,03

0,034

21

0,04

0,042

 

11

0,05

0,049

22

0,03

0,031

 

12

0,06

0,066

23

0,01

0,021

 

13

0,09

0,081

24

0,01

0,014

 

14

0,10

0,093

25

0,01

0,009

 

15

0,11

0,099

 

 

 

 

 

 

 

 

 

 

 

Определим среднее значение числа поступлений за 10 минут, эту величину обозначим через λ:

λ = 5 0.01 + 6 0 + 7 0.01 + ... + 25 0.01 1.6 ,

то есть в среднем за 10 минут приходят 16 рабочих или 1.6 рабочих в минуту. В третьем столбце табл.1 приведены частоты, соответствующие распределению Пуассона, согласно которому

P (t) =

e−λt (λt)n

,

 

n

n!

 

 

 


71

где Pn (t) - вероятность того, что за время t произойдет n поступлений при

λ=1.6 и t=10. Затем с помощью критерия χ2 проверяется приемлемость гипотезы о пуассоновском распределении исследуемой совокупности.

Оказывается, что при числе степеней свободы, равном 19, χ2 =12.

Вероятность того, что экспериментальное распределение является пуассоновским, больше 0.88. Таким образом, можно допустить, что распределение соответствует закону Пуассона с λ=1.6.

Для измерения длительности обслуживания используют счетчик. Счетчик включают в начале обслуживания и выключают в конце. Таким образом регистрируется продолжительность 1000 обслуживаний. Затем вычисляются частоты, соответствующие интервалам 0 – 15, 16 – 30, 31 – 45,

..., 301 – 315 секунд, которые приведены в табл. 2.

Используя значения наблюдаемых частот, вычислим среднее значение длительности обслуживания ν =1.1 минуты. Отсюда следует, что интенсивность обслуживания μ=1/1.1=0.9 требований в минуту. В третьем столбце табл.2 приведены значения статистической функции распределения, которая строится по наблюдаемым частотам. Четвертый столбец табл.2 содержит значения вероятностей, соответствующих экспоненциальному распределению P{ξ < t}=1 e−μt при μ=0.9, где ξ – случайная величина,

длительность обслуживания требования. В данном случае критерий χ2 дает

величину χ2 =1.91 при числе степеней свободы, равном 19. Вероятность того,

что гипотеза об экспоненциальном распределении справедлива, больше 0.99. Таким образом, можно допустить, что в рассматриваемом случае имеет место экспоненциальное распределение с параметром μ=0.9.

 

 

 

Таблица 2.

 

 

 

 

Интервал времени

Наблюдаемая

Статистическая

Экспоненциальная

в секундах

частота

функция

функция

 

 

распределения

распределения

15 (1 – 15)

0,187

0,0187

0,202

30 (16 – 30)

0,161

0,348

0,362

45 (31 – 45)

0,140

0,488

0,491


72

60 (46

– 60)

0,104

0,592

0,593

75 (61

– 75)

0,078

0,670

0,675

90 (76

– 90)

0,069

0,739

0,740

105 (91

– 105)

0,051

0,790

0,793

120 (106 – 120)

0,047

0,837

0,835

135 (121 – 135)

0,038

0,875

0,868

150 (136 – 150)

0,030

0,905

0,895

165 (151 – 165)

0,016

0,921

0,916

180 (166 – 180)

0,017

0,938

0,933

195 (181 – 195)

0,011

0,949

0,946

210 (196 – 210)

0,007

0,956

0,957

225 (211 – 225)

0,009

0,965

0,966

240 (226 – 240)

0,009

0,974

0,973

255 (241 – 255)

0,005

0,979

0,978

270 (256 – 270)

0,004

0,983

0,983

285 (271 – 285)

0,004

0,987

0,986

300 (286 – 300)

0,003

0,990

0,989

315 (301 – 315)

0,010

1

0,999

 

 

 

 

 

В качестве модели рассматриваемой системы будем использовать систему массового обслуживания. Поток рабочих на склад соответствует входящему потоку требований в СМО. Длительность обслуживания рабочего кладовщиком соответствует длительности обслуживания требования прибором обслуживания. Число кладовщиков на складе соответствует числу обслуживающих приборов в СМО.

Как было показано выше, входящий поток требования в систему является пуассоновским с интенсивностью λ, а длительность обслуживания требований на каждом приборе обслуживания представляет собой последовательность независимых одинаково распределенных экспоненциально случайных величин.

Определим функцию потерь R(æ) в единицу времени следующим образом:

R(æ)=6

b

+ 3

g

,

(1)

где æ – число обслуживающих приборов СМО, b - м.о. числа требований в очереди СМО (рабочих), g - м.о. числа свободных приборов (кладовщиков).

Из теории массового обслуживания известно, что


73

 

 

 

 

 

κκψκ+1

 

 

 

 

 

 

b =

P ,

 

 

 

 

κ!(1 − ψ)2

 

 

 

 

 

 

 

0

 

 

 

 

 

 

(κψ)κ

κ−1

(κψ)n 1

 

P0

 

 

 

 

 

+

 

 

,

 

 

 

 

 

=

κ!(1 − ψ)

n!

 

 

 

n=0

 

 

g = (1 − ψ) æ,

где g = (1 − ψ) æ.

Определим значение числа обслуживающих приборов æ, обеспечивающее минимум функции потерь. Можно показать, что функция R(æ) имеет только один минимум. Определим минимум этой функции путем ряда последовательных приближений. Вычислим м.о. числа требований в очереди при æ, равном 2, 3, 4, 5:

 

 

 

 

 

 

 

 

æ =2 , b =6,693;

æ =4, b =0,099;

æ =3 ,

 

 

=0,502;

æ =5,

 

 

=0,021.

b

b

Вычислим м.о. числа свободных приборов при тех же значениях æ:

æ =2 ,

 

=0,222;

æ =4,

 

=2,222;

g

g

æ =3 ,

 

=1,222;

æ =5,

 

=3,222.

g

g

Используя формулу (1), вычислим функцию потерь для различных

значений æ:

 

 

 

 

R(2)=40.824;

R(4)=7.262;

R(3)=6.68;

R(5)=9.795.

Отсюда видно, что минимальное значение функции потерь достигается при æ=3. Таким образом, оптимальное число кладовщиков равно трем.

4.3. Вопросы для самопроверки

1.Что такое обслуживающая система?

2.Что такое требование?