Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

КИС_Лекции / Глава 4

.pdf
Скачиваний:
84
Добавлен:
15.03.2015
Размер:
1.12 Mб
Скачать

Глава 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) =110048.e0.02t +0.0085e10.5t 0.0053e22.5t +0002.e37.8t 0.0004e56.8t +0.00003e82.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 ez 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) =1exp(λ 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 )α eNαt CVx (1)V x (1+ μ /α )x exαt (μ /α)V x =

x=0

=[V ](N V )α eNαt ((1+ μ /α)eαt μ /α)V .

Очевидно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fV (t)dt =1,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

откуда для любого значения V

 

 

 

 

 

 

 

[V] = EN,V (A) =

 

 

1

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(N V)α ((1+μ /α)eαt μ /α )V eNα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+ Ae

N V )V ez dz

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

где

I (N,V , A) = (1 + A e

 

)V ez 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

xV

[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 ! 1A /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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eVμ t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=0

i !

 

 

 

 

 

 

 

Поэтому

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(> t) = [k](Vμ t)

eVμ t =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=V

i=0

i!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

kV

kV

(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! 1AV

 

 

 

 

 

k=V

 

 

 

i=0

 

 

 

V >A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При A>V стационарного процесса обслуживания вызовов нет: длина очереди будет монотонно стремиться к бесконечности, каждую единицу времени занятия очередь будет увеличиваться на (A-V) вызовов.

Функция распределения времени ожидания (ФРВО) начала обслуживания Fw(t) или вероятность окончания времени ожидания начала обслуживания за время t

F

(t) =1P(> t) =1P(> 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 11 a , ФРВО конца обслуживания – F e(t) =1e(μλ) t =1e(1a) μ t ,

среднее время ожидания конца обслуживания – te = μ 1λ ,

вероятность нахождения вызова (пакета) на k-ом месте ожидания – [k] = ak +1 (1a) ,

средняя длина очереди (в пакетах)

n =

 

a2

.

1

a

 

 

 

При обслуживании только одного источника очередь на обслуживание отсутствует,

ФРВО конца обслуживания –

F e(t) =1e(μλ) t λ=0 =1eμ t .

226

Соседние файлы в папке КИС_Лекции