Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи

.pdf
Скачиваний:
21
Добавлен:
23.10.2023
Размер:
7.76 Mб
Скачать

Обозначим для

краткости

 

К

1 + 2 aN0

КІ = 2К*ХХЩ

 

2/1 „Я*

 

4 3 —

ЛГ,

 

В этих обозначениях уравнения (5.131), (5.132) прини •мают вид

X* = - К,л* + Ка

|^(0 cos К ( + T * ) -

2

, (5.133)

9 = _ (/С, +

Л>*) 5 (0 sin к * +

«р*).

(5.134)

Алгоритм (5.133), (5.134) реализуется схемой, представ­ ленной на рис. 5.5. На этой схеме ПГ подстраиваемый

генератор гармонических

колебаний

частоты œo; УЭ —

управляющий

элемент. Оптимальный

приемник состоит

 

 

 

-к,

tftj

X

 

X*(tl

 

 

 

Cosfûj0t

+ (fi*)

2

 

 

 

 

 

1

 

V

 

 

 

X

ПГ

УЭ

 

X

Рис. 5.5.

из двух каналов. Один из них основной, на выходе кото­ рого вырабатывается оценочное значение x*(t). Этот ка­ нал моделирует уравнение (5.133). Второй канал, пост­ роенный в соответствии с (5.134), представляет собой систему фазовой автоподстройки частоты (ФАПЧ), ко-

14-186

201

іорая выраоатывает опорным сигнал на основе приняв. . •ііизации l(t) и оценочного значения x*(t).

2. Фазовая модуляция. Сигнал при фазовой мо , ::..ч,ян записывается в виде

 

 

 

S(x,

0 =

AcosK^-|-0(0],

(5.135)

 

 

в (0 =

» (') +

<Р(0 = ЛГфл: (0 +

(t),

(5.136)

где ' Мф

постоянный

коэффициент,

характеризующий

глубину

ФМ. Случайная

фаза

и

полезное

сообще­

ние x(t)

подчиняются

 

априорным

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

уравнениям:

 

 

 

 

 

 

 

О == - аМфх +

Мфпх

(t) -f. пч (t); х =

-

ах + пх

(t).(5.137)

 

 

 

 

 

 

 

 

-CL

 

 

 

 

_2

 

 

Кх*В

 

 

 

№1

 

 

No

 

 

 

 

x*(t)

X

 

>

 

 

>

 

 

 

 

 

 

 

 

 

•2Мф

А0 Sin(u0t

+ 8

}

S

пг

УЗ

 

 

 

 

Рис. 5.6.

 

Основываясь на

(5.135) — (5.137)

и производя выклад­

ки в той последовательности, в которой они проводились •при анализе амплитудной модуляции, можно получить следующие уравнения фильтрации

определяющие структурную схему оптимального прием­ ника (рис. -5.6). Приемник представляет собой автома­ тическую систему, осуществляющую слежение за фазой принимаемого сигнала.

202

Относительная среднеквадратігческая ошибка филь­ трации при ФМ оказывается равной

х[Vi1+v^f+4?з»~(і+ѵщ)}>

 

 

 

 

 

(5-139)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

•^о

 

2

 

д,2

2

Л>2

ЛГ*

р,

Л

~~ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Величина

D

есть

дисперсия

случайного

набега

фазы

за

время

корреляции

сообщения

% = 1 / а .

 

 

 

 

Результаты

 

 

вычислений

 

 

 

 

 

 

 

по формуле (5:139) для не­

 

 

 

 

 

7—™——.

скольких

значений

аа

 

 

 

 

 

 

 

 

представлены нарнс.5.7—5.9.

 

ft 7

 

 

 

 

 

Ход зависимостей

указывает

 

 

 

 

 

 

на то,

что

лри

заданных q

 

 

 

 

 

 

 

и аь ошибка

фильтрации уве­

 

 

 

 

 

\

i

личивается

с

ростом D^,

по­

 

 

 

 

 

этому

следует

 

добиваться

 

 

 

5

10

30

50c

более

высокой

 

стабильности

 

 

 

 

 

 

 

генератора передатчика. При

 

 

 

Рис.

5.7.

 

 

фиксированных

 

 

q

и

D

 

 

 

 

 

 

 

ошибка

уменьшается

с

увеличением oft, однако

необхо­

димо помнить,

-что

рост

влечет

за собой

расширение

спектра

'радиосигнала.

 

 

 

 

 

 

 

 

 

3.

Частотная

модуляция.

При

частотной

модуляции

полезный радиосигнал представлен в виде

 

 

 

 

 

 

 

 

S(x,

f) =

A,cosKf_|_<j.(f)],

 

(5.140)

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

где ф (/)= ср (/) 4_ Мч

j x

(t>) dt',

<j> =

+

д 7 (/).

 

о

Опуская выкладки, приведем окончательные выра­ жения, определяющие структуру оптимального прием­ ника и качество его работы

= % = - A T ( 0 sin к * + 1 * ) ;

No

H*

203

 

4M'

 

 

 

 

 

 

- / 1

+ ^ + 4 / ^ ( Р І л , + - 2 - ^

 

X

 

X j / H - 2 ^ +

4 j / ^

p V b T ^

)

- (5.142)

 

Al

Р™

— .

Величина

ß V M

характе-

Г Д 6 ? =

2 ^ Ѵ

ризует

индекс

частотной

модуляции.

 

 

^7

 

~JDy>=1-

 

 

 

 

0,5

 

 

 

 

 

 

0,3

 

 

« 3

 

 

 

 

 

 

 

 

 

L

0,1

 

 

a 2

 

 

 

6>=J

 

 

 

 

 

0,1

 

 

a*

V

 

 

 

 

\

 

 

 

 

 

 

 

 

Рис. 5.8.

 

Рис. 5.9.

 

Структурная схема оптимального фильтрующего уст­ ройства, .составленная по уравнениям (5.141), изображе­ на на рис. 5.10. Основная особенность приемника 4 M сигналов, равно как двух других приемников, схемы ко­ торых приведены на рис. 5.5, 5.6, состоит в реализации •квазикогерентной обработки принятого колебания. Ины­ ми словами, оптимальная фильтрация сообщения х(і) осуществляется в том случае, когда, из колебания извле­ кается максимум информации, в том числе и информация о «несущественном» параметре <ç(t). Вычисление в каж­ дый момент значений фазы сигнала позволяет подстраи­ вать колебание опорного генератора с таким расчетом, чтобы получалась когерентная обработка принятой реа-

204

лизации ^(t). Постоянно присутствующий шум п(і) и случайный характер сообщения x(t) не дают возможно­ сти вычислить точное значение фазы сигнала, в резуль­ тате чего оптимальное устройство осуществляет квази­ когерентный прием.

Качество работы оптимального приемника 4 M сиг­ налов характеризуют графики зависимостей, рассчитан­

ные по формуле

(5.142). Они представлены на рис. 5.11.

 

_2_

 

 

 

'No

 

fftj

X

>

>

 

 

 

 

Mit

 

 

фф

Ф *

 

 

 

 

 

 

f

пг

УЭ

 

 

 

Рис.

5.10.

 

Ход зависимостей свидетельствует о том, что при

фиксированных

значениях

q и ß V M

ошибка

фильтрации

минимальна, когда D

=0.

При

заданных

q и

вели­

чина 8цМ тем

меньше,

чем больше $ ч м . Однако при увели­

чении индекса модуляции ßt / y M расширяется спектр сигнала Это обстоятельство следует иметь в виду при решении конкретных задач.

В работах Р. Л. Стратоновича, Н. К. Кульмана, В. И. Тихонова, Ю. В. Саютина и других авторов [15— 24] решены многочисленные задачи по синтезу разнооб­ разных систем связи с непрерывными видами модуля­ ции. Помехоустойчивость импульсных систем передачи непрерывных сообщений исследована в работахМ. С. Яр­

лыков а [25—27].

Это оказалось

возможным

благодаря

использованию

аппарата

теории

условных

марковских

процессов. Рассмотренные

выше

примеры

показывают

205

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

0,7

0,5

0,3

0,2

0,1

0,07

0,05

0,03

*

7 3

5 7 to

20 30 50 70 q

Рис. 5.11.

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

5.8.Оптимальное обнаружение марковских сигналов на фоне белого шума и марковских помех

при дискретном наблюдении

Аппарат теории условных марковских процессов по­ мимо задач фильтрации позволил решить широкий круг задач о б н а р у ж е н и я случайных и детерминированных сигналов на фоне белого шума и разнообразных типов помех. Основополагающие работы по этой тематике при­ надлежат Р. Л. Стратоновичу « Ю. Г. Сосулину [13, 28—32]. Задача любого оптимального приемника заклю­ чается в извлечении из принятой реализации %(t) макси­ мума информации об интерсующем параметре сигнала. При обнаружении сигнала эта информация заключен?

206

в апостериорных вероятностях того, что на интервале наблюдения сигнал есть или сигнала нет. Введем тараметр наличия сигнала

1, когда сигнал есть, О, когда сигнала нет.

Будем полагать, что сигнал S(x,t)

зависит

от одного

случайного марковского сообщения

x(t) и в канале по­

мимо белого шума действует помеха

Ѵ(г\,і),

являющая­

ся детерминированной функцией случайного марковско­ го процесса т](t). Таким образом,

(5.143)

Компоненты l(t), x(t), r\(t) образуют в совокупности трехмерный марковский процесс. По физическому содер­ жанию задачи компонент %(t) известен, следовательно, два остальные компонента образуют условный марков­ ский процесс, который характеризуется соответствующи­ ми апостериорными плотностями вероятности. Введем новые обозначения для апостериорных плотностей веро­ ятностей, удобные при анализе задачи обнаружения:

WpS(xi, Ï](|'Ê'O, =1) —апостериорная плотность веро­ ятности компонентов xt, r\t в момент t при наблюдении реализации 1(4) на интервале [0, t] и при наличии сигнала;

wPs(r\t\Vo,Q—0)— апостериорная плотность вероят­ ности компонента щ в момент t при известной реализа­ ции \(t) и при отсутствии сигнала.

При дискретном наблюдении в момент t = kAt апосте­ риорные плотности вероятности запишем в виде

w(xk, i f c l t Ѳ=1 ) и ш ( % | і £ Ѳ = 0).

(5.144)

Для плотностей вероятностей (5.144) справедливы уравнения нелинейной фильтрации, которые были полу­ чены выше. Однако для принятия решения о наличии или отсутствии сигнала необходимо в конечный момент вре­ мени Т = пМ знать иные апостериорные вероятности, а именно:

w (Ѳ = 1 I £" )

вероятность наличия сигнала

при усло­

вии наблюдения

?" и ш(0 = О|£'ог) — вероятность

отсут­

ствия сигнала при том же условии.

207

Если использовать критерий максимальной апосте­ риорной вероятности, то правило решения формулирует­ ся следующим образом:

п

И Ѳ = 0 | Ц )

ѵ

'

Обычно оптимальные приемники обнаружения работают по другому алгоритму

. ( S i e - . )

ш(і'э МѲ=0)

где Л,і отношение правдоподобия,

Я порог.

Величины In и An связаны между собою простым со­

отношением

 

/„ = -£-Л„.

(5.147)

где р и 7—априорные вероятности наличия и отсутст­ вия сигнала соответственно.

Основная задача состоит в том, чтобы, используятеорию нелинейной фильтрации, вычислить сначала ве­ роятности (5.144), затем на их основе определить вероят­

ности до(£"|Ѳ = 1), ви(£"|Ѳ = 0) иja

заключение реализо­

вать правило

(5.146).

 

 

Отметим,

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

плотности

(5.144) вы­

числяются в каждый момент времени t^[0,

Т]. Следова­

тельно, после

соответствующего

пересчета

отношение

правдоподобия можно также определять с течением вре­

мени по мере поступления

и обработки

реализации

Это означает,

во-первых,

справедливость соотношений

(5.145) — (5.147)

для любого l ^ f e ^ n

и, во-вторых, соз­

дает благоприятные условия для применения двухпорогового решающего устройства, работающего по крите­

рию последовательного

наблюдателя.

Это устройство

принимает следующие решения:

 

Ш если

Л й > Нъ,

то

сигнал^есть;

 

если

Л ь < Я "

, то

сигнала нет;

(5.148)

если • Н" < Л Й <^На

, то наблюдение

продолжается.

Здесь Я® и

Н" — верхний и нижний

пороги соответ­

ственно.

 

 

 

 

 

208

Впервые задача обнаружения марковского сигнала на фоне белого шума при дискретном наблюдении была по­ ставлена и решена 'в работе [28]. Затем в [29] было полу­ чено дифференциальное уравнение для вычисления лога­ рифма отношения правдоподобия Zt = \nAt применитель­ но к задаче обнаружения диффузионного марковского процесса в шуме

 

 

 

 

 

. k =

<F(x,t)>,

 

 

(5.149)

где

<F(x,t)>

 

определяется (5.98).

 

 

 

Наконец,

в

[13] 'было

выведено

уравнение

 

 

 

k

=<F

(х, т], 0 > -

< ^

(х, 0 >

(5.150)

где

<F(xJt\,t)>

 

находится согласно (5.95).

величину

 

Соотношение

(5.150)

позволяет

вычислить

zt в задаче

обнаружения, когда

на вход приемника по­

ступает

смесь

W)=S(x,t)

+

V(i\,t)+n(t),

 

 

 

 

 

 

где

x(t)

и и (t) — одномерные диффузионные

процессы.

 

Уравнения

(5.149) и

(5.150)

были получены с исполь­

зованием понятия ненормированных апостериорных мер [6]. Этот метод хотя в известной степени и сокращает выкладки, но вместе с тем сложен для понимания. По­ этому представляющее основной интерес уравнение

(5.150) ниже выводится иным способом,

основывающем­

ся на результатах работы [28].

 

 

Пусть для простоты S(x,t)=x(t)

и

Ѵ(ц, t) — r\(t).

Разберем сначала случай дискретного наблюдения, что­

бы получить рекуррентное соотношение для

отношения

апостериорных вероятностей Ль.

 

 

 

 

 

 

 

При

дискретном

времени

на

вход

приемника

посту­

пает

последовательность значений

 

=

(£„,

£х ,...,

 

которая

 

непременно

содержит

в

себе

отсчеты

белого

шума

пко

+ 1

(,гй, /г,

пк+і),

помехи

тт*+ 1 =

(т]0,

TJ,,... ,

%+,)

и с

вероятностью

р (р—

1 q)

может

содержать

значения

 

подлежащего

обнаружению

 

сигнала

xko+i

=

==0>

 

хх,...,

хк+1).

 

 

 

 

 

 

 

 

 

В зависимости

от

наличия

(0=1)

или

отсутствия

(0 = 0)

сигнала будем рассматривать

соответственно трех­

мерный

( ^ + 1 ,

xko+1,

7)*)

или

двумерный

 

т,*+1)

мар­

ковские

процессы.

 

 

 

 

 

 

 

 

 

 

2G9

Для обоих процессов запишем плотности вероятности переходов за элементарный .интервал наблюдения. По аналогии с (5.34) имеем

 

 

 

Ѵ(%и

Xi,

 

 

Хі-и T U - l ) =

 

 

 

 

=

V (Xi,

î] i \ Xi-i,

T) i-i)

p (li—Xi—T]

i) ;

(5.151)

 

«(Іь

iliIli-u Щ - І )

=v(r\i\Ï]Ï_0p(îi—m),

(5.152)

где для 'краткости

обозначено

 

 

 

 

Р fa Xi

T),

 

 

 

 

 

 

] ;

(5.153)"

 

P (5» —

 

 

 

 

 

] •

(5.154)

С

учетом

(5.151) — (5.154)

совместная

плотность рас­

пределения

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

последовательностей

х о + 1 ,

\ + [ и

значений

Ѳ р вна:

 

 

 

 

 

 

со (50

, л 0

,

TJ0

, ö) =

 

 

 

P П

P & -

Xi

дао

(•*<>. le) П 0 ^г'+"

 

 

1=0

 

 

 

 

 

і=0

 

 

 

 

 

 

1 } г + 1

I Xi,

Т),:)

при Ѳ = 1;

 

(5.155)

 

ft-r

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я П P fa -

 

Ш О

do) П v ^

>1

'ч*)

X

 

где

8 (л;г-) — дельта-функция. .

 

 

 

 

 

По формуле

условной

вероятности

 

 

 

®(*o

.öl*,,

)

» ( * 5

+ 1 , - ч 5 + 1 . Ё § + , . Ѳ )

(5.156)

 

 

 

 

 

где

tw(^+,)==

2

J

... Jatf(jc*+!,

 

$+i,B)dx0...dxk+1X

 

 

0=0,1

 

 

 

 

 

 

 

 

Xd%

Учитывая,

что

да(^+1)

=

ш

 

| $ )œ(Ç) .

210

Соседние файлы в папке книги из ГПНТБ