КИС_Лекции / Глава 4
.pdfГлава 4 ВЕРОЯТНОСТНО-ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ СЕТЕЙ КОММУТАЦИИ
____________________________________________________________________________________
Тогда вследствие того, что f x (0)=1, |
(4.3.8) можно представить в виде |
|||||
|
|
|
|
x − 1 |
||
|
F x ( z ) = |
λ x + |
1 / ∏ |
f i ( z ) |
||
|
|
|
i = 0 |
|
, |
|
|
x + 1 |
|
||||
|
|
z ∏ |
( z |
+ γ k ) |
||
которое имеет оригинал |
|
k |
= 1 |
|
|
|
x −1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
x +1 λx +1 / ∏ fi (γ k ) |
|
|
|
|
||
Fx (t) = 1 − ∑ |
i =1 |
|
eγ k t . |
(4.3.9) |
||
x +1 |
|
|||||
k =1 γ k ∏ (γ j + γ k ) |
|
|
|
|
i =1 j ≠k
Таким образом, в неявном виде получено аналитическое выражение (4.3.9) для расчета функции распределения промежутков времени между вызовами потерянного потока на (х-1) линиях пучка. Любое аналитическое выражение, полученное в неявном виде, обладает, по крайней мере, двумя недостатками.
Во-первых, оно лишено наглядности для анализа зависимости рассчитываемого параметра от исходных данных, поскольку последние в конечном аналитическом выражении отсутствуют.
Во-вторых, дальнейшее использование рассчитанного параметра в качестве составной части другого сложного аналитического выражения еще более усложняет решение и теряет наглядность. Эти причины заставляют искать эмпирическое
приближение точного решения. |
|
|
В некоторых случаях хорошее эмпирическое |
приближение |
дает |
экспоненциальное распределение в виде |
(4.3.10) |
|
Fx (t) =1 − exp(−λ tE x (A)) |
|
Для сравнения (4.3.9) и (4.3.10) в табл. 4.1 приведены результаты расчета функции распределения промежутков времени между вызовами потерянного потока при следующих исходных данных: λ = 10 час -1, μ = 10 час -1 , x = 5. Точное решение вычислялось по формуле
F5 (t) =1−10048.e−0.02t +0.0085e−10.5t −0.0053e−22.5t +0002.e−37.8t −0.0004e−56.8t +0.00003e−82.9t .
|
|
Табл. 4.1. |
Значение t , час |
Значение Fx(t), вычисленное |
по формуле |
|
(4.3.9) |
(4.3.10) |
0,50 |
0,0074 |
0,0152 |
1,00 |
0,0194 |
0,0302 |
2,00 |
0,0431 |
0,0595 |
5,00 |
0,1107 |
0,1421 |
10,00 |
0,2130 |
0,2641 |
50,00 |
0,7040 |
0,7842 |
100,00 |
0,9128 |
0,9534 |
Как видно из таблицы в диапазоне t>10*λ-1 приближение вполне пригодно для инженерных расчетов.
В процессе обслуживания вызовов сетью связи каждый узел коммутации каналов анализирует адресную часть поступившего сигнального сообщения и устанавливает соединение в соответствии с таблицей маршрутов. Если возникает дефицит ресурсов (каналов) на пути первого выбора, то на альтернативные пути поступает просеянный поток вызовов, не принятый к обслуживанию на пути первого выбора. На
217
Ю.Ф.Кожанов, Колбанев М.О ИНТЕРФЕЙСЫ И ПРОТОКОЛЫ СЕТЕЙ СЛЕДУЮЩЕГО ПОКОЛЕНИЯ
________________________________________________________________________
последующих узлах происходит объединение собственного потока вызовов с просеянным потоком вызовов.
Рассмотрим фрагмент сети связи, состоящей из трех узлов коммутации (A, B, C), изображенной на рис. 4.13, и проследим прохождение потоков при установлении соединений от узла А к узлу В.
Основной поток вызовов пойдет по пути первого выбора (на рис. 4.13 – жирная линия), при недостаточной емкости основного пучка – по пути второго выбора через узел С (на рис. 4.13 – пунктирная линия). Таким образом, в пучке соединительных линий от узла А к узлу В происходит просеивание потока собственных вызовов узла
А с параметром |
ΛА , а на входе узла С – объединение двух потоков: потока |
собственных вызовов с параметром ΛC и просеянного потока абонентов узла А с |
|
параметром λA . |
При расчете характеристик системы коммутации в предыдущих |
разделах было необходимо знать тип и параметры поступающего потока вызовов. Если поток собственных вызовов узла С можно вполне считать простейшим с параметром ΛC , то степень влияния просеянного потока с параметром λA на поток ΛC можно определить, только зная процедуру просеивания потока ΛΑ узлом А. Тогда встают вопросы: как классифицировать суммарный поток вызовов с параметром (ΛC + λA)? В частности, в каких случаях суммарный поток вызовов можно считать простейшим и уверенно пользоваться соответствующими формулами?
ΛΑ |
|
|
|
|
|
|
A |
λA |
|
|
B |
|
|
|
|
C |
|
|
ΛC
Рис.4.13. Объединение и просеивание потоков в сети
Ранее были рассмотрены характеристики двух просеянных потоков: в 4.3.3 – потока Эрланга и в 4.3.4 – потока Пальма. На рис. 4.14 для примера показаны: а)порождающий поток вызовов б)просеянный поток Эрланга 4-го порядка и в)потерянный поток Пальма 4-го порядка.
ts
t
а)
t
б)
t
в)
Рис. 4.14. Потоки вызовов а)порождающий поток вызовов б)просеянный поток Эрланга 4-го порядка в)потерянный поток Пальма 4-го порядка
218
Глава 4 ВЕРОЯТНОСТНО-ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ СЕТЕЙ КОММУТАЦИИ
____________________________________________________________________________________
Как видно из рисунка при одинаковой интенсивности просеянных потоков на интервале среднего времени занятия ts характер потоков существенно разнится, хотя они имеют один и тот же параметр потока.
Просеянный поток Эрланга k-го порядка описывается достаточно простым аналитическим выражением
∞ |
k |
|
i |
|
I (k, λ t) |
|
|||
Fk (t) = ∫ fk (t)dt =1− ∑ |
(λ t) |
|
e−λ t |
=1− |
e−λ t , |
||||
|
|
|
|||||||
0 |
i=0 |
i! |
|
|
|
I (k,0) |
|||
∞ |
|
|
|
|
|
|
|
|
|
где I (k, λ t) = ∫(λ t + z)k e−z dz |
– интеграл |
первого рода, |
λ = Λ/(κ+1) − параметр |
||||||
0 |
|
|
|
|
|
|
|
|
|
просеянного потока, Λ − параметр порождающего потока вызовов. |
|||||||||
Потерянный поток Пальма |
считается только численными методами и не имеет |
||||||||
явного аналитического выражения, хотя |
его |
параметр |
после просеивания |
||||||
полнодоступной схемой известен – Λ E V (Λ). |
Поэтому для аналитического описания |
потерянного потока вызовов на полнодоступном пучке приходится производить подмену потока Пальма потоком Эрланга той же интенсивности.
Пусть Λi – параметр простейшего порождающего i-го потока вызовов, λi – параметр просеянного i-го потока вызовов. Тогда процедура классификации суммарного n-мерного потока вызовов сводится к следующему.
1. Определяется эквивалентный параметр порождающего потока вызовов
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
λi |
|
|
|
λ2 |
|
|
|
λn |
|
|
|
∑λ i Λi |
|
Λ |
s |
= |
Λ |
i |
+ |
Λ |
2 |
+... + |
Λ |
n |
= |
i=1 |
. |
|||
n |
n |
n |
n |
|||||||||||||
|
|
|
∑λi |
|
|
|
∑λi |
|
|
|
∑λi |
|
|
|
∑λ i |
|
|
|
|
i=1 |
|
|
|
i=1 |
|
|
|
i=1 |
|
|
|
i=1 |
|
2. Определяется суммарный параметр просеянного потока вызовов
n
λ s = ∑λ i .
i=1
3.Вычисляется коэффициент просеивания эквивалентного порождающего потока вызовов
k = Λs / λ s – 1.
4. Рекомендуются функции распределения (ФРВО)
|
|
n |
|
|
F (t) =1− exp(−∑λ i t), |
k <1, |
|||
|
I (k, Λst) |
i =1 |
|
|
F(t) =1− |
exp(−Λst), |
k ≥1. |
||
|
||||
|
I (k,0) |
|
ЗАДАЧА 8. Показать, что суммарный поток, полученный объединением простейших потоков, – простейший. Решение.
В этом случае Λi = λ i .
|
n |
|
|
|
|
n |
|
|
|
|
|
∑Λi |
2 |
|
n |
|
∑Λi |
2 |
|
|
|
Λs = |
i=1 |
|
, |
λ s = ∑Λi , |
k = |
i=1 |
|
|
−1 |
< 0. |
n |
|
n |
|
2 |
||||||
|
∑Λi |
i=1 |
|
|
|
|
||||
|
|
|
∑Λi |
|
|
|
|
|||
|
i=1 |
|
|
|
|
i=1 |
|
|
|
|
219
Ю.Ф.Кожанов, Колбанев М.О ИНТЕРФЕЙСЫ И ПРОТОКОЛЫ СЕТЕЙ СЛЕДУЮЩЕГО ПОКОЛЕНИЯ
________________________________________________________________________
ЗАДАЧА 9. Показать, что суммарный поток, простейшему.
Решение.
n
∑λi Λi
Λ |
s |
= |
i=1 |
, |
n |
||||
|
|
|
∑λi |
|
полученный объединением многих просеянных потоков, стремится к
n |
|
λ s = ∑λ i , |
k = Λs / λs −1. |
i=1 |
|
При n Æ ∞ |
i=1 |
среднее значение порождающего потока Λs возрастает гораздо медленнее, чем λ s , что приводит к |
постепенному снижению значения k = Λ s / λ s – 1 , а это и означает, что по своим свойствам суммарный поток вызовов все более будет приближаться к простейшему.
В частности, |
при |
объединении |
n |
одинаково |
просеянных |
потоков |
с параметром λ каждый, среднее значение |
|||||||||||||
параметра порождающего |
потока остается постоянным |
Λs = Λ , |
суммарный |
параметр |
|
просеянного потока λ s = n λ , k = |
||||||||||||||
Λ / nλ – 1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ЗАДАЧА 10. |
Классифицировать |
суммарный поток, |
полученный объединением |
|
простейшего |
потока c параметром |
||||||||||||||
Λ1 = 100 час -1 и потока с параметром Λ2 = 200 час -1 , |
потерянном |
на 100-линейном пучке. |
|
|
||||||||||||||||
Решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Λ1 = 100, |
|
|
|
|
λ1 = 100, |
|
|
|
|
|
|
k1 = Λ1 / λ1 – 1 = 0, |
||||||||
Λ2 = 200, |
|
|
|
|
λ2 = 200 E100(200) =100, |
|
k2 = Λ2 / λ2 – 1= 1, |
|||||||||||||
Λs = 150, |
|
|
|
|
λs = 200, |
|
|
|
|
|
|
k = Λs / λs – 1 < 0. |
||||||||
Суммарный |
поток – простейший с |
ФРВО |
F(t) = 1 - exp (-200 t). |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ЗАДАЧА 11. |
Классифицировать |
суммарный поток, |
полученный объединением |
|
простейшего |
потока c параметром |
||||||||||||||
Λ1 = 100 час -1 и потока с параметром Λ2 = 200 час -1 , |
потерянном |
на 220-линейном пучке. |
|
|
||||||||||||||||
Решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Λ1 = 100, |
|
|
|
|
λ1 = 100, |
|
|
|
|
|
|
k1 = 0, |
|
|||||||
Λ2 = 200, |
|
|
|
|
λ2 = 200 E220(200) =2.2, |
|
k2 = 89, |
|
||||||||||||
Λs = 102.1, |
|
|
|
|
λs = 102.2, |
|
|
|
|
|
|
k < 0. |
|
|
||||||
Суммарный |
поток – простейший с |
ФРВО |
F(t) = 1 - exp (-102.2 t). |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ЗАДАЧА 12. |
Классифицировать суммарный поток, |
|
полученный объединением |
потока c |
параметром Λ1 = 100 час-1 , |
|||||||||||||||
потерянном на 84-линейном пучке и |
потока |
с параметром |
|
Λ2 = 200 час -1 , |
потерянном |
на 153-линейном пучке. |
||||||||||||||
Решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Λ1 = 100, |
|
|
|
|
λ1 = 100 E84(100) =20, |
|
|
|
|
|
k1 = 4, |
|||||||||
Λ2 = 200, |
|
|
|
|
λ2 = 200 E153(200) =50, |
|
|
|
|
k2 = 3, |
||||||||||
Λs = 171.4, |
|
|
|
|
λs = 70, |
|
|
|
|
|
|
|
|
|
k = 1,45. |
|||||
Суммарный |
поток – поток Эрланга 1.45-порядка |
с |
|
ФРВО |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
, |
171 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
I 1 . 45 |
. 4 t |
|
|
|
|
|
|
|||||||
|
|
F ( t ) = |
1 |
− |
|
|
|
|
|
|
|
exp( |
− 171 |
. 4 t ). |
|
|||||
|
|
|
I ( 1 |
. 45 , |
|
0 ) |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ЗАДАЧА 13. |
Классифицировать |
суммарный поток, |
полученный объединением |
|
простейшего |
потока c параметром |
||||||||||||||
Λ1 = 1 час -1 и потока с параметром Λ2 = 200 час -1 , потерянном |
на 220-линейном |
пучке. |
|
|
||||||||||||||||
Решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Λ1 = 1, |
|
|
|
|
|
λ1 = 1, |
|
|
|
|
|
|
|
|
|
k1 = 0, |
||||
Λ2 = 200, |
|
|
|
|
λ2 = 200 E220(200) =2.2, |
|
|
|
k2 = 89, |
|||||||||||
Λs = 137.8, |
|
|
|
|
λs = 3.2, |
|
|
|
|
|
|
|
|
k = 42. |
||||||
Суммарный |
поток – поток Эрланга |
42-порядка |
с |
|
ФРВО |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
I |
|
, |
|
137 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
42 |
|
|
. 8 t |
|
|
|
|
|
|
|||||
|
|
F ( t ) = |
1 |
− |
|
|
|
|
|
|
|
exp( |
|
− 137 |
. 8 t ). |
|
||||
|
|
|
I ( 42 |
, |
0 |
) |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
220
Глава 4 ВЕРОЯТНОСТНО-ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ СЕТЕЙ КОММУТАЦИИ
____________________________________________________________________________________
4.3.5.Модель Энгсета (M/M/V/V/N). Дискретная формула Энгсета
Всистемах распределения информации малой емкости, когда число источников нагрузки (абонентов) мало, параметр потока от одного источника сравним с суммарным потоком вызовов. В этом случае суммарный поток вызовов зависит от числа источников и от состояния {x}, (x=0,1,...V) системы обслуживания. Такой поток
называется |
примитивным |
потоком |
с |
простым |
последействием. Для |
него |
||
λ 0 ≠ λ 1 ≠... ≠ λ V −1 , μx = xμ . |
|
|
|
|
|
|||
Модель |
Энгсета |
для |
расчета |
вероятности |
потерь |
справедлива |
при |
|
предположениях: |
|
|
|
|
|
|
|
|
– параметр поступающего потока вызовов в момент занятости х выходов |
||||||||
пропорционален числу свободных источников, т.е. |
|
|
|
|||||
|
|
λ x = (N − x)α , |
0 ≤ х ≤ V, |
|
|
|||
где N – число источников |
вызовов |
(число входов |
в КП), |
α – интенсивность |
||||
поступления вызова от свободного источника; |
|
|
|
–длительность занятия подчиняется экспоненциальному распределению с параметром μ;
–вызов, не принятый к обслуживанию в момент занятости V линий пучка, теряется, не влияя на моменты поступления последующих вызовов;
–любой из V выходов пучка доступен, когда он свободен, для любого поступающего вызова.
–исходной для расчета является поступающая нагрузка;
–система находится в стационарном режиме.
Подставляя значения параметров λx и μx в формулу (4.3.1) вероятностей
стационарного процесса рождения и гибели, получим |
|
|||||
[ x ] = |
C Nx A x |
|
V |
|
−1 |
(4.3.11) |
|
∑ C Ni |
A i |
, |
|||
|
|
|
i = 0 |
|
|
|
где A = α/μ – пуассоновская нагрузка второго рода. |
|
||
Финальная вероятность |
|
|
−1 |
Pt = [V ] = EN ,V (A) = CNV |
V |
|
|
AV ∑CNi |
A i |
= CNV AV [0] |
|
|
i=0 |
|
|
определяет потери по времени и носит название формулы Энгсета.
Рассмотрим источник без потерь. Для него N=V=1, [1] = a – реальная удельная нагрузка, поэтому из (4.3.11)
|
a = |
|
|
α / μ |
= |
|
|
A |
|
|
= αr ts , |
A = |
α = |
|
|
a |
, |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
1 − a |
|
|
|
|
|
||||||||||
|
1 |
+α / μ |
1 |
+ A |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
μ |
|
|
|
|
|
|
|||||||||||||
где α r – |
реальный параметр |
потока, поступающий |
от |
источника |
в системе без |
||||||||||||||||||||
потерь, ts – |
среднее время занятия. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
При численных расчетах производят замену реальной удельной нагрузки а на |
|||||||||||||||||||||||||
пуассоновскую нагрузку второго рода по формуле А = а/(1-a). |
|
|
|
|
|
|
|
||||||||||||||||||
Параметр потерянного потока вызовов |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
λ π = λ V [V ] = (N −V )α [V ] = (N −V )α EN ,V ( A) = (N −V )α Pt . |
|
|
||||||||||||||||||||||
Потери по вызовамP = |
[V ]λ V |
= |
|
CV AV [0](N −V )α |
= |
CV |
AV |
|
= E |
|
( A). |
||||||||||||||
|
|
|
|
|
|
N |
|
|
|
|
N −1 |
|
|
N −1,V |
|||||||||||
V |
|
|
|
|
V |
|
|
|
V |
|
|
|
|||||||||||||
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
∑[x]λ x |
|
|
|
∑CNx Ax [0](N − x)α |
|
|
∑CNx −1 Ax |
|
|
|
|
|||||||||
|
|
|
|
|
x=0 |
|
|
|
|
|
|
|
x=0 |
|
|
|
|
|
x=0 |
|
|
|
|
|
|
221
Ю.Ф.Кожанов, Колбанев М.О ИНТЕРФЕЙСЫ И ПРОТОКОЛЫ СЕТЕЙ СЛЕДУЮЩЕГО ПОКОЛЕНИЯ
________________________________________________________________________
Потери по нагрузке
|
V |
|
|
|
|
|
|
|
|
|
|
|
Na − ∑x[x] |
|
V V |
|
N −V |
|
|
|
|
|
|
P = |
x=0 |
|
= |
CN −1 A |
= |
E |
N ,V |
( A) = (1 − V )E |
N ,V |
( A). |
|
|
V |
|
|||||||||
y |
Na |
|
|
N |
N |
|
|||||
|
|
∑CNx Ax |
|
|
|
|
|||||
|
|
|
|
x=0 |
|
|
|
|
|
|
|
Среднее число занятых линий (обслуженная нагрузка)
V |
Nα |
|
|
V |
|
|
Vm = ∑x[x] = (N − (N −V )Pt )a = |
|
1 |
− (1 |
− |
|
)EN ,V ( A) . |
|
N |
|||||
x=0 |
α + μ |
|
|
|
Прямой расчет формулы Энгсета во многих практических случаях невозможен из-за переполнения разрядной сетки вычислительного устройства (при больших значениях A и V). Поэтому для ее расчета пользуются рекуррентным соотношением
|
|
E N ,V ( A ) = |
A ( |
N − V + 1) E N ,V −1 ( A ) |
, |
|
|
|
|
|
|
|
|||||
|
|
A ( N |
− V + 1) E N ,V −1 ( A ) + V |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
последовательно вычисляя |
EN ,1 (A), |
EN ,2 (A), . |
. . |
|
EN ,V (A), |
при начальном |
|||||||||||
значении EN ,0 (A) =1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Следует запомнить: |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
EN ,0 ( A) = 1 , |
EN ,V (0) = 0 , |
EN ,0 (0) = 0 , |
|
||||||||||
P |
|
N >V |
= (N −V +1)A −V |
, |
P |
|
N >V = |
(N −V )A −V , |
P |
|
N >V |
= |
|
N −V |
(N −V +1)A −V . |
||
|
|
|
|
|
|||||||||||||
t |
|
A→∞ |
(N −V +1)A |
|
c |
|
A→∞ |
(N −V )A |
y |
|
A→∞ |
|
|
N |
(N −V +1)A |
В инженерной практике при выполнении обеих условий N>>V и N>100 для расчета потерь рекомендуется пользоваться моделью Эрланга, в противном случае следует использовать модель Энгсета.
ЗАДАЧА 14. Определить потери по времени, вызовам и нагрузке для N=5, V=5. Удельная нагрузка a = 0,8 Эрл. Решение.
В нашем случае N=V, поэтому используем модель Энгсета.
A = |
|
a |
|
= |
|
|
0.8 |
= 4 |
Э р |
л. |
|
|
|
|
|
|
|
||
1 − a |
1 |
− 0 .8 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
P = |
|
C 55 4 5 |
|
|
= 0.327 , |
P |
|
= |
C 45 4 5 |
= 0, |
P |
|
= (1 − |
5 |
) P = 0. |
||||
|
∑5 |
|
|
|
|
|
|
|
|
||||||||||
t |
|
C 5i 4 i |
|
|
|
c |
|
∑5 |
C 4i 4 i |
|
|
y |
|
5 |
t |
||||
|
|
i = 0 |
|
|
|
|
|
|
|
|
|
i = 0 |
|
|
|
|
|
|
|
ЗАДАЧА 15. Оператор использует для доступа в Интернет учрежденческую АТС на N=50 номеров через пучок V=30 линий. Измерения потерь по вызовам оказались равными 0.026. Оператор хочет снизить потери до значения не более 0.001. Сколько нужно добавить линий ?
Решение.
Из уравнения Pc = E N-1,V (A) при известных значениях Pc =0.026, N=50, V=30 находим параметр свободного источника
A=0.1 Эрланг.
При известных значениях Pc =0.001, N=50, A=0.1 Эрланг находим, что уравнение Pc =0.1 ≥ E N-1,V (A) удовлетворяется при емкости пучка V=35. Следовательно, нужно добавить 5 линий.
4.3.6. Модель Энгсета (M/M/V/V/N). Интегральная формула Энгсета
Пусть {x} – состояние V-линейного пучка (наличие х установленных соединений), [x] – стационарная вероятность состояния {x}, (0 ≤ x ≤ V). Тогда первый потерянный вызов, поступивший на (V+1)-ю фиктивную линию пучка при условии не освобождения занятых линий, будет принадлежать потоку Энгсета (V-x)-го порядка, поскольку потери наступят только после поступления (V-x+1)-го вызова.
222
Глава 4 ВЕРОЯТНОСТНО-ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ СЕТЕЙ КОММУТАЦИИ
____________________________________________________________________________________
Изображение плотности условно потерянного потока Энгсета (V-k)-го порядка имеет вид
|
|
V |
|
λ i |
|
|
|
f k ( z ) = L [ f k ( t )] = ∏ |
|
|
|
, |
|||
λ i |
+ |
z |
|||||
а сам оригинал |
|
i = k |
|
||||
V |
V |
|
|
|
|
||
V |
|
|
|
|
|||
f k ( t ) = L −1 [ f k ( z )] = ∑ e − λ i t ∏ λ i / |
∏ ( λ j − λ i ) . |
||||||
i = k |
i = k |
j = k , j ≠ k |
|
|
|
Плотность потока условно потерянных вызовов, поступающих на (V+1)-ую a фиктивную линию пучка, находим по формуле полной вероятности
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
V |
x |
−1 |
|
|
λ |
j |
|
V |
|
|
|
|
∏λ j |
||||
V |
[ x] fV −x (t) = ∑ |
∏ |
|
|
|
[0]∑e −λ j t |
|
|
j =k |
|
= |
|||||||||||||
|
μ j +1 |
|
V |
|
|
|
||||||||||||||||||
fV (t) = ∑ |
|
|
|
|
|
|
x =0 j =0 |
|
|
i=x |
|
|
∏(λ j − λi ) |
|||||||||||
x=0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j =k , j ≠i |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.3.12) |
|
|
|
|
|
|
V |
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
V |
|
−λ |
t |
|
∏ μ j |
|
|
|
|
|
|
x |
|
∏ μ j |
|
|
|
|
|
|
||||
= [V ]λV ∑e |
|
j =x +1 |
|
|
|
|
|
1 + ∑ |
|
i = j |
|
|
|
|
|
|
|
|||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
||||||
|
V |
(λ − λ ) |
x −1 |
(λ − |
λ |
|
|
|
|
|||||||||||||||
x =0 |
|
|
|
|
|
|
|
j =1 |
|
|
) |
|
|
|||||||||||
|
|
|
∏ |
|
j |
|
x |
|
|
|
|
|
∏ i |
|
|
x |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
j =x+1 |
|
|
|
|
|
|
|
|
|
i = j −1 |
|
|
|
|
|
|
|
|||
Подставляя в (4.3.12) |
|
λx |
= |
|
(N - x)α , |
μ x = |
x μ , |
получим |
V
fV (t) =[V ](N −V )α e−Nαt ∑CVx (−1)V −x (1+ μ /α )x exαt (μ /α)V −x =
x=0
=[V ](N −V )α e−Nαt ((1+ μ /α)eαt − μ /α)V .
Очевидно,
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ fV (t)dt =1, |
|
|
|
|
|
|||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
откуда для любого значения V |
|
|
|
|
|
|
|
|||||||
[V] = EN,V (A) = |
|
|
1 |
|
|
|
|
= |
|
|
|
|||
∞ |
|
|
|
|
|
|
|
|
|
|
||||
|
|
∫(N −V)α ((1+μ /α)eαt −μ /α )V e−Nαt dt |
|
|
|
|
|
|||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.3.13) |
|
= |
|
(α / μ)V |
|
|
|
= |
|
AV |
|
|
= |
AV |
||
|
|
|
|
|
|
|
|
|
|
, |
||||
∞ |
|
|
|
∞ |
− |
|
z |
|
I(N,V, A) |
|||||
|
∫(1+α / μ−e−αt )V e−(N −V )αt d((N −V)α t) ∫(1+ A−e |
N −V )V e−z dz |
|
|||||||||||
|
|
|
|
|
||||||||||
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|||
|
|
∞ |
− |
z |
|
|
|
|
|
|
|
|||
где |
I (N,V , A) = ∫(1 + A − e |
|
)V e−z dz – интеграл |
второго рода. |
|
|
||||||||
N −V |
|
|
||||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
223
Ю.Ф.Кожанов, Колбанев М.О ИНТЕРФЕЙСЫ И ПРОТОКОЛЫ СЕТЕЙ СЛЕДУЮЩЕГО ПОКОЛЕНИЯ
________________________________________________________________________
Выражение (4.3.13) – |
|
интегральная формула Энгсета, позволяющая |
||||
производить расчет потерь при любом действительном значении V. |
||||||
Для интегральной формулы Эрланга верно рекуррентное соотношение |
||||||
E N ,V ( A ) |
= |
|
A (N − V + 1) E N ,V −1 ( A ) |
|
||
|
A (N − V |
+ 1)E N ,V −1 ( A ) + V |
||||
|
|
|
||||
с начальным значением |
|
|
|
|
|
|
|
|
EN ,v (A) = |
A v |
, |
|
|
|
|
I (N , v, A) |
|
|||
|
|
|
|
|
|
где v = V – INT (V) – дробная часть емкости пучка.
Интегральная формула Энгсета используется при расчете потерь в пучке малой емкости при нецелочисленном значении емкости пучка.
4.4.Система с ожиданием (Blocked Calls Queued, BCQ)
4.4.1.Система M/M/V. Стационарные вероятности
Всистеме с ожиданием при недостатке обслуживающих приборов блокированные вызовы становятся в очередь и обслуживаются в порядке поступления. Такая дисциплина обслуживания применяется в маршрутизаторах при использовании алгоритма “с наибольшими усилиями” (Best Effort).
Диаграмма переходов Марковского процесса системы обслуживания с ожиданием изображена на рис. 4.15.
λ |
|
λ |
|
λ |
λ |
|
λ |
|
|
|
|
||||
0 |
1 |
μ2 |
. . . |
|
X |
. . . |
V |
μ1 |
|
μ |
x |
μx+1 |
|
μV |
|
|
|
|
|
|
μV λ
V+1
μV λ
. . .
Рис. 4.15. Диаграмма переходов системы с ожиданием
Модель справедлива при предположениях:
–вызовы, поступающие на вход системы, образуют пуассоновский поток постоянной интенсивности с параметром λ;
–длительность занятия подчиняется экспоненциальному распределению с
параметром μ;
– любой из V выходов пучка доступен, когда он свободен, для любого поступающего вызова;
224
Глава 4 ВЕРОЯТНОСТНО-ВРЕМЕННЫЕ ХАРАКТЕРИСТИКИ СЕТЕЙ КОММУТАЦИИ
____________________________________________________________________________________
–вызов, не принятый к обслуживанию в момент занятости V линий пучка, поступает в очередь и обслуживается в порядке поступления (дисциплина FIFO);
–исходной для расчета является поступающая нагрузка;
–система находится в стационарном режиме.
Используя уравнения стационарного процесса рождения и гибели при A<V, имеем
|
|
|
(λ / μ )x |
|
|
|
x = 01,,...V , |
|||||||||
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
[0], |
|
|
||||
[x ] = |
|
|
|
x ! |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
(λ / μ )V |
λ x −V |
|
x =V +1,V +2,... |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[0], |
|
|
|
|
|
V ! |
|
|
|
|
||||||||
|
|
|
|
|
|
Vμ |
|
|
||||||||
или, принимая A=λ/μ , |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
Ax |
[0], |
|
|
|
|
|
|
|
x = 0,1,...V , |
|||
|
|
|
|
|
|
|
|
|
|
|||||||
[x] = |
|
|
x! |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
AV A |
x−V |
[0], |
|
x =V +1,V + 2,..., |
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
V! V |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где [0] определяется из нормирующего условия
V |
A x |
|
A V ∞ |
|
A x −V |
|
V |
A x |
|
A V A |
|
|
||||
[0]−1 = ∑ |
|
+ |
|
∑ |
|
|
|
|
= ∑ |
|
+ |
|
|
|
|
. |
|
|
|
|
|
|
|
|
|||||||||
x =0 |
x ! |
|
V ! x =V +1 |
V |
|
A <V x =0 |
x ! |
|
V ! V − |
A |
|
|||||
|
|
|
|
Потери по времени равны вероятности того, что все линии заняты, она же – вероятность ожидания начала обслуживания
∞ |
AV 1 |
|
|
|
E (A) |
|
|
|
|
|
|
||||
Pt =P(>0) = ∑[x]= |
|
|
|
[0]= |
|
|
V |
|
|
. |
|
|
|
(4.4.1) |
|
V ! 1− A /V |
|
A |
( |
|
) |
|
|
|
|||||||
x =V |
1− |
V |
|
|
|
|
|
||||||||
|
|
|
|
V |
|
|
|
|
|
||||||
|
|
|
|
|
1 |
−E (A) |
|
|
|
|
|
||||
(4.4.1) – вторая формула Эрланга. |
|
|
|
|
|
|
|
|
|
||||||
Среднюю длину очереди находим из соотношения |
|
|
|
|
|||||||||||
|
|
∞ |
|
|
|
|
A* EV (A) |
|
|
V |
|
||||
n = ∑(k −V ) *[k] = |
|
|
* |
|
. |
||||||||||
(V − A) + A* EV (A) |
V − A |
||||||||||||||
|
k =V |
|
|
|
|
4.4.2. Система M/M/V. Функция распределения
Обозначим через P(>t) вероятность ожидания начала обслуживания больше, чем t, а через Pk(>t) – условную вероятность ожидания начала обслуживания при поступлении вызова в момент, когда в системе находятся k вызовов, из них V вызовов обслуживаются, а (k-V) вызовов ожидают обслуживания. По формуле полной вероятности
∞
P(> t) = ∑[k]Pk (> t).
k =V
225
Ю.Ф.Кожанов, Колбанев М.О ИНТЕРФЕЙСЫ И ПРОТОКОЛЫ СЕТЕЙ СЛЕДУЮЩЕГО ПОКОЛЕНИЯ
________________________________________________________________________
Вероятность Pk(>t) есть вероятность того, что за время t произойдет не более 1,2,...(k-V) освобождений и вычисляется c использованием распределения Пуассона
|
|
|
|
|
|
|
|
|
|
Pk (> t) = ∑(Vμ t) |
i |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
e−Vμ t |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k −V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i=0 |
i ! |
|
|
|
|
|
|
|
||
Поэтому |
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
P(> t) = ∑[k]∑(Vμ t) |
e−Vμ t = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
∞ |
k−V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
k=V |
i=0 |
i! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
|
|
∞ |
|
k−V |
k−V |
(Vμ t) |
i |
|
|
V |
1 |
|
−(Vμ −λ )t |
|
−(Vμ −λ )t |
|
||||||
|
A |
|
−Vμ t |
|
|
|
|
|
|
A |
|
|
|
|||||||||||
= |
|
[0]e |
|
∑(A/V ) |
|
|
∑ |
|
|
|
= |
|
|
|
|
|
[0]e |
|
= P(> 0)e |
|
. |
|||
V! |
|
|
|
i! |
|
V! 1− AV |
|
|
||||||||||||||||
|
|
|
k=V |
|
|
|
i=0 |
|
|
|
V >A |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
При A>V стационарного процесса обслуживания вызовов нет: длина очереди будет монотонно стремиться к бесконечности, каждую единицу времени занятия очередь будет увеличиваться на (A-V) вызовов.
Функция распределения времени ожидания (ФРВО) начала обслуживания Fw(t) или вероятность окончания времени ожидания начала обслуживания за время t
F |
(t) =1− P(> t) =1− P(> 0)e−(Vμ−λ)t , |
(4.4.2) |
w |
|
|
а ФРВО самого обслуживания |
|
|
|
F s ( t ) = 1 − e − μ t , |
(4.4.3) |
поэтому ФРВО пребывания вызова в системе (ФРВО конца обслуживания вызова)
t |
|
F e (t ) = ∫ F w (τ ) dF s (t − τ ). |
(4.4.4) |
0 |
|
4.4.3. Система M/M/1
Принимая в формулах (4.4.2) – (4.4.4) значение V=1, A=a, получим:
вероятность постановки в очередь – P(> 0) = A = a = λ / μ , |
|
ФРВО начала обслуживания – F (t) =1− λ e−(μ−λ) t , |
|
w |
μ |
|
среднее время ожидания начала обслуживания – tw = λμ μ 1−λ = μa 1−1 a , ФРВО конца обслуживания – F e(t) =1−e−(μ−λ) t =1−e−(1−a) μ t ,
среднее время ожидания конца обслуживания – te = μ 1− λ ,
вероятность нахождения вызова (пакета) на k-ом месте ожидания – [k] = ak +1 (1− a) ,
средняя длина очереди (в пакетах) |
n = |
|
a2 |
. |
||
1 |
− a |
|||||
|
|
|
При обслуживании только одного источника очередь на обслуживание отсутствует,
ФРВО конца обслуживания –
F e(t) =1−e−(μ−λ) t λ=0 =1−e−μ t .
226