Файл: Золотарь, И. А. Экономико-математические методы в дорожном строительстве.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