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

492_Nosov_V._I.__Metody_povyshenija_pomekhoustojchivosti_sistem_radiosvjazi_..

._.pdf
Скачиваний:
33
Добавлен:
12.11.2022
Размер:
6.31 Mб
Скачать

6.7.2 Помехоустойчивость приема

Средняя вероятность ошибки в ресурсном блоке при приёме сигналов

OFDM [66]

1 n

p n n i 1 pi,

где n – число поднесущих в ресурсном блоке; pi – вероятность ошибки на одной поднесущей. В статистически однородных каналах p n pi p.

Вероятность ошибки на одной поднесущей определяется вероятностью невыполнения неравенства (6.90) при передаче r-го сигнала. В системах с OFDM используется квадратурная амплитудно-фазовая модуляция m-QAM, где m = 2, 4, 16, 64. Двухпозиционную (BPSK) и четырёхпозиционную (QPSK) фазовые модуляции можно рассматривать как частный случай m-QAM с m = 2 и 4. Наиболее просто вероятность ошибки вычисляется при m = 2, т.е. для BPSK и при этом наиболее наглядно возможно проанализировать основные особенности приёма с использованием опорных (пилотных) сигналов, в частности, влияние скорости замираний сигналов на качество приёма. Поэтому анализ помехоустойчивости приёма проведём для данного вида модуляции (m = 2). Увеличение позиционности модуляции (m > 2) при заданной символьной скорости сопровождается увеличением битовой скорости в log2m раз ценой энергетического проигрыша в log2m раз.

Для сигналов с m = 2 (относительная фазовая модуляция при z2(t) = - z1(t)) алгоритм (6.90) преобразуется к виду

Q

 

mcq(M ) Xrq msq(M )Yrq 0.

(6.92)

q 1

Здесь согласно (6.91)

(M 1)T

(M 1)T

 

X

rq

 

 

Z

'

t Z

r

t dt,

Y

 

Z

'

t Zr t dt,

(6.93)

 

 

 

q

 

 

rq

 

q

 

 

 

 

 

MT

 

 

 

 

 

 

MT

 

 

 

 

– корреляционные интегралы, вычисленные при передаче первого варианта сигнала. Оценки mcq(M) и msq(M) вычисляются по формулам (6.88), (6.89). При этом

281

следует иметь в виду, что в (6.88), (6.89) при M > 1 формирование оценок должно проводиться по классифицированной обучающей выборке, т.е. со снятием манипуляции. Заметим, что опорные сигналы uq t манипулированы по фазе (BPSK)

псевдослучайными последовательностями Голда (вариант М-последовательности).

Вероятность ошибки найдём как вероятность невыполнения неравенства (6.92) при передаче первого варианта сигнала

0

 

 

p W x dx,

 

(6.94)

 

 

 

где

 

 

Q

 

 

x mcq(M )Xrq msq(M )Yrq

(6.95)

q 1

– при передаче сигнала z1(t).

При независимых релеевских замираниях сигналов в каждой ветви разнесения, левая часть неравенства (6.92) представляет собой квадратичную форму нормальных случайных величин с нулевым средним. Величины mcq(M) и msq(M) , X1q

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

 

m(M)2

m(M) X

 

 

m(M)2

m(M)Y

 

 

Kq

cq

 

cq

rq

sq

sq

rq .

(6.96)

 

 

 

 

 

 

 

 

 

 

 

X 2

 

 

Y2

m(M)X

rq

 

 

m(M)Y

 

 

 

 

cq

 

rq

 

 

sq rq

 

rq

 

 

 

Плотность вероятности квадратичной формы (6.95) определяется известными соотношениями [68]

Q

cq

exp

 

x

 

 

 

 

W x

 

 

,

x 0,

 

2

 

q 1

2

q

 

 

 

 

 

 

 

 

 

 

q

 

 

(6.97)

Q

dq

 

x

 

 

 

W x

 

 

exp

 

 

,

x 0,

2

 

2

 

q 1

q

 

 

 

 

 

 

 

q

 

 

282

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

Q

 

 

k

 

 

Q

k

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

1

 

 

 

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.98)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

k

 

 

Q

 

k

 

 

 

 

 

 

 

 

 

 

 

d

q

 

 

 

1

 

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

k 1

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

q

 

 

q

 

2q ;

 

 

 

 

и

 

 

 

 

 

 

 

 

собственные числа матрицы

 

 

 

1q ;

 

 

 

1q

 

2q

 

 

KqJq ,

где Kq – матрица ковариаций (6.96), Jq 0

1 – матрица квадратичной формы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

(6.95). Решая характеристическое уравнение

 

KqJq qI

 

0, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

(6.99)

 

 

 

 

 

 

 

 

m(M )

X

rq

 

 

m(M )2X2

 

 

 

 

 

 

 

 

 

1,2q

 

 

 

 

cq

 

 

 

 

 

 

 

 

 

cq

 

 

 

 

 

rq

 

 

После подстановки (6.97) в (6.94) и соответствующих преобразований, вероятность ошибки будет определяться формулой [78]

Q

2Q 1

 

Q

 

q

 

 

C2QQ 1 pq,

 

pî ø

 

 

 

 

(6.100)

Q

q k

Q

q k

q 1

 

q 1

 

 

k 1

 

k 1

 

 

 

 

 

 

k q

 

 

 

где pq q q q .

Определим элементы матрицы ковариаций (6.96). Из (6.88), (6.89) после попериодного усреднения произведений быстроменяющихся функций под знаком интеграла, следует

283

 

 

 

 

1

M

kT

 

 

 

 

mcq(N )

k 1

sq(k)(t)uq(t)dt

 

MT

 

 

 

 

 

 

 

 

(k 1)T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.101)

 

1

MT

(k)(t)dt

1

 

MT

(t)u (t)dt.

 

MT

 

MT

 

 

 

cq

 

q q

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

Аналогично из (6.93) и (6.79) получим

 

 

 

 

 

 

 

E (M 1)T

 

 

(M 1)T

 

X1q

 

 

 

 

mcq(t)dt

 

xq (t)zr (t)dt.

(6.102)

 

T

 

 

 

 

 

NT

 

 

 

NT

 

 

 

 

 

 

 

 

 

 

Согласно (6.101) и (6.102) элементы ковариационной матрицы (6.96) определяются выражениями

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(M )2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

h

r

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cq

 

 

 

 

 

2

 

 

 

 

q

 

1q

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

h

r

 

 

 

 

 

 

,

(6.103)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rq

 

 

 

 

2

 

q

 

 

2q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2E

 

r

 

.

 

 

 

 

 

 

 

 

 

 

 

mM X

 

 

h2

 

 

 

 

 

 

 

 

 

 

 

 

cq

 

rq

 

 

2

 

 

 

q

3q

 

 

 

 

 

 

 

 

2 E

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь hq

2 q

 

 

2 cq

 

 

 

 

среднестатистическое значение отношения

2

2

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

энергии принятого сигнала в q -ой ветви разнесения к спектральной плотности мощности шума;

284

 

 

1

 

 

 

MT

 

 

 

 

 

 

 

 

 

 

 

 

r1q

 

 

R1q (t1,t2 )dt1dt2;

 

 

 

 

 

 

(MT)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

(M 1)T

 

 

 

 

 

 

 

 

 

 

 

 

r2q

 

R2q (t1,t2 )dt1dt2;

 

 

 

 

(6.104)

T2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

MT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

NT (M 1)T

 

 

 

 

 

 

 

 

 

 

 

r3q

 

 

 

R3q (t1,t2 )dt1dt2.

 

 

 

 

 

MT2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

MT

 

 

 

 

 

 

 

 

 

 

 

 

В формулах (6.104)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k1)

(t ) (k1)

(t

2

)

 

(k1)

(t ) (k1)

(t

2

)

 

R

(t ,t

2

)

 

cq

1

cq

 

 

 

sq

1

sq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1q

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

q

 

 

 

 

 

корреляционная функция ортогональных составляющих опорного сигнала

вq-ой ветви разнесения;

 

 

 

 

 

(k2)

(t )

(k2)

(t

2

)

 

 

(k2)

(t )

(k2)

(t

2

)

R

(t ,t

2

)

 

cq

1

 

cq

 

 

 

 

sq

1

 

sq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2q

1

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

q

 

 

 

 

– корреляционная функция ортогональных составляющих информационного сигнала;

 

 

 

 

 

(k1)

(t )

(k2)

(t

2

)

 

 

(k1)

(t )

(k2)

(t

2

)

R

(t ,t

2

)

 

cq

1

 

cq

 

 

 

 

sq

1

 

sq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3q

1

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

q

 

 

 

 

– взаимная корреляционная функция ортогональных составляющих опорного и информационного сигналов. Возможна экспоненциальная и гауссовская аппроксимация Rq(t1,t2) [68].

285

= 7,14·10-3.

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

 

 

 

 

 

t ,t

 

 

 

 

 

 

t

t

 

 

 

 

t

 

t

 

 

 

 

 

 

 

 

R

exp

 

 

 

1

2

 

 

1

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

1 2

 

 

 

 

 

 

kq

 

 

 

 

 

kq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выражения (6.104) примут вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

1

MT

;

 

r

1

T

;

 

 

r

 

 

 

1

(M 1)T

.

 

 

 

 

 

 

 

 

 

 

 

1q

 

3 k1

 

 

2q

 

3 k2

 

 

3q

 

 

 

 

 

3 k3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь ki

i 1,2,3 интервал

корреляции изменений параметров сигнала.

В общем случае они могут быть разными для опорного и информационного сигналов. Но при общих замираниях группы поднесущих частот ресурсного блока,

содержащего опорный сигнал k1 k2 k3 k . Например, при скорости движе-

ния абонента = 50 км/час интервал корреляции замираний k составляет величину порядка 10 мс (по данным работ [58, 69]). В сотовой сети LTE длительность OFDM-символа T вместе с префиксом составляет 71,4 мкс, тогда T k

При гауссовской аппроксимации корреляционной функции

 

 

Rq t1,t2 exp

 

 

t

 

t

2

2

 

 

t

t

2

2

 

 

 

 

1

 

 

 

1

1

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 kq

 

 

 

 

 

2 kq

 

 

 

выражения (6.104) примут вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

1 MT 2

;

r

1

T2

;

 

 

r

 

1

T2

 

2M 2 3M 2 .

 

 

 

 

 

 

1q

12 2

 

2q

12 2

 

 

 

 

3q

 

12

2

 

 

 

 

k

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

k

 

 

 

Гауссовская аппроксимация более благоприятна для прогнозирования параметров канала. Многие авторы считают, что в большинстве случаев на линиях связи корреляционная функция близка к гауссовской кривой. Её в дальнейшем и будем использовать при анализе.

286

hq2r3q

Определим собственные числа матрицы KqJq подстановкой (6.103) в (6.104), получим

 

2E

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1,2q

 

 

hq2r3q

hq2r1q

 

hq2r2q 1 .

(6.105)

 

 

 

 

 

2

 

 

M

 

 

 

 

 

 

 

 

 

 

 

После подстановки q 1q и q 2q из (6.105) в (6.104), получим выражение для вероятности ошибки при разнесённом приёме OFDM сигналов

 

 

 

 

 

 

 

 

 

 

 

C2QQ 1

Q

 

pî ø

 

 

1

2

Q

 

 

q 1

 

 

 

 

 

 

. (6.106)

hq2r1q 1 hq2r2q 1

M

Из (6.106) видно, что вероятность ошибки зависит от отношения энергии

символа к спектральной плотности мощности шума hq2 , от размера обучающей

выборки M и от коэффициентов, характеризующих скорость изменения параметров сигнала r1q,r2q,r3q .

На рис. 6.8 приведены зависимости вероятности ошибки p от размера обучающей выборки M при разных Q, вычисленные по формуле (6.106), в предполо-

жении hq2 h2 102 (одинаковая средняя интенсивность сигналов разных ветвей разнесения), при гауссовской аппроксимации Rq t1,t2 при T k = 10-2 (сплошные линии) и при T k = 5·10-2 (штриховые линии). Здесь же для сравнения приведены зависимости потенциальной помехоустойчивости при r1q r2q r3q 1

(штрихпунктирные линии), что практически соответствует значению T k ≈ 10-3 .

Из приведённых зависимостей видно, что при T k = 10-2 заметное уменьшение вероятности ошибки происходит на первых двух-трёх интервалах обучения, затем уменьшение вероятности ошибки замедляется и при M > 7 вероятность ошибки будет возрастать вследствие декорреляции измеренных параметров сигнала. С увеличением M повышение помехоустойчивости приёма происходит более эффективно при большом числе ветвей разнесения. В рассматриваемых условиях значение T k = 10-2 можно считать количественной оценкой (мерой) «медленности» замираний при гауссовской аппроксимации Rq t1,t2 . При

T k = 5·10-2 вероятность ошибки при M = 2 практически не уменьшается по

287

сравнению с M = 1, а при M > 2 вероятность ошибки возрастает вследствие декорреляции измеренных параметров сигнала. Величина T k = 5·10-2 близка к границе, характеризующей быстрые замирания, когда обучение на интервале M > 1 нецелесообразно.

Рис. 6.8. Зависимости вероятности ошибки от размера обучающей выборки (гауссовская аппроксимация Rq t1,t2 )

Монотонное уменьшение вероятности ошибки с увеличением M наблюдается только при r1q r2q r3q 1 (практически при T k ≈ 10-3), однако, и в этих

условиях целесообразно осуществлять обучение на интервале не более двух трёх элементов обучающей выборки, т.к. далее уменьшение вероятности ошибок сильно замедляется.

Заметим, что в системах технологии LTE значение T k = 10-2 будет при скорости движения абонента ≈ 70 км/ч. Технологией LTE предусматривается

288

мобильность до скорости = 350 км/ч, при этом T k = 5·10-2. Значение T k ≈10-3 соответствует движению абонентской станции со скоростью пешехода.

Из выражения (6.106) видно, что в каналах с конечной скоростью изменения

параметров сигнала при неограниченном увеличении величины hq2 вероятность ошибки стремится не к нулю, а к предельному значению

 

 

cQ

Q

 

r

 

 

 

 

p

 

2Q 1

 

 

3q

 

 

 

 

 

 

1

 

 

 

.

 

 

2Q

 

 

 

 

 

ош пр

 

 

r r

 

 

 

 

 

 

 

q 1

 

1q 2q

 

 

Например, при одиночном приёме (Q = 1) при T k

= 10-2 предельная вероят-

ность ошибки при экспоненциальной аппроксимации R

t ,t

равна pср = 3,4·10-3,

 

 

 

 

 

 

 

 

q

1 2

 

а при гауссовской аппроксимации pср = 2,5·10-5. Предельная вероятность ошибки стремится к нулю только в каналах с неизменными параметрами (r1q = r2q = r3q = 1). В таких каналах предельный выигрыш по вероятности ошибки при M >> 1 по сравнения с M = 1 составляет 2Q, реально он меньше.

На рис.

6.9

приведены зависимости вероятности ошибки р от

hq2 h2

при разных Q =1, 2, 3, 4 при обучении на интервале M = 1 при гауссовской

аппроксимации

R

t ,t

при T

k

= 10-2 (сплошные линии) и при T

k

= 5·10-2

 

q

1 2

 

 

 

(штриховые линии).

Штрихпунктирными линиями для сравнения показаны зависимости вероятности ошибки когерентного разнесённого приёма (потенциальная помехоустойчивость) фазоманипулированных сигналов в каналах с релеевскими замираниями. Зависимости p = f (h2) между M = 1 и M >> 1 пройдут на расстоянии (1+1/M)Q , что при M = 2…3 приводит к увеличению вероятности ошибки по сравнению с потенциальной помехоустойчивостью в (1,3…1,5)Q раз.

Из рис. 6.9 видно замедление в уменьшении вероятности ошибки с увеличе-

нием h2 вследствие декорреляции параметров сигнала из-за конечной скорости замираний, однако, адаптивный приём и в таких условиях обеспечивает высокую достоверность связи в системах с OFDM и даже при высокой скорости движения абонента (в LTE до v = 350 км/ч).

Таким образом, приём с обучением по опорным (пилотным) сигналам, в

каналах с гауссовской корреляционной функцией замираний Rq t1,t2 незначительно проигрывает оптимальному приёму при медленных (по сравнению с

289

длительностью OFDM-символа) общих замираниях, которые наблюдаются в каналах функционирования систем WiMAX, Wi-Fi, LTE.

Рис. 6.9. Зависимости вероятности ошибки от отношения сигнал/шум (гауссовская аппроксимация Rq t1,t2 )

Для полноты картины на рис. 6.10 приведены зависимости вероятности

ошибки р от h2 при разных Q =1, 2, 3, 4 при обучении на интервале M = 1 при экспоненциальной аппроксимации Rq t1,t2 , характерной для условий распро-

странения сигналов в офисных помещениях [66]. Зависимости построены при T k = 10-2 (сплошные линии) и при T k = 10-3 (пунктир). Штрихпунктирными линиями показаны зависимости вероятности ошибки когерентного разнесённого приёма (потенциальная помехоустойчивость). Из рис. 6.10 видно, что при экспо-

ненциальной аппроксимации Rq t1,t2 , замедление в уменьшении вероятности ошибки при увеличении h2 происходит гораздо быстрее, чем при гауссовской

290