Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 247
Скачиваний: 3
Общее выражение для корреляционной функции
|
kx (х) == S XiXjP |
{Xi, Xj, t) — |
XiPi^j |
||||
при Xi=x, |
xz=x |
+ a |
с |
помощью |
соотношений (2.68), |
||
(2.75) — (2.78) |
приводится |
к виду |
|
|
|||
|
|
М , ) |
= |
^ е - ^ М . |
(2-79) |
||
Когда |
ХІ = —а, |
хг = а, |
вместо |
(2.79) |
получаем |
||
|
|
М |
х ) |
= |
! ^ е - т |
м . |
(2.80) |
Используя соотношение Винера — Хинчина, находим выражение для спектральной плотности мощности процесса x(t):
S |
X ( , ) = § M , , e ^ x |
= ^ l |
i r L _ . . |
(2.81) |
|
—оо |
|
|
|
Марковский процесс с двумя состояниями |
является |
|||
весьма |
распространенным |
видом |
случайного |
процесса, |
и на практике часто возникает необходимость его мо делирования. Устройство для генерирования такого про цесса весьма просто .можно построить, используя два ра диоактивные элемента в качестве источников пуассоновских случайных потоков [21]. Излучение от одного из ра диоактивных элементов воздействует на специальный де тектор (например, счетчик Гейгера), и импульсы . вы хода последного запускают триггер. Возвращение триг гера в исходное состояние осуществляется импульсами от другого детектора, реагирующего на излучение от второго радиоактивного элемента. Таким образом на триггере создается напряжение, являющееся марковским процессом с двумя состояниями.
В частном случае, |
когда a = ß , |
марковский процесс |
||
с двумя состояниями превращается |
в широко известный |
|||
в радиотехнике т е л е г р а ф н ы й сигнал |
(см., например, |
|||
[3,17]), |
характеристики |
которого |
легко |
находятся из |
(2.80), |
(2.81): |
|
|
|
|
^ ( т ^ а ^ - 2 " 1 * 1 |
, |
(2.82) |
62
Следует подчеркнуть, что выражение (2.82) впервые по лучено С Райеом [22] на основе использования частной закономерности, характерной лишь .для телеграфного сигнала. Применение метода (22] к более общему про цессу (когда «т^р) затруднительно.
Импульсные марковские процессы 'могут служить ма тематической моделью хаотических импульсных помех (ХИН), если при задании пуассоновских потоков поло жить \iî(t)>%u(t)- (і=2, 3, ...,N). Для марковского процесса с двумя состояниями это условие сводится •к выполнению неравенства ß^>a.
Трудности в теории синтеза радиоприемных устройств, находящихся под воздействием белого шума и импульс ных помех объясняются, в частности, сложностью суще ствующих математических моделей импульсных помех* Действительно, даже в простейшем случае необходимо задавать порознь распределения трех основных пара метров:, амплитуды, длительности и времени появления импульсов. В марковской модели задаются только интен сивности пауссоновских потоков. В ряде частных, но важных случаев, число задаваемых параметров может быть сведено к 'минимуму. Так, можно положить <Кц=а, Ati="ß для всех і=2, 3, ..., N. Преимущества марковской модели ХИПособенно проявляются в связи с возмож ностью применения теории условных марковских процес сов к задачам оптимального обнаружения и фильтра: ции сигналов [23] (см. также гл. 5).
t
2.6.Разрывные марковские процессы с непрерывным
множеством состояний
Рассмотрим теперь случай, когда разрывный марков ский процесс имеет непрерывное множество состояний. Одна из возможных реализаций такого процесса изобра жена на рис. 2.3. Поскольку множество состояний не счетно (образует континуум), то интенсивности потоков тина lij(t) задать нельзя и, следовательно, нельзя вос пользоваться формой представления вероятностей пере хода в виде (2.32).
Для описания разрывного марковского процесса с не прерывным множеством состояний воспользуемся мето дикой, •использованной при выводе уравнения (2.28), обобщив ее на непрерывный случай. Для каждого со стояния X определим интенсивность пуассоновского по-
63
тока Х(х, t). Вместо |
матрицы |
Q(i) |
относительных ве |
||||||||
роятностей |
перехода |
<7ij(0 |
зададим |
д в у м е р н у ю |
|||||||
п л о т н о с т ь |
Q(x', |
X, t), |
с помощью которой |
вероят |
|||||||
ность перехода |
из состояния х |
на участок |
(х, |
х+Ах) |
|||||||
п р и у с л о в и и , |
ч т о с к а ч о к п р о и с х о д и т , |
выра |
|||||||||
жается ігзесьма просто |
|
|
|
|
|
|
|
||||
|
|
|
|
|
Q(x', |
X, t)Ax. |
|
• |
|
(2.84) |
|
Функция |
Q (х', |
x, t) |
подчиняется |
условию |
нормировки |
||||||
оо |
|
|
|
|
|
|
|
|
|
|
|
J" Q(x', |
X, t)dx |
— 1. |
|
|
|
|
|
|
|
—оо
Аналогично (2.23), (2.24) запишем вероятности от
сутствия и наличия скачка на интервале (t, |
t+At): |
||
Рі,и(°)-Л |
-l(x',t)M, |
(2.85) |
|
Pt,uW~l(x', |
t)At |
(2.86) |
Тогда с учетом (2.84) — (2.86) вероятности перехода раз-
„рьтвного процесса x(t) из состояния х' за время At вы разятся соотношениями
рх,х, |
(t,t |
+ |
At) = l~X |
(х\ t) At, |
(2.87.) |
|||
Px,,x+Ax{t, |
t]+àt) |
= l{x', |
t)AtQ{x', |
x, t)Ax. |
(2.88) |
|||
Обозначим одномерную плотность распределения со |
||||||||
стояний процесса .через |
w(x, |
t). |
Тогда |
w(x', |
t)Ax' |
есть |
||
вероятность того, что процесс x(t) |
в момент |
t находится |
||||||
на участке (х', |
х'+Ах'). |
Вероятность |
того, |
что процесс |
||||
в момент і находится на участке |
(х, х + Ах) |
и за |
время |
|||||
А^ останется там же, равна |
|
|
|
|
|
|||
|
w(x, |
t)Axll—%(x, |
t)Af\. |
|
(2.89) |
Вероятность того, что процесс в момент t находится на участке (х', х'+Ах') и за время At перейдет на участок (х, х+Ах) определяется соотношением
w(x', t)Ax'h{x', t)AtQ{x\ x, t)Ax. (2.90)
Используя (2.89), (2.90), получим выражение для полной вероятности того, что процесс в момент t+At окажется в состоянии х:
w(x, t-\-At)Ax = w(x, t)Ax[\ —1(х, t)At\-\~
+ S w(x\ t)X{x', t)\àtQ(x\ x, t) AxAx'. |
(2.91) |
46
Сократим (2.91) на Ах, перенесем в левую часть равен ства w(x, і) и .разделим обе части уравнения на At. В результате предельного перехода при Д ^ — У О , АХ'—>0 вместо (2.91) имеем
+ J w (л', I) Я (х\ t) Q (х', x, i) dx! |
(2.92) • |
— o p
Интегро-дифференциальное уравнение (2.92) является основным для разрывных процессов с непрерывным мно жеством состояний. Оно называется уравнением Колмо горова (см., например, [24—26J).
В том случае, когда
|
%(X, |
t) =X(x), |
Q{х', |
X, t)=Q{х', |
X), |
|
|
||
процесс |
становится о д н о р о д н ы м |
(вероятности |
пере |
||||||
хода не зависят от времени) и, следовательно, |
|
|
|||||||
dw (х, |
t) |
|
|
ou |
|
|
|
|
|
|
|
t) + j |
|
|
|
|
|
||
dt |
= |
- |
Я (л-) w (х, |
w {x', |
t) Я (JC) Q (*', |
x) |
dx'. |
||
|
|
|
|
— o o |
|
|
|
|
(2.93) |
|
|
|
|
|
|
|
|
|
|
Решение |
уравнений |
(2.92), |
(2.93) |
даже |
при |
À(x)=À |
в общем виде не получено. При необходимости для ре шения следует пользоваться приближенными методами или рассматривать частные случаи.
В качестве примера построим модель импульсных по мех на основе марковского разрывного процесса с не
прерывным 'множеством состояний. |
|
|
|
|||||
Предположим, что х^О, |
и выберем в качестве опор |
|||||||
ного |
нулевой |
уровень х = 0. |
Двумерную |
функцию |
Q{x', |
|||
x, t) |
нужно |
задать |
специальным |
образом, |
отразив то |
|||
обстоятельство, что |
в состояние |
x(t)=£0 |
за |
время |
Al |
процесс может попасть только из состояния х=0 и не из
какого |
другого; вместе с тем из любого состояния |
x(t) |
||
процесс |
может перейти только в состояние |
х = 0. |
Все |
|
другие |
переходы исключены. Подобный характер функ |
|||
ции |
Q(x', x, t) предопределяет к о н е ч н у ю вероятность |
|||
р ( 0, |
t) |
того, что процесс окажется в нулевом |
состоянии. |
|
Для удобства обозначим интенсивность потока, пере |
||||
водящего процесс из нулевого уровня через |
a(t): |
К (О, |
||
t)=a(t), |
сохранив для остальных состояний |
(x(t)^0) |
5—186 |
65 |
прежнее обозначение интенсивностей ,Х(х, t). Тогда для плотности вероятностей w(x, t) (х>0) и вероятности р(0, t) можно записать следующую систему уравнений:
w{x, |
t + Atf) = w {x, |
t) [1 —X (x, t) M] |
+ |
|
|
'+/7(0, t)a(t)AtQ(0, |
x, t), |
(2.94) |
|
p(0, |
t + bt)=p{0, |
i)[\—a(t)<tâ\ + |
|
|
|
+'At$w(x, |
t)X(x, |
t)dx. |
(2.95) |
После обычных преобразований вместо (2.94), (2.95) получаем
dw(x. t) |
~ - |
X (x, t) w (x, t) + |
p (0, t) a (t) Q (0, x, t); |
(2.96) |
dt |
||||
dp(0, |
t) = |
— a{t)p (0, /) + |
j w (x, t) X (x, t) dx. |
(2.97) |
dt |
|
|
|
|
Очевидно, что для всякого t имеет место условие 'нор мировки
|
|
|
p{Q,t)+ |
lw{x,i)dx=\. |
|
|
|
|
(2.98) |
|||
|
Как |
видно, в систему |
(2.96), |
(2.97) |
естественно вхо |
|||||||
дят |
пуассоновское |
распределение моментов |
появления |
|||||||||
ш*(хл |
|
|
|
|
импульсов |
с |
интенсивностью |
|||||
ф |
* |
|
|
|
|
a(t), |
экспоненциальное рас |
|||||
|
|
|
|
|
|
пределение |
|
длительностей |
||||
|
|
|
|
|
|
помеховых импульсов с пара- |
||||||
|
|
|
|
|
|
метром Х(х, t) и распределе |
||||||
|
|
|
|
|
|
ние Q(0, |
x, |
t)=Q(x, |
|
t). Сле |
||
|
|
|
|
|
|
дует подчеркнуть, что в при |
||||||
|
|
|
|
|
|
кладной |
литературе |
иногда |
||||
|
|
|
|
|
|
под |
распределением |
ампли |
||||
|
|
|
|
|
|
туд |
импульсных |
последова |
||||
|
|
|
Рис. |
2.10. |
|
тельностей |
неверно |
понима |
||||
|
|
|
|
|
|
ется |
функция |
Q(x, t)., кото |
||||
рая «а |
самом деле задает у с л о в н о е |
'распределение ве- - |
||||||||||
роятностей |
п о я в л е н и я |
импульсов |
с |
амплитудой х. |
||||||||
Безусловное |
распределение |
амплитуд |
последовательно |
сти неперекрывающихся случайных импульсов характе
ризуется функциями w(x,t) и ô(x) |
с весом р(0, t). |
||||
Если интересоваться |
стационарным случаем, когда |
||||
a{t)-=a;b(x, |
0 = 4 * ) ; |
Q(*. |
Q=Q'(*)i |
||
dw(x,t) |
dp{0, |
t) |
= |
0, |
|
|
di |
dt |
|
|
|
«6 |
|
|
|
|
|
то 'из (2.96), (2.98) можно .получить соотношения
1-1
Эти соотношения полностью определяют финальное распределение w$(x) процесса x(t):
О»Ф(-АГ) = [ |
1 + А J ^ - |
^ J |
- 1 p W + x T ^ - Q W ] - ( 2 - 9 9 ) |
При Я(лс) = |
Я выражение |
(2.99) упрощается |
|
|
Щ (х) : |
X |
•8(jc)- |
Вид финального распределения [показан на рис. 2.10.
3
НЕПРЕРЫВНЫЕ МАРКОВСКИЕ ПРОЦЕССЫ
3.1. |
Вводные замечания |
|
|||
Рассмотрим основные характеристики непрерывных марковских |
|||||
процессов, основываясь |
на общих |
представлениях |
о непрерывных |
||
случайных процессах. |
|
наиболее |
полно |
описывается многомер |
|
Непрерывный процесс |
|||||
ной (n-мерной) плотностью |
вероятности wn(x\, |
х2, |
h; ...; xn, tn), |
которая определяет вероятность того, что значения случайной функ
ции x(t) в моменты времени U, і2 |
f„ заключены соответственно |
|||||
винтер-валах (хі, Xi+Axi), |
(х2, х2+Ах2), |
|
i(xn, xn+Axn). |
При ма |
||
лых AX{ эта вероятность равна wn (xi it; x2, t2\ ...; xn, |
tn)dxi |
... |
dxn. |
|||
Многомерная плотность вероятности |
wn(xi, |
i\\ ...; xn, |
t„) |
дает |
воз |
|
можность' судить о связи |
между значениями случайной |
функции |
||||
в я моментов времени и характеризует случайный процесс тем де |
||||||
тальнее, чем больше число п. |
|
|
|
|
|
|
По правилу умножения вероятностей |
зависимых |
событий |
|
и>п(*ь |
U; |
хп, |
tnldXidXi, ... . dxn |
= |
|
|
= Wn-i(xi, ti\ |
x2, |
t2\ . . . ; |
X n - i , / n - i ) ^ i . . . |
dXn-іУ. |
|
|
Xv„(x„, |
tn\x,, <ti; xn-i, |
tn-i)dxn. |
' |
(3.1) |
||
Величина vn\xu, t n \x u |
t\\ |
...; xn-i, tn-i)dxn |
есть |
условная |
вероят |
|
ность того, что значение случайного процесса в момент tn' |
окажется |
|||||
5* • |
|
- |
|
|
67 |
|
t