Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 218
Скачиваний: 3
дЧврДх, i) |
|
|
d-\nwps{x, |
t) |
I /ô Intop,(x, I) |
|
||||
lu wps (X, t) = - |
|
- |
i - ln & - |
|
- |
i - In s2 (0 - |
j £ i |
^ ô ~ |
||
Тогда |
из (5.61) |
получим |
|
|
|
|
|
|||
- |
°3 |
(0 + |
|
|
m [i)f + ^ |
à |
. { x _ m {t)Y |
= |
|
|
|
|
= aa2 |
(^) — a,t [л — m (/)] -f- — -f- |
|
|
|||||
|
• |
AL |
[JC - |
m (Ol2 + |
32 |
(0 [25 {t) X - x~], |
(5.70) |
|||
|
|
•(0 |
где точка сверху означает производную по времени. Прежде всего отметим, что правая и левая части дей
ствительно являются многочленами типа а + Ьх + сх2. Это обстоятельство дает возможность, приравнивая коэффи циенты при одинаковых степенях х, получить дифферен циальные уравнения для математического ожидания m(t) и дисперсии а2(() апостериорной плотности вероят ности. Однако прежде чем выписать эти уравнения, обра тим внимание на то, что у нормального распределения мода совпадает с математическим ожиданием. Следова тельно, математическое ожидание в каждый момент времени совпадает с оценкой х* по критерию максималь ной апостериорной вероятности. Дисперсия апостериор ного распределения o2(t) служит оценкой качества филь трации. В дальнейшем параметры апостериорного рас пределения будут отмечаться звездочкой:
m(t) = x*, |
a2(t)=a*2. |
|
|
С учетом новых обозначений из |
(5.70) получаем |
||
X* = ах* + 4Іг X* - |
-£=r х* + |
^ С - Щ, |
(5.71) |
|
— |
ѵ О |
|
Начальными условиями для (5.71), (5.72) служат апри орные значения, определяемые стохастическим диффе ренциальным уравнением процесса х(і).
Для стационарного режима, когда a*2 (t)=0, из (5.72) можно найти среднюю дисперсию апостериорного рас-. пределения
а*2 = ; * 2 ^ - І У Л ^ + Л Ѵ Ѵ , |
(5.73) |
182
•Как видно, дифференциальное уравнение (5.71) яв
ляется линейным, что находится |
в полном |
соответствии |
|
с общими |
представлениями: фильтрация |
нормального |
|
процесса |
из другого нормального |
процесса |
осуществля |
ется линейным устройством. Линейная теория Колмого рова— Винера (§ 5.1) и теория условных марковских процессов приводят в данном случае к одинаковым ре зультатам.
Действительно, если учесть, что в соответствии
с (3.87) при /—>-оо
Nx = Aaa,
жх >
то (5.73) легко представляется в виде (5.9):
s* = - f ( Х - * ) .
где у определяется выражением (5.7).
Рассмотренный вариант решения задачи фильтрации является более простым по сравнению с методом реше ния соответствующего уравнения Винера — Хопфа. Кро ме того, теория условных марковских процеосов имеет еще ряд преимуществ, состоящих в том, что с ее помощью описывается нестационарный режим фильтрации и ре шается задача фильтрации процессов, характеристики которых меняются во времени [6].
5.5.Уравнения фильтрации марковского
сообщения из белого шума и марковской помехи
В § 5.2—5.4 рассматривался наиболее простой случай, когда £ (t) =х (t) + п (t). Этот вариант характерен для задач автоматического управления. В радиотехнике важна задача фильтрации в общей постановке (5.1), из ложению которой посвящен настоящий параграф. Ниже следующие выкладки хотя и более громоздки, но по ха рактеру весьма сходны с теми, которые приводились в §§ 5.2—5.4, поэтому комментарии к ним будут менее подробными.
Итак, на вход приемника поступает реализация (5.1)
W)=S(x,t) + V(i\,t)+n(t).
Необходимо оптимальным образом отфильтровать из этой смеси сообщение x(t).
13* |
183 |
Ё такой постановке задача была впервые решена Р. Л. Стратоновичем и Ю. Г. Сосулиным [13]. Однако авторы проанализировали случай, когда х(і), r\(t)— диффузионные процессы. Подобное ограничение не яв ляется принципиальным, и уравнения оптимального при ема в такой помеховой ситуации можно получить для
дискретных марковских |
процессов x(t), r\(t) [14]. |
|
|||||
Вывод уравнений фильтрации будем производить по |
|||||||
методике [5] *>. |
|
|
|
|
|
|
|
По аналогии |
с § 5.2 |
рассмотрим трехмерный маркое- . |
|||||
скин |
процесс [%(t), x(t), |
r\(t)\ который |
подвергнем |
дис |
|||
кретизации. Дискретный процесс [%т, хт, |
г\т] |
будем |
ха |
||||
рактеризовать |
начальной |
Wo(£,o, Хо, щ) |
и |
переходной |
|||
v(îh, |
Xh, r\k/i,h-i, |
Xh-u *\h-i) |
плотностями |
вероятностей. |
Повторяя выкладки § 5.2, получаем рекуррентное соот ношение, подобное (5.17):
Л S Wp, (хт, Ч[т, •»)
|
|
S |
2 |
S S |
|
|
A ' m + 1 1 m + 1 xm Xi |
||
v |
( S m - H ' Xm + 4 |
| Siw |
|
"4m) |
ffi'ps ( * m . |
f\m, ГЛ.) V ( S m + 1 . * m + l - |
" W l |
I I m . Хт,Ч\„, |
|
Располагая |
двумерной |
плотностью |
wvs{xm+i, |
т + 1), легко найти одномерное распределение для щения
(5.74)
цт+и
сооб
Wps{xm+V |
m + l ) = |
Yi |
Щв{хт+г. |
тіт+ 1 . m + 1 ) . |
|||||
|
|
|
1 ш + 1 |
|
|
|
|
|
|
Упростим |
обозначения: |
S(xJt')=S(x), |
V(-r\Jt) |
= |
V(r\); |
||||
Xh-i, |
Цк-і] |
= [І', x', |
11'], |
[lh, Xh, щ]=[1, |
x, |
T)] |
и заме |
||
тим, что двумерный процесс [x(it), |
r\{t)] |
характеризуется |
|||||||
априорным |
уравнением |
|
|
|
|
|
|
|
|
|
дщгІ»,-п, { ) = Ь |
р г { х 1 |
Шрг |
{ х > Чі |
t), |
|
(5.75) |
||
где Lpr(x, |
t]) — некоторый априорный |
оператор. |
|
||||||
По аналогии с (5.57) плотность вероятности |
перехода |
||||||||
из состояния [|' |
х', г\'] в |
состояние |
[£, |
х,-і\] |
за |
время Ai |
*' В [5] приведен вывод -уравнений для случая, когда f(0 =
=S(x, о+л(0-
184
в предположении независимости процессов x(t) nr\(t) за пишем следующим образом:
|
V(6,х.ц\Ѵ, |
|
X ' , |
т,') = |
[§,,, |
+ |
А а ( X ! х ' ) ] [ 8 ч Ѵ |
+ |
|||||||
+ Д*Я (TJ |
I т)')] |
|
|
ехр { |
|
|
[Ç - 5 (X) - |
V |
} . |
||||||
или в иной |
форме |
|
|
|
|
|
|
|
|
|
|
|
|||
X, |
л', |
7)') = |
/ ^ |
[ |
8 |
^ , |
8^, + |
Д Щ х , т,|х", |
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5.76) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л |
(X, 7, J X', |
у ) |
= |
ехр |
/ |
- [ |
Е |
|
- |
S |
(X) - |
V ( T , ) ] 2 J |
X |
||
|
X [ M ^ : | ^ ) 8 w + A ( 4 | V ) 8 ^ ] + - i r X • |
||||||||||||||
Х {ехр ( - ^ - [5 - |
Six) |
- |
V (T,)]2) - |
1} Ьхх, 8 ^ . |
(5.77) |
||||||||||
Подставим |
(5.76) в рекуррентное соотношение (5.74). |
||||||||||||||
Обозначив для краткости знаменатель через |
С, с учетом |
||||||||||||||
изменившихся |
обозначений |
получим |
|
|
|
|
|||||||||
|
|
|
|
|
wps(x, |
ц, t-\-M) |
= |
|
|
|
|||||
|
Е |
SWPS{Х''F , I ) V { L |
*'711%'!Х''Ѵ ) |
= |
|
||||||||||
|
|
X' |
г,' |
|
|
|
|
|
|
|
|
|
|
|
|
=~^Е YàV^h[5-' SII'+Л'Л{Х>711Х ''Ѵ ) 1ШР5 |
^0 ^ |
||||||||||||||
X' |
-Г)' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ Ä T S SЛ{Х'711 |
Y )% S |
|
^ |
|
(5'78) |
|||||||||
|
|
|
X' |
у |
|
|
|
|
|
|
|
|
|
|
|
Нетрудно заметить [см. (5.74)], что для вычисления |
|||||||||||||||
знаменателя |
С следует |
числитель |
(5.74) |
просуммировать |
|||||||||||
по X, ц, т. е. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
1) |
|
|
|
|
|
|
|
|
|
|
|
+ Л ' Е Е Л {Х~711 |
|
Ѵ ) |
^ |
|
|
^ / ) ] = |
|
|
||||||
|
|
X' |
ц> |
|
|
|
|
|
|
|
|
|
|
|
|
18§
в |
Разложив величину |
С - 1 |
в |
ряд и подставив |
результат |
||||||||||||||||
(5.78), лолучим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
wvs |
{х, |
т,, t-{-àt)=wps |
{х, |
-т], |
О+ЛгЕ ЦЛ {х,і) 1 X',t)') |
Wps |
{X', |
7)', |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Л ' |
1)' |
|
|
|
|
|
|
|
|
О - а>Р 8 (-*Л1.')А* |
S |
S |
Л (л, Tj|Jc',^')>pe (Jf', |
т)', 9. |
(5.79) |
||||||||||||||||
В соотношении (5.79) опущено слагаемое, имеющее по |
|||||||||||||||||||||
рядок |
малости (АО2 - |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Перенесем |
wps(x,r\,t) |
|
|
влево, |
разделим |
обе |
части |
|||||||||||||
уравнения |
(5.79) |
на |
At |
|
и |
перейдем |
к |
пределу |
при |
||||||||||||
At—ѵО: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
l i m |
wps |
{x, |
У), |
t+At) |
|
— wvs |
(x, |
7), |
/ ) |
|
|
|
|
|
||||
|
|
= J ™ |
J ] |
J ] A |
(•*• |
|
ч |
I |
|
ч') ° Ѵ (*'. "4', 0 |
- |
|
|
||||||||
|
- Шпшр, (x, |
TJ, 0 J |
] |
J |
] |
Л (x, |
7j |
I x', |
7)') W p |
a |
(x', |
n',t). |
(5.80) |
||||||||
|
Предел |
величины |
Л(х, ті|х', г\') |
находится |
просто |
и |
|||||||||||||||
выражается соотношением |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
lira Л (л-, т) |
I x', |
у ) = |
Я (л-1 |
8 |
|
+ |
|
|
|
|||||||||
|
+ я (ѴѴ) 8Д,. - |
|
|
[5 - |
|
S (x) - |
V (г,)]3 Ьхх, |
8 ч ѵ . |
(5.81) |
||||||||||||
Используя |
(5.81), |
преобразуем |
формулу |
(5.80) к |
виду |
|
|||||||||||||||
|
|
^ ( * . ^ 0 а = 2 ] я ( ^ | ^ > р . ( ^ , , , / ) + |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
Л:' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
wps |
(х.-ц.і) |
£ j |
^ |
Л |
|
У |
' |
|
0 |
1 5 - 5 |
(x) |
- |
V (i)f. |
(5.82) |
|||||
|
|
|
|
|
* |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При записи (5.82) учтено, что
£ а (*1 * . ' ) = £ Мч h ' ) = Q.
*1
186