Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 219
Скачиваний: 3
ми); если же заявка застает свободным хотя бы один
канал, то она обслуживается любым из |
свободных |
каналов. |
|
' П о д со с то я н и е м с и с т е м ы Хк будем |
понимать |
ч и с л о з а н я т ы х о б с л у ж и в а н и е м к а н а л о в не-
заівисимо от |
их |
порядковых номеров. |
В |
таком |
случае |
||
состояние х0 |
означает, что все каналы |
свободны; состоя |
|||||
ние Хи—-занято |
k |
каналов. В состоянии |
Xh (k=l, |
2, |
|||
Я |
|
|
|
|
|
|
|
|
|
|
Чкф |
miß |
Wf? |
|
|
|
|
|
Рис. 2.8. |
|
|
|
|
N—1) на систему воздействует два потока: поток |
заявок |
||||||
с интенсивностью |
%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