Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.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