492_Nosov_V._I.__Metody_povyshenija_pomekhoustojchivosti_sistem_radiosvjazi_..
._.pdf6.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
При экспоненциальной аппроксимации корреляционной функции
|
|
|
|
|
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
сравнению с 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