книги из ГПНТБ / Шабалин Н.Н. Оптимизация процесса переработки вагонов на станциях
.pdfное (экспоненциальное), эрланговские, гамма-распреде ление, гнперэ'кспоненциалы-юе, вейбулловское и нор мальное. Часто, однако, для .решения практических за дач достаточно иметь только числовые характеристики интервалов в потоках, такие как среднее значение (ма тематическое ожидание), среднее квадрэтическое откло нение, коэффициент вариации.
Коэффициент вариации интервалов характеризует степень неравномерности потока. Он равен отношению среднего квадратнческого отклонения к математическо му ожиданию и выражает относительное отклонение интервалов от среднего значения. При регулярном по токе, когда все интервалы равны, коэффициент вариа ции равен нулю, так как нет отклонения интервалов от среднего значения. С увеличением неравномерности по тока возрастает и коэффициент вариации, ухудшаются условия использования устройств станции и увеличива
ется потребность в резервах мощности для |
смягчения |
|
влияния неравномерности. |
|
|
П о к а з а т е л ь н о е |
(экспоненциальное) |
распреде |
ление определяет закономерности появления тех или иных интервалов в простейшем потоке однородных со бытий. Простейший поток, как известно, обладает свой ствами стационарности, ординарности и отсутствия по следействия.
Под стационарностью подразумевается свойство по
тока, при котором вероятность поступления |
определен |
||||
ного числа |
требований |
(поездов, |
вагонов) |
в течение |
|
некоторого |
отрезка |
времени |
зависит |
только |
от |
продолжительности этого отрезка времени, числа требо ваний и не зависит от положения этого отрезка на оси времени.
Ординарность потока означает практическую невоз
можность появления двух |
и |
более требований в |
тече |
ние малого отрезка времени dt. Это утверждение |
того, |
||
что требования поступают |
по |
одному, а не по два |
или |
больше. |
|
|
|
Отсутствие последействия заключается в том, что ве роятность поступления определенного числа требований за принятый отрезок времени не зависит от того, сколь ко требований поступило до этого. Отсутствие последей ствия определяет взаимную независимость событий в по токе.
10
С простейшим потоком |
связало |
р а с п р е д е л е н и е |
||||
П у а с с о н а |
|
|
|
|
|
|
|
|
Pa(t)-SWLe-\ |
' |
|
(1) |
|
где Pn(t) |
— |
вероятность появления |
п событий |
(требо |
||
п — |
ваний) за период t; |
|
|
|
||
целочисленное количество событий, появля |
||||||
X — |
ющихся за период времени t {п=0, |
1, 2,...); |
||||
интенсивность |
потока |
событий |
(среднее |
|||
|
|
количество требований |
потока, |
поступа |
||
|
|
ющих за единицу времени). |
|
|
||
Зависимость (1) определяет вероятность того, что |
за |
|||||
время [0, |
t] |
поступит п требований |
потока, например |
за |
||
время t поступит п поездов |
или завершится накопление |
|||||
п составов |
в сортировочном |
парке. |
|
|
|
Плотность распределения интервалов между собы тиями (дифференциальная функция распределения) в этом случае определится по формуле
/ (t) = |
le~xt; |
(t>0), |
f(t) = 0, для t > |
0)- |
(2) |
Известно |
[ 2 ] , |
что для |
пуассоновского |
потока |
интен |
сивности % интервалы между двумя соседними события ми потока независимы и распределены показательно с плотностью распределения (2). Дискретное пуассоновское распределение числа событий за определенное время t по формуле (1) и непрерывное показательное распреде ление независимых интервалов между событиями по фор муле (2) выражают один и тот же простейший поток событий.
Среднее значение (математическое ожидание) ин тервалов t определяется как величина, обратная интен сивности .потока,
Дисперсия интервалов будет равна
а среднее квадратическое отклонение как корень квадрат ный из дисперсии
=J _
и
Коэффициент вариации интервалов в простейшем по токе равен единице. Это значит, что интервалы отклоня
ются |
от среднего значения на 100%. |
Так, если |
на стан |
|||
цию поступают поезда с интенсивностью |
i = 3 поезда |
в |
||||
час, |
то средняя величина интервалов между ними f = 4 - |
ч |
||||
или 20 мин. |
|
|
|
О |
|
|
При показательном распределении |
интерва |
|||||
лов |
среднее |
квадратическое отклонение |
от |
среднего |
||
значения также будет равно 20 мин. |
На рис. 1 представ |
лен график непрерывного распределения интервалов по
показательному закону между прибывающими |
поездами |
|
при Я = 3. На этом графике по |
оси абсцисс |
отмечено |
среднее значение интервалов. |
|
|
В ряде случаев может иметь |
место так называемое |
сдвинутое показательное распределение, когда интерва лы не могут быть меньше определенной величины ta, на пример при поступлении поездов только с одного участ ка, где минимальный интервал определится условиями пропускной способности. В этом случае плотность рас пределения .интервалов
( 1е->- с-'.) |
при |
/ > * 0 |
|
|
/ " Н о . |
• • • |
п Р И |
, « , • |
<2 а > |
Средняя величина интервалов |
при этом будет |
|
||
? |
- ' ° + |
т - |
|
|
График плотности 'вероятности сдвинутого показа тельного распределения интервалов приведен на рис. 2 при / 0 = 6 мин (0,1 ч). Площадь под кривой плотности распределения, как известно, равна единице.
Интегральная функция показательного распределе ния
|
|
|
t |
|
|
|
F(t) = |
Р(Т |
<t)= |
l\e-udt=\ |
— e-* |
( * > 0 ) , |
(3) |
|
|
|
о |
|
|
|
а для сдвинутого |
показательного распределения |
|
||||
|
F(t)= |
1 — е - М ' - ' . ) , |
|
|
(За) |
|
Интегральная функция распределения интервалов по |
||||||
показательному |
закону приведена |
на рис. 3 для случая, |
||||
когда Я = 3 |
поезда в час, а по сдвинутому |
показательно |
||||
му закону |
при сдвижке на * о = 0 , 1 |
ч — на |
рис. 4. |
|
12
Рис. 1. Непрерыв 5!If ное распределение a: l интервалов меж I
ду поступлением поездов в парк прибытия по по казательному за кону
Рис. 2. Плотность 6 вероятности сдви §
нутого показа тельного распре деления интерва лов поступления поездов в . парк
прибытия
Рис. 3. Интеграль ная функция рас пределения интер валов по показа тельному закону
№
Рис. 4. Интеграль ная функция рас пределения интер валов по сдвину
тому показатель ОД
ному закону
-
t
0,1 0,2 Щ 0,6 0,8 1,0 1,г Величина интердалод поступления поездов t,v
0,1 0,2
— -=
I i
J
г
т
0,1 0,2 |
Ofi |
0,6 |
Ofi |
1,0 |
t,4- |
Между функцией плотности распределения и инте гральной функцией распределения оуществует тесная связь. Производная от функции распределения для непре рывных случайных величин равна плотности распреде ления. В связи с этим ордината точки t кривой плотно сти распределения равна тангенсу угла наклона каса тельной в точке t к интегральной кривой. Чем ниже точ ка на кривой плотности распределения, тем меньше тан
генс наклона |
касательной к интегральной кривой, тем |
медленнее она |
растет. |
Показательное распределение является наиболее |
лростым, зависящим только -от одного параметра К. Для потоков с показательным распределением интервалов получены аналитические зависимости, характеризующие показатели функционирования систем массового обслу живания. Кроме того, при помощи показательных рас пределений можно сформировать потоки с различной степенью неравномерности, т. е. с различными коэффи циентами вариации интервалов.
Р а с п р е д е л е н и е Э р л а н г а
Каждый интервал в потоке Эрланга состоит из К подын тервалов, распределенных по юказательному закону с ин тенсивностью X/. На рис. 5 по оси П нанесен поток собы тий с показательным распределением интервалов (пуассоновский поток) и интенсивностью X/. На оси Э показан по ток с эрланговским распределением порядка К = 2 и интен сивностью X, у которого каждый интервал состоит из двух подынтервалов. Математическое ожидание интервалов в по токе Эрланга равно — , а математическое ожидание подын-
х
тервалов равно .
11 |
|
|
X т |
. |
\ II |
|
|
|
/ |
|
|
1 |
1 |
1 |
^ ' |
J |
JL |
|
|
|
|
|
|
Рис. 5. Схема образования потока Эрланга: |
|||||
Л — пауссоновский |
поток; |
|
Э — эрланговский поток второго по |
||
|
|
|
рядка |
|
|
14
Следовательно, среднее значение интервалов в эрланговском потоке
т |
1 |
1 . 1 |
2 |
. |
X' |
,X'
или в оощем виде л = — .
Для потока Эрланга порядка К'- плотность распределения интервалов
т |
= |
(Л" — 1)! |
(4) |
у w |
|
v 7 |
дисперсия интервалов
/с
среднее квадратическое отклонение
а - № •
коэффициент вариации
t 1к
Формула (4) выражает плотность распределения интер валов эрланговского потока при фиксации интенсивности X' исходного пуассоновского потока и поэтому с увеличением порядка К увеличивается среднее значение интервалов
{t = —-^, т. е. поток как бы прореживается, делается бо лее редким. Для устранения этого применяется нормиро ванное распределение Эрланга, которое можно получить, подставив в зависимость (4) интенсивность эрланговского
потока |
л = |
вместо |
интенсивности |
исходного |
потока |
|
У = |
\к, |
|
W W - * |
- |
щ |
|
|
|
|
( 5 ) |
|||
|
|
•> w |
{к — 1) |
|
|
|
В этом случае с увеличением параметра К поток сохраняет свою интенсивность неизменной, а будет ме няться лишь величина рассеивания интервалов в потоке.
15
Числовые характеристики нормированного распреде ления Эрланга примут вид:
среднее значение интервалов
дисперсия
1 .
2 •среднее квадратическое отклонение
1
° IVк '
коэффициент вариации
мода, т. е. интервал, имеющий наибольшее значение плот ности распределения,
|
ЛГ„< |
к — 1 |
|
|
|
|
|
|
К\ |
|
|
При К=1 эрланговское |
распределение |
принимает |
|||
вид показательного распределения |
с коэффициентом |
||||
вариации, |
равным единице, |
а при /(-»-оо поток |
Эрланга |
||
сходится |
к регулярному |
потоку с одинаковыми интерва |
|||
лами и коэффициентом |
вариации, |
равным нулю. Зави |
симость степени неравномерности эрланговских потоков
от порядка |
К приведена в следующем ряду: |
|
|
||||||||
А"... |
1; |
2; |
3; |
4; |
5; |
6; |
7; |
8; |
9; |
10; |
.... «>• |
V . . . |
1; |
0,71; |
0,58; |
0,5; |
0,45; |
0,41; |
0,38; |
0,35; |
0,33; |
0,31; |
0 |
На рис. 6 приведены кривые плотности распределе ния нормированных потоков Эрланга при различных значениях порядка К и интенсивности %=4, а значит среднем значении интервалов t=0,25. Для примера по следовательность расчета при /С=2 приведена в табл. 1. Как видно из рис. 6, с увеличением порядка К уменьша ется рассеивание интервалов относительно среднего зна-
.чения, уменьшается положительная правосторонняя (хвост справа) асимметрия, а мода перемещается от ну-
16
Рис. 6. Зависи мость плотности распределения от параметра К при нормированном за
коне Эрланга
ля до среднего значения. Так, при К=\, М0=0 — пока
зательное |
распределение: |
||||
К |
= |
2 |
ЛГ0 |
= |
0,125 |
К = |
5 |
М0 |
= |
0,2 |
|
К = 25 |
М0 |
= |
0,24 |
||
ЛГ -»- оо |
Ж 0 |
= |
^ = 0,25 — регулярный поток. |
При большом значении К эрланговское распределе ние приближается к нормальному (Гаусса) и может быть заменено последним (при К^9).
Сумма всех элементов вероятностей (вероятностей попадания случайной величины в данный разряд), выра жающая площадь под кривой плотности распределения, должна быть равна единице (итог последней графы табл. 1).
17
|
|
|
|
|
|
Т а б л и ц а 1 |
Рас чет значений вероятности |
тех или иных |
интервалов |
||||
поступления требований в эрланговском |
потоке |
|||||
|
|
при А = 4, |
К = 2 |
|
|
|
Границы |
Середина интерва |
|
/(') |
|
Элемент вероят |
|
интервалов, ч |
ла |
группировки |
по формуле |
(5) |
ности Ц1) At |
|
0 —0,1 |
|
0,05 |
|
2,14 |
|
2,145-0,1=0,214 |
0,1—0,2 |
|
0,15 |
|
2,80 |
|
0,280 |
0,2—0,3 |
|
0,25 |
|
2,16 |
|
0,216 |
0,3—0,4 |
|
0,35 |
|
1,30 |
|
0,130 |
0,4-0,5 |
|
0,45 |
|
0,78 |
|
0,078 |
0,5-0,6 |
|
0,55 |
|
0,40 |
|
0,040 |
0,6—0,7 |
|
0,65 |
|
0,22 |
|
0,022 |
0,7—0,8 |
|
0,75 |
|
0,11 |
|
0,011 |
0,8—0,9 |
|
0,85 |
|
0,05 |
|
0,005 |
0,9—1,0 |
|
0,95 |
|
0,03 |
|
0,003 |
1,0—1,1 |
|
1,05 |
|
0,01 |
|
0,001 |
1,1—1,2 |
|
1,15 |
|
0,00 |
|
0,000 |
Эрланговским |
распределением с |
соответствующим |
значением порядка К можно аппроксимировать потоки с различной степенью неравномерности с коэффициентом вариации интервалов от 0 до 1. Однако при переходе от К—\ до К=2 имеется значительный скачок значений коэффициента 'вариации от 1 до 0,71. Для потоков, име ющих коэффициент вариации в границах от 0,71 до I , плотность .распределения интервалов может быть пред ставлена в виде о б о б щ е н н о г о з а к о н а Э р л а н г а второго порядка, который является композицией законов с показательным распределением, но с различной ин тенсивностью. В этом случае интервалы обобщенного по тока состоят из подынтервалов, распределенных по пока
зательному закону и имеющих разные |
интенсивности. |
Плотность распределения интервалов для обобщенно- |
|
то закона Эрланга второго порядка |
|
/ ( 0 |
(6) |
Aj — A i |
|
Здесь интенсивность составляющих потоков определяет ся в зависимости от коэффициента вариации интервалов по формуле [20]:
, |
X ± X V l —2(1 —У») |
, f i , |
18
Если |
коэффициент |
вариации |
интервалов |
|
прибытия |
по |
||||||||||
ездов |
на |
станцию |
|
составляет |
V = 0,8 |
при |
интенсивности |
|||||||||
h=4 |
поезда в |
час, |
то для |
аппроксимации |
"этого |
потока |
||||||||||
•применим |
обобщенный |
закон |
Эрланга |
порядка |
|
К=2. |
||||||||||
По формуле (6а) найдем интенсивности |
|
составляющих |
||||||||||||||
исходных |
потоков: |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
4 — 4У1 — 2(1 — 0,82 ) |
= |
5,24; |
|
|
|
|||||||
|
|
|
|
|
|
|
1 — |
0,82 |
|
|
|
|
|
|
|
|
|
|
|
|
|
4 + 4JT— |
2(1 — 0,8; ) |
= |
17. |
|
|
|
|||||
|
|
|
|
|
|
|
1 — |
0,8' |
|
|
|
|
|
|
|
|
Среднее |
значение |
|
первого |
подынтервала |
^=1/5,24 = |
|||||||||||
= 0,191 |
ч, а второго |
1/17=0,059 |
ч, и средняя |
величи |
||||||||||||
на интервала в потоке Эрланга составит |
0,25 |
ч. |
|
|
||||||||||||
Плотность распределения интервалов для обобщен |
||||||||||||||||
ного потока Эрланга определяем по формуле |
(6), |
под |
||||||||||||||
ставив полученные значения hi и h%, |
|
|
|
|
|
|
|
|||||||||
/ ( Q - |
|
5 ' 2 4 - |
1 7 |
(е~5'™ |
е ~ ш ) = 7,56 (е -5,141 |
|
|
|
||||||||
|
|
17 — |
5,24 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Показанная на рис. 6 пунктиром кривая распределе |
||||||||||||||||
ния по обобщенному закону для h—4 |
и |
V = 0 , 8 |
заняла |
|||||||||||||
промежуточное (положение |
между |
кривыми |
с К= |
1 и |
||||||||||||
/С=2, /имеющими |
коэффициент вариации |
соответствен |
||||||||||||||
но 1 и 0,71. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
С приближением коэффициента вариации к единице |
||||||||||||||||
увеличивается |
различие |
параметров |
исходных |
|
потоков |
и средних значений интервалов. Один подынтервал стре мится к среднему интервалу t, а второй к нулю. С приб лижением коэффициента вариации к 0,71 параметры ис
ходных |
потоков |
и средних |
значений |
подынтервалов |
вы |
||
равниваются. Это хорошо видно из данных, |
приведенных |
||||||
в табл. 1а. |
|
|
|
|
|
|
|
|
|
|
|
1 |
*. |
Т а б л и ц а |
1а |
V |
\ |
X, |
|
|
|
||
|
|
|
|
|
|||
0,71 |
6 |
6 |
10 |
|
10 |
20 |
|
0,75 |
9,2 |
4,45 |
6,5 |
|
13,5 |
20 |
|
0,8 |
12,7 |
3,93 |
4,7 |
|
15,3 |
20 |
|
0,85 |
18,0 |
3,6 |
3,3 |
|
16,7 |
20 |
|
0,9 |
28,22 |
3,36 |
2,1 |
|
17,9 |
. 20 |
|
0,95 |
56,8 |
3,2 |
1,1 |
|
18,9 |
20 |
|
Г
!