Опубликован: 20.04.2011 | Уровень: для всех | Доступ: платный
Лекция 9:

Полнодоступные системы с потерями

< Лекция 8 || Лекция 9: 123 || Лекция 10 >

Распределение Энгсета

Единственная разница по сравнению с материалами секции 8.2 - то, что число источников S теперь больше или равно числу пучков каналов (каналам), S \ge n . Поэтому, попытки вызова могут быть потеряны.

 Диаграмма переходов состояний для случая распределения Энгсета с S>n, где S - число источников и n - число каналов.

Рис. 8.4. Диаграмма переходов состояний для случая распределения Энгсета с S>n, где S - число источников и n - число каналов.
Вероятности состояния

Уравнения сечения идентичны (8.1), но они существуют только для 0 \le i \le n (рис.8.4). Уравнение нормализации (8.2):

1=p(0)*\left \{ 1+ {S\choose 1}*{\gamma \choose \mu}+\dots +{S\choose n}*{\gamma \choose \mu}^n \right \}.

Из него мы получаем p(0) и, подставляя \beta = \frac{\gamma}{\mu} , получаем вероятности состояния, которые равны:

p(i)=\frac{{S\choose i}*\beta^i}{\sum_{j=0}^n {S\choose j}* \beta^j}. ( 8.24)

Тем же самым способом, который мы применяли выше, используя (8.8), мы можем переписать это выражение в форме, которая является аналогом (8.4):

p(i)=\frac{{S\choose i}*a^i*(1-a)^{S-i}}{\sum_{j=0}^n {S\choose j}*a^j*(1-a)^{S-j}}, 0 \le i \le n, ( 8.25)
Характеристики нагрузки в модели Энгсета

Распределение Энгсета сопровождается более сложными вычислениями, чем Эрланговская система с потерями. Основная проблема в том, что нужно понять, как найти критерии качества работы непосредственно из вероятностей состояния, используя определения. Энгсетовская система характеризуется следующими параметрами: \beta = \frac{\gamma}{\mu} - предложенная нагрузка на свободный источник, S - число источников и n - число каналов.

Потери по времени Е, по определению, пропорциональны времени блокирования системы для новых попыток вызова, то есть р(n) (8.24):

E_{n,S}(\beta}=p(n)=\frac{{S\choose n}* \beta^n}{\sum_{j=0}^n {S\choose j}* \beta^j}, S \ge n. ( 8.26)

Потери по вызовам В, по определению, пропорциональны числу потерянных попыток вызова. Блокируются только попытки вызова, достигающие системы в состоянии п. В течение одной единицы времени мы получаем следующее отношение между числом блокированных попыток вызова и общим количеством попыток вызова:

B_{n,S}(\beta)=\frac{p(n)*(S-n) \gamma}{\sum_{j=0}^n p(i)*(S-j) \gamma}\\
=\frac{{S\choose n}* \beta^n *(S-n) \gamma}{\sum_{j=0}^n {S\choose j}* \beta^j *(S-j) \gamma}.

Используя

{S\choose i}*\frac{S-i}{S}={S-1\choose i}

мы имеем:

B_{n,S}(\beta)=\frac{{S-1\choose n}*\beta^n}{\sum_{j=0}^n{S-1\choose j}*\beta^j},\\
B_{n,S}(\beta)=E_{n,S-1}(\beta), S \ge n. ( 8.27)

Этот результат можно интерпретировать следующим образом. Вероятность того, что попытка вызова от случайного источника (абонента) будет отклонена, равна вероятности того, что остальные (S-1) источники заняли все п каналов. Это называется теоремой поступления, и можно показать, что она справедлива для систем с явными потерями и для систем с ожиданием и ограниченным числом источников. Результат основан на вычислении произведения среди источников и свертывании источников. Поскольку Е увеличивается, когда увеличивается S, мы имеем:

B_{n,S}(\beta)=E_{n,S-1}(\beta) < E_{n,S}(\beta).

Теорема 8.1. Теорема поступления: Для всех систем с ограниченным числом источников, случайный источник при поступление заявки будет наблюдать состояние системы так, как будто сам источник не принадлежит этой системе.

Свойство PASTA включено в этот случай, потому что бесконечное число источников минус один есть бесконечное число.

Обслуженная нагрузка: применяя уравнение сечения между состоянием [i - 1] и состоянием [i], мы получаем:

Y=\sum_{i=1}^n i*p(i) ( 8.28)
=\sum_{i=1}^n \frac{\gamma}{\mu}*(S-i+1)*p(i-1)\\
=\sum_{i=0}^{n-1} \beta * (S-i)*p(i) ( 8.29)
=\sum_{i=0}^n \beta *(S-i)*p(i)-\beta *(S-n)*p(n),\\
Y=\beta * (S-Y)- \beta *(S-n)*E, ( 8.30)

Поскольку Е= Е_п, S (\beta) =р(п) . Последнее уравнение решается относительно Y:

Y=\frac{\beta}{1+\beta}*\{S-(S-n)*E\}. ( 8.31)

Потери по нагрузке С = C_{n,S} (А) . Это самая важная характеристика потерь. Предложенная нагрузка дается в (8.20), и мы получаем:

C=\frac{A-Y}{A}\\
=\frac{\frac{S \beta}{1 +\beta}-\frac{\beta}{1 + \beta}*\{S-(S-n)*E\}}{\frac{S \beta}{1+ \beta}},\\
C=\frac{S-n}{S}*E. ( 8.32)

Мы можем также найти обслуженную нагрузку, если знаем потери по вызовам В. Число попыток принятия вызовов от источника, который находится в свободном состоянии, в среднем \frac{1}{\gamma} в единицу времени прежде, чем источник сгенерирует одну попытку вызова - 1(1-В), и каждый принятый вызов имеет среднюю продолжительность 1/\mu - Таким образом, обслуженная нагрузка на один источник есть соотношение времени: когда источник является занятым, она будет:

y=\frac{(1-B)/\mu}{1/\mu +(1-B)/\mu}.

Полная обслуженная нагрузка будет:

Y=S*y=S*\frac{\beta (1-B)}{1+\beta (1-B)}. ( 8.33)

Приравнивая эти два выражения для обслуженной нагрузки (8.31) и (8.33), мы получаем следующее отношение между Е и В:

E=\frac{S}{S-n}*\frac{B}{1+ \beta (1-B)}. ( 8.34)

Число попыток вызова в единицу времени:

\Lambda = \sum_{i=0}^n p(i)*(S-i) \gamma\\
\Lambda = (S-Y)*\gamma, ( 8.35)

где Y - обслуженная нагрузка (8.28). Таким образом, из уравнения (S-Y) среднее число свободных источников является очевидным. Исторически, полная предложенная нагрузка была определена как \Lambda / \mu. Это, однако, вводит в заблуждение, потому что мы не можем принять, что каждая повторная попытка вызова имеет среднее время пребывания в системе, равное 1/ \mu -Также это определение создает большое неудобство, потому что предложенная нагрузка по этому определению зависит от состояния системы (числа занятых каналов). Также возможно, что немногие из доступных обслуживающих много попыток вызова устройств блокированы, а свободные источники с более высоким средним временем поступления вызовов генерируют больше попыток вызова в единицу времени.

Потерянная нагрузка:

A_l=A*C\\
=S\frac{\beta}{1+ \beta}*\frac{S-n}{S} E\\
= \frac{(S-n) \beta}{1+ \beta}*E. ( 8.36)

Продолжительность состояния i Оно является экспоненциально распределенным с интенсивностью:

\gamma (i)=(S-i) * \gamma +i* \mu, \qquad 0 \le i \le n,\\
\gamma (n) = n \mu, \qquad i=n ( 8.37)

Функция увеличения:

F_{n,S}(A)=Y_{n+1}-Y_n ( 8.38)

Выше мы определили вероятности состояния p(i), согласно предположению о статистическом равновесии, как часть времени, которое система находится в состоянии i, то есть как математическое ожидание времени. Мы можем также исследовать, как выглядит система, когда выполняется равновесие между поступлением и возвратом вызова отправляющим источником (пользователем) (математическое ожидание вызова). Если мы рассматриваем одну единицу времени то, в среднем, в системе будет (S-i)\gamma p(i) источников в состоянии [i] перед моментом поступления вызова, и если вызов принят, то он переведёт систему в состояние [i+1].

Источники, которые наблюдают систему в состоянии п, блокированы или остаются свободными. Поэтому источники, от которых поступает вызов прибытия, наблюдают систему в состоянии [i] с вероятностью:

\pi_{n,S, \beta}(i)=\frac{(S-i) \gamma * p(i)}{\sum_{j=0}^n (S-j) \gamma * p(j)}, i=0,1, \dots , n. ( 8.39)

Используя метод аналогового дифференцирования выражения (8.27), мы можем показать, что в соответствии с теоремой поступления (Теорема 8.1) мы имеем:

\pi_{n,S, \beta}(i)=p_{n,S-1, \beta}(i-1), i=0,1, \dots , n. ( 8.40)

Когда источник оставляет систему, система находится в состоянии [i-1] с вероятностью:

\Psi_{n,S, \beta}(i-1)=\frac{i \mu *p(i)}{\sum_{i=1}^n j \mu * p(j)}, i=1,2, \dots, n. ( 8.41)

Применяя уравнения сечения, мы немедленно убеждаемся, что это соотношение идентично (8.39), если включить в рассмотрение блокированных абонентов. Таким образом, источники покидают систему в том же самом состоянии, в котором они в нее прибывают. Процесс будет, обратим и нечувствителен к распределению времени обслуживания. Если бы мы непрерывно наблюдали состояние системы, то не могли бы определить, выполняется ли время прямое действие (вызов) или обратный (выход из системы).

Отношения между Е, В и С

Из (8.34) мы получаем следующее отношение между E=E_{n,S}(\beta) и B=B_{n,S}(\beta)=E_{n,S-1}(\beta)

E=\frac{S}{S-n}*\frac{B}{1+ \beta (1-B)} \mbox{ или } \frac 1E=\frac{S-n}{S}\left \{(1+\beta)*\frac 1B - \beta \right \}, ( 8.42)
B=\frac{(S-n)*E*(1+ \beta)}{S+(S-n)*E* \beta} \mbox{ или } \frac 1B=\frac{1}{1+ \beta} \left \{ \frac{S}{S-n}* \frac 1E + \beta \right \} ( 8.43)

Выражения с правой стороны линейны по отношению к вероятностям блокировки. В (8.32) мы получили следующее простое отношение между С и Е:

C=\frac{S-n}{S}*E, ( 8.44)
E=\frac{S}{S-n}*C. ( 8.45)

Если мы в (8.44) подставим Е, выраженное через (8.42), то мы получаем С, выраженное через В:

C=\frac{B}{1 + \beta *(1-B)}, ( 8.46)
B=\frac{(1+ \beta) C}{1+ \beta C}. ( 8.47)

Это отношение между В и С является общим и может также быть получено следующим образом. Обслуженная нагрузка Y соответствует (У* \mu ) принятых попыток вызова в единицу времени. Среднее число свободных источников - (S-Y), так что среднее число попыток вызова в единицу времени - (S-Y) \gamma (8.35).

Потери по вызовам тогда равны:

B=\frac{(S-Y) \gamma - Y * \mu}{(S-Y) \gamma}\\
=\frac{(S-Y) \beta - Y}{(S-Y) \beta}.

По определению Y= (1 - C), и из (8.20) имеем S=А( 1 + \beta)/ \beta. Подставляя это выражение, мы имеем:

B=\frac{A(1+\beta)-A(1-C) \beta - A(1-C)}{A(1+\bets) - A(1-C) \beta}\\
B=\frac{(1+ \beta)C}{1+ \beta C}

Из последнего уравнения видно, что для малых значений потерь по вызовам В потери по нагрузке Z в несколько раз больше, чем потери по вызовам:

C \approx \frac{B}{1+ \beta}=Z*B ( 8.48)

Расчеты по формуле Энгсета

Если мы пробуем вычислить числовые значения непосредственно по формуле Энгсета из (8.26) (потери по времени Е ), то возникают проблемы расчета для больших значений S и п. Ниже мы получим различные рекурсивные формулы для Е и для его обратной величины I=1/E. Когда потери по времени Е известны, можно просто получить потери по вызовам В и потери по нагрузке С, используя формулы (8.43) и (8.44). В числовой форме также просто найти любой из этих четырех параметров S, \beta, n, Е, когда мы знаем три из них. Математически мы можем предположить, что п и S не являются целыми числами.

Рекурсивная формула для п

Из общей, рекурсивной формулы (7.27) для я, используя \lambda_х = (S-x) \gamma и \beta = \frac{\gamma}{\mu}, мы получаем

E_{x,S}(\beta)=\frac{\frac{\gamma_{x-1}}{x \mu}*E_{x-1,S}(\beta)}{1+\frac{\gamma_{x-1}}{x \mu}*E_{x-1,S}(\beta)},\\
E_{x,S}(\beta)=\frac{(S-x+1) \beta * E_{x-1, S}(\beta)}{x+(S-x+1) \beta * E_{x-1, S}(\beta)}, E_{0,S}(\beta)=1. ( 8.50)

Используя обратную формулу потерь по времени I_{n,S}(\beta)=\frac{1}{E_{n,S}(\beta)}, мы находим рекурсивную формулу:

I_{x,S}(\beta)=1+\frac{x}{(S-x+1) \beta} I_{x-1, S}(\beta), I_{0,S}(\beta)=1. ( 8.51)

Число итераций равно п. Обе формулы (8.50) и (8.51) аналитически точны и представляют собой устойчивые при приближенных вычислениях и точные рекурсии для растущих значений х. Однако для уменьшающихся значений х числовые ошибки накапливаются и рекурсии недостоверны.

Рекурсивная формула по S

Обозначим нормализованные вероятности состояния системы с п каналами и S-1 источниками p_{n, S-1}(i). Мы получаем вероятности состояния системы с S источниками и п каналами, сочетая эти вероятности состояния с вероятностями состояния одиночного источника, которые обозначаются \{р_{1,1}(0) = 1 - а, р_{1,1} (1) = а_g \}. Мы тогда получим состояние от нуля до п + 1, ограничиваем пространство состояний до n и нормализуем вероятности состояния (см. Пример 3.2.1) (принимая, что р(х) = 0, когда x < 0 ):

q_{n,S}(i)=(1-a)*p_{n,S-1}(i)+a*p_{n,S-1}(i-1), i=0,1, \dots, n. ( 8.52)

Полученные вероятности состояния q_{n,S}(i) не нормализованы, потому что мы ограничиваем состояния числом [n] и исключаем последние элементы, начиная с состояния [n+1], q_{n,S}(n+1)=a*р_{m,S-1}(n) . Нормализованная вероятность состояния р_{n,S}(i) для системы с S источниками и п каналами, таким образом, получается из нормализованной вероятности состояния p_{n,S-1}(i) для системы с S-1 источником:

p_{n,S}(i)=\frac{q_{n,S}(i)}{1-a*p_{n,S-1}(n)}, i=0,1, \dots, n. ( 8.53)

Потери по времени Е_{n,S}(\beta) для системы с S источниками могут быть выражены потерями по времени Е_{n,S-1}(\beta) для системы с S-1 источниками. Подставляя (8.52) в (8.53), получаем:

E_{n.S}(\beta)=p_{n,S}(n)\\
=\frac{(1-a)*p_{n,S-1}(n)+a*p_{n,S-1}(n-1)}{1-a*p_{n,S-1}(n)}\\
=\frac{(1-a)*E_{n,S-1}(\beta)+a*\frac{n \mu}{(S-n) \gamma} E_{n,S-1}(\beta)}{1-a*E_{n,S-1}(\beta)},

где мы использовали уравнение равновесия между состоянием [n-1, S - 1] и состоянием [n, S-1]. Заменяя а на(8.8), мы получаем:

E_{n,S}(\beta)=\frac{E_{n,S-1}(\beta)+\frac{n}{S-n} E_{n,S-1}(\beta)}{1+\beta-\beta E_{n,S-1}(\beta)}.

Таким образом, получаем следующую рекурсивную формулу:

E_{n,S}(\beta)=\frac{S}{S-n}*\frac{E_{n,S-1}(\beta)}{1+\beta\{1-E_{n,S-1}(\beta)\}}, S > n, E_{n,n}(\beta)=a^n. ( 8.55)

Начальное значение получено от (8.12). Находим обратное значение вероятности блокировки I= 1/E:

I_{n,S}(\beta)=\frac{S-n}{S(1-a)}*\{I_{n,S-1}(\beta)-a\}, S > n, I_{n,n}(\beta)=a^{-n}. ( 8.56)

Для того чтобы увеличить S, нужно, чтобы число итераций было S-n. Однако числовые ошибки накапливаются из-за умножения с (S/(S-n)) таких умножений больше чем одно, и применимость этой формулы ограничена. Поэтому рекомендуют использовать рекурсию (8.58), приведенную в следующей секции для того, чтобы увеличить S. Чтобы уменьшить S, вышеупомянутая формула аналитически точна и в числовой форме устойчива. Однако начальное значение должно быть известно заранее.

Рекурсивная формула и по n и по S

Если подставить (8.50) в (8.55), соответственно (8.51) в (8.56), мы находим:

E_{n,S}(\beta)=\frac{Sa \cdot E_{n-1, S-1}(\beta)}{n+(S-n)a \cdot E_{n-1,S-1}(\beta)},\ E_{0,S-n}(\beta)=1, ( 8.57)
I_{n,S}(\beta)=\frac{n}{Sa}*I_{n-1, S-1}(\beta)+\frac{S-n}{S},\ I_{0,S-n}(\beta)=1, ( 8.58)

которые являются рекурсивными и относительно числа обслуживающих приборов, и относительно числа источников. Обе из этих рекурсий в числовой форме точны при процессе увеличения показателей и числа итераций n (Joys, 1967 [54]).

Из материалов, рассмотренных выше, мы можем сделать следующие заключения для рекурсивной формулы Энегсета. При вычислении с увеличением значения параметра рекурсивные формулы (8.50) и (8.51) очень точны, и формулы (8.57), и (8.58) почти хороши. Рекурсивные формулы (8.55) и (8.56) неточны при увеличении значения параметра, но в отличие от других, устойчивы при уменьшении значения. Вообще, можно заметить, что рекурсия, которая является точной в одном направлении, будет неточна в противоположном направлении.

Пример 8.5.1: система с потерями Энгсета

Мы рассматриваем систему с потерями Энгсета, имеющую n=3 канала и S= 4 источника. Скорость поступления вызовов от свободного источника - \gamma= 1/3 вызовов в единицу времени, и среднее время обслуживания ( 1/\mu ) равно 1 в единицу времени. Мы находим следующие параметры:

\beta=\frac{\gamma}{\mu}=\frac 13 Эрл. ( предложенная нагрузка на один свободный источник);
a=\frac{\beta}{1+\beta}=\frac 14 Эрл. ( предложенная нагрузка на один свободный источник);
Z=1- \frac AS=\frac 34 (пиковость).

Из диаграммы переходов состояний мы получаем следующую таблицу:

i \gamma(i) \mu(i) q(i) p(i) i*p(i) \gamma(i)*p(i)
0 4/3 0 1.0000 0.3176 0.0000 0.4235
1 3/3 1 1.3333 0.4235 0.4235 0.4235
2 2/3 2 0.6667 0.2118 0.4235 0.1412
3 1/3 3 0.1481 0.0471 0.1412 0.0157
Total 3.1481 1.0000 0.9882 1.0039

Мы находим следующие вероятности блокировки.

Потери по времени:

E_{3,4} \left ( \frac 13 \right )=p(3)=0.0471,

Потери по нагрузке:

C_{3,4} \left ( \frac 13 \right )=\frac{A-Y}{A}=\frac{1-0.9882}{1}=0.0118,

Потери по вызовам:

B_{3,4} \left (\frac 13 \right )=\{\gamma (3)*p(3)\}/\left \{ \sum_{i=0}^3 \gamma(i)*p(i)\right \}=\frac{0.0157}{1.0039}=0.0156.

Заметим, что Е>В>С, и это является общим результатом для модели Энгсета (8.49) и (рис.8.6).

Применяя рекурсивную формулу (8.51) мы, конечно, получаем те же самые результаты:

E_{0,4}\left (\frac 13 \right)=1\\
E_{1,4}\left (\frac 13 \right)=\frac{(4-1+1)*\frac 13*1}{1+(4-1+1)*\frac 13*1}=\frac 47,\\
E_{2,4}\left (\frac 13 \right)=\frac{(4-2+1)*\frac 13*\frac 47}{2+(4-3+1)*\fac 13 *\frac 47}=\frac 29,\\
E_{3,4}\left (\frac 13 \right)=\frac{(4-3+1)*\frac 13*\frac 29}{3+(4-3+1)*\frac 13 *\frac 29}=\frac{4}{85}=0.0471
< Лекция 8 || Лекция 9: 123 || Лекция 10 >
Нияз Сабиров
Нияз Сабиров

Здравствуйте. А уточните, пожалуйста, по какой причине стоимость изменилась? Была стоимость в 1 рубль, стала в 9900 рублей.

Елена Сапегова
Елена Сапегова

для получения диплома нужно ли кроме теоретической части еще и практическую делать? написание самого диплома требуется?

Иван Бузмаков
Иван Бузмаков
Россия, Сарапул
Никита Сомов
Никита Сомов
Россия, Удмуртская республика