
- •Отношения между е, в и с
- •Расчеты по формуле Энгсета
- •Рекурсивная формула для п
- •Рекурсивная формула по s
- •Рекурсивная формула и по п и по s
- •Пример 8.5.1: система с потерями Энгсета
- •Пример 8.5.2: Ограниченное число источников
- •Распределение Паскаля (отрицательное биноминально е распределение)
- •Усеченное распределение Паскаля
•р
п
адр) =
(8.27)Этот результат можно интерпретировать следующим образом. Вероятность того, что попытка вызова от случайного источника (абонента) будет отклонена, равна вероятности того, что остальные (S-\) источники заняли все п каналов. Это называется теоремой поступления, и можно показать, что она справедлива для систем с явными потерями и для систем с ожиданием и ограниченным числом источников. Результат основан на вычислении произведения среди источников и свертывании источников. Поскольку Е увеличивается, когда увеличивается S, мы имеем:
5„,ЛР)= £Л)5-1(Р)<£я>ИР).
Теорема 8.1. Теорема поступления: Для всех систем с ограниченным числом источников, случайный источник при поступление заявки будет наблюдать состояние системы так, как будто сам источник не принадлежит этой системе.
Свойство PASTA включено в этот случай, потому что бесконечное число источников минус один есть бесконечное число.
П
Y
П
(8.28)
/=0
п
(8.29)
Yj p-(j- i)-p(о - р-(*у- п)-р(п),
(8.30)
= p. (S- Y)-р-(S-п) • Е,
поскольку Е = Еп, S (Р) = р(п). Последнее уравнение решается относительно Y:
Y= j^-{S-(S-n)-E) . (8.31)
Потери по нагрузке С - CnS(A). Это самая важная характеристика потерь. Предложенная нагрузка дается в (8.20), и мы получаем:
С
=
S/i (3
■{S- (S-п) • Е)
1 + р 1 + р
.Sp 1 + Р
С = --Е. (8.32)
Мы можем также найти обслуженную нагрузку, если знаем потери по вызовам В. Число попыток принятия вызовов от источника, который находится в свободном состоянии, в среднем Уув единицу времени прежде, чем источник сгенерирует одну попытку вызова - 1(1-6), и каждый принятый вызов имеет среднюю продолжительность 1/ц- Таким образом, обслуженная нагрузка на один источник есть соотношение времени: когда источник является занятым, она будет:
у= V-Wp
1/7 + (1 - В)/)х '
Полная обслуженная нагрузка будет:
Приравнивая эти два выражения для обслуженной нагрузки (8.31) и (8.33), мы получаем следующее отношение между ЕиВ:
Е = — . (8.34)
S-n 1+0(1 -В) v ;
Число попыток вызова в единицу времени:
П
А = ^p(i)-(S- /) у
i=0
А = (S - Y)-y, (8-35)
где Y — обслуженная нагрузка (8.28). Таким образом, из уравнения (S- Y) среднее число свободных источников является очевидным. Исторически, полная предложенная нагрузка была определена как Это, однако, вводит в заблуждение, потому что мы не можем принять, что каждая повторная попытка вызова имеет среднее время пребывания в системе, равное 1/ju. • Это определение создает большое неудобство, потому что предложенная нагрузка по этому определению зависит от состояния системы (числа занятых каналов). Также возможно, что немногие из доступных обслуживающих много попыток вызова устройств блокированы, а свободные источники с более высоким средним временем поступления вызовов генерируют больше попыток вызова в единицу времени.
Потерянная нагрузка:
Ае = А-С
Р S-n Г
S
= ft?" .Е (8.36)
1+ р
Продолжительность состояния /. Оно является экспоненциально распределенным с интенсивностью:
у (/) = (S - /') • у + i • ц, 0 < / < п, 1
\ (8.37)
у (л) = лц, i = п. J
Функция увеличения:
f„,s(A)= Y„+i - Y„. (8 38)
Выше мы определили вероятности состоянияp(i), согласно предположению о статистическом равновесии, как часть времени, которое система находится в состоянии г, то есть как математическое ожидание времени. Мы можем также исследовать, как выглядит система, когда выполняется равновесие между поступлением и возвратом вызова отправляющим
источником (пользователем) (математическое ожидание вызова). Если мы рассматриваем одну единицу времени то, в среднем, в системе будет (S-i) jp(i) источников в состоянии [/] перед моментом поступления вызова, и если вызов принят, то он переведёт систему в состояние [/+1].
Источники, которые наблюдают систему в состоянии п, блокированы или остаются свободными. Поэтому источники, от которых поступает вызов прибытия, наблюдают систему в состоянии [/'] с вероятностью:
n„Atm= . / = 0,1...... <8.39)
Ё (s~J')y ’р0)
;=о
Используя метод аналогового дифференцирования выражения
, мы можем показать, что в соответствии с теоремой поступления (теорема 8.1) мы имеем:
ли,5',р(0 = Pn,s-\,p(i~ 1), / = 0,1,...,и. (8.40)
Когда источник оставляет систему, система находится в состоянии [/-1] с вероятностью:
У„АР(/-1)
= п1Ц‘Р(1)
, /=1,2,...,/!. (8.41)
YjJv 'р0)
/=1
Применяя уравнения сечения, мы немедленно убеждаемся, что это соотношение идентично (8.39), если включить в рассмотрение блокированных абонентов. Таким образом, источники покидают систему в том же самом состоянии, в котором они в нее прибывают. Процесс будет обратим и нечувствителен к распределению времени обслуживания. Если бы мы непрерывно наблюдали состояние системы, то не могли бы определить, выполняется ли время прямое действие (вызов) или обратный (выход из системы).
Отношения между е, в и с
Из (8.34) мы получаем следующее отношение между Е ~ EnS (р) и
в= B„,s№= En,s-№-
Е- sh ■ .Тр^гт, ! = ^{<1+»-Н' <8-42>
„ {S- п)-Е- (1+ (3) 1 1 Г S 1 , „1 /0
в=
SHS-nl-E-р
ИЛИ
5= (8'43>
Выражения с правой стороны линейны по отношению к вероятностям блокировки. В (8.32) мы получили следующее простое отношение между С и Е:
С = ‘ Е, (8.44)
Е = — ■ С. (8.45)
S - п
Если мы в (8.44) подставим Е, выраженное через (8.42), то мы получаем С, выраженное через В\
В (8.46)
1 + р-(1- ву
В = <1±Ш. (8.47)
1 + р С
Это отношение между В и С является общим и может также быть получено следующим образом. Обслуженная нагрузка У соответствует (У*|х) принятых попыток вызова в единицу времени. Среднее число свободных источников - (Л--К), так что среднее число попыток вызова в единицу времени - (S-Y)j (8.35).
Потери по вызовам тогда равны:
в = (S- Y)у- У-ц
(S- Y)у
(S- У)Р - У (S- Y) р •
По определению У= (1- С), и из (8.20) имеем 5=А(1 + р)/р. Подставляя это выражение, мы имеем:
= Л(1 + Р)-Л(1- С)Р-Л(1- С)
Л(1 + р)-Л(1-С)р
о _ а + Р)с 8 ~ TTW
Из последнего уравнения видно, что для малых значений потерь по вызовам В потери по нагрузке Z в несколько раз больше, чем потери по вызовам:
C*j!j = Z.5. (8.48)
Расчеты по формуле Энгсета
Если мы пробуем вычислить числовые значения непосредственно по формуле Энгсета из (8.26) (потери по времени Е), то возникают проблемы расчета для больших значений Sun. Ниже мы получим различные рекурсивные формулы для Е и для его обратной величины I = '/£. Когда потери по времени Е известны, можно просто получить потери по вызовам В и потери по нагрузке С, используя формулы (8.43) и (8.44). В числовой форме также просто найти любой из этих четырех параметров S, (3, п, Е, когда мы знаем три из них. Математически мы можем предположить, что п и S не являются целыми числами.
Рекурсивная формула для п
Из общей, рекурсивной формулы (7.27) для п, используя Хх = (S-x)y и Р = V\i, мы получаем
£*,ИР)
—
1 + Р)'
F m - 0?-* + 1)Р-Дг-1„у(Р) р .ост
£«*<(» - * + (j-*+i)p..e,_uo)- <8-50>
Используя обратную формулу потерь по времени 7 S(P) = Уе 5(3), мы находим рекурсивную формулу:
4,5(Р)
= 1 + (ty—/»-ЫР),
4j(P)
= 1. (8-51)
Число итераций равно п. Обе формулы (8.50) и (8.51) аналитически точны и представляют собой устойчивые при приближенных вычислениях и точные рекурсии для растущих значений х Однако для уменьшающихся значений х числовые ошибки накапливаются и рекурсии недостоверны.
Рекурсивная формула по s
Обозначим нормализованные вероятности состояния системы с я каналами и 5—1 источниками рп Мы получаем вероятности состояния системы с 5 источниками и я каналами, сочетая эти вероятности состояния с вероятностями состояния одиночного источника, которые обозначаются {/?, ,(0) = 1 - a, pt ,(1) = а}. Мы тогда получим состояние от нуля до и + 1, ограничиваем пространство состояний до и и нормализуем вероятности состояния (см. пример 3.2.1) (принимая, чтор(х) = 0, когда х< 0):
qn,s(i) = (l-a)'Pn,s-i(i) + a’pn,s-i(i-l), /=0,1,...,п. (g 52)
Полученные вероятности состояния qn s(i) не нормализованы, потому что мы ограничиваем состояния числом [я] и исключаем последние элементы, начиная с состояния [n+l],qn s(n+l)=a 'рп 5_,(я). Нормализованная вероятность состояния рп s(i) для системы с 5 источниками и я каналами, таким образом, получается из нормализованной вероятности состояния p„ s^(i) для системы с 5-1 источником:
Р*М0
=
-п——-,
i
= Q,\,...,n. (8.53)
- а • p„iS-i(n)
Потери по времени Еп S(P) для системы с 5 источниками могут быть выражены потерями по времени Еп w (Р) для системы с 5-1 источниками. Подставляя (8.52) в (8.53), получаем:
En,s( Р) = p„,s(n)
= (1- а) •p„,s-\(n) + q -pnrS-\(n-1)
1- a-pntS-\in)
(l-e)-^-,(P) + fl-^5^y^(J-i(P)
1- а-Е^чф) ’
где мы использовали уравнение равновесия между состоянием [я -1, 5-1 ] и состоянием [я, 5-1]. Заменяя а на(8.8), мы получаем:
En,s-xW)+ c^£„,S-i(P)
*•*<» i
+ n-SuU)
•
Таким образом, получаем следующую рекурсивную формулу.
Д*№>-^Т+»'п^(го-
■S>"'
£'"№)=“" (8'55>
Начальное значение получено от (8.12). Находим обратное значение вероятности блокировки I = Ур
Л,ИР)
= 5^)•
{V-i(P)-«}> 1п,пФ)
=
а~п
(8.56)
Для того чтобы увеличить 5, нужно, чтобы число итераций было S'-и. Однако числовые ошибки накапливаются из-за умножения с (S/(S-n))> таких умножений больше чем одно, и применимость этой формулы ограничена. Поэтому рекомендуют использовать рекурсию (8.58), приведенную в следующей секции, для того, чтобы увеличить S. Чтобы уменьшить S, вышеупомянутая формула аналитически точна и в числовой форме устойчива. Однако начальное значение должно быть известно заранее.