492_Nosov_V._I.__Metody_povyshenija_pomekhoustojchivosti_sistem_radiosvjazi_..
._.pdfгде
|
|
ˆcik |
|
1 |
|
kT tn |
zk t zr t ti dt, |
|||
|
|
|
|
|
||||||
|
E |
|||||||||
|
|
|
|
|
|
i |
k 1 T tn |
|
||
|
|
ˆsik |
1 |
|
kT tn |
(6.39) |
||||
|
|
|
zk t zr t ti dt |
|||||||
|
|
|
|
|||||||
|
E |
|
||||||||
|
|
|
|
|
|
i k 1 T tn |
||||
– максимально-правдоподобные оценки параметров сигнала на k-м этапе |
||||||||||
наблюдения в i-м луче; ˆi2 |
|
2 |
|
– дисперсия этих оценок. |
||||||
2Ei |
||||||||||
|
|
|
|
|
|
|
|
Оценки параметров сигнала согласно (6.38) определяются сверткой принятой смеси и образца сигнала с задержкой на время запаздывания i-го луча. При этом необходимо снятие манипуляции с рабочих посылок сигнала в реализа-
ции (6.36).
Апостериорное распределение вероятностей параметров сигнала с учетом полного разделения лучей можно представить в виде
n |
|
W(μc,μs /z ) W( ci, si /z ). |
(6.40) |
i 1 |
|
Будем полагать, что в обучающей выборке (6.36) последовательность пара- |
метров ci1 , , ciN образуют цепь Маркова с плотностью вероятностей перехода
s s
|
k 1 |
k |
|
W |
ci |
/ ci |
. |
|
s |
s |
|
|
1 |
|
Полагая априорное распределение W |
ci |
равномерным по формуле Байе- |
|
s |
|
са находим апостериорное распределение параметров сигнала с учетом марковских свойств последних [46]
|
|
|
|
|
, ,zNi |
N 1 |
|
|
/ ci(k) |
|
. (6.41) |
WN ci(1), , ci(N) /z1i |
, ,zNi |
constWN z1i |
/ ci W ci(k 1) |
|
|||||||
|
s |
s |
|
|
|
s k 1 |
|
s |
s |
|
|
262
Искомая апостериорная плотность вероятностей на N-м этапе наблюдения
получается из (6.40) (N-1)-кратным интегрированием по ci1 |
, , ciN 1 |
|
||||||
|
|
|
|
|
|
s |
s |
|
|
|
|
|
|
|
|
|
|
WN |
ci(N) /z1i, ,zNi |
W ci(1), , ci(N) /z1i, ,zNi |
d ci(1), ,d ci(N 1). (6.42) |
|||||
|
s |
|
|
s |
s |
|
s |
s |
При поступлении нового наблюдения z N 1 i плотность вероятностей (6.42)
преобразуется следующим образом
|
(N 1) |
|
, ,z |
|
|
|
|
(N 1) |
|
W(N 1) |
ci |
/z1i |
N 1 i |
constW(N 1) z |
N 1 i |
/ ci |
|
||
|
s |
|
|
|
|
|
|
s |
|
(6.43)
|
|
/ ci(N) |
|
|
|
|
|
W ci(N 1) |
WN ci(N) /z1i, ,zNi d ci(N). |
||||||
|
s |
s |
|
|
s |
|
s |
Это выражение определяет основной рекуррентный алгоритм для нахождения апостериорной плотности вероятностей параметров сигнала. Подставляя (6.39) с учетом (6.40) и (6.36) в (6.31) и интегрируя, можно получить адаптивный алгоритм работы приемника с обучением по сигналу в каждом конкретном случае. Для практических целей в условиях высокой точности измерения параметров сигнала достаточно ограничиться гауссовским приближением (6.43). Тогда целесообразна аппроксимация переходной плотности вероятностей формулой [65]
|
|
|
|
|
|
|
|
|
|
(N 1) |
|
(N)R2 |
|
|
|
|
|
|
1 |
|
|
|
|
|
ci |
|
ci |
i |
|
|
|
|
|
|
|
|
s |
|
s |
|
|
||||||
W ci(N 1) |
/ ciN |
|
|
|
|
exp |
|
|
|
|
|
|
|
. |
(6.44) |
|
|
|
|
2 |
|
2 |
) |
||||||||
|
|
|
|||||||||||||
s |
s |
|
2 i2(1 Ri2) |
|
|
|
2 i |
(1 Ri |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
После подстановки в (6.43) выражений (6.38) и (6.44), получим апостериорную плотность вероятностей на N-ом этапе [7]
N |
N |
|
1 |
|
|
1 |
|
N |
|
2 |
N |
|
2 |
|
|
|
|
|
|
|
|
|
|
||||||||
WN ci |
, si |
/z1i, ,zNi |
|
|
exp |
|
|
ci |
mci |
|
|
si msi |
|
|
, (6.45) |
2 |
2 |
2 |
2 |
|
|
||||||||||
|
|
|
Ni |
|
Ni |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
263
параметры которой определяются рекуррентными соотношениями
|
m |
(N 1) |
|
2 |
R |
ˆ |
(N) |
|
|
2 |
|
|
|
R |
2 |
|
2 |
(1 R |
2 |
) |
|
|
||||||||||
|
ci |
|
ˆ |
i |
|
ci |
|
|
|
|
|
|
|
i |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
i |
|
|
|
|
(N 1)i i |
|
|
|
|
|
|
i |
|
|
||||||||||||
mci(N) |
s |
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
(6.46) |
|
|
|
|
|
ˆ2 |
2 |
|
|
R2 |
2(1 R2 ) |
|
|
|
|
|||||||||||||||||||
s |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
i |
(N 1)i |
|
i |
|
|
|
i |
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
ˆ |
2 |
2 |
R |
2 |
|
2 |
(1 R |
2 |
) |
|
|
|
|
|
|
|
||||||||||||
2 |
|
|
i |
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
(N 1)i i |
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
||||||||||||
Ni |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
(6.47) |
||||||||||||||||
ˆ2 |
|
2 |
1)i |
R2 |
|
2 |
(1 R2) |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
i |
|
|
|
(N |
|
|
i |
|
|
i |
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
Здесь i2 и Ri – параметры (дисперсия и коэффициент корреляции) пере-
ходной плотности вероятностей величин ci, si .
Начальные условия
|
mci(1) |
ˆci(k), |
12i |
ˆi2. |
|
s |
s |
|
|
ˆci(k) |
– определяются выражениями (6.39). |
|
||
s |
|
|
|
|
Из (6.46) и (6.47) следуют частные случаи.
1.В каналах с постоянными параметрами сигнала ( Ri 1)
|
1 |
N |
|
2 |
|
|
|
mci(N) |
ˆci(k), |
Ni2 |
|
i |
. |
(6.48) |
|
|
|
||||||
s |
N k 1 s |
|
|
N |
|
По мере увеличения размера обучающей выборки параметры определяются все точнее.
2.В каналах с быстрыми изменениями параметров сигнала ( Ri 0)
|
ˆ(N)h2 |
|
|
|
|
|
|
|
|
|||
|
|
2 |
2 |
|
|
|
||||||
|
ci |
i |
|
|
|
|
||||||
m(N) |
s |
|
|
, |
2 |
|
ˆi |
|
hi |
|
, |
(6.49) |
|
|
|
|
|
|
|
||||||
ci |
1 hi2 |
Ni |
1 hi2 |
|||||||||
s |
|
|
264
|
|
|
2 |
|
2 2E |
|
|
|
h2 |
|
|||||
где |
|
|
i |
|
i i |
– среднестатистическое значение отношения энергии |
|
|
|
|
|||||
i |
ˆi2 |
2 |
|||||
|
|
|
принятого сигнала в i-м луче к спектральной плотности шума.
При hi2 1 соотношения (6.49) примут вид
mciN |
ˆci(N), |
Ni2 |
ˆi2, |
(6.50) |
s |
s |
|
|
|
что соответствует (6.49) при N = 1.
Из (6.49) и (6.50) видно, что при быстрых изменениях параметров сигнала формирование оценок происходит независимо на каждом этапе наблюдения. В реальных каналах имеют место достаточно медленные изменения параметров сигнала ( Ri 0,999). Это обеспечивает возможность формирования оценок по правилу (6.38) для достижения высокой достоверности приема.
После подстановки в (6.31) выражений (6.36) и (6.45) с учетом (6.40) и интегрирования получим адаптивный алгоритм приема широкополосных сигналов в n-лучевом канале в общем виде
n |
|
|
|
Yri mciN |
2 |
|
2 |
Eri |
|
|
|
|
2 |
|
|
|
|
||
2 mciN Xri msiN |
msiN |
|
Xri2 Yri2 |
Ni |
|
|
|
||||||||||||
|
2 |
|
|
|
|||||||||||||||
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆi |
Eri |
|
(6.51) |
||
|
n |
|
|
|
|
2 |
|
|
2 |
|
|
|
|
|
|
2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
2 mciN Xli msiN Yli mciN |
msiN |
|
Eli Xli2 Yli2 |
Ni |
, |
r l. |
||||||||||||
|
2 |
||||||||||||||||||
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆi |
Eli |
|
|||
Здесь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
(N 1)T tn |
z t zr t ti dt, |
|
|
|
|
|
(N 1)T Δtn |
z t zr t ti dt. (6.52) |
||||||||
|
|
Xri |
|
|
Yri |
|
|
||||||||||||
|
|
|
NT tn |
|
|
|
|
|
|
|
|
NT Δtn |
|
|
|
|
|
|
|
При использовании сигналов с одинаковыми энергиями, правило приема (6.51) существенно упрощается
n |
n |
|
|
mci(N) Xri |
msi(N)Yri mci(N)Xli msi(N)Yli , |
r l, |
(6.53) |
i 1 |
i 1 |
|
|
265
6.6Применение шумоподобных сигналов с нулевой зоной корреляции
вмобильных системах MIMO. Анализ алгоритмов
Проведем анализ помехоустойчивости приема сигнала (6.28) по алгоритму (6.53) с учетом мешающих факторов традиционной системы CDMA, т.е. при влиянии интерференционных помех. Из полученного выражения для вероятности ошибки помехоустойчивость системы LAS-CDMA будет вытекать как частный случай.
В рамках концепции MIMO возможны виды модуляции: однократная относительная фазовая (BPSK), квадратурная относительная фазовая (QPSK), с минимальным сдвигом (GMSK), квадратурная амплитудная (m-QAM) и их разновидности.
Вначале проанализируем помехоустойчивость приема двоичных противоположных сигналов z2 t z1 t , т.е. сигналов с относительной фазовой манипу-
ляцией. Для таких сигналов правило приема (6.53) преобразуется к виду
|
|
n |
msi(N ) Xri 0. |
|
||
|
|
mci(N ) Xri |
(6.55) |
|||
|
|
i 1 |
|
|
|
|
Здесь |
|
|
|
|
|
|
|
(N 1)T tn |
z t zr t ti dt, |
|
(N 1)T tn |
z t zr t ti dt. |
|
Xri |
|
Yri |
|
(6.56) |
||
|
NT tn |
|
|
NT tn |
|
|
В неравенство (6.41) входят mciN и msiN , которые с учётом формулы (6.40), равны
(N) |
|
1 |
N |
k |
|
(N) |
|
1 |
N |
k |
|
|
mci |
|
|
|
ˆci |
; |
msi |
|
|
|
ˆsi |
, |
(6.57) |
|
|
|||||||||||
|
|
N k 1 |
|
|
|
|
N k 1 |
|
|
|
267
в которых
ˆcik |
|
1 |
kT tn |
zk t z1 t ti dt, |
|
|
|
||||
E |
|||||
|
|
i k 1 T tn |
(6.58) |
||
|
|
|
kT tn |
||
ˆsik |
|
1 |
zk t z1 t ti dt. |
||
|
|
||||
E |
|||||
|
|
i k 1 T tn |
|
n
zk ciz1 t ti siz1 t ti
i 1
(6.59)
M1
ciqzq t ti siqzq t ti t .
q 1
В(6.59), n – общее число лучей в многолучевом канале, воздействующих на вход приемника. Оно может отличаться от n в (6.55), где n – число используемых для обработки лучей.
Вканалах с переменными случайными параметрами прием противоположных сигналов, как известно, можно реализовать системами с измерением параметров канала на предыдущем элементе или на N предыдущих элементах. Измерение параметров канала при N 1 должно проводиться по классифицированной обу-
чающей выборке, т.е. со снятием манипуляции, поэтому mc(1Ni) mc(N2i) mci(N) .
s s s
Вероятность ошибки найдем как вероятность невыполнения неравенства (6.55) при передаче варианта сигнала
0 |
|
p W x dx, |
(6.60) |
|
|
где |
|
n |
|
x mciN Xli msiN Yli |
(6.61) |
i 1
– при передаче сигнала z1 t .
268
При независимых релеевских замираниях сигналов в лучах, при полном их разделении, левая часть неравенства (6.55) представляет собой квадратичную форму нормальных случайных величин с нулевым средним. Величины mci(N) и
msi(N), X1i и Y1i попарно независимы с матрицей ковариаций
K |
|
|
m N 2 |
|
|
m N 2X |
li |
|
|
m N 2 |
|
m N 2Y |
|
|
|||||
i |
|
|
ci |
|
|
|
ci |
|
|
|
si |
|
|
si |
li |
. |
(6.62) |
||
|
N |
|
|
X 2 |
|
|
N |
Y2 |
|
||||||||||
|
|
m |
X |
|
|
|
|
|
m |
|
|
|
|
||||||
|
|
|
|
li |
|
li |
|
|
|
si |
Y |
|
li |
|
|
|
|||
|
|
ci |
|
|
|
|
|
li |
|
|
|
|
Плотность вероятности квадратичной формы (6.61) определяется известными соотношениями [28]
n |
ci |
|
|
|
x |
|
|
|
|
|
W x |
|
exp |
|
|
|
, |
x 0, |
|
||
2 i |
|
|
||||||||
i 1 |
|
|
|
2 i |
|
, |
(6.63) |
|||
n |
di |
|
|
|
x |
|
|
|
||
|
|
|
|
|
|
|||||
W x |
|
exp |
|
|
|
, |
x 0, |
|
||
2 i |
|
|
|
|
||||||
i 1 |
|
|
|
2 i |
|
|
|
где
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
n |
|
|
k |
|
|
n |
|
k |
|
|
||||
ci 1 |
|
|
|
|
1 |
|
|
|
, |
||||||
|
|
|
|
|
|
||||||||||
k 1 |
i k 1 |
i |
|||||||||||||
|
|
|
|
|
|
|
|
|
k i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(6.64) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
n |
|
|
k |
|
|
n |
|
|
k |
|
|
|||
di 1 |
|
|
|
|
|
1 |
|
|
|
|
. |
||||
|
|
|
|
|
|
|
|
|
|
||||||
k 1 |
|
i k 1 |
|
i |
|||||||||||
|
|
|
|
|
|
|
|
|
k i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
269