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

книги из ГПНТБ / Гинзбург, В. В. Теория синхронизации демодуляторов

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

где I (to, t*) — функция потерь, рассматриваемая как плата за вынесение решения t* о положении границы посылок (в общем случае — оцениваемого параметра), в то время как в действи­ тельности граница попала в момент t0. При функции потерь I(t0, t*)=p(t*t0) из (7.2) получаем (7.1).

Величина Рош R) зависит от вида оператора t*(x). Задача синтеза оптимального УС сводится к выбору оператора, достав­ ляющего минимум Р0ш-

Как видно из гл. 6 , условная вероятность ошибки p(t*10) яв­ ляется довольно сложной функцией момента t* (или фазы синхро­ сигнала). Приходится поэтому заменить p(t*to), например, так

называемой простой функцией потерь

(7.3)

I(t0, t*) = C - 8 ( t * - t 0),

где С — постоянная, вводящаяся для удобства решения

задачи

•синтеза. Функция потерь (7.3) означает, что любая ошибка при

оценке параметра

t0 одинаково

опасна и

лишь точное

решение

t* — to считается выгодным.

функция

поте|рь (7.3)

близка

С точ1ки зрения

синтеза УС

к логарифму условной плотности вероятности, еши последняя имеет острый минимум в окрестности решений t*, близких к t0. По-ви­ димому, такая ситуация довольно типична при оптимальной об­ работке сигнала в демодуляторе, особенно при больших отноше­ ниях сигнал/помеха.

Часто рассматривается также квадратичная функция потерь

I(t0, t*) = ( t * - t 0)\ (7-4)

которая в ряде случаев позволяет упростить синтез. Иногда, к то­ му же, оказывается, что синтезированные с функциями потерь (7.3) и (7.4) оптимальные оценки совпадают [82, 136].

Методы синтеза и структурные схемы оптимальных УС. При синтезе оптимальных УС удобно выразить плотность вероятности

wx (x) в (7.1) через условную плотность вероятности вектора па­ раметров сигнала х при данном значении оцениваемого парамет­

ра to, называемую функцией правдоподобия и обозначаемую как

•—►

F(x\t0) или через апостериорную плотность вероятности оценивае­

мого параметра Wt(to\x). В соответствии с формулами полной ве­ роятности и умножения вероятностей

«М *)= f F(x\to)w0{to)dt0.

(7.5)

—6,57*

Решение задачи синтеза при функциях потерь (7.3) и (7.4) приведено, например, в [15, 85, 88, 125, 136]. Оптимальный опера­ тор при квадратичной функции потерь, определяющий так назы­ ваемую эффективную оценку [88], равен

°-5г

t* (х) =

Г twt (11х) dt.

(7.6)

-0,5Г

170

Оптимальный оператор три простой 'футкции потерь, опреде­ ляющий оценку с наибольшей апостериорной .вероятностью

t* (х) — argmaxtt>f(/|je).

(7.7)

Алгоритмы обоих оптимальных УС содержат операцию вычис­ ления функции апостериорной плотности вероятности (фапв) по

реализации сигнала х. Соответствующее устройство в оптималь­ ных УС (рис. 7.1) будем называть вычислителем апостериорной

Рис. 7.1. Структурные схемы оптимальных УС при функциях потерь:

а) простой; б) квадратичной

плотности вероятности (ВАП). Первое УС содержит, помимо уст­ ройства ВАП, фиксатор максимумов (ФМ). Во втором устройстве изменяющаяся во времени фапв подается на перемножитель, на другом входе которого действует линейно изменяющееся на ин­ тервале (—Г/2, 7/2) напряжение. Интегратор (И) к концу инте­ грирования выдает на выходе напряжение, пропорциональное оп­ тимальной оценке положения границы между посылками относи­ тельно, например, момента, когда линейно изменяющееся напря­ жение равно 0. В отличие от устройства рис. 7.1а, здесь может потребоваться еще преобразователь «напряжение—время», не по­ казанный на рис. 7.16.

О вычислении апостериорных вероятностей. Непосредственное вычисление фапв обычно сложно и при синтезе ее представляют через функцию правдоподобия (фп). Выразим алгоритмы (7.6) и (7.7) через фп. По теореме умножения вероятностей

wt{t\x)

F (x \t) w0 (t)

(7.8)

 

wx (x)

Если априорная информация о положении границ посылок от­ сутствует, то их априорное распределение можно считать равно­ мерным на интервале (—7/2, 7/2), т. е. w0(t) = \!T. С учетом это­ го алгоритм (7.6) при подстановке (7.8) принимает вид

t*{x)=g(x) f

tF(x\t)dt,

(7.9)

-0 .5 Т

 

 

где

0,5Г

_

 

g(x) = [Twx (x)]~i =

 

j

F{x\t)dt

 

 

0,5Г

 

 

— известная функция вектора параметров сигнала х.

7*

171

При подстановке (7.8) в (7.7) учтем, что умножение функции на 'положительное число не меняет положения ее максимума. По­ этому вместо (7.7) получим

t* (х) = arg шах F (я 11).

(7.10)

i

 

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

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

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

 

F{x\t) = Y. PiF (x\t, i).

(7.11)

i

 

Ниже используется в основном второй подход.

Приведенные алгоритмы и схемы рис. 7.1 решают задачу син­ теза оптимальных УС. Такое решение, однако, не всегда пригод­ но для практического использования, так как простота схем 7.1 является кажущейся. Реализационно ВАП (или используемый вместо него вычислитель фп) оказывается чрезвычайно сложным устройством, алгоритм которого, к тому же, изменяется во време­ ни. Материал последующих параграфов представляет собой, по существу, попытку упростить структуру ВАП применительно к не­ которым частным задачам синтеза.

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

Методика синтеза оптимальных УС зависит от наложенных /при синтезе ограничений «а память УС и на то, какие из парамет­ ров сигнала, за исключением границы между посылками и номе­ ра варианта сигнала, неизвестны. Смысл термина «память» опре­ делен ниже.

7.2. Синтез УС с ограниченной памятью

Исходные соотношения. Задача синтеза оптимальной оценки неизвестного параметра сигнала Ао обычно рассматривается для

сигнала вида суммы

A0) + n{t), t ' < t < t ' + т,

(7.12)

'x (0 = s(/,

где s(t, Ао) — полезный

сигнал; n(t) — нормальный

шум; (?,

^ 4 -т) — интервал наблюдения.

172

Применительно к этой сумме фп

параметра

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

в виде [85, 125, 136]

 

f'+т Г+т

 

 

 

 

F (х | Л) =

g' exp

 

 

 

 

 

f

f

X (^i) x (^г) 0 (^1,

д л д

-

f'+т

 

f

t'

 

r+ x

 

 

 

-----—Гs(t, A)u(t,

 

 

 

 

 

 

A)dt -f

Г

x(t)u(t,

A)dt

(7.13)

2

,)

 

 

 

 

,)

 

 

 

 

<'

 

 

 

 

r

 

 

 

Здесь g' — постоянная величина; 0(71, /2)

— решение интегрально­

го уравнения:

 

t’+x

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.14)

 

б & - * ,) =

j

0(g,

t ) R ( t - Q d i -

 

 

 

 

t'

 

 

 

 

 

 

 

u(f,

f'+ T

 

A)Q(t, tjdt*

 

(7.15)

 

Л) =

[ s(tv

 

R(t) — корреляционная функция шума n(t), а под вектором x по­ нимается значение суммы x(t) на интервале наблюдения. Функ­ ция u(t, А) может быть определена также как решение инте­ грального уравнения:

<'4-т

 

J R ( t - U ) u ( t , A)dt = s(t„ А).

(7.16)

v

 

Обычно спектры сигнала и помехи сосредоточены в одной и той же полосе частот и, кроме того, частотная характеристика канала может считаться прямоугольной с полосой пропускания от <о0— £2 до (оо+£2> где Q= jtA/\ AF — ширина полосы. Тогда корреляцион­ ная функция помехи при действии на входе канала белого шума равна

п 1 v

 

Sill ОТ

tm 1 7 ,

R (т) =

—------------cos со0т.

(7.17)

 

я

От

 

Подставив (7.17) в (7.16) и полагая, что s(t) вне интервала наблюдения равна нулю, заменим пределы интегрирования в (7.16) на бесконечные

ГSmo Q..{t ~ [г) cos(о0 (t tt) и (t, A)dt = s(tv А). (7.18)

ЛJ Q (t — t2)

со

Если спектр сигнала s(t) сосредоточен в полосе частот (о)0—Q, (oo+Q)1), то из (7.18) следует, что

_________

и(/, A) = (2/oJ)s(l, А),

(7.19)

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

173

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

Первое слагаемое в показателе экспоненты в (7.13) определяет энергетические характеристики сигнала. Во многих задачах, в ча­ стности в задачах синхронизации, влияние оцениваемого пара­ метра на энергетические характеристики незначительно и им мож­ но пренебречь. Второе слагаемое не зависит от сигнала и, кроме того, как показано ниже, часто не зависит от Л. Поэтому оба сла­ гаемые изменяют лишь масштаб фп и могут быть учтены измене­ нием постоянного множителя g' в (7.13). Тогда, приняв во внима­ ние (7.19), запишем фп в виде

 

 

 

 

 

t'-И

 

 

A)dt J,

(7.20)

 

*4*1 Л) = gyexp

\

if

X(t)s(t,

 

 

 

 

i'+xt'+x

 

 

 

 

где

g = g' exp

у

J

J

x(t1)x(t2)Q(tv

t%)* 1 , ^ 2>

 

 

 

V

t'

 

 

 

 

 

 

 

 

 

t'+T

 

 

 

(7.21)

 

у =

exp

2

j*

s(t, A)u(t,

A)dt

 

 

 

 

 

 

 

 

 

Если под Л понимать неизвестное положение границы посы­

лок,

то подстановка (7.20)

в

(7.9) и (7.10) определяет алгоритмы

оптимальных УС по реализации сигнала x(t)

с ограниченной дли­

тельностью (0 , т).

по одной

посылке

при

полностью

известном

Синхронизация

сигнале. Допустим, что интервал наблюдений (?, f-\-x) заведомо перекрывает интервал посылки (toТ, to), т. е. t'<Cto—Т <.t0< t ' + т, причем полезный сигнал вне посылки равен нулю, а внутри ее

принимает с вероятностью pi(i = >1 ,..., I) одно из ‘своих возможных значений

 

s(t, Л) = *(*-*„).

 

(7.22)

Если бы номер переданного варианта был известен, то фп име­

ла бы вид

 

t'+x

 

 

 

 

 

 

d/x

(7.23)

F(x\t,

i) = gYiexP

4 -

f x V J s ^ - t )

 

 

A

i

 

 

 

где определяется ф-лой (7.21). Усредняя (7.23)

по всем вариан­

там сигнала, находим в соответствии

с (7.11) фп:

 

 

м

г

г+т

 

 

(7.24)

F (x \t) =

g \]p ;Y iexP

- у

f

*(*i)*i(*i—if )*!

 

t i

LCT/

r

 

 

 

174

Так как Si(t)=0 вне интервала (О, Т), то пределы интегриро­ вания в (7.24) можно заменить па (tТ, t), если tT>t', t< t'+ т.

Тогда интеграл в

(7.24) можно вычислить с помощью согласован­

ного с /-м вариантом сигнала

фильтра с импульсной

реакцией

s°i(t)=Si(Ti), откуда

 

 

 

 

м

 

 

(7.25)

^ ( * =1

g0 ^ A Y iex p

- т

U

 

£=1

°f

Л

 

Интегралу в показателе экспоненты можно придать более при­

вычный вид, если, приняв во внимание, что (0 = 0 ПРИ t~>T, за­ менить нижний предел интегрирования на нуль. При этом

м

г

t

 

 

(7.26)

’(*10= % Pi Y<exP

4

j *(*i)s?(* — tjdti

?=\

L af

о

 

 

 

'В частности, три однократной

модуляции

(М =12)

шротивопо-

ложными равновероятными сигналами из (7.26) получаем

 

t

 

 

 

(7.27)

F(x\t) = gych - у

 

dtt

J

°f

j

 

 

 

Итак, оптимальное УС по одной посылке содержит набор со­ гласованных с вариантами сигнала фильтров, набор нелинейных преобразователей с экспоненциальными характеристиками, уст­ ройство сложения с весами, а также преобразователи с алгорит­ мами (7.9) или (7.10). Существенно отметить, что на входы согла­ сованных фильтров сигнал x(t) поступает непрерывно, а сделан­ ное выше допущение об ограниченности сигнала оказалось излиш­ ним и свелось к ограниченной длительности импульсной реакции согласованного фильтра. Эта длительность и определяет память УС.

Синхронизация по одной посылке при неопределенной фазе си­ гнала. Согласованный фильтр в синтезированном выше УС может быть составной частью РУ когерентного приемника. Аналогичное УС можно синтезировать для оптимального некогерентного при­ емника, обеспечивающего наименьшую вероятность ошибки при равномерном распределении начальной фазы сигнала [56, 132, 137].

Считая для общности систему связи многоканальной, г-й ва­

риант сигнала можно записать в виде

 

L

(7.28)

s{(0 = V ал cos (о/(- + фу, + ф/0).

mm

 

/= 1

Здесь фу0 — случайные независимые величины, равномерно рас­ пределенные в (—я, я). Полагая, что спектр помехи равномерен в полосе частот, занимаемой сигналом, подставим (7.28) с учетом (7.19) в (7.24). Проинтегрировав полученное выражение по пе­ ременным фуо и заменив с учетом конечной длительности сигнала

175

пределы интегрирования, как это было сделано при известной фор­ ме сигнала, видим, что фп границ посылок

М

я

 

я

г

t

L

 

* 4 * 1 0 = g j ] PiYij

• •

j exP

-7

J

 

+

<=!

о

 

0

L af t—т

/= 1

 

 

L

 

М

T

 

 

+ ф/t + ф/о)

П ^Ф/0 = £

 

Pi Yi П

А) Г ~

(0 ’ (7.29)

где

y=l

 

(=1

/= 1

L°f

 

 

 

 

 

 

 

 

 

Z%(t) = X*t (t) + Y%(ty,

 

 

 

 

 

i

 

 

 

 

x ii (0 =

a/i

J * (*1)cos (©</ A +

Ф/i) * 1;

 

 

 

 

t - T

 

 

 

 

 

 

 

t

 

 

 

(7.30)

Уa(t) =

a/i

J * &) sin (Ш/< +

ф ц )^1 ,

<-r

70(X) — функция Бесселя мнимого аргумента нулевого порядка. Подставляя (7.29) и (7.9) в (7.10), получим алгоритм оптималь­ ных УС, определяющих границы посылок по отрезку сигнала с длительностью Г и с неизвестной начальной фазой.

Синхронизация по нескольким посылкам. Увеличив длитель­ ность анализа сигнала в измерителе параметра, можно улучшить оценку параметра сигнала. Случай, когда интервал наблюдения ограничен (7+1)-й посылкой, отличается от рассмотренного ранее только тем, что сигнал представляет собой последовательность

/.....

, (A) = «(ft), S /ft— Т), • • -, S?f t — /Г),

и по аналогии с ф-лой (7.26) при известной форме сигнала

-

м

 

 

 

 

F,(x\t) = g

5 1

h

’ ' ■ * ^ Y ( ‘ .

U ■ ■ ■•?) X

 

i.

i.....?=i

 

 

 

 

X exp

 

 

ff( ' - ^

i

(7.31)

где s°j, j , .... q(t) — импульсная

реакция фильтра, согласованного с

'последовательностью

сигналов

s it

j ....q(t), t^t

(/+ 1 ) T. (При

пере­

даче независимых сигналов вероятность p(i, j, .... q), коэффициент y(i, j, .... q) и экспоненциальный сомножитель в (7.31) представ­ ляются в виде произведений, и фп границ посылок

F ' W - U P o i C x V -

i n

(7.32)

/=о

 

 

где фп Foj(x\tjT) вычисляется по ф-ле

(7.26)

для (/+ 1 —/')-й

посылки.

 

 

176

Представление фп в форме (7.32) обладает определенными до­ стоинствами. Например, устройство, вычисляющее логарифм фпт

скоторым часто удобнее работать, чем с самой фп, в соответствии

с(7.32) реализуется в виде вычислителя ln^of*]/) и накопителя, состоящего из линии задержки на время IT с отводами и сумма­

тора.

При неопределенной начальной фазе сигнала фп относитель­ но совокупности посылок представляется в виде произведения фп относительно отдельных посылок, если начальные фазы <pj0 на со­ седних посылках независимы. Если же начальные фазы зависимы, то фп более сложна [66].

7.3. Характеристики УС с ограниченной памятью

Характеристики УС, оптимальных при полностью известной форме сигнала. Будем считать, что фп в рассматриваемом случае определяется (7.32) и симметрична относительно положения гра­ ниц посылок. Тогда из (7.9) и (7.10) следует, что обе оптималь­ ные оценки совпадают и можно ограничиться исследованием оцен­ ки наибольшего правдоподобия. Так как максимумы фп и ее ло­ гарифма совпадают, то оптимальная оценка

 

t* =

arg шах 2

In Fol (х \ t — jT),

(7.33)

откуда

 

 

*

/-о

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.34)

 

 

%

± \ n F

o l C x \ F - j T ) =

0,

 

 

/*=о

 

 

как

моменты пересечений

т. e. оценки t* можно

рассматривать

вниз нулевого уровня процессом

 

 

 

 

 

 

 

/

 

 

 

 

(7.35)

 

У(*) “ S

Л- ln F°l (* I* _

/Т)'

 

 

 

 

 

/=о

 

 

 

 

 

Рассмотрим интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.36)

через

который в

 

a f

0

с (7.26)

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

процесс

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

F(x\t).

Представив

аналогично

(6.3)

входной сигнал

суммой

x(t) = z(t)-\-en(t), где

коэффициент е введен для удобства после­

дующих рассуждений,

и обозначив

 

 

 

Z l (t)=Zi (t, rt_ v r/, r u+) = -^-j‘Z( 0

1о

^ ( / - ^ ) ^ 1

(7.37)

t

Ni (0

af {

177

видим, что

 

Ji{t) = Zt (t, r(_ v rh гг+1) +sNi(t).

(7.38)

Здесь величины Z£ и Л/£ с точностью до постоянного коэффициен­ та совпадают с аналогичными величинами в гл. 6, а ги , ri, rt+l — информационные символы, переданные на трех соседних посыл­ ках. С учетом (7.26), (7.36), (7.38) ф-ла (7.35) принимает вид

мм

у^^Si lnSPiу‘ехр[Ziу~^ +8Ni^~уТ)] =

/=о г=1

/ 2 P /7 /exp[Z1(< - /T ) + 8/Vl ( / - / T ) l [ 2 i ( / - / T ) + e / / | ( < - / T ) ]

2 Pi Y« ехР lzi (< — /т ) + 8 (* — /Т)]

i=o

(7.39)

При больших отношениях сигнал/помеха можно считать пара­ метр е малой величиной и разложить (7.39) в ряд Тейлора по сте­ пеням е. Ограничиваясь первыми двумя членами разложения, имеем

 

y (t) = y{t,

e) = y(t,

0 ) +

в 0 ',(/, 0 ),

(7.40)

где

 

 

м

 

 

 

 

 

 

 

 

 

 

 

1

2

P m Z iV — iT)ex?iz i(t — iT)]

(7.41)

 

 

°)=%

*

4 i---------------------------

 

 

 

 

 

/=0

 

'2xPni**9[zt { t - m \

 

 

м

м

 

i=i

—2 Г М

 

 

 

 

 

у ' Л 0 ) = V

^ p {Y1 exp(Zi (^ — /Т))

\ т)

i=i

 

 

 

 

1 = 0

 

м

 

1=1

 

 

 

 

 

 

м

+ Wi (f — ;Т) Z, (/ — )Т)] ехр £ рк ykехр [Z* {t — /Т)] — £

Pi Y< X

 

 

 

 

*=1

 

i=i

 

 

 

 

м

 

 

\

x N i (t -

}Т) ехр [Z£(f - jT)) V р* у, Z* (* -

/Т) ехр [Z*(/ - /Т)] . (7.42)

 

 

 

 

А = 1

 

 

)

Если

n(i) — нормальный процесс, то при фиксированных гц ,

гг, гг+1 процессы Ni(t—jT) и

0) тоже нормальные. Известно

(см., например,

[85]),

что корреляционная функция

случайного

процесса N.-t(t) близка по форме к функции Z£(7) и монотонно убывает с ростом (tto). Нетрудно показать, что при этом кор­ реляционная функция процесса y'(t, 0 ) также монотонно убывает с ростом (t10).

178

Разложив (7.40) в ряд Тейлора в окрестности точки t0, прирав­ няв первые члены ряда нулю

y(i)zny (t 0) + гу'е (t0) + (/„) + е у\ (*„)] (/ — U

(7 -43)

где точка сверху означает производную по t, и (решив получен­ ное уравнение относительно t to, имеем

У (to) +

8 у е (to)

 

t * — tn = —

 

 

У (to) +

е у е (to)

 

Если у (to) “ 0, то

 

 

ry'e (to)

ey'g(to)

(7-44)

t* to —

У (<o)

У (to) + e у г (to)

 

Приняв e = l, находим математическое ожидание и дисперсию оценки')

 

<Уе Vo)> .

(7.45)

 

У (to)

 

 

< и ; м а >

 

(7.46)

У* (to)

 

 

 

Заметим, что здесь, как и в {85], можно уточнять статистиче­ ские характеристики оценок вторым и последующими приближе­

ниями.

Формулы (7.45) и (7.46) определяют характеристики условных оценок положения границ посылок, соответствующих конкретным реализациям сигнала. Если необходимо найти усредненные харак­ теристики оценок, то в ф-лы (7.45) и (7.46) следует подставлять среднее всех возможных реализаций сигнала.

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

видно из (7.42), математические ожидания процессов

у

е’

( t ) ,

Ni(t) равны нулю. Следовательно, оценки

границ посылок

с

по­

мощью синтезированных УС являются несмещенными

относи­

тельно точки, в которой y ( t , 0 ) = 0 .

этой точкой

является

При выводе (7.44) предполагалось, что

t0, чего в действительности может и не быть, как видно из рас­ смотренного ниже примера. Итак, оптимальное УС дает несме­ щенные оценки границ посылок только для таких сигналов, для

которых равенство y ( t ,

0 ) = 0

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

t,

равном точной

границе посылок.

в

силу

(7.46) обратно

пропорциональна

Дисперсия оценки

величине y * ( t o ) - Для

разных

реализаций сигнала

эти величины

<) Формулы (7.45) и (7.46) ораведливы, если t—/о < т к, где тк — время корре­ ляции процесса y(t).

179

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