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

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

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

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

Добавлен: 28.06.2024

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

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

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

ней дисперсий апостериорного распределения вероятно­ стей. Процесс улучшения качества фильтрации прекра­

щается

с

того момента,

когда

апостериорные

сведения

в силу

действия белого

шума

уже не

могут

уточнить

априорных

сведений,

выработанных

к предыдущему

шагу. Иными словами, в стационарном режиме устанав­ ливается равновесие между двумя противоборствующими факторами — вырабатываемыми априорными сведения­ ми, которые расширяют апостериорное распределение вследствие априори возможных переходов процесса, и извлекаемой из реализации %{kAt) информации, которая уменьшает дисперсию апостериорной плотности вероят­ ности.

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

Перейдем

теперь

от

рекуррентных

соотношений

(5.35) — (5.37)

к

дифференциальным

уравнениям

филь­

трации.

 

 

 

 

 

 

 

 

 

Преобразуем выражение (5.35). С учетом того, что

ехр(—zAt) = l—zAt+...,

 

пренебрегая

членами

второго

порядка малости, получаем

 

 

 

 

 

W p s (0, / +

Al) =

У

 

i j j g y l [ W p s

(0, f)

-

I

 

 

 

 

 

 

 

 

 

Сомножитель

 

 

соотношения

(5,45) также уп-

ростам, используя

(5.37)

и

разложение

(1 —

zAt)'1^

= l + z A * + ...

 

 

 

 

 

 

 

 

 

w (lm+l)

 

 

L

Ü2 (t + AQ

At-

 

 

 

 

r,N0

^

Л',

 

 

-^rwps([,t)At№(t

 

+

At)-l)

(5.46)

173


Подставляя (5.46)

в

(5.45),

имеем

 

wpe{0,t

+

M) =

wpa{0.t)(l-àÀt)-\-

 

+twp, (1,0 р д / - w p

s (0, t) wps

( 1, t) -àjjt \2\{t-\-M) — 1].

(5.47)

То обстоятельство,

что

 

 

 

ps ( i , 0 +

Œ)p ,(0,o = i .

(5.48)

позволяет, во-первых, исключить из (5.47) шР і ,(1, 0 и, во-вторых, дает возможность не рассматривать уравне­ ния (5.36).

Если перенести wps(Q, t) в левую часть уравнения (5.47), разделить обе части уравнения на At и затем перейти к пределу, то из (5.47) следует н е л и н е й н о е дифференциальное уравнение для апостериорной плотно­ сти вероятности

dw-, (0.0

_

•atüpe (0.0 + ß[l

-wps(0,t)]

dt

 

 

 

 

No •wps(0,t)[l-wps(0,t)][2t(t)-l].

(5.49)

Схема оптимального приемника фильтрации, реали­ зующая алгоритм (5.49) и дополненная пороговым устройством, изображена на рис. 5.1. На вход приемника

 

 

-а.

 

 

 

1

<

 

 

 

 

Wps(Q,t)

 

 

No

 

Реше­

 

 

 

>

> 3HZ

Порога

Зое уст.

ние

X

Рис. 5.1.

поступает реализация £(0, которая обрабатывается в со­ ответствии с правилом (5.49). В результате обработки [интегрирования дифференциального уравнения (5.49)] на вход 'порогового устройства подается вычисленное зна­ чение апостериорной плотности вероятности wps(0, t),

174


Которое

сравнивается

 

с

порогом. На

основе

еравиёния

выносится решение о приеме того или

иного

значения

x(t).

Важно

подчеркнуть,

что как вычисление

плотности

вероятности, так и сравнение в пороговом

устройстве

осуществляется

н е п р е р ы в н о

во

времени.

Благодаря

этому получаемое оценочное значение переданной реали­

зации также непрерывно во времени.

 

 

 

 

 

 

 

 

 

Специально отметим, что точное аналитическое

реше­

ние уравнений нелинейной фильтрации, не считая вы­

рожденного лнненногослучая

(см. § 5.4), возможно

край­

не

редко. В

частности,

 

известно

[6] решение

уравнения

(5.49).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Основная ценность уравнений нелинейной фильтра­

ции состоит в том, что по ним составляются схемы

опти­

мальных приемных устройств. В условиях, когда спра­

ведливо

гауссово

приближение

(см. § 5.6),

эти

урав­

нения позволяют также оценить качество фильтрации.

 

 

Физический

смысл

членов

уравнения

(5.49)

сводится

к следующему [8]. При большом

шуме извлечь

информа­

цию из реализации £(і) практически невозможно. Дейст­

вительно, при N0-+oo уравнение (5.4-9)

переходит в (5.'20),

так что первые два члена уравнения

(5.49)

выражают

изменение

апостериорной

вероятности

в

соответствии

с изменением априорных

сведений.

 

 

 

 

 

 

 

 

 

Если a=iß 0,

то

в

правой

части

уравнения

(5.49)

остается только третий

член

 

 

 

 

 

 

 

 

 

 

 

(О,

i) =

_

 

( 0 . t )

[

l

_

^

( 0 | {

) ] _ 1 _

 

{ t ) _ l

^

( 5 5

0 )

 

 

Этот случай соответствует отсутствию переходов из

нуля в единицу и наоборот, т. е. на

интервале

времени

[О,

t] передается

либо

 

x(t)==\,

 

либо

x(t)=0.

 

Следова­

тельно,

 

 

 

wp.(0,t)=w(0\Ut)),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Wps(0\l(t))

 

определяется

формулой

(5.38).

 

 

 

Как известно из теории оптимального

радиоприема

(см., например, [3]),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ш (0 | Ü (0) =

kwvr

(0) exp j -

-

i -

j

V (f) dt' J,

(5.51)

w(l

I б (0) =

Ät»p r

(І) exp j

j

"

 

(^) — Ifdt'

J.

(5.52)

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

175,


Ёеря

отношение

апостериорных

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

(5.51),

(5.52), приходим к алгоритму

 

 

 

 

 

 

 

«(оIi(0)

=

^ ѵ В

 

(

 

L

ГШ

п

-

№'

\ $

1

 

 

 

 

 

 

 

 

 

à

 

 

 

'

(5.53)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Этот же самый

результат

получается

и

из

уравнения

(5.50). Действительно, умножив (5.50)

на (1—ayp s (0, / ) ] ~ 2

и обозначив

 

 

 

(о, о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-Q(t),

 

 

 

 

 

 

 

 

[І-а-р, (0, О]2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вместо

(5.50)

имеем

 

 

 

 

 

 

 

 

 

 

 

 

dQ (t)

_

- О ( 0 т ^ [ 2 £ ( 0 - 1 ] .

 

 

(5.54)

 

 

d t "

 

 

Решение (5.54)

с начальным

условием

 

 

 

 

 

 

 

 

 

д ( 0

) = В!гЖ

 

 

 

 

 

 

 

 

 

 

ѵ

'

а>рг([)

 

 

 

 

 

полностью совпадает с (5.53).

 

 

 

 

 

 

 

Таким образом, рассмотрение частного случая a = ß =

= 0 оказалось

полезным

с двух точек зрения: во-первых,

•выяснено, что третий

член

правой

 

части

уравнения

(5.49)

отражает

влияние

на

апостериорную вероятность

наблюдаемых

данных,

и, во-вторых,

установлена

связь

нелинейной фильтрации с линейной фильтрацией детер­

минированного

сигнала конечной длительности.

 

Отметим, что несмотря на то, что в данном

случае

S(x,

t)=x(t),

уравнение (5.49) для

плотности

wps(0,t)

оказалось -нелинейным. Это обстоятельство

объясняется

иегауссовскнм

характером

процесса

x(t).

Если

x(t) —

нормальный процесс, то нелинейная фильтрация

вырож­

дается в линейную, осуществляемую

по критерию мини­

мума

среднеквадратической

ошибки

(см.

следующий

параграф).

5.4.Фильтрация марковского процесса

сдискретными состояниями и гауссовского марковского процесса в белом шуме

Пусть по-прежнему

S(x, t)=x(t),

однако

теперь x(t)

представляет собой разрывный марковский процесс с N

состояниями, который

рассматривался в §

2.3. Вероят-

176


ность

перехода

из состояния ХІ

в'состояние

х-}

{\Ф'\)

за

малое время Ы равна

 

 

 

 

 

 

 

 

 

Pij(t,

t+'M)=öi}+bi]M.

 

 

(5.55)

Матрица

перехода

для такого

процесса

записывается

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

P = (Pii(t,t + M)) =

 

 

 

 

 

 

1 +

XwAt

(\0,d^ . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.56)

На входе приемника наблюдается реализация

%(t) =

—x(it)

+n(t).

Рассмотрим

двумерный марковский про­

цесс

ll(t),

 

x(t)],

выражение для плотности вероятности

перехода

которого [8]

 

 

,

'

 

 

ѵ(%,

ЛЙИЙ - , .

^ f e - i ) =

( 8 f t - i . f t + A ^ f t - i , f e ) X

 

 

 

 

x / ^ r e x p f -

 

 

 

 

(5.57)

можно получить, если повторить все рассуждения, ис­ пользованные при выводе формулы (5.34). Вновь вос­ пользуемся рекуррентным соотношением (5.17), для чего сначала выпишем и преобразуем его числитель

^аУрв (ХІ, t) V (^, Xj I le, ХІ) =

i

 

 

= S » , . ( X i . №

+ AH«)

e x p [ - ( g j - g ) 2 A / ] -

= Wps(Xj,

t) j /

^ [ l

i

 

 

(5.58)

Здесь учтено (5.57), экспонента разложена в ряд и от­ брошены члены, имеющие порядок (М)2 и выше.

12—186

177