Программные продукты и системы

Многоантенные системы связи, в которых участники общения используют несколько передающих и принимающих антенн, позволяют передавать данные с высокой скоростью и обеспечивать их надежный прием. Помехоустойчивое пространственное кодирование, применяемое в таких системах связи в качестве канального кодирования, базируется на внесении избыточности во временную, частотную и пространственную области.

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

В литературе уже описаны алгебраические слоеные пространственные коды, использующие, помимо слоения, свойства алгебраических чисел [3, 4], показано, что для некоторых из таких кодов существуют порождающие матрицы [5].

В ряде работ были представлены позиционирующие функции, позволяющие конструировать кодовые слова слоеных пространственных кодов [6, 7].

В данной статье представлен слоеный пространственный код, удовлетворяющий критериям для проектирования и ориентированный на каналы с медленными, быстрыми или блоковыми рэлеевскими замираниями. Особенностью кода является то, что в нем отсутствует символьная избыточность: число элементов в кодовом слове совпадает с числом элементов в кодируемом им векторе.

Модели каналов связи

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

где iÎ{1, 2, …, NRx} – номер принимающей антенны; jÎ{1, 2, …, NTx} – номер передающей антенны; mÎ{1, 2, …, Nt}  – номер посылки; hi,j "(i, j) – независимые комплексные гауссовы случайные величины с математическим ожиданием, равным нулю, и дисперсией, равной единице; cj(m) – элемент кодового слова – комплексное число; wi(m) " (i, m) – независимые комплексные гауссовы случайные величины с математическим ожиданием, равным нулю, и дисперсией, равной единице; r – коэффициент, управляющий отношением сигнал/шум в модели канала. Нормирующий коэффициент NTx обеспечивает одинаковое количество мощности, излучаемой каждой передающей антенной.

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

где hi,j(m) " (i, j, m) – независимые комплексные гауссовы случайные величины с математическим ожиданием, равным нулю, и дисперсией, равной единице.

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

где kÎ{1, 2, …, Nc} – номер частотного подканала;

где ; L – число поступающих на приемник рассеянных реплик переданного сигнала; ai,j,l(m) " (i, j, l, m) – независимые комплексные гауссовы случайные величины с математическим ожиданием, равным нулю, и дисперсией, равной , при этом ; Ts – длительность передачи одной посылки; tl – задержка распространения сигнала.

Пространственное кодирование

Организуя передачу сигналов в многоантенных системах связи, необходимо указывать номер передающей антенны, номер частотного подканала (при многочастотной передаче) и номер посылки для каждого передаваемого сигнала. Таким образом, если нет многочастотной передачи, каждый передаваемый сигнал увязывается с номером передающей антенны и номером посылки. Если же применяется многочастотная передача, передаваемый сигнал увязывается с номерами передающей антенны, частотного подканала и посылки. Эти пары и тройки номеров называются позиционными парами и тройками соответственно.

Чтобы задать пространственное кодовое слово, требуется указать множество элементов кодового слова и позиционных пар (троек), а также правила соответствия элементов этих множеств друг другу. Пространственные кодовые слова представляются матрицами или векторами. В этих матричных или векторных формах каждый элемент – это элемент кодового слова. Ячейка в такой форме увязана с номерами передающей антенны, частотного подканала и посылки. Элемент кодового слова, заняв ячейку, увязывается с этими номерами, иными словами – с позиционной парой (или тройкой).

Порождающая матрица пространственного кода выполняет линейное преобразование кодируемого вектора в кодовое слово, тем самым выполняя сразу два действия для создания пространственного кодового слова – получение элементов кодового слова и увязывание их с позиционными парами (тройками). Порождающая матрица пространственного кода – это произведение из двух матриц: матрицы для преобразования и матрицы для выборки. Матрица для преобразования выполняет линейное преобразование кодируемого вектора в вектор из элементов кодового слова; таким образом осуществляется первое действие, требуемое для создания пространственного кодового слова, – получение элементов кодового слова. Матрица для выборки составляется из нулей и единиц так, что на каждой строке только одна единица, остальные элементы – нули. Каждой строке матрицы для выборки сопоставлены номера передающей антенны, частотного подканала и посылки. Линейное преобразование, выполняемое над вектором из элементов кодового слова, – это второе и завершающее действие, необходимое для конструирования пространственного кодового слова, – увязывание элементов кодового слова с номерами антенн, подканалов и посылок.

Слоеный пространственный блоковый код

Предположим, что NTx – число передающих антенн, Nc – число частотных подканалов (в случае, когда применяется многочастотная передача). Требуется за Nt посылок передать вектор = =(s1 s2 … )T, где si (iÎ{1, 2, …, NTxNcNt}) – комплексные числа, мнимая и вещественная части которых – целые числа; (·)T – оператор транспонирования. Одна посылка – это время передачи NTxNc элементов, которое может включать заградительный интервал, как, например, c использованием OFDM-технологии.

Представленный в данной статье слоеный пространственный код имеет порождающую матрицу. Для задания порождающей матрицы используется несколько численных параметров, принимающих значения из множества целых положительных чисел. Эти параметры обозначены символами J, Nq, NL, Nb. Требуется, чтобы между величинами NTx, Nc, Nt и параметрами J, Nq, NL, Nb выполнялись следующие соотношения. Прежде всего необходимо, чтобы Nq³NTx. Если Nc=1 (нет многочастотной передачи), параметрам требуется задать такие целые положительные численные значения, при которых выполнялось бы равенство JNqNLNb= Nt; код при таком значении Nc называется пространственно-временным (ПВ). Если же Nc>1, то параметрам задаются такие значения, которые удовлетворяют равенствам JNqNL=Nc и Nb=Nt; код при таком значении Nc называется пространственно-частотно-временным (ПЧВ).

Матрица для преобразования обозначена символом T и представляет собой блоково-диагональ­ную матрицу

где f – алгебраический элемент; Ä – оператор кронекерова произведения; IJ – единичная матрица размером J´J;

где N=NTxNLNb; если N=2n, то

где (·)H – оператор эрмитова сопряжения; F(N) – матрица прямого преобразования Фурье без нормирующего коэффициента, размер которой N´N; если же N¹2n, то в качестве Q(N) используется первая главная подматрица размера N´N из матрицы  где écù – ближайшее целое, которое больше или равно c.

Матрицы Q(N) называют комплексными матрицами вращения, в [8] описаны варианты таких матриц. Эти матрицы совместно с f выполняют преобразование части кодируемого вектора в слои кодового слова. Величиной N задается число элементов в слое, а произведением JNq – число слоев в кодовом слове.

Величина f – это алгебраический элемент степени как минимум NLNqNb над полем F, которое представляет собой такое расширение поля  (поле рациональных чисел), что содержит все элементы матрицы Q(N), информационный алфавит  (комплексные числа, мнимая и вещественная части которых – целые числа) и величины exp(–I×2ptl/Ts). Одним из возможных вариантов значений для этой величины является  где p – целое число, большее или равное NLNqNb, а n – наибольшее общее кратное элементов множества {4N, s1, s2, …, sL}, где sl (lÎ{1, 2, …, L}) определяются из пропорции  при условии, что  и что наибольший общий делитель rl и sl равен 1. Кроме того, в качестве f можно использовать трансцендентное над полем F число. В [9] приведены примеры вычисления этой величины.

Матрица для выборки обозначена символом P и составляется с применением позиционирующих функций x и z:

где

где

и

где iÎ{1, 2, …, NTxNcNt }; [a] – целая часть числа a;   и  где ëcû – ближайшее целое, меньшее или равное c.

Порождающая матрица слоеного пространственного кода – это произведение из матриц T и P:

G(J, Nq, NL, Nb, NTx, f)=P×T.

Чтобы закодировать вектор  слоеным кодом, достаточно выполнить линейное преобразование этого вектора порождающей матрицей G:

=G(J, Nq, NL, Nb, NTx, f)×,

где – кодовое слово слоеного пространственного кода. Если Nc=1, то

где cj(m) – элемент ПВ кодового слова, передаваемый j-й передающей антенной в течение m-й посылки. Если же Nc>1, то

где cj,k(m) – элемент ПЧВ кодового слова, передаваемый j-й передающей антенной по k-му частотному подканалу в течение m-й посылки.

Слоеный пространственный код, в котором NTx=Nq=A, J=1, NL=1, Nb=B, по свойствам подобен коду, предложенному в [3], а код, в котором J=A, NLNq=2B, Nb=C, – предложенному в [4].

Позиционирующие функции

Представим функции, через которые определяются позиционирующие функции.

Функция  где ; écù – ближайшее целое, которое больше или равно c; [c] – целая часть (вместе со знаком) действительного числа c. Значение функции goe(a, b) либо равно 1, если ½a½³½b½, либо равно 0, если ½a½<½b½.

Функция

где

Функция

где  и

Функция

где  и

Функция

где

Функция  где ;  – элемент единичной матрицы, расположенный на a-й строке и в b-м столбце, то есть  и

Функция  

где  и

Через эти функции вычисляются позиционирующие функции, которые определены для случаев, когда Nq³NTx. Число NL может быть положительной рациональной дробью, однако знаменатель этой дроби должен быть делителем и Nq, и Nb. Позиционирующие функции – это пара функций x и z:

и

где

и

где

Примеры оценок помехоустойчивости слоеного пространственного блокового кода

Пространственные коды – это относительно новый инструмент, обещающий высокую помехоустойчивость беспроводной связи. Далее приведены несколько примеров слоеных пространственных кодов и графики зависимости вероятности ошибки декодирования по минимуму евклидова расстояния от отношения сигнал/шум на бит. Для декодирования кодов использовался сферический декодер [10, 11].

Для примера сравним помехоустойчивость двух слоеных пространственных кодов, а именно ПВ-кодов. Первый ПВ-код используется для случая, когда NTx=2, Nc=1, Nt=2 и кодируются векторы  где . Параметры кода: J=1, Nq=2, NL=1, Nb=1. Порождающая матрица первого кода G(1, 2, 1, 1, 2, f)=P×T, где ,  

и

где

Кодовые слова первого кода: =G(1, 2, 1, 1, 2, exp(I×p/8))×, где =(c1(1)  c2(1)  c1(2)  c2(2))T, где cj(m) – элемент кодового слова, передаваемый j-й передающей антенной в течение m-й посылки.

Второй ПВ-код используется, когда NTx=2, Nc=1, Nt=4 и кодируются векторы =(s1  s2  … s8)T, где . Параметры кода: J=1, Nq=2, NL=1, Nb=2. Порождающая матрица второго кода G(1, 2, 1, 2, 2, f)=P×T, где .

и

где

Кодовые слова второго кода: =G(1, 2, 1, 2, 2, exp(I×p/32))×, где =(c1(1)  c2(1)  c1(2)  c2(2)  c1(3)  c2(3)  c1(4)  c2(4))T, где cj(m) – элемент ПВ кодового слова, передаваемый j-й передающей антенной в течение m-й посылки.

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

Графики зависимости вероятности ошибки декодирования от отношения сигнал/шум в канале с медленными замираниями для первого и второго ПВ-кодов представлены на рисунке 1. Система связи 1´1 (NRx=NTx=1), не использующая кодиро­вание, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 38 дБ. Система связи 2´2 (NRx=NTx=2), не использующая пространственное кодирование, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 26 дБ, то есть выигрыш от увеличения числа передающих и принимающих антенн по сравнению с системой связи 1´1 составляет 12 дБ. Система связи 2´2, использующая первый ПВ-код, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 24 дБ, то есть выигрыш за счет применения первого ПВ-кода по сравнению с системой связи 2´2, не использующей пространственное кодирование, составляет 2 дБ. Такой же выигрыш обеспечивает и система связи 2´2, использующая второй ПВ-код, то есть в этом канале связи нет преимущества друг над другом у первого и второго ПВ-кодов. А система связи 1´1, в канале которой есть только аддитивный белый гауссовский шум и которая не использует кодирование, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 12,5 дБ.

Графики зависимости вероятности ошибки декодирования от отношения сигнал/шум в канале с быстрыми замираниями для первого и второго ПВ-кодов представлены на рисунке 2. Система связи 2´2, использующая первый ПВ-код, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 20 дБ, то есть выигрыш за счет применения первого ПВ-кода по сравнению с системой связи 2´2, не использующей пространственное кодирование, составляет 6 дБ. А система связи 2´2, использующая второй ПВ-код, обеспечивает вероятность ошибки 10-4 при отношении сигнал/шум на бит, приблизительно равном 16 дБ, то есть выигрыш за счет применения второго ПВ-кода по сравнению с системой связи 2´2, использующей первый ПВ-код, составляет 4 дБ.

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

Литература

1.     Tarokh V., Seshadri N., Calderbank A.R. Space-time codes for high data rate wireless communication: performance criterion and code construction, IEEE Transactions on Information Theory. 1998. Vol. 44, no. 2, pp. 744–765.

2.     Weifeng Su, Safar Z., Olfat M., Liu K.J.R. Obtaining full-diversity space-frequency codes from space-time codes via mapping, IEEE Transactions on Signal Processing, 2003, Vol. 51, no. 11, pp. 2905–2916.

3.     El Gamal H., Damen M.O. Universal Space–Time Coding, IEEE Transactions on Information Theory. 2003. Vol. 49, no. 5, pp. 1097–1119.

4.     Zhang W., Xia X.-G., Ching P.C. High-Rate Full-Diversity Space–Time–Frequency Codes for Broadband MIMO Block-Fading Channels, IEEE Transactions on Communications. 2007. Vol. 55, no. 1, pp. 25–34.

5.     Гофман М.В. Построение кодовых слов пространственно-частотно-временных кодов // Программные продукты и системы. 2010. № 3. C. 149–151.

6.     Гофман М.В. Метод построения алгебраических пространственно-частотно-временных кодов // Изв. ПГУПС. 2010. № 4. С. 88–98.

7.     Гофман М.В. Алгебраический пространственно-частотно-временной код // Информационно-управляющие системы. 2011. № 3. С. 39–46.

8.     Su W., Safar Z., Liu K.J.R. Full-rate full-diversity space-frequency codes with optimum coding advantage, IEEE Transac­tions on Information Theory. 2005. Vol. 51, no. 1, pp. 229–249.

9.     Kiran T., Rajan B.S. A systematic design of high-rate full-diversity space-frequency codes for MIMO-OFDM systems, International Symposium on Information Theory, ISIT-2005. Proc., 2005, pp. 2075–2079.

10.  Agrell E., Eriksson T., Vardy A., Zeger K. Closest point search in lattices, IEEE Transactions on Information Theory. 2002. Vol. 48, no. 8, pp. 2201–2214.

11.  Viterbo E., Boutros J. A universal lattice code decoder for fading channels, IEEE Transactions on Information Theory. 1999. Vol. 45, no. 5, pp. 1639–1642.