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

книги из ГПНТБ / Филиппов, Л. И. Основы теории радиоприема дискретных сигналов (синтез и анализ) [монография]

.pdf
Скачиваний:
20
Добавлен:
19.10.2023
Размер:
7.31 Mб
Скачать

если строить схему согласно формуле (4. 17), то она будет опти­ мальной для обоих режимов работы.

Рассмотрим режим установления. Для этого рассмотрим такой интервал оценки Т, в течение которого а (£) постоянно и равно а0, т. е. принимаемая реализация

 

 

Y (t)^ a 0S(t)-\-n{t),

— Г < г < 0.

Апостериорное распределение параметра а0 имеет вид

„ r„ t v

и м __Р (“о) Р (О I “о]----

 

р Lao I 1

р [У (t)1

=

 

= eip{ - s ; !

+

 

(5-15)

где

г

о

О

L = [S *{t)dt; й0 =

1 (

Y(t)S{t)dt.

- т

- т

Формулы (4. 15) и (5. 15) для апостериорного распределения

параметра а0 совпадают. Поэтому bq будут определяться по фор­

муле, аналогичной (4. 17). Однако теперь величина Ьа0 будет

результатом интегрирования произведения Y

(t) S (t):

 

 

 

 

о

 

 

 

 

 

 

 

La0=

j

Y{t)S(t)dt,

 

(5.16)

в то время как в (4. 17) La0

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

низкочастотного

фильтра

при

подаче

на

него произведения

Y (t) S (t).

Но

так как

в

течение

интервала

Т величина a (£)

постоянна и, следовательно, Т

« 0,

то идеальный узкополосный

фильтр с

переходной характеристикой

h(t) =

sin [( ф а) t] выпол­

 

 

 

 

 

 

 

 

 

{ Ф а ) *

нит это интегрирование достаточно точно. На рис. 5.2 приведена схема с использованием в качестве S (t) суммы сигналов s (t, ma)-(- + s (if, тпг), которая является оптимальной как для неустановивщегося, так и для установившегося режима.

Аналогичные доказательства можно произвести и для всех остальных случаев (узкополосные и широкополосные сигналы). В схеме рис. 5.2 и аналогичных ей для узкополосных и широко­ полосных сигналов режим будет меняться от некогерентного до когерентного. Вероятность ошибки также будет изменяться, начиная с величины, равной вероятности ошибки при некогерент­ ном приеме 1 и, кончая величиной, равной вероятности ошибки при когерентном приеме.

Ш

Рис. 5.2. Блок-схема приемника, оптимального в режиме установления

ив установившемся режиме

2.На примере простейшего случая, когда сигнальная состав­ ляющая У (£) имеет вид a (t) S (if), рассмотрим требования по­

стоянства огибающей сигнала S (t) или его квазипериодичности (повторяемости последовательности сигналов s (t) длительностью

Т0 < \)-

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

3. В каналах с переменными параметрами искажение узко­ полосного сигнала сводится к флуктуациям его амплитуды и фазы. Этот эффект называется замиранием сигнала. Для борьбы с зами­ раниями сигналов применяются разнесенные методы приема: разнесение по пространству, по частоте, по углу прихода. При этом на вход всей приемной системы в целом поступает колебание

У(*) = 2

Уь (*) = 2 К (t, т ) -f- пк (t)].

(5.17)

к

к

 

Узкополосные сигналы sk {t, mq) отличаются своей задержкой, амплитудой, начальной фазой и частотой (при разнесении по ча­ стоте). Каждая составляющая в (5. 17) принимается отдельно,

112

а затем производятся различные методы сложения, в ре­ зультате чего достоверность приема увеличивается.

При применении широкополосных сигналов входной сигнал представляется в виде

У(*) = 2 h (t, т ) -f- п (*),

(5. 18)

к

 

где сигналы sk (t, mt) имеют задержку друг относительно друга, кратную 2п/Асое, а также отличаются своими амплитудами и на­ чальными фазами. Выражение (5. 18) аналогично (5. 17). Следо­ вательно, если сигнальные составляющие в (5. 18) разделимы, то метод приема широкополосного сигнала будет подобен разнесен­ ному методу приема. Выше было показано, что требование прямоугольности энергетического спектра сигналов s (t, и s (t, т2) является необходимым для разделимости сигнальных составляю­ щих в (5. 18). Более общим требованием является не прямоугольность энергетического спектра сигналов, а равенство (или близость) нулю автокорреляционной функции сигналов s (t, щ ) и s (t, т2) в точках, кратных 2д/Д и>0. Прямоугольный энергетический спектр сигналов идеально удовлетворяет этому требованию.

Рассмотренные в § 4.5 и 4.6 оптимальные методы приема широко­ полосных сигналов являются методами разнесения по времени с оптимальными когерентным и некогерентным сложением. Можно показать, что если в формуле (5. 17) случайные амплитуды сиг­ нальных составляющих распределены по закону Релея с разными дисперсиями, случайные начальные фазы — равномерно, а по­ мехи пк (t) представляют стационарные белые гауссовы шумы со спектральной плотностью N0, то рассмотренные методы когерент­ ного и некогерентного сложения будут оптимальными при любом способе разнесения.

4.При синтезе оптимального приемника широкополосных

сигналов сначала накладывалось требование Та гкаа, и лишь после дополнительного анализа полученных формул был сделан

переход к более мягкому требованию Т0 ^ тК|Ш.

Возможность

применения сигналов длительностью Тв ^ \ ап, т. е.

достаточной,

чтобы избежать перекрытия сигналов, можно показать и другим более простым способом. Действительно, пусть излучается сигнал с прямоугольным энергетическим спектром, полученный в резуль­

тате возбуждения 8-импульсом

некоторого линейного фильтра

с коэффициентом передачи

F( ш),

абсолютная величина которого

постоянна в полосе А шс.

Длительность этого сигнала Т0 пусть

удовлетворяет требованию

Та

так что за длительность сиг­

нала коэффициент передачи канала можно считать постоянным:

It (ш, *) = £(«>), 0 < t < T c.

На входе приемного устройства поставим фильтр, коэффи­ циент передачи которого F* (ш) является комплексно-сопряжен-

8 Л. И. Филиппов из

ной к F ( со) функцией. Методом, аналогичным методу сведения небелого шума к белому [1], можно показать, что приемник, оп­ тимальный к выходному напряжению фильтра F* (и>), в сово­ купности с этим фильтром будет оптимальным к входному на­ пряжению приемника. Шум на выходе фильтра F* ( ш) останется белым. Легко видеть, что вследствие линейности фильтров R (ш), F ( ш) и F* ( си) последний можно перенести, включив между фильт­

ром F ( (о) и каналом. Помеховые составляющие в обеих схемах также статистически эквивалентны. Последовательное соединение фильтров F ( со) и F* ( ш) является идеальным полосовым фильт­ ром. Поэтому сигнальная составляющая на выходе новой схемы

равна импульсному

отклику

(у, t), длительность которого

тгап+1/А /0 или при

ткав 1/Д/с ее можно считать равной ттоп.

Из сказанного следует, что сигнальную составляющую на входе приемника можно предварительно сжать до длительности тгап. Поэтому, если Тс ^ тгап, то перекрытия сигналов не будет. При Тс < тЕап этого избежать невозможно. В схемах рис. 4.6 и';4.8 предварительное сжатие осуществляется при помощи пропускания входной реализации через линию задержки с последующим выде­ лением на каждом отводе составляющей, имеющей суммарную задержку, равную тган.

В предыдущем изложении мы основывались на том, что дли­ тельность импульсного отклика канала конечна и равна ткап. Однако это не всегда может выполняться из-за случайности ка­ нала. Будем говорить поэтому, что имеется интервал длитель­ ностью ткап, в течение которого переходная характеристика прини­ мает наиболее существенные, в среднеквадратичном смысле, значения и имеются «хвосты», вес которых мал, но длительность которых может быть достаточно большой. Так же как и раньше, будем рассматривать диффузионный канал. При этих предполо­ жениях на каждом интервале решения, кроме основной сигнальной составляющей, будет присутствовать добавка в виде суммы хвостов от предыдущих сигналов. Они будут создавать дополнительную помеховую составляющую. Так как канал диффузионный, то эта составляющая будет нормальным случайным процессом. Преды­ дущие сигналы статистически не зависят от сигнала на данном интервале решения, поэтому помехи, создаваемые ими, являются независимыми от сигнальной составляющей. При применении сигналов с прямоугольным энергетическим спектром энергети­ ческий спектр помехи, создаваемой предыдущими сигналами, будет постоянным в полосе сигнала. Из этого следует, что общая помеха будет иметь вид нормального белого шума с увеличенной спек­ тральной плотностью. Поэтому структурная схема оптимального приемника при сделанных предположениях останется прежней.

114

5.3. Обобщение результатов синтеза

•-оптимальных приемников. Широкополосный прием и разнесенный прием

1. Если сигналы sx (t, mq) проходят канал связи с импульсным откликом h (t, v), то согласно (2.8) колебания (сигналы) на выходе представляются (без учета помех) в виде

(*, mt) = ^ Н (*) si ( * - /c£ o’ mt)

k

(5.19)

Таким образом, в точку приема поступает несколько сигналов sx (t, mq), получивших в процессе прохождения через канал вре­ менное запаздывание, изменение амплитуд и начальных фаз. Это явление можно назвать разнесением сигналов. Его можно ис­ пользовать для повышения помехоустойчивости радиоприема. Если свойства канала таковы, что флуктуации величин ак (t), ак (t) (при разных к) происходят не одновременно (иногда говорят «не дружно»), то, разделяя составляющие в сигналах (5. 19) и обрабатывая их порознь с последующим суммированием, можно получить выигрыш в вероятности ошибок.

В рассмотренном нами до сих пор случае широкополосных сиг­ налов разнесение явилось естественным следствием свойств канала и оказалось разнесением по времени прихода сигналов. Однако часто целесообразно создавать разнесение искусственно, путем направления в точку приема нескольких разных сигналов sk(t,mq), несущих одну и ту же информацию тпд. Такой радиоприем назы­ вается разнесенным или многоканальным г. Сигналы sk (t, mq) могут отличаться средними (несущими) частотами. При этом говорят 0 разнесении по частоте. Возможно, очевидно, разнесение и по пространственным координатам, по поляризациям векторов элек­ тромагнитного поля, по начальным фазам и др. Широкополосный прием, рассмотренный в § 4.5, 4.6, является частным случаем разнесения по времени прихода сигналов sk (t, mq). Приемник в § 5.1 можно рассматривать как частный случай разнесения по частоте (когда все сигналы расположены так, что их спектры при­ мыкают один к другому).

2. Задачу о синтезе оптимального радиоприемника разнесен­ ных сигналов можно рассматривать в общем виде [25]. Пусть на передающем конце излучается непрерывная последовательность равновероятных бинарных сигналов s (t, mq) (<7= 0, или 1) длитель­ ностью Та, с единичной удельной энергией и коэффициентом взаим­ ной корреляции х. Предположим, что задержки и центральные

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

8 *

115

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

(«)} = {У1 («),

У , (f),

. . ., У,

(0) будет иметь вид

 

Yk(0 = Н (*) 2

(*— iTc> mq) —

 

 

 

(f) 2

Sk {t — tTa, niq) - f nk(t).

(5. 20)

Как и ранее, будем рассматривать посимвольный прием, когда на каждом интервале длительностью Таприемник должен выносить решение о переданном сигнале. Необходимую для этого стати­ стику амплитуд квадратурных компонент (неизвестных параметров) будем получать путем непрерывного вычисления их апостериорной вероятности на основании анализа предыдущих наблюдений. Как и ранее, будем предполагать, что <хк (t) и &к (t) постоянны на интер­ вале Т Гс, который будем использовать в качестве интервала оценки. Ниже укажем, как влияет на схему приемника непостоян­ ство этих функций.

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

Предполагая функции ак (t) и (t) постоянными на интер­ вале (О, Т) и равными ак и ак соответственно, найдем апостериор­ ное распределение параметров {ак, ак}. Предположим сначала, что на интервале (О, Т) переданная последовательность символов тп

известна, т. е. известны колебания Sk(t) = 2 s* — iT0, mg). Тогда,

*

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

p IK> I {*"*(*)}» {5fc(0}l = CiPKa*. a*}lx

(5.21)

где

T

0

116

1.

(*)•?*(*)*;

(5. 22)

« * = — ^rS

Cj — несущественный коэффициент пропорциональности. Множителем при рЦяд., я/с}] в формуле (5. 21) является функ­

ция, представляющая нормальный закон распределения пара­ метров со средними ак и ак и дисперсиями N0/m. При большом т N0/m может оказаться малым, в результате функция, представ­ ляющая нормальный закон распределения, будет играть решаю­ щую роль в произведении. Поэтому априорным распределением рЦосд., а.к) ] можно пренебречь. Нетрудно убедиться, что в случае

нормально

распределенных независимых

параметров

{aft, a^.)

с

дисперсиями

о|

необходимым

условием

для этого

является

a|

NJm,

что

при качественном приеме

всегда выполняется.

Пренебрегая априорным распределением, получаем

 

р [{“»• «,)|<г»<‘» .

{J , (0}1 —

 

 

 

 

= ( к У еч>'

т

— §1; 2

<5* — “ »)"I ■

<5-23>

 

г; 2 <“ * —

 

 

 

 

к=1

к=1

 

Это выражение нужно усреднить по всем возможным гипотезам относительно последовательностей равновероятных символов тд. С последовательностью символов тд формулами (5. 22) однозначно связан вектор параметров {ак, ак}, поэтому усреднение по гипоте­ зам эквивалентно усреднению по параметрам {а&, ак). Так как число гипотез фиксировано (2т), то распределение р [{aJ:, ак) ] дискретно. В таком виде оно неудобно для построения приемника, однако его можно аппроксимировать непрерывным распределением.

Введем обозначения:

 

(••+1)т0

 

 

sk(t — 1Т0, 1)1 dt,

Чы=

Т

S

(*)[»*(*— iT* °) +

 

(»'+l)T0

 

 

h (t — iT°> l)]dt,

4w =

y

S

Yk(t) lh {t — iT^ °) +

 

(*'+l)Tc

 

 

 

 

(5. 24)

 

 

 

 

 

1)] dt,

 

 

\

^ k ( 0 LSfc @

^ o > 0 )

s k ( t

t T 0 ,

 

(>'+l)T0

^ 0»0)

 

 

 

%lci =

~2~

1

^k Wf^fc

S,. (t

iT0,

1)] dt.

^0

117

Тогда

(5. 25)

Из формулы (5. 25) следует, что

(5. 26)

Принятые реализации Yk (t), а следовательно, и величины £/с< могут быть любыми. Поэтому интервал (5. 26) заранее не фикси­ рован. Из свойства эргодичности и закона больших чисел следует,

что при m -м ю он сходится по вероятности к 2 |£ь. j = 2 1£к|, где усреднение взято по всем возможным реализациям Yk(£). При

большом т интервал (5. 26)

в вероятностном смысле будет мало

отличаться от

2 ||.

 

Отношение

т | |/2iV0

определяет пределы изменения пока­

зателя экспоненты в формуле (5. 23) при рассмотрении всех воз­ можных ак. Оно увеличивается прямо пропорционально тп, а число значений параметра ак увеличивается экспоненциально (2т). Это говорит о возможности замены дискретного распределения параметров {ак, ак) (суммирования при усреднении) при большом тп непрерывным (интегрированием). Вид распределения легко определить из формулы (5. 25). Сумма большого числа случайных величин (—I)5,- £к{ (случайны д{) согласно центральной предельной теореме нормализуется, поэтому распределение р [{afc, ак) ] можно считать нормальным и

i=0

тт-1

=

\

Yk Q) 2

Is* (* ~ tT°’ °) + sk (*— 1Tv 1)1 dt>

 

0

 

(5.27)

 

 

 

 

T

 

 

 

S

2

tT» ° ) + §к(*— 1То- i ) i * -

 

0

 

 

118

Корреляционная матрица имеет вид

7 7 1 -1

(ч — ч) ( ч— «,)=^5 2

i-0

7 7 1 -1

( « * - a fc)(5r- a r) = ^ 2

(5- 28)

7 = 0

7 7 1 -1

(ак Ч) (Ч Ч) = ^2 2 t=0

Усредняя распределение (5. 23) по случайным параметрам {аЛ., аА), находим

Р [{<**. Ч ) I { Y k- ( О Д = c i е х р

[

— j

1 —

* * . Ч —

Ч If X

Х ^

+

^

Г К - * * ’

( 5 ' 2 9 )

где I — единичная матрица; индекс t

означает

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

(вектор с этим индексом

означает

вектор-строку); Съ — несу­

щественный коэффициент

пропорциональности.

Из формул (5. 27) и (5.

28)

следует, что

элементы вектора

{ак, ак] вычисляются при помощи интеграторов на длительность Т, а’ элементы корреляционной матрицы Ка — при помощи идеаль­ ных интеграторов на длительность Те, умножителей и сумматоров. Нетрудно показать, что с учетом непостоянства функций ak (t) и S.k (t), но при условии, что на интервале Тс их можно считать постоянными, распределение параметров останется прежним, однако вектор (aft, ак) и матрица Ка будут вычисляться иначе. Суммирование в (5.27) и (5.28) будет происходить с некоторыми весами р{ и р% а вместо 1/гъ будет другой коэффициент С. Весо­ вое суммирование эквивалентно замене идеального интегратора узкополосным фильтром низкой частоты. Определение весов р{ (формы переходной характеристики фильтра низкой частоты) и коэффициента С эквивалентно задаче линейной фильтрации про­ цесса ak (t), находящегося в смеси с белой помехой.

Распределение (5. 29) используем для построения оптималь­ ного алгоритма распознавания на интервале (Т, Т-\-Тс), для чего на основании принятой на этом интервале совокупности реализа­ ций {ук (t)} необходимо найти апостериорные вероятности илифункции правдоподобия символов тд (д=0 или 1), предполагаемых равновероятными. После обычных вычислений получаем следую­

щую

формулу

для

функций правдоподобия:

 

 

Р

К г / * ( * ) } I

" V

{ Y k ( О Д =

с з е х р

1|

% к II* X

,

X

1 II

\ q k II

2/V q II

^ к '

Ч f

X

119

т •*■)

а А> %Тс “ fcll}>

( 5 . 3 0 )

где

 

 

г+гс

 

 

^ = J Ук(t)Sk (t -

тпТ0, mt) dt,

 

L r .

 

(5- 31)

%тс=— \ УА*)§Л* — тТ* mq)dt.

т

Как и ранее, приемнику достаточно вычислять и сравнивать ве­ личины Ъд, равные показателю экспоненты выражения (5. 30).

3. Осуществить приемник, соответствующий этой формуле, сложно. Используем медленность изменения параметров для ее упрощения. Члены матрицы Ка (5. 28) при те, —> то сходятся по

вероятности к величинам 1/те. £г, которые обратно пропорцио­ нальны те. Удобно поэтому ввести новую матрицу Къ—тКа. При большом те члены матрицы (Kb/KJ/m могут оказаться ма­ лыми. Тогда, разлагая матрицу в ряд по степеням 1/те, можно ограничиться двумя членами:

[1+ т г ( 1+ £ ) Г “ , - и ( 1+ £ ) -

(5' 32)

В соответствии с этим приемник должен вычислять и сравнивать с нулем величину

= &0 -

й 1 =

2 ||; -0А-

 

?0i

if I

fl a k ’

II ~ h

 

+

W

i f ( / ^ - ^ ’ i

+

I )ll^ofc>

iofcll-

 

 

 

Ы ‘ ( ^ ^

+

I ) l l ^ .

u

-

 

 

2

lit

6

 

 

+

 

8 J .

( 5 . 3 3 )

т

1

 

 

 

 

 

 

 

 

 

 

 

При простом методе кодирования, когда каждому символу бинарной информации тд сопоставляется свой сигнал s (t, те}), первый член в выражении (5. 33) намного больше остальных, т. е. схема приемника в этом случае определяется формулой

= ^ |

£lft, %йк |[Ь |f 11| dk, ak|=

= 2 2 tdk (t0k ilft) -f- dk(l0fc

(5. 34)

k = l

 

Физическая трактовка получившегося алгоритма распознавания (5. 34) следующая. Сначала происходит оценка величин <хк и ак, а затем — разнесенный прием по алгоритму точно известных па­ раметров канала (5. 34), который называется когерентным разне­ сенным приемом. Для оценки неизвестных параметров канала

120

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