эффективный расчет весовых коэффициентов фильтра для системы mimo

Формула изобретения

1. Устройство для получения матрицы пространственного фильтра, содержащее:

первый процессор, который во время работы получает матрицу отклика канала; и

второй процессор, который во время работы итерационно получает первую матрицу на основе матрицы отклика канала и выводит матрицу пространственного фильтра на основе первой матрицы и матрицы отклика канала, при этом второй процессор рассчитывает обращение матрицы для матрицы, полученной от матрицы отклика канала, путем итерационного получения первой матрицы.

2. Устройство по п.1, в котором второй процессор во время работы инициализирует первую матрицу до единичной матрицы.

3. Устройство по п.1, в котором второй процессор во время работы для каждой из множества итераций получает промежуточный вектор строки на основе первой матрицы и вектора строки отклика канала, соответствующего строке матрицы отклика канала, для получения скалярной величины на основе промежуточного вектора строки и вектора строки отклика канала для получения промежуточной матрицы на основе промежуточного вектора строки и для обновления первой матрицы на основе скалярной величины и этой промежуточной матрицы.

4. Устройство по п.1, в котором первая матрица предназначена для получения матрицы пространственного фильтра с минимальной среднеквадратической ошибкой (MMSE).

5. Устройство по п.1, в котором второй процессор во время работы получает первую матрицу на основе следующего уравнения:



где P_i представляет собой первую матрицу для i-й итерации, h_i представляет собой i-ую строку матрицы отклика канала, r_i представляет собой скалярную величину, полученную на основе h_i и P_i-1, и "^H" представляет собой сопряженную перестановку.

6. Устройство по п.1, в котором второй процессор во время работы получает первую матрицу на основе следующих уравнений:





где P_i представляет собой первую матрицу для i-й итерации, h_i представляет собой i-ю строку матрицы отклика канала, a_i представляет собой промежуточный вектор строки для i-й итерации, C _i представляет собой промежуточную матрицу для i-й итерации, r_i представляет собой скалярную величину для i-й итерации, ² _n представляет собой дисперсию шума, и " ^H" представляет собой сопряженную перестановку.

7. Устройство по п.1, в котором второй процессор во время работы получает матрицу пространственного фильтра на основе следующего уравнения:



где М представляет собой матрицу пространственного фильтра, P представляет собой первую матрицу, H представляет собой матрицу отклика канала, и "^H" представляет собой сопряженную перестановку.

8. Способ получения матрицы пространственного фильтра, содержащий этапы, на которых:

итерационно получают первую матрицу на основе матрицы отклика канала, при этом обращение матрицы для матрицы, полученной от матрицы отклика канала, рассчитывают путем итерационного получения первой матрицы; и

получают матрицу пространственного фильтра на основе первой матрицы и матрицы отклика канала.

9. Способ по п.8, содержащий также этап, на котором инициализируют первую матрицу до единичной матрицы.

10. Способ по п.8, в котором получение первой матрицы содержит, для каждой из множества итераций этапы, на которых

получают промежуточный вектор строки на основе первой матрицы и вектора строки отклика канала, соответствующего строке матрицы отклика канала,

получают скалярную величину на основе вектора промежуточной строки и вектора строки отклика канала,

получают промежуточную матрицу на основе промежуточного вектора строки, и

обновляют первую матрицу на основе скалярной величины и промежуточной матрицы.

11. Устройство для получения матрицы пространственного фильтра, содержащее:

средство итерационного получения первой матрицы на основе матрицы отклика канала, при этом обращение матрицы для матрицы, полученной от матрицы отклика канала, рассчитывают путем итерационного получения первой матрицы; и

средство получения матрицы пространственного фильтра на основе первой матрицы и матрицы отклика канала.

12. Устройство по п.11, которое также содержит средство инициализации первой матрицы до единичной матрицы.

13. Устройство по п.11, в котором средство получения первой матрицы содержит для каждой из множества итераций:

средство получения промежуточного вектора строки на основе первой матрицы и вектора строки отклика канала, соответствующего строке матрицы отклика канала,

средство получения скалярной величины на основе промежуточного вектора строки и вектора строки отклика канала,

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

средство обновления первой матрицы на основе скалярной величины и промежуточной матрицы.

Описание изобретения к патенту

Область техники, к которой относится изобретение

Изобретение, в общем, относится к области связи и, более конкретно, к методикам расчета весовых коэффициентов фильтров в системе связи.

Уровень техники

В системе связи с множеством входов и множеством выходов (MIMO, МВМВ) для передачи данных используется множество (T) передающих антенн передающей станции и множество (R) приемных антенн приемной станции. Канал MIMO, формируемый T передающими антеннами и R приемными антеннами, может быть разложен на S пространственных каналов, где S min {T, R}. S пространственных каналов можно использовать для передачи данных таким образом, чтобы достичь большей общей пропускной способности и/или более высокой надежности.

Передающая станция может одновременно передавать T потоков данных через T передающих антенн. В этих потоках данных возникают искажения в соответствии с откликом канала MIMO, и их качество дополнительно ухудшается в результате воздействия шумов и взаимных помех. Приемная станция принимает передаваемые потоки данных через R приемных антенн. Принимаемый сигнал от каждой приемной антенны содержит масштабированную версию T потоков данных, переданных передающей станцией. Переданные потоки данных, таким образом, диспергированы среди R сигналов, принятых через R приемных антенн. Приемная станция затем выполняет пространственную обработку приемника для R принятых сигналов, используя матрицу пространственного фильтра, для восстановления переданных потоков данных.

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

Таким образом, в данной области техники требуется разработка методики для эффективного расчета весовых коэффициентов фильтра.

Сущность изобретения

Здесь описаны методики эффективного расчета весовых коэффициентов матрицы пространственного фильтра. Эти методики позволяют исключить непосредственный расчет обращения матрицы.

В первом варианте воплощения для получения матрицы М пространственного фильтра Эрмитову матрицу P итерационно получают на основе матрицы H отклика канала, и обращение матрицы опосредованно рассчитывают путем итерационного получения Эрмитовой матрицы. Эрмитова матрица может быть инициализирована до единичной матрицы. Одну итерацию затем выполняют для каждой строки матрицы отклика канала, и эффективную последовательность расчетов выполняют для каждой итерации. Для i-й итерации получают промежуточный вектор a_i строки на основе вектора h_i отклика канала, который представляет собой i-ю строку матрицы отклика канала. Скалярное значение r _i получают на основе промежуточного вектора строки и вектора строки отклика канала. Промежуточную матрицу C _i также получают на основе промежуточного вектора строки. Эрмитову матрицу затем обновляют на основе скалярного значения и промежуточной матрицы. После окончания всех итераций, получают матрицу пространственного фильтра на основе Эрмитовой матрицы и матрицы отклика канала.

Во втором варианте воплощения выполняют множество поворотов для итерационного получения первый матрицы P^1/2 и второй матрицы B для псевдообращенной матрицы отклика канала. Одну итерацию выполняют для каждой строки матрицы отклика канала. Для каждой итерации формируют матрицу Y, содержащую первую и вторую матрицы из предыдущей итерации. Множество поворотов Гивенса затем выполняют для матрицы Y для обнуления элементов в первой строке матрицы, для получения обновленных первой и второй матриц для следующей итерации. После того как все итерации будут закончены, получают матрицу пространственного фильтра на основе первой и второй матриц.

В третьем варианте воплощения формируют матрицу X на основе матрицы отклика канала и разлагают (например, с использованием разложения по собственным значениям), для получения унитарной матрицы V и диагональной матрицы . Разложение может быть получено в результате итерационного выполнения поворотов Якоби для матрицы X. Матрицу пространственного фильтра затем получают на основе унитарной матрицы, диагональной матрицы и матрицы отклика канала.

Различные аспекты и варианты выполнения изобретения более подробно описаны ниже.

Краткое описание чертежей

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

На фиг.1, 2 и 3 представлены процессы, выполняемые для расчета матрицы пространственного фильтра MMSE (МСКО, минимальная среднеквадратическая ошибка), на основе первого, второго и третьего вариантов воплощения, соответственно.

На фиг.4 показана блок-схема точки доступа и терминала пользователя.

Подробное описание изобретения

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

Описанные здесь методики расчета весовых коэффициентов фильтра можно использовать для системы MIMO с одной несущей и системы MIMO с множеством несущих. Множество несущих могут быть получены с использованием мультиплексирования с ортогональным частотным разделением сигналов (OFDM), при множественном доступе с разделением частот с перемежением (IFDMA), локализованном множественном доступе с разделением частот (LFDMA), или некоторой другой методики модуляции. OFDM, IFDMA и LFDMA эффективно разделяют общую полосу пропускания системы на множество (K) ортогональных частотных подполос, которые также называются тонами, поднесущими, элементами сигнала и частотными каналами. Каждую подполосу ассоциируют с соответствующей поднесущей, которая может быть модулирована данными. В системе OFDM символы модуляции передают в области частот для всех или поднабора K подполос. В IFDMA передают символы модуляции в области времени по подполосам, которые равномерно распределены по K подполосам. В LFDMA передают символы модуляции в области времени и обычно в соседних подполосах. Для ясности большая часть следующего описания направлена на систему MIMO с одной несущей, в которой используется одна поднесущая.

Канал MIMO, сформированный множеством (T) передающих антенн в передающей станции и множеством ® приемных антенн в приемной станции, может характеризоваться матрицей H отклика канала размером RxT, которая может быть задана как:

Уравнение (1)

где h _i,j, для i=1,...,R и j=1,...,Tобозначает усиление связи или комплексное усиление канала между передающей антенной jи приемной антенной i; и

h_i представляет собой вектор строки отклика канала размером 1×T для приемной антенны i, который представляет собой i-ю строку матрицы H.

Для простоты в следующем описании предполагается, что канал MIMO имеет полный ранг и, что количество пространственных каналов (S) задается как: S=T R.

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

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

Уравнение (2)

где s представляет собой вектор размером T×1, в котором T символов модуляции переданы передающей станцией;

r представляет собой вектор размером R×1, где R принимаемых символов получают в приемной станции через R приемных антенн; и

n представляет собой вектор шумов размером R×1.

Для простоты можно предположить, что шумы представляют собой аддитивный белый Гауссов шум (AWGN) с нулевым средним вектором и матрицей ковариации ² _n×I, где ² _nпредставляетсобой дисперсию шума и I представляет собой единичную матрицу.

В приемной станции могут использоваться различные методики пространственной обработки для восстановления символов модуляции, переданных передающей станцией. Например, приемная станция может выполнять пространственную обработку приемника с минимальной среднеквадратичной ошибкой (MMSE), следующим образом:

Уравнение (3)

где M представляет собой матрицу пространственного фильтра MMSE размером T×R;

P представляет собой Эрмитову матрицу ковариации размером T×T для ошибки оценки s- ;

представляет собой вектор размером T×1, который представляет собой оценку s; и

"^H " обозначает сопряженную перестановку.

Матрица P ковариации может быть задана как P=E[( s- )×(s- )^H], где E[] представляет собой операцию математического ожидания. P также представляет собой Эрмитову матрицу, недиагональные элементы которой имеют следующие свойства p _i,j=p*_i,j, где "*" обозначает комплексно сопряженное число.

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

В первом варианте воплощения расчета матрицы М пространственного фильтра MMSE рассчитывают Эрмитову матрицу P на основе уравнении Риккати. Эрмитова матрица P может быть выражена следующим образом:

Уравнение (4)

Эрмитова матрица P_i размером TxT может быть определена как:

Уравнение (5)

Лемма обращения матрицы может быть применена к уравнению (5) для получения следующего:

Уравнение (6)

где r_iпредставляет собой скалярную действительную величину. Уравнение (6) называется уравнением Риккати. Матрица P_i может быть инициализирована как P ₀= I. После выполнения R итераций уравнения (6), для i=1,...,R, получают матрицу P_R как матрицу P, или P=P_R.

Уравнение (6) может быть перемножено на определенные коэффициенты для получения следующего:

Уравнение (7)

где матрицу P_i инициализируют как P₀= I, и матрицу P получают как P= P_R. Уравнения (6) и (7) отличаются от решения уравнения (5). Для простоты использовали те же переменные P_i и r_i для обоих уравнений (6) и (7), даже, несмотря на то, что эти переменные имеют разные значения в двух уравнениях. Конечные результаты, полученные по уравнениям (6) и (7), то есть, P_R для уравнения (6) и P_R для уравнения(7), эквивалентны.Однако расчеты для первой итерации уравнения (7) упрощаются, благодаря использованию P₀ в качестве единичной матрицы.

Каждая итерация уравнения (7) может быть выполнена следующим образом:

Уравнение (8a)

Уравнение (8b)

Уравнение (8c)

Уравнение (8d)

где a _i представляет собой вектор промежуточной строки размером 1×T элементов с комплексным значением; и

C представляет собой промежуточную Эрмитову матрицу размером T×T.

В системе (8) уравнений последовательность операций структурирована для эффективного расчета с помощью аппаратных средств. Скалярное значение r_i рассчитывают перед матрицей C_i. Разделение по r_i в уравнении (7) достигается с помощью обращения и умножения. Обращение r_i может быть выполнено параллельно с расчетом C_i. Обращение r_i достигается со сдвигом для нормализации r_i и с применением справочной таблицы для получения обращенного значения r_i. Нормализация r_i может быть компенсирована умножением на C _i.

Матрицу P_i инициализируют как Эрмитову матрицу или P₀=I, и она остается Эрмитовой матрицей во всех следующих итерациях. Следовательно, только верхнюю (или нижнюю) диагональную матрицу требуется рассчитывать для каждой итерации. После окончания R итераций получают матрицу P, как P= P_R. Матрица пространственного фильтра MMSE может затем быть рассчитана следующим образом:

Уравнение (9)

На фиг.1 показан процесс 100 расчета матрицы M пространственного фильтра MMSE, основанной на первом варианте воплощения. Матрица P _i инициализирована как P₀=1 (блок 112), и индекс i используется для обозначения номера итерации, инициализируется как i=1 (блок 114). Затем выполняют R итераций уравнения Риккати.

Каждую итерацию уравнения Риккати выполняют с помощью блока 120. Для i-й итерации вектор a_i промежуточной строки рассчитывают на основе заявленного вектора h_i отклика канала и Эрмитовой матрицы P_i+1 предшествующей итерации, как показано в уравнении (8a) (блок 122). Скалярное значение r_i рассчитываютна основе дисперсии ² _n шумов, вектора a_i промежуточной строки и вектора h_i строки отклика канала, как показано в уравнении (8b) (блок 124). Скалярная величина r_i после этого становится обращенной (блок 126). Промежуточную матрицу C_i рассчитывают на основе a_i промежуточной строки, как показано в уравнении (8c) (блок 128). Матрицу P_i затем обновляют на основе инвертированной скалярной величины r _iи промежуточной матрицы C_i , как показано в уравнении (8d) (блок 130).

Затем определяют, были ли выполнены все R итераций (блок 132). В случае отрицательного ответа, выполняют последовательное приращение индекса i(блок 134), и процесс возвращается к блоку 122 для выполнения другой итерации. В противном случае, если все R итераций были выполнены, рассчитывают матрицу М пространственного фильтра MMSE на основе Эрмитовой матрицы P_R для последней итерации, матрицы H отклика каналов и дисперсии ² _n шумов,как показано в уравнении (9) (блок 136). Матрицу М затем можно использовать для пространственной обработки приемника,как показано в уравнении (3).

Во втором варианте воплощения расчета матрицы М пространственного фильтра MMSE определяют Эрмитову матрицу P путем получения квадратного корня P, который представляет собой P^1/2, на основе итерационной процедуры. Пространственная обработка приемника в уравнении (3) может быть выражена следующим образом:

Уравнение (10)

где U= представляет собой дополненную матрицу канала размером (R+T)×T;

U^p представляет собой псевдообращенную матрицу размером T×(R+T), полученную в результате операции обращения или псевдообращения Мура-Пенроуза для U, или U^p=(U^H U)^-1 U^H;

0_Tx1 представляет собой вектор размером T×1, содержащий одни нули; и

H^p _n представляет собой подматрицу размером T×R, содержащую первые R столбцов U^p.

Разложение QR может быть выполнено для матрицы с дополненным каналом следующим образом:

Уравнение (11)

где Q представляет собой матрицу размером (R+T)×T с ортонормированными столбцами;

R представляет собой матрицу размером T×T, которая не является единичной матрицей;

B представляет собой матрицу размером R×T, содержащую первые R строк матрицы Q; и

Q₂ представляет собой матрицу размером T×T, содержащую последние T строк матрицы Q.

QR (КО, квазиобращенное) разложение в уравнении (11) разлагает матрицу дополненного канала на ортогональную матрицу Q и на не единичную матрицу R. Ортогональная матрица Q имеет следующее свойство: Q^H Q=I, что означает, что столбцы ортогональной матрицы являются ортогональными относительно друг друга, и каждый столбец имеет единичную степень. Не единичная матрица представляет собой матрицу, для которой может быть рассчитана обращенная матрица.

Эрмитова матрица P может быть затем выражена как:

Уравнение (12)

R представляет собой разложение Колецкого или квадратный корень матрицы P^-1. Следовательно, P^1/2 равно R^-1 и называется квадратным корнем матрицы P.

Псевдообращенная матрица в уравнении (10) затем может быть выражена как:

Уравнение (13)

Подматрица H^P _n, которая также представляет собой матрицу пространственного фильтра MMSE, затем может быть выражена как:

Уравнение (14)

Уравнение (10) затем может быть выражено как:

Уравнение (15)

Матрицы P^1/2 и B могут быть рассчитаны итерационно следующим образом:

или Уравнение (16)

Уравнение (17)

где Y_i представляет собой матрицу размером (T+R+1)×(T+1), содержащую элементы, полученные на основе P^1/2 _i-1, B_i-1 и h_i ;

_i представляет собой унитарную матрицу преобразования размером (T+1)×(T+1);

Z_i представляет собой преобразованную матрицу размером (T+R+1)×(T+1), содержащую элементы для P_i ^1/2, B_i и r_i ;

e_i представляет собой вектор размером R×1 с единицей (1,0) в качестве i-го элемента и с остальными нулевыми элементами; и

k _i представляет собой вектор размером T×1 и I _i представляет собой вектор R×1, причем оба они являются несущественными.

Матрицы P ^1/2 и B инициализируют как P₀ ^1/2= I и B₀=0_RxT.

Преобразование в уравнение (17) может быть выполнено итерационно, как описано ниже. Для ясности каждую итерацию уравнения (17) называют внешней итерацией. R внешних итераций уравнения (17) выполняют для R векторов h_i отклика канала для i=1,...,R. Для каждой внешней итерации унитарная матрица _i преобразования в уравнении (17) преобразуется в трансформированную матрицу Z_i, содержащую все нули в первой строке, за исключением первого элемента. Первый столбец преобразованной матрицы Z_i содержит r_i ^1/2, k_i и I_i . Последние T столбцов Z_i содержат обновленные P_i ^1/2 и B_i. Первый столбец Z_i не требуется рассчитывать, поскольку только P_i ^1/2 и B_i используются в следующей итерации. P_i ^1/2 представляет собой верхнюю треугольную матрицу. После окончания R внешних итераций получают P_R ^1/2 как P^1/2, и B _R получают как B. Матрица М пространственного фильтра MMSE может быть затем рассчитана на основе P ^1/2 и B, как представлено в уравнении (14).

Для каждой внешней итерации i преобразование по уравнению (17) может выполняться путем последовательного обнуления одного элемента в первой строке Y_i одновременно с 2×2 поворотами Гивенса. T внутренних итераций поворота Гивенса могут быть выполнены для обнуления последних T элементов в первой строке Y_i.

Для каждой внешней итерации i, матрица Y_i,j может быть инициализирована как Y_i1=Yi. Для каждой внутренней итерации j для j=1,...,T, внешней итерации i, первоначально формируют подматрицу Y _i,j размером (T+R+1)×2, содержащую первый и (j+1)-й столбцы Y_i,j. Затем выполняют поворот Гивенса для подматрицы Y _i,j для генерирования подматрицы Y" _i,j размером (T+R+1)×2, содержащей ноль во втором элементе в первой строке. Поворот Гивенса может быть выражен как:

Уравнение (18)

где G_i,j представляет собой матрицу поворота Гивенса размером 2×2 для j-й внутренней итерации i-й внешней итерации, которая описана ниже. Матрицу Y_i,j+1 затем формируют вначале путем установки Y_i,j+1= Y_i,j, затем замены первого столбца Y _i,j+1 первым столбцом Y"_i,j и затем замены (j+1)-го столбца матрицы Y_i,j+1 вторым столбцом Y"_i,j. Поворот Гивенса, таким образом, модифицирует только два столбца Y_i,j j-й внутренней итерации для получения Y_i,j+1 для следующей внутренней итерации. Поворот Гивенса может быть выполнен на месте двух столбцов Y_i для каждой внутренней итерации, в результате чего промежуточные матрицы Y_i,j,Y _i,j,Y"_i,j иY_i,j+1 не нужны и описаны выше для ясности изложения.

Для j-й внутренней итерации i-й внешней итерации матрица G_i,j поворота Гивенса определяется на основе первого элемента (который всегда представляет собой действительную величину) и (j+1)-го элемента в первой строке матрицы Y_i,j. Первый элемент может быть обозначен как , и (j+1)-й элемент может быть обозначен как b-e^j . Матрица G_i,j поворота Гивенса затем может быть получена следующим образом:

Уравнение (19)

где c= и s= для уравнения (19).

На фиг.2 показан процесс 200, предназначенный для расчета матрицы М пространственного фильтра MMSE на основе второго варианта воплощения. Матрицу P_i ^1/2 инициализируют как P₀ ^1/2= I, и матрицу B_i инициализируют как B₀=0 (блок 212). Индекс i для обозначения номера внешней итерации инициализируют как i=1, и индекс j, используемый для обозначения номера внутренней итерации, инициализируют как j=1 (блок 214). Затем выполняют R внешних итераций унитарного преобразования в соответствии с уравнением (17) (блок 220).

Для i-й внешней итерации вначале формируют матрицу Y_i с вектором h_i строки отклика канала и матрицы P_i-1 ^1/2 и B_i-1, как показано в уравнении (17) (блок 222). Матрица Y_i затем обозначается как матрица Y_i,j для внутренних итераций (блок 224). T внутренних итераций поворота Гивенса затем выполняют для матрицы Y_i,j (блоке 230).

Для j-й внутренней итерации получают матрицу G_i,j поворота Гивенса на основе первого и (j+1)-го элементов в первой строке Y_i,j, как показано в уравнении (19) (блок 232). Матрицу G_i,j поворота Гивенса затем применяют для первого и (j+1)-го столбцов Y _i,j для получения Y_i,j+1, как показано в уравнении (18) (блок 234). Затем определяют, были ли выполнены все T внутренние итерации (блок 236). Если ответ представляет собой "Нет", тогда индекс j увеличивают на единицу (блок 238), и обработка возвращается к блоку 232 для выполнения другой внутренней итерации.

Если все T внутренние итерации были выполнены для текущей внешней итерации, и ответ представляет собой "Да" для блока 236, тогда последняя Y_i,j+1 равна Z_i в уравнении (17). Обновленные матрицы P_i ^1/2 и B_i получают из последней Y_i,j+1 (блок 240). Затем определяют, были ли выполнены все R внешних итераций (блок 242). Если будет получен ответ "Нет", тогда индекс i увеличивают на единицу, и индекс j повторно инициализируют как j=1 (блок 244). Обработка затем возвращается в блок 222 для выполнения другой внешней итерации с P_i ^1/2 и B_i. В противном случае, если все R внешние итерации были выполнены, и ответ представляет собой "Да" для блока 242, тогда рассчитывают матрицу М пространственного фильтра MMSE, на основе P _i ^1/2 и B_i, как показано в уравнении (14) (блок 246). Матрицу М можно затем использовать для пространственной обработки приемника, как показано в уравнении (15).

В третьем варианте воплощения расчета матрицы М пространственного фильтра MMSE выполняют разложение по собственным значениям P^-1 следующим образом:

Уравнение (20)

где V представляет собой унитарную матрицу T×T собственных векторов; и

представляет собой диагональную матрицу размером T×T с действительнымисобственными значениями вдоль диагонали.

Разложение по собственным значениям Эрмитовой матрицы X_2x2 размером 2×2 может быть получено с использованием разных методик. В варианте воплощения разложение по собственным значениям X_2x2 получают путем выполнения сложного поворота Якоби для X_2x2 , для получения матрицы V_2x2 размером 2×2 собственных векторов X_2x2. Элементы X _2x2 и V_2x2 могут быть заданы как:

Уравнение (21)

Элементы V_2x2 могут быть рассчитаны непосредственно из элементов X_2x2 следующим образом:

Уравнение (22a)

Уравнение (22b)

Уравнение (22c)

Уравнение (22d)

Уравнение (22e)

Уравнение (22f)

Уравнение(22g)

Уравнение (22h)

Уравнение (22i)

Уравнение (22j)

Уравнение (22k)

Разложение по собственным значениям Эрмитовой матрицы X размером T×T, которая больше чем 2×2, может быть выполнено в итеративном процессе. В таком итеративном процессе неоднократно используется поворот Якоби для обнуления недиагональных элементов в X. Для итеративного процесса индекс i обозначает номер итерации и инициализируется как i=1. X представляет собой Эрмитову матрицу размером T×T, которая должна быть разложена, и устанавливается как X=P^-1. Матрица D_i представляет собой аппроксимацию диагональной матрицы в уравнении (20) и инициализируется как D ₀=X. Матрица V представляет собой аппроксимацию унитарной матрицы V в уравнении (20) и инициализируется как V₀=I.

Одиночная итерация поворота Якоби для обновления матриц D_i и V_i может быть выполнена следующим образом. Вначале Эрмитову матрицу D_pq размером 2×2 формируют на основе текущей матрицы D_i следующим образом:

Уравнение (23)

где d _p,q представляет собой элемент в местоположении (p,q) матрицы D_i, p {1,...,T}, g {1,...,T} и p q. D_pq представляетсобой подматрицу D_i размером 2×2, и четыре элемента D _pq представляют собой четыре элемента в местоположениях (p, p), (p, q), (q, p) и (q, q) матрицы D_i. Индексы p и q могут быть выбраны, как описано ниже.

Затем выполняют разложение по собственным значениям D _pq, как показано в уравнении (22), для получения унитарной матрицы V_pq размером 2×2 собственных векторов D_pq. Для разложения по собственным значениям D_pq,X_2x2 в уравнении (21) заменяют на D_pq, и V_2x2 из уравнения (22j) или (22k) представляют как V_pq .

Матрица T_pq сложного поворота Якоби размером TxT затем формируется с V_pq,T _pq и представляет собой единичную матрицу с четырьмя элементами в местоположениях (p, p), (p, q), (q, p) и (q, q), которые заменены элементами v_1,1, v_1,2, v_2,1 и v_2,2, соответственно, в матрице Y_pq.

Матрица D_i затем обновляется следующим образом:

Уравнение (24)

Уравнение (24) обнуляет два недиагональных элемента в местоположениях (p, q) и (q, p)в матрице D_i. Расчет может изменять значения других недиагональных элементов в D _i.

Матрицу V_i также обновляют следующим образом:

Уравнение (25)

V _i можно рассматривать как матрицу кумулятивного преобразования, которая содержит все матрицы T_pq поворота Якоби, использовавшиеся для D_i.

Каждая итерация с поворотом Якоби обнуляет два недиагональных элемента матрицы Di. Множество итераций поворота Якоби могут быть выполнены для разных значений индексов p и q, для обнуления всех недиагональных элементов D_i. Один проход по всем возможным значениям индексов p и q может быть выполнен следующим образом. Индекс p последовательно изменяется от 1 до T-1 с приращением единица. Для каждого значения p индекс q последовательно изменяется от p+1 до T с приращением единица. Поворот Якоби выполняют для каждой разной комбинации значений p и q. Множество проходов могут быть выполнены до тех пор, пока D_i и V _i не станут представлять достаточно точные оценки и V, соответственно.

Уравнение (20) может быть переписано следующим образом:

Уравнение (26)

где ^-1 представляет собой диагональную матрицу, элементы которой представляют собой обращенные значения соответствующих элементов . Разложение по собственным величинам X=P ^-1 представляет оценки и V. может быть инвертирована для получения ^-1.

Матрица пространственного фильтра MMSE может быть затем рассчитана следующим образом:

Уравнение (27)

На фиг.3 представлен процесс 300, предназначенный для расчета матрицы М пространственного фильтра MMSE, на основе третьего варианта воплощения. Эрмитову матрицу P^-1 первоначально получают на основе матрицы H отклика канала, как показано в уравнении (20) (блок 312). Затем выполняют разложение по собственным величинам P^-1 для получения унитарной матрицы V и диагональной матрицы , как также показано в уравнении (20) (блок 314). Разложение по собственным величинам может итерационно выполняться с множеством поворотов Якоби, как описано выше. Матрицу М пространственного фильтра MMSE затем получают на основе унитарной матрицы V , диагональной матрицы и матрицы H отклика канала, как показано в уравнении (27) (блок 316).

Матрица М пространственного фильтра MMSE, полученная на основе каждого из описанных выше вариантов воплощения, представляет собой смещенное решение MMSE. Смещенную матрицу М пространственного фильтра можно масштабировать с помощью диагональной матрицы D_mmse для получения несмещенной матрицы MMSE M _mmse пространственного фильтра. Матрица D _mmse может быть получена как D_mmse=[diag[ M H]]^-1, где diag[М H] представляет собой диагональную матрицу, содержащую диагональные элементы М H.

Описанные выше расчеты также можно использовать для получения матриц пространственного фильтра для методики с обращением в нуль не значащих коэффициентов (ZF) (которая также называется методикой обращения матрицы с корреляцией канала (CCMI, ОМКК)), методики комбинирования максимального отношения (MRC, КМО) и так далее. Например, приемная станция может выполнять пространственную обработку приемника с обращением в нуль не значащих коэффициентов и MRC, следующим образом:

Уравнение (28)

Уравнение (29)

где M_zf представляет собой матрицу пространственного фильтра размером T×R с обращенными в ноль не значащими коэффициентами;

M_mrc представляет собой матрицу пространственного фильтра MRC с размером T×R;

P_zf=(H^H H)^-1 представляет собой Эрмитову матрицу размером T×T; и

[diag(P_zf)] представляет собой диагональную матрицу размером T×T, содержащую диагональные элементы P_zf.

Обращение матрицы необходимо для непосредственного расчета P_zf . P_zf могут быть рассчитаны с использованием вариантов воплощения, описанных выше для матрицы пространственного фильтра MMSE.

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

Уравнение (30)

где x представляет собой вектор размером T×1 с T символами передачи, которые должны быть переданы через T передающих антенн; и

W представляет собой матрицу передачи размером T×S. Матрица W передачи может представлять собой (1) матрицу с правыми единичными векторами, полученными путем выполнения разложения единичного значения H, (2) матрицу собственных векторов, полученную путем выполнения разложения по собственным величинам H^HH, или (3) управляющую матрицу, выбранную для пространственного распределения символов модуляции по S пространственным каналам в канале MIMO. Матрица H_eff эффективного отклика канала, наблюдаемого по символам модуляции, затем может быть задана как H _eff =H W. Описанная выше комбинация может выполняться на основе H_eff вместо H.

Для ясности приведенного выше описания, в нем приведена система MIMO с одной несущей, с одной подполосой. Для системы MIMO с множеством несущих матрица H(k) отклика канала может быть получена для каждой подполосы k, представляющей интерес. Матрица М(k) пространственного фильтра затем может быть получена для каждой подполосы k на основе матрицы H(k) отклика канала для этой подполосы.

Описанные выше расчеты для матрицы пространственного фильтра могут быть выполнены с использованием процессоров различного типа, таких как процессор с плавающей запятой, процессор с фиксированной запятой, процессор цифрового компьютера поворота координат (CORDIC), таблица поиска и так далее, или их комбинации. Процессор CORDIC воплощает итеративный алгоритм, который обеспечивает возможность выполнения быстрого расчета с использованием аппаратных средств тригонометрических функций, таких как синус, косинус, магнитуда и фаза, с использованием простых аппаратных средств сдвига и сложения/вычитания. Процессор CORDIC может итерационно рассчитывать каждую из переменных r, c₁ и s₁ системы (22) уравнений с большим количеством итераций, что позволяет достичь более высокой точности для переменной.

На фиг.4 показана блок-схема точки 410 доступа и терминала 450 пользователя в системе 400 MIMO. Точка 410 доступа оборудована N_ap антеннами _, и терминал 450 пользователя оборудован N_ut антеннами, где N_ap>1 и N_ut>1. По нисходящему каналу передачи, в точке 410 доступа процессор 414 передаваемых (TX) данных принимает данные трафика из источника 412 данных и другие данные из контроллера/процессора 430. Процессор 414 TX данных форматирует, кодирует, выполняет перемежение, модулирует данные и генерирует символы данных, которые представляют собой символы модуляции для данных. Пространственный процессор 420 TX мультиплексирует символы данных с пилотными символами, выполняет пространственную обработку с матрицей W передачи, если это применимо, и предоставляет N_ap потоков передаваемых символов. Каждый модуль 422 передатчиков (TMTR) обрабатывает соответствующий поток передаваемых символов и генерирует модулированный сигнал для нисходящего канала передачи. Модулированные сигналы для нисходящего канала передачи N_ap из модулей 422a-422ap передатчика передают через антенны 424a-424ap, соответственно.

В терминале 450 пользователя N_ut антенн 452a-452ut принимают переданные модулированные сигналы нисходящего канала передачи, и каждая антенна передает принятый сигнал в соответствующий модуль (RCVR) 454 приемника. Каждый модуль 454 приемника выполняет обработку, взаимодополняющую обработку, выполнявшуюся модулями 422 передачи, и предоставляет принятые пилотные символы и принятые символы данных. Блок/процессор 478 оценки канала обрабатывает принятые пилотные символы и предоставляет оценку отклика H _dn канала нисходящего канала передачи. Процессор 480 получает матрицу M_dn пространственного фильтра нисходящего канала передачи на основе H_dn, и используя любой из вариантов воплощения, описанных выше. Пространственный процессор 460 приемника (RX) выполняет пространственную обработку приемника (или пространственную согласованную фильтрацию) для принятых символов данных из всех N_ut модулей 454a-454ut приемникас матрицей M_dn пространственного фильтра нисходящего канала передачи и предоставляет детектируемые символы данных, которые являются оценками символов данных, переданных точкой 410 доступа. Процессор 470 приемника обрабатывает (например, выполняет обратное отображение символа, устраняет перемежение и декодирует) детектированные символы данных и предоставляет декодированные данные в приемник 472 данных и/или контроллер 480.

Обработка для восходящего канала передачи может быть такой же или может отличаться от обработки для нисходящего канала передачи. Данные из источника 486 данных и сигналы из контроллера 480 обрабатывают (например, кодируют, выполняют перемежение и модулируют) с помощью процессора 488 TX данных, мультиплексированных с пилотными символами и, возможно, пространственно обработанных пространственным процессором 490 TX. Символы передачи из пространственного процессора 490 TX дополнительно обрабатывают с помощью модулей 454a-454ut передатчика для генерирования N_ut модулированных сигналов восходящего канала передачи, которые передают через антенны 452a-452ut.

В точке 410 доступа модулированные сигналы восходящего канала передачи принимают с помощью антенн 424a-424ap и обрабатывают с помощью модулей 422a-422ap приемника для генерирования принятых пилотных символов и принятых символов данных, для передачи по восходящему каналу передачи. Блок/процессор 428 оценки канала обрабатывает принимаемые пилотные символы и предоставляет оценку отклика H_up канала восходящей передачи. Процессор 430 получает матрицу M_up пространственного фильтра восходящего канала передачи, и используя один из вариантов выполнения, описанных выше. Пространственный процессор 440 RX выполняет пространственную обработку приемника для принимаемых символов данных с матрицей M_up пространственного фильтра восходящего канала передачи и предоставляет детектированные символы данных. Процессор 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). Модуль памяти может быть воплощен внутри процессора или вне процессора.

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