книги из ГПНТБ / Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи
.pdfОбозначим для |
краткости |
|
К |
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