Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Шабалин Н.Н. Оптимизация процесса переработки вагонов на станциях

.pdf
Скачиваний:
9
Добавлен:
23.10.2023
Размер:
6.13 Mб
Скачать

ное (экспоненциальное), эрланговские, гамма-распреде­ ление, гнперэ'кспоненциалы-юе, вейбулловское и нор­ мальное. Часто, однако, для .решения практических за­ дач достаточно иметь только числовые характеристики интервалов в потоках, такие как среднее значение (ма­ тематическое ожидание), среднее квадрэтическое откло­ нение, коэффициент вариации.

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

ется потребность в резервах мощности для

смягчения

влияния неравномерности.

 

 

П о к а з а т е л ь н о е

(экспоненциальное)

распреде­

ление определяет закономерности появления тех или иных интервалов в простейшем потоке однородных со­ бытий. Простейший поток, как известно, обладает свой­ ствами стационарности, ординарности и отсутствия по­ следействия.

Под стационарностью подразумевается свойство по­

тока, при котором вероятность поступления

определен­

ного числа

требований

(поездов,

вагонов)

в течение

некоторого

отрезка

времени

зависит

только

от

продолжительности этого отрезка времени, числа требо­ ваний и не зависит от положения этого отрезка на оси времени.

Ординарность потока означает практическую невоз­

можность появления двух

и

более требований в

тече­

ние малого отрезка времени 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

 

Г

!

Соседние файлы в папке книги из ГПНТБ