Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 219

Скачиваний: 3

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ми); если же заявка застает свободным хотя бы один

канал, то она обслуживается любым из

свободных

каналов.

 

' П о д со с то я н и е м с и с т е м ы Хк будем

понимать

ч и с л о з а н я т ы х о б с л у ж и в а н и е м к а н а л о в не-

заівисимо от

их

порядковых номеров.

В

таком

случае

состояние х0

означает, что все каналы

свободны; состоя­

ние Хи—-занято

k

каналов. В состоянии

Xh (k=l,

2,

Я

 

 

 

 

 

 

 

 

 

 

Чкф

miß

Wf?

 

 

 

 

Рис. 2.8.

 

 

 

 

N1) на систему воздействует два потока: поток

заявок

с интенсивностью

%h—X, который

. переводит

систему

в состояние Xh+u и поток освобождений всех k занятых обслуживанием каналов, стремящийся перевести систе­ му в состояние Хи-і. Интенсивность потока освобождений

равна

 

Работа

подобной системы массового

обслуживания

иллюстрируется графом состояний (рис

2.8). Используя

мнемоническое правило, выпишем систему дифференци­ альных уравнений для вероятностей состояний

% ^ = - * / > . ( ' ) + w \ W ;

dp*, (t) = -

(я +

kv)рк (t)+

хРк_,

(t)

+

+ (H-i)w>*+ 1 tf)

 

 

 

( £ = 1 , 2 ,

...,N-l),

 

 

 

dpN (t)fdt =

- NwN (t) + lpN_x

(t).

 

В стационарном режиме (2.52) превращается в си­

стему алгебраических

уравнений

 

 

 

0 = - ( * +

/e!i.)pft_j_Apft_I

+

 

 

+ ( / г + 1 ) № ^

( Ä = 1 , 2

Л^— 1),

{ ' '

57


Ее следует решать совместно с условием нормировки (2.47).

•По аналогии с (2.48) введем обозначения

 

Я Й = — %ph-i+k\iph

 

(Is^kg^N).

 

(2.54)

С учетом

(2.54) вместо

(2.53)

получаем

 

 

я.і=0, îtft-M—nf t = 0,

jtjv = 0

( й = 1 , 2, .. , N—l).

(2.55)

Система

(2.55)

имеет

решение

J T / I = Ö

(iß—1,

2,

..., N).

•Следовательно,

 

 

 

у 1

 

 

 

 

 

X

/

 

 

 

 

P>=(-j)'

wPo

р Л = ( т ) 6

-W Po(k = 0,\,2,

...,N).

 

 

 

 

 

 

 

 

(2.56)

Из условия нормировки (2.47) с

помощью

(2.56)

находим

 

 

 

 

 

 

 

 

 

 

*-(S тдгГ-

 

 

( 2 '5 7 )

так что финальные вероятности определяются

соотно­

шением

 

 

 

 

 

 

 

 

 

*~ (т)Ч (É (т)ЧГ •

 

<2'58)

Располагая

формулами для

вероятностей

состояний

ph, можно получить ряд других

важных характеристик

системы обслуживания. Например, вероятность обслу­ живания заявки р о б е л равна вреояпности того, что.заяв­ ка застанет свободным хотя бы один канал р о б с л = 1 — P N - Среднее число занятых каналов < & > ' также вычисляет­ ся с помощью вероятностей (2.58) :

N

<é> = 2 kPk

Нетрудно также найти соотношения для среднего времени 'простоя канала, для вероятностей занятости и полной загрузки системы {5, 8].

®


2.5.Импульсные марковские процессы с дискретными

состояниями

В статистической радиотехнике, теории автоматиче­ ского регулирования можно указать целый ряд задач, в которых используются марковские разрывные процес­ сы и м п у л ь с н о г о хар актер а [3, 17—20].

Характерной особенностью большинства импульсных процессов является наличие опорного (в частном случае, нулевого) уровня, на котором начинаются и заканчива­ ются отдельные импульсы. Это обстоятельство непосред­ ственным образом отражается на матрице вероятностей перехода (2.33), которая в данном случае приобретает вид

P ( U + A0 =

- 1 + Х И ( / ) Д <

 

\lt(t)àt

Х„(0А'

•••

(0Д<

Х„(0Д*

 

1 +Хн(0Д<

о

... о

Х„

0

 

i+K„(t)àt

...

о

| _ Ѵ № "

о

 

о

...

i+xNN(t)&t

(2.59)

Здесь в качестве опорного уровня выбран уровень (со­ стояние) с номером 1. Ненулевыми в матрице (2.59) яв­ ляются элементы, которые соответствуют переходам из состояния ХІ во все остальные; из состояний ХІ (іфі) переходы возможны только в опорное состояние ХІ и в самоё себя, остальные переходы запрещены. Одна из возможных реализаций импульсного процесса, описывае­

мого

матрицей (2.59), представлена на рис. 2.9.

Для марковского процесса потоки, переводящие си­

стему

.из состояния в состояние, являются

пуассонов-

скими

и, следовательно, р а с п р е д е л е н и е

д л и т е л ь ­

н о с т е й импульсов (равно как и пауз на первом уров­ не) является э к с п о н е н ц и а л ь н ы м с соответствую­ щей интенсивностью. Составляя по матрице (2.59) граф состояний и используя мнемоническое правило (§ 2.3),

получим систему уравнений (2.40)

для вероятностей со-

' стояний иміпульоного процесса в виде

 

Чг=

- S х« (*) л W + £

(0 * (0.

(£-6 °)

 

 

і=2

 

59



*ej[P = -liAt)Pi(t)

+ ki(t)PAt).

/ = 2,3

N.

В частном случае,

когда N = 2,

рассматриваемый

процесс вырождается

в марковский процесс с двумя со­

стояниями, возможная реализация которого изображена на рис. 2.1. Тогда из (2.60) получаем

 

=

- аІ 2 (0 Рг (t) + я„ (0

Р2 (*).

(2.61)

^

= -

я и (0 р2 (0 + *„ (0 л

(0.

(2-62)

Поскольку для каждого момента времени имеет место соотношение Pi(i)+pz(t) = l, то вместо двух уравнений (2.61), (2.62) можно рассматривать одно из них.

N

3

1

Рис. 2.9.

Исследуем подробнее однородный процесс с двумя состояниями, для которого

' Xa(t)=3^2=<iJ

•ÄaiOO=^ai='ß.

(2.63)

Тогда вместо (2.61) имеем

^Р- = - aPl'(t) + ßp2 (0 = - YP, (0 + ß, (2.64)

где

(2.65)

Решение уравнения (2.64) при начальном условии Рі(іо) имеет вид

P , ( f o . 0 = P 1 ( g ^ ( ^ ) + J r ( l - e - ^ > )

(2.66)

60

Или

 

 

 

p 2 ( C ) = P 2 ( ' o ) e - 7 ( ' - ' ° 4 - f - е _ Т < ^ 0 ) ) -

( 2 - 6 ? )

Для стационарного режима

(t—>-оо) из

(2.66),

(2.67)

получаем

 

 

 

р , {to. * ) = Р . = - у >

Рг (to, 0 = р 2 =

-у.

(2.68)

Найдем теперь корреляционную функцию импульсно­ го марковского процесса с двумя состояниями, для чего

вычислим совместные

вероятности p(Xi, Xj, т)

состояний

Хі и х% разделенных

временным интервалом

и

 

Составим согласно

 

(2/.37)

систему уравнений для ве­

роятностей перехода

 

 

 

 

 

 

 

 

Щ і £

=

-

aPil

(*„, t) +

ßp*. (t0, t),

 

(2.69)

dPi*V°J)

=

aPil

(t0, t) - p P f t (t0, t).

 

(2.70)

Начальное условие

для нее задается

соотношением

(2.34). Решение системы

(2.69), (2.70) имеет вид

 

А, ('..') = -$-(!

- е - ^ - ' ° ' ) + е - ^ - ' ° \

(2.71)

Р-АК. t) =

~ { \

- е ^ < ' - « )

+ е - т

<'-'•»,

(2.72)

AJ^ . ,

0 =

^ ( 1 - е - 7

" - ' " ' ) ,

 

 

(2.73)

А, (^. 0 =

^ (

1 - е -

7 ( М о ) ) .

 

 

(2.74)

Поскольку при a=const, ß=const марковский про­ цесс однороден, то вероятности перехода (2.71)—>(2.74) зависят только от разности аргументов t—U=x. Для стационарного режима совместные вероятности р(хі, Xj, %) на основании (2.68), (2.71) — (2.74) записываются сле­ дующим образом:

р ( х 1 Д ^ ) = № І «

=

^

( І -

е

- 1

Н

{ ^ ,

(2.75)

р

хѵ

х) =

p l

P l 2

(*) =

І

( і _

е - т ^

(2.76)

P(x 2 ,x a ,x) =

p2 pä 2 (x) =

^

( l -

е

- ^ +

^ - е - ^ ,

(2.77)

р (*2,

X i , т) =р 2 р2і (т . ) = р (х ь

х2,

т).

(2.78)

61