
книги из ГПНТБ / Филиппов, Л. И. Основы теории радиоприема дискретных сигналов (синтез и анализ) [монография]
.pdfИз последних двух соотношений получаем энергетический спектр фазово-манипулированных колебаний
|G,MI! = ( ' + ! ) / v V 2 |
‘ ( - " Й |
- |
р= —00 |
ЛЧ„ 2 |
|
|
< 3 - 3 7 > |
Энергетический спектр на выходе идеального фильтра получаем по общему правилу, перемножая (3. 37) на квадрат коэффициента передачи фильтра:
У-1
2
G (4P = ( i + t ) tvV 2 N-1 |
|
р=- |
|
ЛЧ0 |
<338> |
И —£)• |
Автокорреляционная функция выходного колебания есть об ратное преобразование Фурье для (3. 38):
У—1 |
|
|
|
2 |
1 |
|
|
N J /Ут0 2 |
(3. 39) |
||
2тсДЧп " |
|||
р=— -У—1 |
|
|
Просуммировав конечную сумму геометрической прогрессии в (3.39), получаем
niM — — |
( ‘\ I |
1 ^ |
1 sin (^foo) |
1 |
1 \ — 2% |
\ ' |
N ) |
А'т0 sin (ти/ЛЧ0) |
2tzN,zq ‘ |
При т = О |
|
|
|
|
1 |
1 |
|
1 |
1 |
Ра(0): 2и |
/VV0+4)А - 2тсЛЧо |
2titq |
Нормированная автокорреляционная функция выходного колебания равна
1 ( л , |
1 \ sin (ит/т0) |
1 |
(3.40) |
РW — ж ( 1 + ж ) sin (m/N'Zo) |
N |
||
В точках, кратных Ат0, она равна единице, |
а в точках, |
кратных т0, |
|
но не кратных Ат0, |
|
|
|
Р(/” о ) = ^ 4 ‘ * |
* = •■•— 2, |
— 1 , О, 1 , |
2 , .. . |
При достаточно большом А |
эту величину можно считать равной |
71
нулю |
и соответственно функцию автокорреляции в точках т = |
= Рхо |
(Р ¥= Nk) считать равной нулю. |
Аналогично изложенному можно показать, что если исходной последовательностью возбуждать идеальный полосовой фильтр с центральной частотой ш0, то функция автокорреляции выход ного колебания будет иметь огибающую (3. 40) с косинусоидаль ным заполнением, так что
Условие (3. 32) при большом N будет выполнено.
Ч а с т ь II
МАТЕМАТИЧЕСКИЙ СИНТЕЗ РАДИОПРИЕМНИКОВ ДИСКРЕТНЫХ СИГНАЛОВ
Г л а в а 4
СИНТЕЗ ОПТИМАЛЬНЫХ РАДИОПРИЕМНИКОВ НА ОСНОВЕ ВРЕМЕННОГО ПОДХОДА
4.1. Вводные замечания. Синтез и анализ
Концепция синтеза оптимальных структур в радиотехнике воз никла сравнительно недавно в связи с разработкой общей (или статистической) теории решений В. А. Котельниковым, Н. Вине ром, Ф. Вудвортом, А. Вальдом и рядом других ученых.
Понятие «найти оптимальное радиоприемное устройство» озна чает определить характер и последовательность операций, кото рые необходимо выполнить над входным процессом, т. е. над сово купностью сигналов и помех, чтобы обеспечить наименьшую возможную в заданных условиях вероятность ошибок.
В течение первого периода развития радиотехники возникнове ние предложений об улучшении качества радиоприемников носило в определенном смысле случайный характер и было связано с усовершенствованием понимания путей осуществления основных операций (усиления, детектирования и др.). Если в процессе этого изучения обнаруживалось несколько возможностей выпол нения операции, то путем сопоставления преимуществ и недостат ков выбирался лучший вариант. Понимание слова «лучший» но сило цри этом часто качественный полуинтуитивный характер. Уверенности в том, что нельзя отыскать еще более совершенный способ при таком подходе, не было.
Лишь с развитием статистического подхода к системе передачи сигналов появилась возможность теоретически отыскать структуру «наилучшего» радиоприемника. При этом термину «наилучший» был придан вполне определенный количественный смысл, основан ный на математической концепции критерия. Это направление в радиотехнике получило название синтеза оптимальных прием ников.
Теория не может, конечно, указать способ конкретного выпол нения той или иной операции, например, перемножения, фильтра ции, ограничения и др. Эти действия могут, в частности, выпол
73
няться при помощи схем аналоговых операций, т. е. схем, состоя щих из электронных приборов и линейных элементов. Но они могут выполняться и при помощи цифровых вычислительных ма шин после представления входного процесса в цифровой форме и обработки данных по программе. Таким образом, теория син теза указывает структуру радиоприемного устройства или про грамму математической обработки принимаемых колебаний.
В реальных условиях результаты синтеза обычно не могут рассматриваться как окончательное решение вопроса о способе построения приемника, так как синтез, во-первых, проводится на упрощенных моделях реальных условий работы и, во-вторых, полученная структура оптимального приемника может оказаться весьма сложной и даже физически нереализуемой. Путем незна чительной потери качества часто удается существенно упростить устройство приемника. Получаемые структуры при этом являются подоптимальными. Степень цотери качества может быть оценена лишь после сравнительного анализа вероятности ошибок, созда ваемых оптимальным и подоптимальным приемниками. Таким об разом, математический синтез указывает пути, по которым сле дует идти. Совокупность же синтеза с последующим анализом позволяет найти разумный компромисс теории с требованиями практики.
4.2.Математическая формулировка задачи
осинтезе оптимального приемника
1.В дискретной системе связи любое сообщение (например последовательность букв, цифр, фраз — команд и т. д.) передается путем сопоставления каждому сообщению отдельного сигнала или их комбинации из конечного набора М сигналов. В наиболее
широко распространенной бинарной системе связи М —2, так что
в канал |
связи посылается один из двух сигналов s (t, та) или |
s (t, т2) |
одинаковой длительности Таи полосы частот А ш0. Струк |
тура (форма, спектры) сигналов может быть самой различной. Существен в данном случае тот факт, что на вход приемного устройства может поступить на каждом интервале Таодин из двух (в общем случае М) сигналов, неискаженные формы которых из вестны («запасены») в точке приема.
Сигналы s (t, mx), s (t, т2), . . ., s (t, тм) подвергаются иска жениям в канале связи, кроме того, к ним добавляются помехи п (t). В итоге, в соответствии с (2. 17) колебание на входе прием ника на интервале 0 ^ t ^ Тс, которое будем называть принятой реализацией, представляется как
74
где q равно 1 , 2 , . . Ж, в зависимости того, который иэ сигна лов передавался на рассматриваемом интервале \
Задала радиоприемного устройства — решить на основании принятой реализации (4. 1), какой именно из М сигналов содер жится в реализации у (t). Из-за налитая помех при вынесении решений неизбежны ошибки, состоящие в «перепутывании» сиг налов.
Разные по устройству приемники будут иметь разные вероят ности ошибок при приеме сигналов. Приемник, который будет давать наименьшую вероятность ошибок, назовем оптимальным. Возможны и другие критерии оптимальности. Однако мы будем пользоваться указанным. Критерий минимальных ошибок оправ дан логически, если мы не имеем основания предполагать,, что опасность последствия принять, например, i-й сигнал за к-й по чему-либо больше, чем опасность принять к-й сигнал за i-й.
В дальнейшем мы будем предполагать, что сигналов может быть только два, т. е. будем рассматривать бинарные системы связи. В рамках теории синтеза это, как будет видно, не является существенным ограничением.
2. Аналогично тому, что было отмечено при рассмотрении оценки параметров в § 2.4, наиболее полным описанием резуль тата радиоприема бинарных сигналов служит определение апосте риорных вероятностей передачи сигналов s (t, тх) и s (t, т2), кото рые даются теоремой Байеса. Таким образом, мы будем рассматри вать передачи сигналов s (t, тщ) и s (£, m2) как случайные собы тия А жВ, происходящие с вероятностями Р (тх) и Р (т2) соот ветственно. Появление реализации у (t) также является случай ным событием В, которое может произойти совместно с любым из
событий А х или А 2. Согласно теории вероятностей |
вероятность |
совмещения событий |
|
Р (АВ) = Р {А) РА (В) = Р (В) РВ{А), |
(4. 2) |
где А — любое из событий А х или А 2, Р (В) — вероятность появ ления события В (реализации у (t)) независимо от того, имело место событие А х или А 2; РА(В) — вероятность события В, если имело место определенное событие А, т. е. А х или А 2 (посылка первого или второго сигналов); Рв (А) — вероятность того, что имело место определенное событие А , т. е. А х или А 2, если из вестно событие В (принятая реализация). Эта величина и является апостериорной вероятностью события А.
Так как событие В может произойти совместно с любым из событий А, то
Р(В) = |
Р (Ах) PAl(В) + Р (А2) Р^В).1 |
1 Здесь мы игнорируем |
время распространения сигналов от передатчика |
к приемнику. Начало отсчета времени совмещено с моментом прихода сигнала к приемнику. Практически это достигается путем предваритель ного поиска момента прихода сигналов.
75
В терминах передачи сигналов из (4. 2) получаем
P„[s(t> »»,)! = |
Р К ) Ртд [V (01 |
Р (mq) Ртд№(0] |
|
(4. 3) |
|
р [г/ (*)] |
—1 |
|
* |
||
|
|
2 |
Р (т *) Ртд[У (01 |
|
|
|
|
q = l |
|
|
|
где р,Пд [у (<)] — плотность вероятности |
конкретной |
реализации |
|||
на входе приемника, |
если передавался сигнал s (t, mq). Эти вели |
||||
чины называются правдоподобиями сигналов s (t, тх) |
и s (t, |
тг). |
Величины Ру [s (t, mq) ] являются апостериорными вероятностями передачи обоих сигналов.
В теории статистических решений доказывается, что, вычислив апостериорные вероятности первого и второго сигналов, опти мальный приемник должен принимать решение в пользу того сиг нала, апостериорная вероятность которого больше. Обычно с вы сокой степенью достоверности можно полагать, что априорные вероятности обоих сигналов равны, так что Р [/rea]= P [?7i2] = 1/2. Кроме того, энергии этих сигналов одинаковы. В этом случае вероятность приема реализации у (t), p,„q [у (t) ] одинакова для обоих сигналов. Поэтому вместо вычисления и сравнения апосте риорных вероятностей сигналов s (t, mq), Ру [s (t, mq)] можно вычислять и сравнивать правдоподобия сигналов p„lq [у it) ] — плотности вероятностей приема реализации у (t) в предположе нии, что передавались первый или второй сигналы.
Таким образом, оптимальный приемник принимает решение, что передавался сигнал s (t, m j, если
Рт1\У (t)]> РтгI*/ (<)l
или если |
|
* = рЖ11у (*)1 — |
(4- 4) |
Реализацию у (t) в (4. 1) перепишем так: |
|
y(t) = s(t, mq, а0, otj, а21 . . .) + тг (t). |
(4.5) |
Эта форма записи подчеркивает, что составляющая сигнала в у (t) зависит от времени, номера переданного сигнала (первого или второго) и совокупности случайных параметров а0, ах, а2, . . ., зависящих от свойств канала связи. Составляющая помехи п (t) может иметь самые различные свойства. Мы будем предполагать, что это стационарный нормальный шум с равномерной спектраль ной плотностью N„.
Плотность вероятности реализации такого шума выражается так [4, 6]:
(4.6)
76
где N = k fN 0 — мощность |
шума; |
ге==2Д/Г; Д/ — полоса частот, |
в которой сосредоточен |
спектр |
реализации; Т — длительность |
реализации; п{ — выборочные значения п (t).
Пусть в выражении (4. 5) параметры а0, ах, а2, . . . приняли некоторые постоянные значения. Тогда вследствие аддитивности сигналов и помех вероятность данной реализации у (t) равна, очевидно, вероятности того, что помеха примет значение, равное
разности г/ (if) — s (t, |
mq, a0, a1) a2, . . .). |
Следовательно, |
можно |
утверждать, что |
|
|
|
P„lq[y(t)] = |
P„[y(t) — s{t, mf, a0, |
a: , a2, ...)], |
(4.7) |
где pn — функция, описывающая плотность вероятности реализа ции шума (4. 6). Но параметры а0, аг, а2, . . . в (4. 7) — случай ные величины, которые при приеме сигналов могут принимать различные значения. Поэтому для получения вероятности реали зации у (<) в среднем необходимо произвести усреднение получен ного значения (4. 7) по всем возможным значениям случайных параметров. Пусть функция р (а0, а1г а2, . . .) описывает сов местную плотность вероятности величин к0, я2, . . . Тогда в соответствии с правилом усреднения
P„t \У(0 1 = S••' SP’»v (01Р К> *j. a2. ■■•) da0daLda2. . . |
(4. 8) |
Именно эти величины (вычисленные в предположении, что q = t и 2) сравнивает оптимальный приемник для принятия решения о том, какой из сигналов передавался на данном отрезке вре мени Та.
Рассмотрим далее решение задачи о синтезе оптимальных, приемников для ряда последовательно усложняющихся случаев.
4.3.Простейший случай приема без оценки
ис оценкой параметра1
1.Рассмотрим задачу синтеза оптимального приемника би нарных сигналов, прошедших канал с медленной мультиплика тивной и быстрой аддитивной помехами. На входе приемника та
кой сигнал на интервале 0 ^ t ^ Тв имеет вид
y(t) = a(t)s(t, mq) + n(t). |
(4.9) |
Это наиболее простой случай общего выражения (4. 1), когда во всей сумме остается только один член. Мультипликативной поме хой иногда называют случайный множитель a (t). Такой случай редко может встретиться на практике. Однако рассмотрение его имеет большое значение для ясного понимания более реальных ситуаций.
Аддитивная помеха п (<) по условию является стационарным шумом с равномерной спектральной плотностью и нормальным законом распределения. Плотность распределения вероятности
77
реализации шума п (t) дается выражением (4. 6). Мультиплика тивная помеха а (t) пусть также является стационарным процес сом с нормальным законом распределения. В соответствии с (2. 10) его можно представить рядом Котельникова:
а ( * ) = 2 а ^ |
sin 1(*/тв) (t — /т„)] |
sin [(к/ха) (t — /т„)] |
(4.10) |
|
(те/та) (t — j%) |
(*K) {* — Ю ' |
|||
|
||||
3 |
|
|
|
|
где ха=1/Д /а ~ |
время корреляции процесса a (t); a .= a(jza) — |
|||
выборка из процесса а (t) в момент времени t= j тв. |
|
|||
Пусть случайные величины |
независимы между собой и рас |
пределены по нормальному закону. Как и ранее, будем считать, что процесс « (t) изменяется медленно по сравнению со скоростью
передачи |
сигналов, так что та |
Тс. При этом можно считать, что |
|
за время |
длительности сигнала величина a (t) не изменяется, и |
||
в формуле (4. 9) вместо функции а (t) можно поставить |
выборку |
||
из нее, взятую в момент времени t = 0. Процессы п (t) и |
а (t) бу |
||
дем также считать независимыми между собой. |
|
2. Рассмотрим сначала случай, когда величина а0 постоянна и известна в точке приема. Этот случай можно назвать случаем приема «точно известных сигналов». Апостериорная плотность
вероятности реализации у (t) |
в соответствии с |
общей форму |
|||
лой (4. 7) равна |
|
|
|
|
|
Pmq[y(t)]=P„\y(t) — |
mt)]. |
|
|||
Учитывая (4. 6), |
получаем |
То |
|
|
|
|
|
|
|
|
|
Р ^ У ^ = { - ш ) Фе^ |
J |
[?(*) — |
|
m9W dt |
|
/ 1 у /2exp |
To |
Nо о |
Tc |
|
|
|
|
|
|||
[ y2(t)dt + |
-]^- j y(t)s{t, |
mg) dt — |
|||
\2%N ) |
No о |
|
0 |
|
|
|
|
|
|
s2(*> |
(4.11) |
|
Эта формула позволяет синтезировать структурную схему прием ника бинарных сигналов с известными параметрами. Он должен сопоставлять (сравнивать) величины (4. И) для q= 1 и 2. Но, учи тывая монотонность экспоненциальной функции, достаточно вы числять и сравнивать величины, входящие в показатель степени экспоненты. Рассмотрим три слагаемых, определяющих этот по казатель.
Первое слагаемое (интеграл) не зависит от индекса q (номера сигнала) и при сравнении несущественно. Если энергии сигналов
равны, так что
То Т0
Е12= | s2(t, m j d t ^ j s2 (t, m2)dt,
оо
78
Рис. 4.1. Блок-схема оптимального приемника точно известных сигналов а — а„ — постоянная величина; б — а0 — неизвестный параметр
то последнее слагаемое также не зависит от q. Таким образом,
‘приемник будет оптимальным, если он будет вычислять и сравни вать величины
То |
|
q — 1,2 . |
(4.12) |
ьд = 1 |
mq)dt, |
||
о |
|
|
|
Соотношения (4. 12) можно назвать корреляционными интегра лами. Для получения величин Ьдследует проинтегрировать произ ведение реализации у (t) на первый и на второй сигналы, а затем сравнить результат. Решение принимается в пользу того сигнала, для которого корреляционный интеграл будет больше. Соответ ствующая структурная схема приведена на рис. 4.1. Из приведен ного анализа видна необходимость иметь в приемнике копии пере даваемых сигналов s (t, тп^) и s (£, тпг).
3. Рассмотрим далее случай, когда а0 — случайная величина (постоянная за время длительности сигнала), распределение кото рой нормальное, с заранее известными параметрами (средним и дисперсией). В этом случае согласно (4. 8) необходимо при вы числении правдоподобия сигналов произвести усреднение его по
79
параметру а0: |
|
PmqГг/ (г)1= \ р п Гг/ (*) —v (*>mq)\p К)йао- |
(4-13) |
В дальнейшем, как и ранее, множители в p„,q [у (i) ], |
одинаковые |
для q= 1 и 2 , будем опускать, так как при сравнении они не будут играть роли. При этом вместо знака равенства будем писать знак пропорциональности.
Подставляя в (4. 13) выражения для плотностей вероятностей Р„ Гр (/)] и р (а0) (при а = 0 и дисперсии а^), получаем
»Г
P„lq[Р (*)] = \ exp |
|
— Щ $ |
1р (*) — (i. я»,)1 dt— |
da.. |
—со |
l |
0 |
*> |
|
Пусть сигналы имеют единичную энергию, так что |
|
|||
|
|
\ S2(t, |
mq) d t ~ i . |
|
|
|
о |
|
|
Тогда интеграл, входящий в показатель степени экспоненты, равен
То |
тв |
То |
S IР (*) — V (г. mg) f dt = |
j |
Р2 ( * ) — 2а0J р (f) s (t, mq) dt -J- a*. |
0 |
0 |
0 |
Первое слагаемое в правой части этого выражения не зависит ни от номера q, ни от а0, и оно, следовательно, не будет иметь зна чения при сравнении правдоподобий первого и второго сигналов. Поэтому
Pmqly(t)}= $ |
ex p fe -S y(t)s(t, mq) |
d% |
—со |
l 0 |
“ > |
Дополним показатель экспоненты до полного квадрата и пред ставим выражение так:
То Г
•*о
J v(t)s (i, mq) dt |
J у (t) s (t , mq) dt |
A^o \/l//V0H- l/2o2 _ |
/V|j(l//V0+l/2e*) |
Сделаем замену переменной интегрирования по формуле
____ |
То |
JУ (t) s (t, mq) dt |
|
_ лГ i , i |
yV0\/l//V0+ l/2o2 |
0 v Nn+ 2 ol |
80