Файл: Золотарь, И. А. Экономико-математические методы в дорожном строительстве.pdf

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

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

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

Добавлен: 19.10.2024

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

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

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

При решении конкретных задач методами ТМО необходима про­ верка гипотезы о том, что поток заявок на обслуживание является пуассоновским. Методика такой проверки была рассмотрена в гл. III.

§16. Метод получения расчетных зависимостей теории массового обслуживания

Для многих практически важных случаев соответствующие рас­ четные зависимости получены и приводятся в литературе по ТМО. Однако практика ставит новые задачи, для решения которых тре­ буется владение методом получения выражения для рп (х). Рас­ смотрим основы этого метода на простейших видах систем МО.

Функция рп (т) находится из решения системы дифференциаль­ ных уравнений вероятностей состояния системы. Обычно диффе­ ренциальные уравнения для процессов, развивающихся во вре­ мени, характеризуют закономерности изучаемого процесса при­ менительно к бесконечно малому промежутку времени dx. В после­ дующем на основе их решения удается выявить общие свойства процесса. Это полностью относится к системам МО. Для написа­ ния дифференциальных уравнений, описывающих изменения состо­ яния системы МО за бесконечно малый промежуток времени dx, необходимо предварительно установить два основных соотношения, а именно: для вероятности поступления в систему одной заявки на обслуживание за время dx, т. е. для величины p x(dx); для вероятно­ сти убыли из системы (вследствие завершения обслуживания) за время dx одного требования, т. е. для величины p - X(dx).

Соотношение для

р х(dx) легко получить

из общего

выражен

для пуассоновского

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

рп(т ) =

\%пе~и

а именно

------------,

 

 

 

п\

 

(Xrft)ie

xd-z

p x(dx) =

Mx.

(VII.9)

p x(dx) = -------------------, что дает

Аналогично, вероятность окончания обслуживания одного требо­ вания за время dx выразится

1 (dx) = pdx, (VII. 10)

где р — пропускная способность (номинальная производительность

прибора обслуживания в единицу времени).

основе

теоремы

Зависимости (VII.9)

и (VII. 10)

позволяют на

сложения вероятностей

записать

сразу еще два

важных

соотно­

шения:

 

 

 

 

1) 1 — Ых — вероятность того, что за время dx в систему МО не поступит требование;

2) 1 — p it — вероятность того, что за время dx систему не поки­ нет обслуженное требование.

114


Если в системе на момент времени т нет требований, то на мо­ мент времени x + dx возможны два состояния системы МО: •

1)Ро-^Ро, т. е. в систему не поступит требование на обслужива­

ние;

2)ро-э-рь т. е. в систему за время dx поступит одно требование. На основе (VII.9) и (VII. 10), используя теоремы сложения и

умножения вероятностей (см. гл. III), получим следующие соотно­ шения для вероятности указанных переходов системы из одного со­ стояния в другое:

Ро—» Ро= 1

(VII. 11)

Ро- ” Pi~^dx.

(VII.12)

Если на момент времени т в системе было п требований, то воз­ можны ее переходы в следующие три состояния:

Р п * Рл> Р п * Р п + 1’ Р п * Р п — 1-

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

Рп~* р„= (1 Ых){ \ ■—y.dx)A-'kdx-\xdx.

Первое произведение в правой части дает вероятность того, что за время d\t в систему не поступит и из системы не убудет требо­ вание. Второе — вероятность того, что в систему за время dx посту­ пит одно требование и систему покинет одно обслуженное требо­ вание. Других вариантов сохранения в системе п требований, оче­ видно, быть не может.

Опуская бесконечно малые величины второго порядка, получим:

р„ — Рл^ 1 — (X-fix)rft. (VII. 13)

На основе подобных же рассуждений можно получить:

Рп—-, Рл + 1 =

(1 — p d t ) M t

или р„—» p„4.! = >-filt;

(VII. 14)

рл—* p„_i(l —M%)y-dx

или р„—* Рл-^tAflfT.

(VII. 15)

Составим теперь

с помощью

соотношений (VII.11) — (VII. 15)

матрицу вероятностей перехода системы МО за время dx в возмож­ ные состояния:

p{x-\-dx)

Состояние системы на момент t + d z

 

0

1

2

3

4

5

0

1—~kdz

\d z

0

0

0

0

1

\xdz

1—(\+ \x jd z

~kdz

0

0

0

2

0

txdz

1— (X-J-jj.)t/x

Adz

0

0

3

0

0

fulz

1— (X+p.)rfT:

Idz

0

115


Рассматривая матрицу, можно заметить, что, кроме первой, все строки матрицы в записи идентичны. То же можно сказать и о столбцах. Тогда можно записать следующую систему дифферен­ циальных уравнений:

pQ{x-{-dx) =

p0{x)(\—'kix)-\-pl {x)^dX

для п = 0;

 

 

рп-\-dx) =

р п_ х(т) Ых -4- рп(т) [ 1 (X-ф р.)dx\-f pn+1 (т) p-dx для «•> 1.

Эта система может содержать п + 1 уравнение и в принципе мо­

жет быть бесконечной.

 

 

 

Производя упрощения, получим:

 

 

 

p Q{xJr dx) — р 0{х)=

p Q{x)\dx-\-pl (x)]i-dx

для n =

0;

Pn (x -f- dx) -

Рп ( Г ) =

Pn-i М Xdx — рп( t ) +

р.) dx - f р п+1 М

V-dx

 

р0 (т + dx) — /7р (т)

 

 

для

1;

Ра {х)'>-Аг Pi СО V-

 

для п = 0;

dx

 

 

 

 

 

 

 

рп (т 4- dx) рп (т)

Рп--1 (t) л — рп(t) (X +

[*)+ рп+j (Т) ц

для п ^>\.

 

 

dx

Таким образом,

*Ро dx

d p n

dx

P n - l (x j X - P n ( * ) ( X + Iх ) + P n + l ( t ) I*

 

для n = 0;

для 1.

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

Тогда, принимая dpn - = 0 и

dPn --—О, получим следующую си-

d (т)

d (т)

 

 

стему обыкновенных уравнений:

 

 

р<^—Р№

 

для я = 0 ;

+

+

 

Для п > {-

Примем во втором уравнении п = 1. Тогда

 

Pok^Ptf

 

для п = 0;

P i^ + V-)= P ^ + P ^

для п — 1.

Разделим теперь оба уравнения на ц. Получим

 

 

1-

 

для п = 0;

Ро— = P l

 

 

V-

 

 

Pi ( --------b

^ ) ----Ро---------V Р2

для п = 1.

\ Р

I

н-

 

116


Так как — = <р, то будем иметь: н-

J Ро^= Рх или Р! = Р0%

Внося первое уравнение во второе, получим

 

 

Рч= Р ^

и аналогично

 

 

Рп= РаТ-

 

(VII. 16)

Так как сумма вероятностей всех состояний системы МО равна

единице, то можно записать:

 

 

 

 

оо

 

 

 

 

2

Л ,= 1-

 

(VII. 17)

 

/1=0

 

 

Внося (VII.16) в (VI 1.17), получим

 

 

ОО

 

далее

со

(VII. 18)

22 /7оФ"= 1 и

p QV фл = 1.

л = 0

 

 

л = 0

 

оо

представляет

собой сумму

членов геомет-

Известно, что V

/т=0

 

 

 

 

рической прогрессии,

равную

----- 1------. Тогда

 

 

 

1 — О/

 

 

Ро ——Ц — = 1 , ИЛИ p 0= l —ty.

(VII. 19)

Соотношение (VII. 19) можно было бы записать и из чисто логи­ ческих, или, как часто говорят, эвристических соображений. В са­ мом деле, как отмечалось в главе выше, величина а]з представляет собой математическое ожидание требований, находящихся в систе­ ме, или, что то же самое, показывает, какую часть своего времени прибор занят. Тогда 1—^ соответствует той части времени, когда прибор свободен, что и дает вероятность ро.

Внося (VII.19) в (VII.16), получим

Рп= Г { 1 —ф).

(VII.20)

Определим теперь в соответствии с зависимостью (VII.2) мате­ матическое ожидание числа требований в системе in:

__

ОО

оо

00

m =

v

прп = ^

whn{\ — <!>)=( 1 —<|>)2 п-Т-

 

/2=0

л = 0

/1=0

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

на ------— . Тогда (1—Ф)2

т = -----*-----

.

(VII.21)

J —

<]/

 

117


На основе (VI 1.5) получим

57

— I

Ф

I

V

Ф2

(VII.22)

Y =

т — ф = —

— ф, т. е. Y =

1 ~ |

 

 

1 - ф

 

 

 

Далее tf —

1

ф2

 

 

X

то

 

-----;----- .

Так как Ф— — ,

 

 

X

(1 — ф)Х

 

 

 

 

 

 

 

 

1 — ф

 

 

(VII.23)

 

 

 

 

 

 

 

Наконец,

ts = t f -|-— ,

в результате

чего

получаем

 

 

 

V-

 

 

 

 

 

 

 

t . = -

1

 

 

(VII.24)

 

 

\ ! —Ф

 

 

 

 

 

К-

 

 

 

Если время обслуживания постоянно, т. е. если каждое обслу­

живание продолжается точно в течение времени — , то по теории

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

Ф

(VII.25)

2 •(1 — ф)

_ Кроме того, среднее число требований в системе т и в очереди F определяются в этом случае по формулам:

т — ф-}

Ф2

(VII.26)

 

2 (1 - Ф)

Y

ф2

(VII.27)

2-(1 — ф)

 

 

§17. Разомкнутые системы массового обслуживания с одним

инесколькими приборами

Рассмотрим одну задачу, решение которой можно будет полу­ чить с помощью зависимостей (VII. 19) и (VII.22).

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

В качестве критерия оптимальности примем минимум суммар­ ных потерь (в стоимостном выражении) от простоя транспортных средств и средств погрузки. Будем рассматривать разомкнутую СМО с ожиданием, с одним прибором обслуживания, без приори­ тета требований, для которой и справедливы формулы (V.19) и

118