ЭФФЕКТИВНЫЙ РАСЧЕТ ВЕСОВЫХ КОЭФФИЦИЕНТОВ ФИЛЬТРА ДЛЯ СИСТЕМЫ MIMO
EFFICIENT FILTER WEIGHT COMPUTATION FOR MIMO SYSTEM

Область техники, к которой относится изобретение
Изобретение, в общем, относится к области связи и, более конкретно, к методикам расчета весовых коэффициентов фильтров в системе связи.
Уровень техники
В системе связи с множеством входов и множеством выходов (MIMO, МВМВ) для передачи данных используется множество (T) передающих антенн передающей станции и множество (R) приемных антенн приемной станции. Канал MIMO, формируемый T передающими антеннами и R приемными антеннами, может быть разложен на S пространственных каналов, где S≤min {T, R}. S пространственных каналов можно использовать для передачи данных таким образом, чтобы достичь большей общей пропускной способности и/или более высокой надежности.
Передающая станция может одновременно передавать T потоков данных через T передающих антенн. В этих потоках данных возникают искажения в соответствии с откликом канала MIMO, и их качество дополнительно ухудшается в результате воздействия шумов и взаимных помех. Приемная станция принимает передаваемые потоки данных через R приемных антенн. Принимаемый сигнал от каждой приемной антенны содержит масштабированную версию T потоков данных, переданных передающей станцией. Переданные потоки данных, таким образом, диспергированы среди R сигналов, принятых через R приемных антенн. Приемная станция затем выполняет пространственную обработку приемника для R принятых сигналов, используя матрицу пространственного фильтра, для восстановления переданных потоков данных.
Определение весовых коэффициентов матрицы пространственного фильтра требует объемной обработки. Это связано с тем, что матрицу пространственного фильтра обычно получают на основе функции, которая содержит обращение матрицы, и прямые расчеты обращения матрицы требуют объемных вычислений.
Таким образом, в данной области техники требуется разработка методики для эффективного расчета весовых коэффициентов фильтра.
Сущность изобретения
Здесь описаны методики эффективного расчета весовых коэффициентов матрицы пространственного фильтра. Эти методики позволяют исключить непосредственный расчет обращения матрицы.
М
H
a
h
C
i
i
i
i
r
В первом варианте воплощения для получения матрицы пространственного фильтра Эрмитову матрицу P итерационно получают на основе матрицы  отклика канала, и обращение матрицы опосредованно рассчитывают путем итерационного получения Эрмитовой матрицы. Эрмитова матрица может быть инициализирована до единичной матрицы. Одну итерацию затем выполняют для каждой строки матрицы отклика канала, и эффективную последовательность расчетов выполняют для каждой итерации. Для i-й итерации получают промежуточный вектор  строки на основе вектора  отклика канала, который представляет собой i-ю строку матрицы отклика канала. Скалярное значение получают на основе промежуточного вектора строки и вектора строки отклика канала. Промежуточную матрицу  также получают на основе промежуточного вектора строки. Эрмитову матрицу затем обновляют на основе скалярного значения и промежуточной матрицы. После окончания всех итераций, получают матрицу пространственного фильтра на основе Эрмитовой матрицы и матрицы отклика канала.
P
B
Y
Y
1/2
Во втором варианте воплощения выполняют множество поворотов для итерационного получения первый матрицы  и второй матрицы  для псевдообращенной матрицы отклика канала. Одну итерацию выполняют для каждой строки матрицы отклика канала. Для каждой итерации формируют матрицу , содержащую первую и вторую матрицы из предыдущей итерации. Множество поворотов Гивенса затем выполняют для матрицы  для обнуления элементов в первой строке матрицы, для получения обновленных первой и второй матриц для следующей итерации. После того как все итерации будут закончены, получают матрицу пространственного фильтра на основе первой и второй матриц.
X
V
Λ
X
В третьем варианте воплощения формируют матрицу  на основе матрицы отклика канала и разлагают (например, с использованием разложения по собственным значениям), для получения унитарной матрицы  и диагональной матрицы . Разложение может быть получено в результате итерационного выполнения поворотов Якоби для матрицы . Матрицу пространственного фильтра затем получают на основе унитарной матрицы, диагональной матрицы и матрицы отклика канала.
Различные аспекты и варианты выполнения изобретения более подробно описаны ниже.
Краткое описание чертежей
Свойства и сущность настоящего изобретения будут более понятны из подробного описания, приведенного ниже, которое следует рассматривать совместно с чертежами, на которых одинаковыми номерами ссылочных позиций обозначены соответствующие элементы на всех чертежах.
На фиг.1, 2 и 3 представлены процессы, выполняемые для расчета матрицы пространственного фильтра MMSE (МСКО, минимальная среднеквадратическая ошибка), на основе первого, второго и третьего вариантов воплощения, соответственно.
На фиг.4 показана блок-схема точки доступа и терминала пользователя.
Подробное описание изобретения
Слово "примерный", используемое здесь, означает "используемый в качестве примера, случая или иллюстрации". Любой вариант выполнения или конструкция, описанный здесь как "примерный", необязательно следует рассматривать как предпочтительный или преимущественный по сравнению с другими вариантами выполнения или конструкциями.
Описанные здесь методики расчета весовых коэффициентов фильтра можно использовать для системы MIMO с одной несущей и системы MIMO с множеством несущих. Множество несущих могут быть получены с использованием мультиплексирования с ортогональным частотным разделением сигналов (OFDM), при множественном доступе с разделением частот с перемежением (IFDMA), локализованном множественном доступе с разделением частот (LFDMA), или некоторой другой методики модуляции. OFDM, IFDMA и LFDMA эффективно разделяют общую полосу пропускания системы на множество (K) ортогональных частотных подполос, которые также называются тонами, поднесущими, элементами сигнала и частотными каналами. Каждую подполосу ассоциируют с соответствующей поднесущей, которая может быть модулирована данными. В системе OFDM символы модуляции передают в области частот для всех или поднабора K подполос. В IFDMA передают символы модуляции в области времени по подполосам, которые равномерно распределены по K подполосам. В LFDMA передают символы модуляции в области времени и обычно в соседних подполосах. Для ясности большая часть следующего описания направлена на систему MIMO с одной несущей, в которой используется одна поднесущая.
H
Канал MIMO, сформированный множеством (T) передающих антенн в передающей станции и множеством  приемных антенн в приемной станции, может характеризоваться матрицей  отклика канала размером RxT, которая может быть задана как:
Уравнение (1)
i,j
,
;
где hдля i=1,...,R и j=1,...,Tобозначает усиление связи или комплексное усиление канала между передающей антенной jи приемной антенной i и
h
H
i
,
представляет собой вектор строки отклика канала размером 1×T для приемной антенны iкоторый представляет собой i-ю строку матрицы .
Для простоты в следующем описании предполагается, что канал MIMO имеет полный ранг и, что количество пространственных каналов (S) задается как: S=T≤R.
Передающая станция может передавать T символов модуляции одновременно через T передающих антенн в каждый период символа. Передающая станция может выполнять или может не выполнять пространственную обработку для символов модуляции перед передачей. Для простоты в следующем описании предполагается, что каждый символ модуляции передают через передающую антенну без какой-либо пространственной обработки.
Приемная станция получает R принятых символов из R приемных антенн в каждый период символа. Принятые символы могут быть выражены как:
Уравнение (2)
s
где  представляет собой вектор размером T×1, в котором T символов модуляции переданы передающей станцией;
r
представляет собой вектор размером R×1, где R принимаемых символов получают в приемной станции через R приемных антенн; и
n
представляет собой вектор шумов размером R×1.
2
2
n
n
I
I
Для простоты можно предположить, что шумы представляют собой аддитивный белый Гауссов шум (AWGN) с нулевым средним вектором и матрицей ковариации δ ×, где δ представляетсобой дисперсию шума и  представляет собой единичную матрицу.
В приемной станции могут использоваться различные методики пространственной обработки для восстановления символов модуляции, переданных передающей станцией. Например, приемная станция может выполнять пространственную обработку приемника с минимальной среднеквадратичной ошибкой (MMSE), следующим образом:
Уравнение (3)
M
где  представляет собой матрицу пространственного фильтра MMSE размером T×R;
P
s
представляет собой Эрмитову матрицу ковариации размером T×T для ошибки оценки -;
представляет собой вектор размером T×1, который представляет собой оценку s; и
"
"
H
обозначает сопряженную перестановку.
P
P
s
s
P
H
i,j
i,j
Матрица  ковариации может быть задана как =E[(-)×(-)], где E[] представляет собой операцию математического ожидания.  также представляет собой Эрмитову матрицу, недиагональные элементы которой имеют следующие свойства p=p*, где "*" обозначает комплексно сопряженное число.
М
Как показано в уравнении (3), матрица пространственного фильтра MMSE имеет расчет обращенной матрицы. Прямой расчет обращения матрицы требует выполнения большого объема компьютерных операций. Матрица пространственного фильтра MMSE может быть более эффективно получена на основе вариантов выполнения, описанных ниже, которые позволяют опосредованно рассчитывать обращение матрицы с использованием итеративного процесса, вместо непосредственного расчета обращения матрицы.
М
P
В первом варианте воплощения расчета матрицы пространственного фильтра MMSE рассчитывают Эрмитову матрицу  на основе уравнении Риккати. Эрмитова матрица P может быть выражена следующим образом:
Уравнение (4)
P
i
Эрмитова матрица  размером TxT может быть определена как:
Уравнение (5)
Лемма обращения матрицы может быть применена к уравнению (5) для получения следующего:
Уравнение (6)
r
=
i
i
0
R
R
P
P
I
P
P
P
P
где представляет собой скалярную действительную величину. Уравнение (6) называется уравнением Риккати. Матрица  может быть инициализирована как =•. После выполнения R итераций уравнения (6), для i1,...,R, получают матрицу  как матрицу , или =.
Уравнение (6) может быть перемножено на определенные коэффициенты для получения следующего:
Уравнение (7)
P
P
I
P
P
P
P
P
P
P
i
0
R.
i
i
R
R
0
где матрицу  инициализируют как =, и матрицу  получают как =•Уравнения (6) и (7) отличаются от решения уравнения (5). Для простоты использовали те же переменные  и rдля обоих уравнений (6) и (7), даже, несмотря на то, что эти переменные имеют разные значения в двух уравнениях. Конечные результаты, полученные по уравнениям (6) и (7), то есть, для уравнения (6) и • для уравнения(7),эквивалентны.Однако расчеты для первой итерации уравнения (7) упрощаются, благодаря использованию  в качестве единичной матрицы.
Каждая итерация уравнения (7) может быть выполнена следующим образом:
Уравнение (8a)
Уравнение (8b)
Уравнение (8c)
Уравнение (8d)
a
i
где  представляет собой вектор промежуточной строки размером 1×T элементов с комплексным значением; и
C
представляет собой промежуточную Эрмитову матрицу размером T×T.
i
i
i
i
i
i
i
i
i
i
C
C
C
В системе (8) уравнений последовательность операций структурирована для эффективного расчета с помощью аппаратных средств. Скалярное значение r рассчитывают перед матрицей . Разделение по rв уравнении (7) достигается с помощью обращения и умножения. Обращение r может быть выполнено параллельно с расчетом . Обращение r достигается со сдвигом для нормализации r и с применением справочной таблицы для получения обращенного значения r. Нормализация r может быть компенсирована умножением на .
P
P
I
P
P
P
i
0
R.
Матрицу  инициализируют как Эрмитову матрицу или =, и она остается Эрмитовой матрицей во всех следующих итерациях. Следовательно, только верхнюю (или нижнюю) диагональную матрицу требуется рассчитывать для каждой итерации. После окончания R итераций получают матрицу , как =• Матрица пространственного фильтра MMSE может затем быть рассчитана следующим образом:
Уравнение (9)
M
P
P
i
0
На фиг.1 показан процесс 100 расчета матрицы  пространственного фильтра MMSE, основанной на первом варианте воплощения. Матрица  инициализирована как =1 (блок 112), и индекс i используется для обозначения номера итерации, инициализируется как i=1 (блок 114). Затем выполняют R итераций уравнения Риккати.
i-
δ
r
r
a
h
P
a
h
C
a
P
C
i
i
i+1
i
n
i
i
i
i
i
i
i
i
2
Каждую итерацию уравнения Риккати выполняют с помощью блока 120. Для й итерации вектор  промежуточной строки рассчитывают на основе заявленного вектора  отклика канала и Эрмитовой матрицы  предшествующей итерации, как показано в уравнении (8a) (блок 122). Скалярное значение rрассчитываютна основе дисперсии   шумов, вектора  промежуточной строки и вектора  строки отклика канала, как показано в уравнении (8b) (блок 124). Скалярная величина после этого становится обращенной (блок 126). Промежуточную матрицу рассчитывают на основе  промежуточной строки, как показано в уравнении (8c) (блок 128). Матрицу  затем обновляют на основе инвертированной скалярной величины и промежуточной матрицы , как показано в уравнении (8d) (блок 130).
δ
P
H
М
R
n
2
Затем определяют, были ли выполнены все R итераций (блок 132). В случае отрицательного ответа, выполняют последовательное приращение индекса i(блок 134), и процесс возвращается к блоку 122 для выполнения другой итерации. В противном случае, если все R итераций были выполнены, рассчитывают матрицу М пространственного фильтра MMSE на основе Эрмитовой матрицы  для последней итерации, матрицы  отклика каналов и дисперсии   шумов,как показано в уравнении (9) (блок 136). Матрицу  затем можно использовать для пространственной обработки приемника,как показано в уравнении (3).
М
P
P
P
1/2
Во втором варианте воплощения расчета матрицы пространственного фильтра MMSE определяют Эрмитову матрицу  путем получения квадратного корня , который представляет собой , на основе итерационной процедуры. Пространственная обработка приемника в уравнении (3) может быть выражена следующим образом:
Уравнение (10)
U
где = представляет собой дополненную матрицу канала размером (R+T)×T;
U
U
U
U
U
U
p
p
H
-1
H
представляет собой псевдообращенную матрицу размером T×(R+T), полученную в результате операции обращения или псевдообращения Мура-Пенроуза для , или =( ∙ ) ∙ ;
0
Tx1
представляет собой вектор размером T×1, содержащий одни нули; и
H
p
p
δn
U
представляет собой подматрицу размером T×R, содержащую первые R столбцов .
Разложение QR может быть выполнено для матрицы с дополненным каналом следующим образом:
Уравнение (11)
Q
где  представляет собой матрицу размером (R+T)×T с ортонормированными столбцами;
R
представляет собой матрицу размером T×T, которая не является единичной матрицей;
B
Q
представляет собой матрицу размером R×T, содержащую первые R строк матрицы ; и
Q
Q
2
представляет собой матрицу размером T×T, содержащую последние T строк матрицы .
Q
R
Q
Q
Q
I
H
QR (КО, квазиобращенное) разложение в уравнении (11) разлагает матрицу дополненного канала на ортогональную матрицу  и на не единичную матрицу . Ортогональная матрица  имеет следующее свойство: •=, что означает, что столбцы ортогональной матрицы являются ортогональными относительно друг друга, и каждый столбец имеет единичную степень. Не единичная матрица представляет собой матрицу, для которой может быть рассчитана обращенная матрица.
P
Эрмитова матрица  может быть затем выражена как:
Уравнение (12)
R
P
R
P
-1
1/2
-1
представляет собой разложение Колецкого или квадратный корень матрицы P. Следовательно,  равно  и называется квадратным корнем матрицы .
Псевдообращенная матрица в уравнении (10) затем может быть выражена как:
Уравнение (13)
H
P
δn
Подматрица  , которая также представляет собой матрицу пространственного фильтра MMSE, затем может быть выражена как:
Уравнение (14)
Уравнение (10) затем может быть выражено как:
Уравнение (15)
P
B
1/2
Матрицы  и  могут быть рассчитаны итерационно следующим образом:
или Уравнение (16)
Уравнение (17)
Y
P
B
h
i
i-1
i-1
i
1/2
где  представляет собой матрицу размером (T+R+1)×(T+1), содержащую элементы, полученные на основе  ,  и ;
θ
i
представляет собой унитарную матрицу преобразования размером (T+1)×(T+1);
Z
B
i
i
i
i
1/2
r
представляет собой преобразованную матрицу размером (T+R+1)×(T+1), содержащую элементы для P ,  и ;
e
i
представляет собой вектор размером R×1 с единицей (1,0) в качестве i-го элемента и с остальными нулевыми элементами; и
k
I
i
i
представляет собой вектор размером T×1 и  представляет собой вектор R×1, причем оба они являются несущественными.
P
B
P
I
B
0
1/2
1/2
0
0
RxT
Матрицы  и  инициализируют как  =• и =.
h
δ
Z
Z
k
I
Z
P
B
Z
P
B
P
P
P
B
B
М
P
B
i
i
i
i
i
i
i
i
i
i
i
i
i
i
R
R
1/2
1/2
1/2
1/2
1/2
1/2
1/2
Преобразование в уравнение (17) может быть выполнено итерационно, как описано ниже. Для ясности каждую итерацию уравнения (17) называют внешней итерацией. R внешних итераций уравнения (17) выполняют для R векторов  отклика канала для i=1,...,R. Для каждой внешней итерации унитарная матрица  преобразования в уравнении (17) преобразуется в трансформированную матрицу , содержащую все нули в первой строке, за исключением первого элемента. Первый столбец преобразованной матрицы  содержит r ,  и . Последние T столбцов  содержат обновленные   и . Первый столбец  не требуется рассчитывать, поскольку только   и  используются в следующей итерации.   представляет собой верхнюю треугольную матрицу. После окончания R внешних итераций получают   как , и  получают как . Матрица  пространственного фильтра MMSE может быть затем рассчитана на основе  и , как представлено в уравнении (14).
i
Y
Y
i
i
Для каждой внешней итерации преобразование по уравнению (17) может выполняться путем последовательного обнуления одного элемента в первой строке  одновременно с 2×2 поворотами Гивенса. T внутренних итераций поворота Гивенса могут быть выполнены для обнуления последних T элементов в первой строке .
Y
Y
Y
Y
Y
Y
Y
i,j
i1
i,j
i,j
i,j
i,j
Для каждой внешней итерации i, матрица может быть инициализирована как =i. Для каждой внутренней итерации j для j=1,...,T, внешней итерации i, первоначально формируют подматрицу ′размером (T+R+1)×2, содержащую первый и (j+1)-й столбцы . Затем выполняют поворот Гивенса для подматрицы ′для генерирования подматрицы "размером (T+R+1)×2, содержащей ноль во втором элементе в первой строке. Поворот Гивенса может быть выражен как:
Уравнение (18)
G
Y
Y
Y
Y
Y
Y
Y
Y
Y
Y
Y
Y
Y
Y
i,j
i,j+1
i,j+1
i,j
i,j+1
i,j
i,j+1
i,j
i,j
i,j+1
i
i,j
i,j
i,j
i,j+1
где представляет собой матрицу поворота Гивенса размером 2×2 для j-й внутренней итерации i-й внешней итерации, которая описана ниже. Матрицу  затем формируют вначале путем установки =, затем замены первого столбца  первым столбцом "и затем замены (j+1)-го столбца матрицы  вторым столбцом ". Поворот Гивенса, таким образом, модифицирует только два столбца  j-й внутренней итерации для получения  для следующей внутренней итерации. Поворот Гивенса может быть выполнен на месте двух столбцов  для каждой внутренней итерации, в результате чего промежуточные матрицы ,′,"и не нужны и описаны выше для ясности изложения.
G
Y
G
i,j
i,j
i,j
jθ
Для j-й внутренней итерации i-й внешней итерации матрица  поворота Гивенса определяется на основе первого элемента (который всегда представляет собой действительную величину) и (j+1)-го элемента в первой строке матрицы . Первый элемент может быть обозначен как α, и (j+1)-й элемент может быть обозначен как b-e. Матрица поворотаГивенса затем может быть получена следующим образом:
Уравнение (19)
где c= и s= для уравнения (19).
М
P
I
B
B
0
i
0
i
0
1/2
1/2
На фиг.2 показан процесс 200, предназначенный для расчета матрицы  пространственного фильтра MMSE на основе второго варианта воплощения. Матрицу   инициализируют как P =•, и матрицу  инициализируют как = (блок 212). Индекс i для обозначения номера внешней итерации инициализируют как i=1, и индекс j, используемый для обозначения номера внутренней итерации, инициализируют как j=1 (блок 214). Затем выполняют R внешних итераций унитарного преобразования в соответствии с уравнением (17) (блок 220).
Y
h
P
B
Y
Y
Y
i
i
i-1
i-1
i
i,j
i,j
1/2
Для i-й внешней итерации вначале формируют матрицу  с вектором  строки отклика канала и матрицы   и , как показано в уравнении (17) (блок 222). Матрица  затем обозначается как матрица  для внутренних итераций (блок 224). T внутренних итераций поворота Гивенса затем выполняют для матрицы  (блоке 230).
G
Y
G
Y
Y
i,j
i,j
i,j
i,j
i,j+1
Для j-й внутренней итерации получают матрицу поворота Гивенса на основе первого и (j+1)-го элементов в первой строке , как показано в уравнении (19) (блок 232). Матрицу поворота Гивенса затем применяют для первого и (j+1)-го столбцов  для получения , как показано в уравнении (18) (блок 234). Затем определяют, были ли выполнены все T внутренние итерации (блок 236). Если ответ представляет собой "Нет", тогда индекс j увеличивают на единицу (блок 238), и обработка возвращается к блоку 232 для выполнения другой внутренней итерации.
Y
Z
P
B
Y
P
B
М
P
B
i,j+1
i
i
i
i,j+1
i
i.
i
i
1/2
1/2
1/2
Если все T внутренние итерации были выполнены для текущей внешней итерации, и ответ представляет собой "Да" для блока 236, тогда последняя  равна  в уравнении (17). Обновленные матрицы   и  получают из последней  (блок 240). Затем определяют, были ли выполнены все R внешних итераций (блок 242). Если будет получен ответ "Нет", тогда индекс i увеличивают на единицу, и индекс j повторно инициализируют как j=1 (блок 244). Обработка затем возвращается в блок 222 для выполнения другой внешней итерации с   и  В противном случае, если все R внешние итерации были выполнены, и ответ представляет собой "Да" для блока 242, тогда рассчитывают матрицу  пространственного фильтра MMSE, на основе   и , как показано в уравнении (14) (блок 246). Матрицу М можно затем использовать для пространственной обработки приемника, как показано в уравнении (15).
P
-1
В третьем варианте воплощения расчета матрицы М пространственного фильтра MMSE выполняют разложение по собственным значениям  следующим образом:
Уравнение (20)
V
где  представляет собой унитарную матрицу T×T собственных векторов; и
Λ
представляет собой диагональную матрицу размером T×T с действительнымисобственными значениями вдоль диагонали.
X
X
X
V
X
X
V
2x2
2x2
2x2
2x2
2x2
2x2
2x2
Разложение по собственным значениям Эрмитовой матрицы  размером 2×2 может быть получено с использованием разных методик. В варианте воплощения разложение по собственным значениям  получают путем выполнения сложного поворота Якоби для , для получения матрицы  размером 2×2 собственных векторов . Элементы  и могут быть заданы как:
Уравнение (21)
V
X
2x2
2x2
Элементы  могут быть рассчитаны непосредственно из элементов  следующим образом:
Уравнение (22a)
Уравнение (22b)
Уравнение (22c)
Уравнение (22d)
Уравнение (22e)
Уравнение (22f)
Уравнение(22g)
Уравнение (22h)
Уравнение (22i)
Уравнение (22j)
Уравнение (22k)
X
X
X
X
P
D
Λ
D
X
V
V
V
I
-1
i
0
0
Разложение по собственным значениям Эрмитовой матрицы  размером T×T, которая больше чем 2×2, может быть выполнено в итеративном процессе. В таком итеративном процессе неоднократно используется поворот Якоби для обнуления недиагональных элементов в . Для итеративного процесса индекс i обозначает номер итерации и инициализируется как i=1.  представляет собой Эрмитову матрицу размером T×T, которая должна быть разложена, и устанавливается как =. Матрица  представляет собой аппроксимацию диагональной матрицы  в уравнении (20) и инициализируется как =. Матрица  представляет собой аппроксимацию унитарной матрицы  в уравнении (20) и инициализируется как =.
D
V
D
D
i
i
pq
i
Одиночная итерация поворота Якоби для обновления матриц  и  может быть выполнена следующим образом. Вначале Эрмитову матрицу  размером 2×2 формируют на основе текущей матрицы  следующим образом:
Уравнение (23)
p,q
i
pq
i
pq
i.
D
D
D
D
где d представляет собой элемент в местоположении (p,q) матрицы , p{1,...,T}, g{1,...,T} и p≠q. представляетсобой подматрицу  размером 2×2, и четыре элемента Dпредставляют собой четыре элемента в местоположениях (p, p), (p, q), (q, p) и (q, q) матрицы  Индексы p и q могут быть выбраны, как описано ниже.
D
V
D
D
X
D
V
V
pq
pq
pq
pq,
2x2
pq
2x2
pq
Затем выполняют разложение по собственным значениям , как показано в уравнении (22), для получения унитарной матрицы  размером 2×2 собственных векторов . Для разложения по собственным значениям  в уравнении (21) заменяют на , и  из уравнения (22j) или (22k) представляют как .
pq
pq,
pq
1,1
1,2
2,1
2,2
pq
V
T
Y
Матрица T сложного поворота Якоби размером TxT затем формируется с и представляет собой единичную матрицу с четырьмя элементами в местоположениях (p, p), (p, q), (q, p) и (q, q), которые заменены элементами v, v, v и v, соответственно, в матрице .
D
i
Матрица  затем обновляется следующим образом:
Уравнение (24)
D
D
i
i
Уравнение (24) обнуляет два недиагональных элемента в местоположениях (p, q) и (q, p)в матрице . Расчет может изменять значения других недиагональных элементов в .
V
i
Матрицу  также обновляют следующим образом:
Уравнение (25)
V
T
D
i
pq
i
можно рассматривать как матрицу кумулятивного преобразования, которая содержит все матрицы поворота Якоби, использовавшиеся для .
D
D
D
V
Λ
V
i
i
i
Каждая итерация с поворотом Якоби обнуляет два недиагональных элемента матрицы i. Множество итераций поворота Якоби могут быть выполнены для разных значений индексов p и q, для обнуления всех недиагональных элементов . Один проход по всем возможным значениям индексов p и q может быть выполнен следующим образом. Индекс p последовательно изменяется от 1 до T-1 с приращением единица. Для каждого значения p индекс q последовательно изменяется от p+1 до T с приращением единица. Поворот Якоби выполняют для каждой разной комбинации значений p и q. Множество проходов могут быть выполнены до тех пор, пока  и не станут представлять достаточно точные оценки  и , соответственно.
Уравнение (20) может быть переписано следующим образом:
Уравнение (26)
Λ
Λ
X
P
Λ
V
Λ
Λ
-1
-1
-1
где  представляет собой диагональную матрицу, элементы которой представляют собой обращенные значения соответствующих элементов . Разложение по собственным величинам = представляет оценки  и .  может быть инвертирована для получения .
Матрица пространственного фильтра MMSE может быть затем рассчитана следующим образом:
Уравнение (27)
М
P
H
P
V
Λ
М
V
Λ
H
-1
-1
На фиг.3 представлен процесс 300, предназначенный для расчета матрицы  пространственного фильтра MMSE, на основе третьего варианта воплощения. Эрмитову матрицу первоначально получают на основе матрицы  отклика канала, как показано в уравнении (20) (блок 312). Затем выполняют разложение по собственным величинам  для получения унитарной матрицы  и диагональной матрицы , как также показано в уравнении (20) (блок 314). Разложение по собственным величинам может итерационно выполняться с множеством поворотов Якоби, как описано выше. Матрицу  пространственного фильтра MMSE затем получают на основе унитарной матрицы , диагональной матрицы  и матрицы  отклика канала, как показано в уравнении (27) (блок 316).
М
М
D
M
H
М
H
М
H
mmse
mmse
mmse
mmse
M
.
D
-1
Матрица  пространственного фильтра MMSE, полученная на основе каждого из описанных выше вариантов воплощения, представляет собой смещенное решение MMSE. Смещенную матрицу  пространственного фильтра можно масштабировать с помощью диагональной матрицы дляполучения несмещенной матрицы MMSE пространственного фильтраМатрица может быть получена как D=[diag[•]], где diag[•] представляет собой диагональную матрицу, содержащую диагональные элементы •.
Описанные выше расчеты также можно использовать для получения матриц пространственного фильтра для методики с обращением в нуль не значащих коэффициентов (ZF) (которая также называется методикой обращения матрицы с корреляцией канала (CCMI, ОМКК)), методики комбинирования максимального отношения (MRC, КМО) и так далее. Например, приемная станция может выполнять пространственную обработку приемника с обращением в нуль не значащих коэффициентов и MRC, следующим образом:
Уравнение (28)
Уравнение (29)
M
zf
где  представляет собой матрицу пространственного фильтра размером T×R с обращенными в ноль не значащими коэффициентами;
M
mrc
представляет собой матрицу пространственного фильтра MRC с размером T×R;
P
zf
H
-1
=(H•H) представляет собой Эрмитову матрицу размером T×T; и
P
P
zf
zf
[diag()] представляет собой диагональную матрицу размером T×T, содержащую диагональные элементы .
P
P
zf
zf
Обращение матрицы необходимо для непосредственного расчета .  могут быть рассчитаны с использованием вариантов воплощения, описанных выше для матрицы пространственного фильтра MMSE.
В приведенном выше описании предполагается, что T символов модуляции передают одновременно через T передающие антенны, без какой-либо пространственной обработки. Передающая станция может выполнять пространственную обработку перед передачей следующим образом:
Уравнение (30)
x
где  представляет собой вектор размером T×1 с T символами передачи, которые должны быть переданы через T передающих антенн; и
W
W
H
H
H,
H
H
H
W
H
H
H
eff
eff
eff
представляет собой матрицу передачи размером T×S. Матрица  передачи может представлять собой (1) матрицу с правыми единичными векторами, полученными путем выполнения разложения единичного значения , (2) матрицу собственных векторов, полученную путем выполнения разложения по собственным величинам  или (3) управляющую матрицу, выбранную для пространственного распределения символов модуляции по S пространственным каналам в канале MIMO. Матрица эффективного отклика канала, наблюдаемого по символам модуляции, затем может быть задана как =•. Описанная выше комбинация может выполняться на основе  вместо .
H
М
H
Для ясности приведенного выше описания, в нем приведена система MIMO с одной несущей, с одной подполосой. Для системы MIMO с множеством несущих матрица (k) отклика канала может быть получена для каждой подполосы k, представляющей интерес. Матрица (k) пространственного фильтра затем может быть получена для каждой подполосы k на основе матрицы (k) отклика канала для этой подполосы.
1
1
Описанные выше расчеты для матрицы пространственного фильтра могут быть выполнены с использованием процессоров различного типа, таких как процессор с плавающей запятой, процессор с фиксированной запятой, процессор цифрового компьютера поворота координат (CORDIC), таблица поиска и так далее, или их комбинации. Процессор CORDIC воплощает итеративный алгоритм, который обеспечивает возможность выполнения быстрого расчета с использованием аппаратных средств тригонометрических функций, таких как синус, косинус, магнитуда и фаза, с использованием простых аппаратных средств сдвига и сложения/вычитания. Процессор CORDIC может итерационно рассчитывать каждую из переменных r, c и s системы (22) уравнений с большим количеством итераций, что позволяет достичь более высокой точности для переменной.
ap
,
ut
ap
ut
ap
ap
W
На фиг.4 показана блок-схема точки 410 доступа и терминала 450 пользователя в системе 400 MIMO. Точка 410 доступа оборудована N антеннами и терминал 450 пользователя оборудован N антеннами, где N>1 и N>1. По нисходящему каналу передачи, в точке 410 доступа процессор 414 передаваемых (TX) данных принимает данные трафика из источника 412 данных и другие данные из контроллера/процессора 430. Процессор 414 TX данных форматирует, кодирует, выполняет перемежение, модулирует данные и генерирует символы данных, которые представляют собой символы модуляции для данных. Пространственный процессор 420 TX мультиплексирует символы данных с пилотными символами, выполняет пространственную обработку с матрицей  передачи, если это применимо, и предоставляет N потоков передаваемых символов. Каждый модуль 422 передатчиков (TMTR) обрабатывает соответствующий поток передаваемых символов и генерирует модулированный сигнал для нисходящего канала передачи. Модулированные сигналы для нисходящего канала передачи N из модулей 422a-422ap передатчика передают через антенны 424a-424ap, соответственно.
ut
dn
dn
dn
ut
dn
M
H
M
В терминале 450 пользователя N антенн 452a-452ut принимают переданные модулированные сигналы нисходящего канала передачи, и каждая антенна передает принятый сигнал в соответствующий модуль (RCVR) 454 приемника. Каждый модуль 454 приемника выполняет обработку, взаимодополняющую обработку, выполнявшуюся модулями 422 передачи, и предоставляет принятые пилотные символы и принятые символы данных. Блок/процессор 478 оценки канала обрабатывает принятые пилотные символы и предоставляет оценку отклика H канала нисходящего канала передачи. Процессор 480 получает матрицу  пространственного фильтра нисходящего канала передачи на основе , и используя любой из вариантов воплощения, описанных выше. Пространственный процессор 460 приемника (RX) выполняет пространственную обработку приемника (или пространственную согласованную фильтрацию) для принятых символов данных из всех Nмодулей454a-454ut приемникас матрицей  пространственного фильтра нисходящего канала передачи и предоставляет детектируемые символы данных, которые являются оценками символов данных, переданных точкой 410 доступа. Процессор 470 приемника обрабатывает (например, выполняет обратное отображение символа, устраняет перемежение и декодирует) детектированные символы данных и предоставляет декодированные данные в приемник 472 данных и/или контроллер 480.
ut
Обработка для восходящего канала передачи может быть такой же или может отличаться от обработки для нисходящего канала передачи. Данные из источника 486 данных и сигналы из контроллера 480 обрабатывают (например, кодируют, выполняют перемежение и модулируют) с помощью процессора 488 TX данных, мультиплексированных с пилотными символами и, возможно, пространственно обработанных пространственным процессором 490 TX. Символы передачи из пространственного процессора 490 TX дополнительно обрабатывают с помощью модулей 454a-454ut передатчика для генерирования Nмодулированных сигналов восходящего канала передачи, которые передают через антенны 452a-452ut.
H
M
M
up
up
up
В точке 410 доступа модулированные сигналы восходящего канала передачи принимают с помощью антенн 424a-424ap и обрабатывают с помощью модулей 422a-422ap приемника для генерирования принятых пилотных символов и принятых символов данных, для передачи по восходящему каналу передачи. Блок/процессор 428 оценки канала обрабатывает принимаемые пилотные символы и предоставляет оценку отклика  канала восходящей передачи. Процессор 430 получает матрицу  пространственного фильтра восходящего канала передачи, и используя один из вариантов выполнения, описанных выше. Пространственный процессор 440 RX выполняет пространственную обработку приемника для принимаемых символов данных с матрицей  пространственного фильтра восходящего канала передачи и предоставляет детектированные символы данных. Процессор 442 данных RX дополнительно обрабатывает детектированные символы данных и предоставляет декодированные данные в приемник 444 данных и/или контроллер 430.
Контроллеры 430 и 480 управляют операциями в точке 410 доступа и в терминале 450 пользователя, соответственно. В модулях 432 и 482 сохраняются данные и программные коды, используемые контроллерами 430 и 480, соответственно.
Блоки, показанные на фиг.1-4, представляют функциональные блоки, которые могут быть воплощены в виде аппаратных средств (одно или больше устройств), встроенных программ (одно или больше устройств), программного обеспечения (один или больше модулей) или их комбинации. Например, описанные здесь методики расчета весовых коэффициентов фильтра могут быть воплощены в виде аппаратных средств, встроенных программ, программного обеспечения или их комбинации. При выполнении в виде аппаратных средств модули обработки, используемые для расчета весовых коэффициентов фильтра, могут быть воплощены в одной или больше специализированных интегральных схемах (ASIC), цифровых сигнальных процессорах (DSP), устройствах обработки цифровых сигналов (DSPD), программируемых логических устройств (PLD), программируемых пользователем вентильных матриц (FPGA), процессоров, контроллеров, микроконтроллеров, микропроцессоров, электронных устройств, других электронных модулей, разработанных для выполнения описанных здесь функций, или их комбинации. Различные процессоры в точке 410 доступа на фиг.4 также могут быть воплощены с использованием одного или больше аппаратных процессоров. Аналогично, различные процессоры в терминале 450 пользователя могут быть воплощены с одними или больше аппаратных процессоров.
Для варианта воплощения с использованием встроенного программного обеспечения или программного обеспечения, методики расчета весового коэффициента фильтра могут быть воплощены с модулями (например, процедурами, функциями и так далее), которые выполняют описанные здесь функции. Программные коды могут быть сохранены в модуле памяти (например, модуле 432 или 482 памяти по фиг.4) и могут выполняться процессором (например, процессором 430 или 480). Модуль памяти может быть воплощен внутри процессора или вне процессора.
Приведенное выше описание раскрытых вариантов выполнения предоставлено для обеспечения возможности для специалиста в данной области техники использовать настоящее изобретение. Различные модификации этих вариантов воплощения будут понятны для специалиста в данной области техники, и общие принципы, определенные здесь, могут применяться для других вариантов воплощения, без выхода за пределы сущности или объема изобретения. Таким образом, настоящее изобретение не предназначено для ограничения вариантами воплощения, описанными здесь, но его следует понимать в самом широком объеме, который соответствует раскрытым здесь принципам и новым свойствам.