СПОСОБ ИТЕРАТИВНОГО ДЕТЕКТИРОВАНИЯ И ДЕКОДИРОВАНИЯ СИГНАЛА В СИСТЕМАХ СВЯЗИ С MIMO КАНАЛОМ Российский патент 2014 года по МПК H04L27/04 

Изобретение относится к телекоммуникационным технологиям, а более конкретно - к способам приема сигнала в системе передачи, использующей технологию на основе MIMO (Multiple-In Multiple-Out).

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

Одним из наиболее распространенных способов передачи сигналов по MIMO каналам является пространственное мультиплексирование, известное также как V-BLAST архитектура (Vertical Bell Laboratories Layered Space Time architecture) [1], при котором сигнал разбивается на несколько параллельных потоков, передающихся по нескольким пространственным каналам, образующимся между множеством передающих и приемных антенн. MIMO система в частотной области может быть описана матричным уравнением (1)

$$Y=HX+\eta \left(1\right)$$

где Y - M-мерный вектор отсчетов комплексного сигнала на входе MIMO детектора, который можно считать вектором выходных отчетов MIMO канала, X - N-мерный вектор переданных модулированных QAM символов в передатчике, который можно считать вектором входных отсчетов MIMO канала, η - M-мерный вектор отсчетов комплексного шума, H - комплексная канальная матрица MIMO канала размера M×N.

При передаче данных с использованием пространственного мультиплексирования возникает взаимная интерференция между различными потоками данных и поэтому возникает задача разработки эффективных методов детектирования, демодуляции и декодирования сигнала на фоне внешних шумов и взаимной интерференции. Наиболее качественным методом детектирования сигнала считается метод максимального правдоподобия ML (Maximum Likelihood), при котором оценивается вероятность передачи всех возможных комбинаций символов. Главным недостатком данного метода является его высокая сложность, которая растет экспоненциально с ростом количества параллельно передающихся потоков. Существует множество альтернативных методов детектирования сигнала за счет подавления взаимной интерференции, наиболее известными из которых являются метод Zero Forcing (минимизация усредненной мощности интерференции на приеме) и MMSE (Minimum Mean Square Error - минимизация средне-квадратичной ошибки оценивания). При использовании метода MMSE уравнение (1) преобразуется к виду (2):

$$Y=G_{MMSE}^{-1}{\tilde{X}}_{MMSE}$$

где $$G_{MMSE}={\left(H^{H}H+\frac{{\sigma }_{\eta }^2}{{\sigma }_s^2}{I}_N\right)}^{-1},$$ H^H - матрица линейного преобразования MMSE фильтра, $${\sigma }_{\eta }^2$$ - дисперсия шума, $${\sigma }_s^2$$ - среднее значение мощности сигнала, I_N - единичная матрица размера N×N, $${\tilde{X}}_{MMSE}$$ - вектор MMSE оценок переданных QAM символов (MMSE решение).

Методы Zero Forcing и MMSE обладают сравнительно низкой сложностью реализации, но при этом значительно проигрывают методу ML по точности оценки принятых символов и, в конечном итоге, по вероятности приема неверных бит в сообщении (bit error rate).

Помимо пространственного мультиплексирования современные цифровые системы связи широко применяют также другие методы передачи и приема сигналов, улучшающие эффективность и качество передачи данных. В частности, в них широко применяются методы канального кодирования цифровых данных в передающем устройстве и соответственно их декодирования в приемнике. Помимо этого применяют также скремблирование данных и их перемежение, обозначаемое термином «интерливинг» (interleaving). Интерливинг позволяет более эффективно распределить транслируемый сигнал в пространственно-частотно-временном континууме, что улучшает качество передачи данных.

Известны способы итеративного детектирования и декодирования сигнала, в которых детектирование и декодирование производится совместно в рамках итерационного процесса. Такие способы принято называть также турбо (turbo) обработкой. Практически во всех современных системах связи цифровые данные подвергаются процедуре кодирования в канальном кодере, где исходная входная последовательность информационных бит преобразуется в выходную последовательность кодированных бит, содержащую избыточные проверочные биты. Известно множество типов канальных кодов, применяемых в современных системах цифровой связи. Наиболее широко применяемые на практике коды - сверточные, турбо коды, коды с низкой плотностью проверки на четность (LDPC, Low Density Parity Check Codes). После кодирования данные могут также подвергаться процедуре интерливинга, модуляции, пространственного демультиплексирования, преобразованию цифрового сигнала в аналоговый, преобразованию частоты, усилению и излучению в пространство посредством нескольких передающих антенн. Типовая схема генерации MIMO сигнала на передающей стороне (transmit system) представлена на Фиг.1, где входные информационные биты преобразуются в канальном кодере 101, образованные кодированные биты подвергаются интерливингу в модуле 103, модулируются в модуле 105, демультиплексируются в модуле 107 для разделения на несколько пространственных потоков и, после преобразования в аналоговый сигнал, преобразования частоты и усиления в передатчиках 109, излучаются через несколько передающих антенн. Соответственно, на приемной стороне (receiver system), типовая схема которой представлена на Фиг.2, излученные сигналы принимаются посредством нескольких приемных антенн, усиливаются, преобразуются по частоте, фильтруются и преобразуются в цифровой сигнал в соответствующих приемниках 201, после чего поступают на MIMO детектор 203, где осуществляется предварительное MIMO детектирование, после чего сигналы демодулируются в модуле 205, проходят через деинтерливер 207 и декодируются в модуле 209.

Схема приемной системы, использующей метод итеративного детектирования и декодирования, приведена на Фиг.3. Многомерный входной сигнал детектируется в MIMO детекторе 303, затем демодулируется в модуле 305, где производится первая оценка вероятности кодированных бит, которая, как правило, выражается посредством логарифмического отношения правдоподобия (LLR) или логарифмического отношения апостериорных вероятностей. Детектор и демодулятор (его также часто именуют в англоязычной литературе Demapper, соответственно в передатчике модулятор называют Mapper) могут быть объединены в один модуль 304, куда наряду с входными сигналами может поступать априорная информация о кодированных битах (либо QAM символах, содержащих кодированные биты). Оценку апостериорной вероятности кодированных бит (или логарифмического отношения апостериорных вероятностей) принято также называть мягкими решениями. После деинтерливинга в модуле 307, мягкие решения направляют в канальный декодер 309, который в соответствии с параметрами используемого кода производит оценку апостериорной вероятности (логарифмического отношения апостериорных вероятностей) информационных бит, формируемых на первом (основном) выходе декодера, а также уточняет значения вероятностей кодированных бит (формирует уточненные мягкие оценки кодированных бит) на втором выходе декодера. Мягкие оценки информационных бит с первого выхода декодера 309 преобразуются в жесткие оценки модулем 313, выход которого является общим выходом системы связи с MIMO каналом. Уточненные мягкие оценки кодированных бит со второго выхода декодера 309 после удаления из них информации, поступающей на вход декодера в виде мягких оценок кодированных бит, подвергают процедуре интерливинга в модуле 311 и направляют на дополнительный вход детектора/демодулятора 304 в качестве априорной информации для следующей итерации. В представленной схеме детектор и демодулятор объединены в один общий модуль, поскольку в ряде применений они оказываются тесно переплетены. Например, при применении метода сферического детектирования [5] априорная информация поступает на вход MIMO детектора для уменьшения количества символов-кандидатов, участвующих в расчете апостериорной информации. В других случаях (например, при ML MIMO детекторе, см. [3]), априорная информация поступает на вход демодулятора. В следующей итерации детектор и демодулятор 304 производят новую оценку апостериорной вероятности кодированных бит (логарифмического отношения апостериорных вероятностей), которая после удаления из нее априорной информации (логарифмического отношения априорных вероятностей) вновь подвергается деинтерливингу и декодированию. Повторение итеративной процедуры детектирования - декодирования приводит к повышению надежности оценок кодированных бит, так что окончательная последовательность жестких оценок информационных бит совпадает с переданной последовательностью с большей вероятностью.

Известны различные способы реализации итеративного детектирования и декодирования для систем с MIMO каналом [3-5]. Наиболее близким к заявленному изобретению является способ, предложенный в патентной заявке [3]. Данный способ итеративного детектирования и декодирования выбран в качестве прототипа заявленного изобретения.

Схема приема и обработки сигнала в [3] соответствует схеме на Фиг.3. При этом оценка апостериорной вероятности после MIMO детектора производится согласно алгоритму ML:

$$L\left(b_{n,k}\right)=\ln \left(\frac{{\displaystyle \sum _{x:f\left(b_{n,k}=1\right)}{\displaystyle \prod _{m=1}^M\exp \left({\beta }_m\right)}}{\displaystyle \prod _{p=1}^N{\displaystyle \prod _{\underset{t\ne kwhenp=n}{t=1}}^K\mathrm{Pr}\left(b_{p,t}\right)}}}{{\displaystyle \sum _{x:f\left(b_{n,k}=-1\right)}{\displaystyle \prod _{m=1}^M\exp \left({\beta }_m\right)}}{\displaystyle \prod _{p=1}^N{\displaystyle \prod _{\underset{t\ne kwhenp=n}{t=1}}^K\mathrm{Pr}\left(b_{p,t}\right)}}}\right)+\ln \left(\frac{\mathrm{Pr}\left(b_{n,k}\right)=+1}{\mathrm{Pr}\left(b_{n,k}\right)=-1}\right)\left(3\right)$$

где L(b_n,k) - значение логарифмического отношения правдоподобия, n=1, 2,… N, k=1, 2,… K, K - число бит в QAM символе, определяемым созвездием x:f(b₁,…b_K), суммирование проводится по всем возможным комбинациям символов, включающих бит, равный +1 или -1 (+1 и -1 соответствует изначальным значениям бита 0 и 1), β_m - значение евклидовых расстояний по каждой приемной антенне, определяемого выражением (4):

$${\beta }_m=\frac{1}{2{\sigma }_{\eta }^2}{\Vert Y_m-H_{m}X\Vert }^2\left(4\right)$$

Второе слагаемое в (3) характеризует логарифмическое отношение априорных вероятностей.

Число возможных символов, участвующих в подсчете L(b_n,k), экспоненциально растет с ростом числа передающих антенн 2^NK, что делает такой подход трудно реализуемым. В [3] предложена также упрощенная процедура MIMO детектирования с использованием процедуры подавления интерференции. Для этого вектор Y умножают слева на матрицу Q_m размера (M-N+1)×M, которая зануляет все переданные символы (символы от всех передающих антенн) кроме одного. В результате для подсчета LLR необходимо вычислять только (M-N+1)K евклидовых расстояний. Нетрудно видеть, что при M=N (число передающих антенн равно числу приемных), данный подход эквивалентен Zero Forcing, который демонстрирует значительную деградацию по сравнению с ML и который серьезно проигрывает другому относительно простому методу линейной фильтрации входного сигнала MMSE.

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

Технический результат достигается за счет разработки нового итеративного способа детектирования и декодирования сигнала в системах связи с МIМО каналом, который включает в себя выполнение следующих операций:

- принимают сигналы, модулированные переданными QAM (Quadrature Amplitude Modulation) символами через несколько приемных антенн;

- осуществляют частотное преобразование принятых сигналов на нулевую несущую частоту;

- формируют отсчеты принятых и преобразованных сигналов путем квантования и дискретизации, которые рассматриваются как выходные отсчеты MIMO канала и описываются формулой:

Y=HX+η,

где Y - M-мерный вектор выходных отсчетов MIMO канала, Х - N-мерный вектор входных отсчетов MIMO канала, η - M-мерный вектор отсчетов шума, H - канальная матрица MIMO канала размера М×N, при этом полагают, что частота дискретизации равна частоте следования QAM символов, поэтому вектор входных отсчетов MIMO канала X является также вектором переданных QAM символов,

- по полученным выходным отсчетам MIMO канала оценивают

канальную матрицу MIMO канала H, дисперсию шума $${\sigma }_{\eta }^2$$ и среднюю

мощность сигнала $${\sigma }_s^2$$ ;

- детектируют QAM символы на основе выходных отсчетов MIMO канала, производя оценку апостериорной вероятности (логарифмического отношения апостериорных вероятностей) кодированных бит;

- декодируют информационные биты на основе оценки апостериорной вероятности (логарифмического отношения апостериорных вероятностей) кодированных бит, полученной в детекторе QAM символов, и, в соответствии с параметрами канального кода, производят оценку апостериорной вероятности информационных бит и уточняют оценки апостериорной вероятности (логарифмического отношения апостериорных вероятностей) кодированных бит;

- на основе информации о вероятности кодированных бит на выходе канального декодера формируют оценки априорной вероятности о принятых QAM символах и направляют их на дополнительный вход MIMO детектора в качестве априорной информации для формирования новых оценок апостериорных вероятностей (логарифмических отношений апостериорных вероятностей) кодированных и информационных бит на следующем цикле итеративного процесса;

- после завершения всех итераций на основе оценки апостериорной вероятности информационных бит на основном выходе канального декодера формируют окончательную последовательность жестких оценок информационных бит, восстанавливая оригинальные данные;

отличающегося тем, что

- детектирование QAM символов осуществляют в два этапа, причем

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

$$\begin{array}{c}{\tilde{X}}_{MMSE}={\overline{X}}_{pr}+G_{MMSE}\left(Y-H{\overline{X}}_{pr}\right)\\ G_{MMSE}=V_{pr}H\text{'}{\left(H{V}_{pr}H\text{'}+\frac{{\sigma }_{\eta }^2}{{\sigma }_s^2}I\right)}^{-1}\end{array}\left(5\right)$$

$${\overline{X}}_{pr}$$ - N-мерный вектор априорных математических ожиданий переданных QAM символов, V_pr - диагональная матрица, характеризующая априорную дисперсию переданных QAM символов, нормированную средней мощностью сигнала $${\sigma }_s^2$$ ,

для каждого компонента $${\tilde{x}}_{MMSE,n}$$ (n=-1, 2,…N) вектора MMSE оценок $${\tilde{X}}_{MMSE}$$ вычисляют дисперсию ошибки оценивания $${\sigma }_{{}_{MMSE,n}}^2$$ , которая совпадает с соответствующим диагональным элементом корреляционной матрицы Кд/ж, вычисляемой по формуле

$$V_{MMSE}=V_{pr}-G_{MMSE}H{V}_{pr}\left(6\right)$$ ,

на втором этапе путем линейного преобразование каждого n-го (где n=1, 2,…N) компонента $${\tilde{x}}_{MMSE,n}$$ вектора MMSE оценок $${\tilde{X}}_{MMSE}$$ и соответствующего ему априорного математического ожидания $$\overline{x}pr,n$$ осуществляют коррекцию и находят новую линейную оценку $${\tilde{x}}_{LMS,n}$$ переданного QAM символа и дисперсию ошибки для этой оценки $${\sigma }_{{}_{LMS,n}}^2$$

$$\begin{array}{l}{\tilde{x}}_{LMS,n}=-\frac{{\sigma }_{{}_{MMSE,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}{\overline{x}}_{pr,n}+\frac{{\sigma }_{{}_{pr,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}{\tilde{x}}_{MMSE,n}\left(7\right)\\ {\tilde{\sigma }}_{{}_{LMS,n}}^2=\frac{{\sigma }_{{}_{pr,n}}^2{\sigma }_{{}_{MMSE,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}\end{array}$$

где $${\sigma }_{{}_{pr,n}}^2$$ - априорная дисперсия n-го переданного QAM символа,

- осуществляют демодуляцию оценок QAM символов, полученных в детекторе, и вычисляют логарифмическое отношение апостериорных вероятностей кодированных бит в соответствии с выражением:

$${\lambda }_{n,k}=\ln \left(\frac{{\displaystyle \sum _{x:f\left(b_{n,k}=1\right)}\exp }\left(-\frac{1}{{\tilde{\sigma }}_{LMS,n}^2}{\left|{\tilde{x}}_{LMS,n}-x\left(b_{n\mathrm{,1}},b_{n\mathrm{,2}}\dots ,b_{n,K}\right)\right|}^2\right){\displaystyle \prod _{t=1}^K\mathrm{Pr}\left(b_{n,t}\right)}}{{\displaystyle \sum _{x:f\left(b_{n,k}=-1\right)}\exp }\left(-\frac{1}{{\tilde{\sigma }}_{LMS,n}^2}{\left|{\tilde{x}}_{LMS,n}-x\left(b_{n\mathrm{,1}},b_{n\mathrm{,2}}\dots b_{n,K}\right)\right|}^2\right){\displaystyle \prod _{t=1}^K\mathrm{Pr}\left(b_{n,t}\right)}}\right)\left(9\right)$$

где x(b₁,… b_K) - табличная функция с числом состояний 2, описывающая сигнальное созвездие переданного QAM символа в зависимости от комбинации кодированных бит b_n,k∈{-1;1}, k=1, 2,…K, n=1, 2,…N, Pr(b_n,k) - априорная вероятность k-го кодированного бита в n-м QAM символе, полученная на предыдущей итерации, λ(b_n,k) - логарифмическое отношение апостериорных вероятностей для k-го бита в n-м QAM символе, K - число бит в одном QAM символе;

- удаляют из вычисленного в демодуляторе логарифмического отношения апостериорных вероятностей кодированных бит априорную информацию для каждого кодированного бита P_pr(b_n,k) путем вычитания логарифмического отношения априорных вероятностей

$${\stackrel{⌢}{\lambda }}_{n,k}={\lambda }_{n,k}-{\lambda }_{n,k,pr}\left(10\right)$$

где $${\lambda }_{n,k,pr}=\ln \left(\frac{\mathrm{Pr}\left(b_{n,k}\right)=+1}{\mathrm{Pr}\left(b_{n,k}\right)=-1}\right)$$ и на основании этой информации осуществляют декодирование информационных бит в соответствии с параметрами канального кода;

- выделяют из уточненной апостериорной вероятности кодированных бит, полученную после декодирования в виде уточненного логарифмического отношения вероятностей $${\tilde{\lambda }}_{n,k}$$ , внешнюю информацию (вероятность), удаляя из нее информацию, поступающую на вход декодера в виде мягких оценок кодированных бит (логарифмического отношения апостериорных вероятностей) $${\stackrel{⌢}{\lambda }}_{n,k}$$ , и на основании этой информации формируют оценку априорной вероятности, поступающую в QAM демодулятор на следующей итерации

$$P\left(b_{n,k}\right)=\frac{e^{b_{n,k}\frac{{\lambda }_{n,k,pr}}{2}}}{e^{\frac{{\lambda }_{n,k,pr}}{2}}+e^{\frac{{\lambda }_{n,k,pr}}{2}}}\left(11\right)$$

$${\lambda }_{n,k,pr}={\tilde{\lambda }}_{n,k}-{\stackrel{⌢}{\lambda }}_{n,k}$$

- на основе уточненных оценок апостериорной вероятности кодированных бит на выходе канального декодера после их интерливинга формируют оценки переданных QAM символов путем весового сложения значений сигнального QAM созвездия, где веса определяются уточненными апостериорными вероятностями кодированных бит после интерливинга в соответствии с выражением:

где

и вычисляют дисперсии этих оценок

- сравнивают полученные дисперсии $${\sigma }_n^2$$ с оценками дисперсий после детектирования QAM символов $${\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , полученными на предыдущей итерации, и если новая оценка дисперсии оказывается меньше дисперсии детектирования, т.е. $${\sigma }_n^2<{\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , то осуществляют коррекцию априорного математического ожидания и априорной дисперсии

$$\begin{array}{l}{\tilde{x}}_{pr,n}=-\frac{{\sigma }_n^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}{\tilde{x}}_{LMS,n}+\frac{{\tilde{\sigma }}_{{}_{LMS,n}}^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}{\widehat{x}}_n\left(14\right)\\ {\tilde{\sigma }}_{{}_{pr,n}}^2=\frac{{\sigma }_n^2{\tilde{\sigma }}_{{}_{LMS,n}}^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}\end{array}$$

если дисперсия $${\sigma }_n^2$$ больше либо равна дисперсии $${\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , т.е. $${\sigma }_n^2\ge {\tilde{\sigma }}_{{}_{LMS,n}}^2,$$ то априорные математических ожидания и дисперсию оставляют без изменений

$$\begin{array}{c}{\tilde{x}}_{pr,n}={\overline{x}}_{pr,n}\\ {\tilde{\sigma }}_{{}_{pr,n}}^2={\sigma }_{{}_{pr,n}}^2\end{array}$$

- априорные математические ожидания $${\tilde{x}}_{pr,n}$$ и дисперсии $${\tilde{\sigma }}_{{}_{pr,n}}^2$$ , для всех n=1, 2,…N, объединенные в вектор новых априорных математических ожиданий $${\tilde{X}}_{pr}$$ и в новую априорную диагональную корреляционную матрицу $${\tilde{V}}_{pr}$$ , являются входными параметрами MIMO детектора, используемыми на следующей итерации.

Сопоставительный анализ заявляемого способа с другими техническими решениями в данной области техники не позволил выявить признаков, заявленных в отличительной части формулы изобретения. Это позволяет утверждать, что заявляемый способ детектирования сигнала в системах связи с MIMO каналом отвечает критериям новизны и неочевидности.

Для лучшего понимания настоящего изобретения далее приводится его подробное описание с соответствующими графическими материалами.

Фиг.1. Типовая схема генерации сигнала в системе с MIMO каналом (уровень техники).

Фиг.2. Типовая схема приема сигнала в системе с MIMO каналом (уровень техники).

Фиг.3. Типовая схема итеративного детектирования и декодирования сигнала в системе с MIMO каналом (уровень техники).

Фиг.4. Заявляемая схема итеративного детектирования и декодирования сигнала в системе с MIMO каналом.

Фиг.5. Зависимость битовой вероятности ошибки от отношения сигнал/шум для MIMO системы стандарта LTE, модуляция 16 QAM, скорость кодирования ½, канал ЕРА-5. Три типа детекторов/демодуляторов: стандартный последовательный с MMSE детектором (обозначение MMSE), стандартный последовательный с ME детектором (обозначение ME), итеративный детектор/демодулятор, согласно заявляемому способу (обозначение ТР FB).

Фиг.6. Зависимость битовой вероятности ошибки от отношения сигнал/шум для MIMO системы стандарта ЕТЕ, модуляция 16 QAM, скорость кодирования 3/4, канал ЕРА-5. Три типа детекторов/демодуляторов: стандартный последовательный с MMSE детектором (обозначение MMSE), стандартный последовательный с ML детектором (обозначение ML), итеративный детектор/демодулятор, согласно заявляемому способу (обозначение ТР FB).

Способ итеративного детектирования и декодирования сигнала по заявляемому изобретению используется в системах связи с MIMO каналом. Типовые схемы генерации и приема сигнала в системе с MIMO каналом представлены на Фиг.1 и Фиг.2. Их описание приведено выше. На Фиг.4 представлена блок-диаграмма, иллюстрирующая предлагаемый в настоящем изобретении способ итеративного детектирования и декодирования сигнала в системе с MIMO каналом.

На входы MIMO детектора 401 поступают усиленные и преобразованные сигналы, принятые несколькими приемными антеннами, на другие входы MIMO детектора 401 подается также априорная информация о принятом символьном векторе, а именно, математическое ожидание и дисперсия его компонентов, полученная на предыдущей итерации. Детектированный символ и оценка дисперсии ошибки оценивания подаются на вход демодулятора 403, куда также подаются оценки априорной вероятности кодированных бит, полученные на предыдущей итерации. Демодулятор 403 производит оценку апостериорной вероятности (логарифмического отношения апостериорных вероятностей) кодированных бит, из которой в модуле 405 удаляют априорную вероятность и тем самым выделяют внешнюю информацию, полученную в детекторе. Мягкие оценки кодированных бит (логарифмическое отношение вероятностей) подвергают деинтерливингу в модуле 407 и направляют на вход канального декодера 409. Канальный декодер на основании входной вероятности кодированных бит и в соответствии с параметрами канального кода осуществляет оценку апостериорной вероятности информационных бит, направляемую на первый выход канального декодера, а также формирует уточненные мягкие оценки кодированных бит, направляемые на второй выход канального декодера, используемый в линии обратной связи. Мягкие оценки кодированных бит с выхода канального декодера 409 в линии обратной связи подвергаются интерливингу в модуле 411 и подаются на вход модуля 413, который удаляет из них информацию, поступающую на вход декодера в виде мягких оценок кодированных бит, после чего они поступают на вход демодулятора 403, где используются в новом цикле итеративного процесса. Мягкие оценки кодированных бит после интерливера 411 по второй ветви линии обратной связи поступают на вход ремодулятора 415, который производит оценку принятого символьного вектора, а именно его математического ожидания и дисперсии. Полученная оценка поступает на вход модуля 417 коррекции априорных данных, куда также поступает оценка символьного вектора, полученная в MIMO детекторе 401 на предыдущем цикле итеративного процесса. После коррекции априорная оценка символьного вектора поступает на вход MIMO детектора 401 для использования на следующем цикле итеративного процесса. После выполнения нескольких циклов итеративной процедуры детектирования-декодирования мягкие оценки информационных бит с первого выхода канального декодера направляют в модуль 419, который осуществляет финальную жесткую оценку информационных бит. На первом цикле итерации в качестве априорной оценки принятого символьного вектора используют вектор с нулевым математическим ожиданием и единичной дисперсией.

Повторение итеративной процедуры детектирования-декодирования приводит к повышению надежности принятых данных, так что окончательная жесткая оценка информационных бит совпадает с переданной последовательностью (с большей вероятностью).

В процессе декодирования-детектирования принятый сигнал, а именно его цифровое представление проходит преобразования, описываемые рядом математических уравнений. MIMO детектор 401 производит оценку принятого символа, состоящую из двух этапов. На первом этапе формируют MMSE оценку переданного символа $${\tilde{X}}_{MMES}$$ и оценку дисперсии ошибки оценивания с компонентами $${\sigma }_{MMES,n}^2$$ . Для этого используют линейный MMSE фильтр описываемого уравнениями (5). На вход детектора подают вектор принятых по каждой антенне сигналов Y, компоненты канальной матрицы Н, а также априорные данные об оцениваемом векторе, а именно его мат. ожидание и дисперсию. На втором этапе путем линейного преобразование каждого n-го (где N=1, 2,…N) компонента $${\tilde{x}}_{MMES,n}$$ , вектора MMSE оценок $${\tilde{X}}_{MMES}$$ и соответствующего ему априорного математического ожидания $${\overline{x}}_{pr,n}$$ осуществляют коррекцию и находят новую линейную оценку $${\tilde{x}}_{LMS,n}$$ переданного QAM символа и дисперсию ошибки для этой оценки $${\tilde{\sigma }}_{LMS,n}^2$$ согласно формулам (7). Тем самым производят выделение внешней информации, получаемой в процессе детектирования. Полученные в детекторе оценки $${\tilde{x}}_{LMS,n}$$ и $${\tilde{\sigma }}_{LMS,n}^2$$ направляют в демодулятор 403, куда также подаются оценки априорной вероятности кодированных бит, полученные на предыдущей итерации. Демодулятор 403 определяет апостериорные вероятности кодированных бит в соответствии с выражением (9). Из оценок апостериорной вероятности кодированных бит в модуле 405 удаляют априорную информацию, подаваемую в демодулятор, в соответствии с выражением (10), и тем самым выделяют внешнюю информацию, полученную в результате детектирования и демодуляции. Аналогичным образом в модуле 413 производят выделение внешней информации, подученной при декодировании, в соответствии с выражениями (11)

Декодирование производится на основании оценки вероятности кодированных бит, поступающей с выхода деинтерливера 407, и в соответствии с параметрами канального кода. Существуют различные методы декодирования, которые могут использоваться в канальных декодерах. Наиболее известные и широко применяемые на практике: метод Витерби, метод максимума апостериорной вероятности (MAP и Log-MAP декодеры), турбо-декодирование, метод распространения доверия (Believe Propagation) и другие. Применение того или иного метода декодирования обусловлено рядом факторов, главными из которых являются тип используемого кода и эффективность используемого алгоритма (повышение помехозащищенности при приемлемой сложности реализации). Применительно к данному методу детектирования-декодирования выбор того или иного способа декодирования не является принципиальным. Любой способ декодирования, позволяющий улучшить (т.е. сделать более точной) оценку вероятности информационных, а также кодированных бит, является приемлемым. Декодер 409 осуществляет оценку апостериорной вероятности информационных бит, а также уточненную оценку апостериорной вероятности кодированных бит. Следует отметить, что оценка апостериорной информации кодированных бит после декодера может быть произведена альтернативным методом, а именно за счет рекодирования, когда оценка вероятности кодированных бит производится на основании оценки апостериорной вероятности информационных бит после декодера. Ремодулятор 415 производит оценку принятого символьного вектора и дисперсии его компонентов на основе апостериорной вероятности кодированных бит, полученной в канальном декодере, в соответствии с выражениями (12), (13). Модуль коррекции априорных данных производит сравнение дисперсии (13), полученной в ремодуляторе, с дисперсией (7), полученной в MIMO детекторе на предыдущем цикле итеративного процесса, и, в случае, если новая оценка дисперсии оказывается меньше предыдущей, корректирует оценку априорных данных на вход MIMO детектора в соответствии с выражениями (14). Заметам, что аналогичная коррекция используется в патенте [6], однако в нашем случае апостериорная вероятность берется с выхода канального декодера, в то время как в [6] - непосредственно с выхода MIMO детектора.

На Фиг.5 и Фиг.6 представлены результаты сравнения работы заявляемого способа итеративного детектирования и декодирования с традиционным способом последовательного детектирования, декодирования, в котором в качестве MIMO детектора используется ML, а также MMSE детекторы. Результаты получены для модели стандарта LTE, конфигурации MIMO 4×4, модуляции 16QAM и модели канала ЕРА-5. Заметим, что при совместном детектировании и декодировании использованы только два внешних цикла итеративного детектирования-декодирования, при этом канальный декодер производил декодирование с использованием своих внутренних итеративных циклов (декодировался турбо-код) таким образом, что после первого цикла работы MIMO детектора декодер производил два внутренних цикла декодирования, после чего производился еще один цикл детектирования с использованием априорной информации, полученной из декодера, после второго цикла детектирования декодер производил еще четыре цикла внутреннего декодирования. Таким образом, общее количество внутренних циклов работы канального декодера равнялось шести. Точно такое же количество циклов работы декодера использовалось в традиционных схемах последовательного детектирования, декодирования, т.е. прирост сложности и задержки происходил лишь за счет одного дополнительного цикла работы MIMO детектора, что можно считать незначительным ввиду относительной простоты базирующегося на MMSE фильтрации метода MIMO детектирования. Как видно из графиков, предложенный метод значительно превосходит традиционный метод с использованием MMSE детектора, а в ряде случаев даже превосходит традиционный метод, использующий ML детектор. Это позволяет сделать вывод, что предложенный способ, обладая невысокой сложностью реализации, приближается по своим характеристикам к традиционным методам детектирования, декодирования, использующим ML детектор, и в ряде случаев даже превосходит его, что делает предложенный способ весьма перспективным кандидатом для реализации в существующих, а также будущих системах связи.

Ссылки

1. P.W.Wolniansky, G.J.Foschini, G.D.Golden, R.A.Valenzuela "V-BLAST: An Architecture for Realizing Very High Data Rates Over the Rich-Scaterring Wireless Channel," Proc. ISSE, Pisa, Italy, September 1998.

2. Dirk Wubben, Ronald Bonke, Volker Kuhn, Karl-Dirk Kammeyer, "MMSE Extension of V-BLAST based on Sorted QR Decomposition, Vehicular Technology Conference, 2003, VTC 2003-Fall, 2003, IEEE 58^th, V.1, Pages 508-512.

3. "Iterative Detecting and Decoding for a M1MO-OFDM System", Patent Application Publication No. US 2005/0157811 Al Jul. 21, 2005

4. "Signal Detection Using Sphere Decoding Technique", Patent No. US 7,443,928 B2.

5. "Iterative Detecting and Decoding Apparatus and Method in MIMO System", Patent No. US 8,027,401.

6. "Способ детектирования сигнала в системах связи с MIMO каналом", патент Российской Федерации № 2444846.

Изобретение относится к способам приема сигнала в системе передачи, использующей технологию на основе MIMO (Multiple-In Multiple-Out).

Технический результат заключается в улучшении эффективности и качества передачи данных. Способ итеративного детектирования и декодирования сигнала в системах связи с MIMO каналом содержит два этапа: на первом этапе формируют MMSE (Minimum Mean Square Error) оценку переданного символа с помощью линейного MMSE фильтра, на втором этапе путем линейного преобразование каждого n-ого (где n=1, 2, (N) компонента $${\tilde{x}}_{MMSE,n}$$ вектора MMSE оценок $${\tilde{X}}_{MMSE}$$ и соответствующего ему априорного математического ожидания осуществляют коррекцию и находят новую линейную оценку $${\tilde{x}}_{LMS,n}$$ переданного QAM символа и дисперсию ошибки для этой оценки $${\sigma }_{{}_{LMS,n}}^2$$ .

4 з.п. ф-лы, 6 ил.

1. Способ итеративного детектирования и декодирования сигнала в системах связи с MIMO (Multiple-In Multiple-Out) каналом, включающий в себя выполнение следующих операций:
- принимают сигналы, модулированные переданными QAM (Quadrature Amplitude Modulation) символами, через несколько приемных антенн;
- осуществляют частотное преобразование принятых сигналов на нулевую несущую частоту;
- формируют отсчеты принятых и преобразованных сигналов путем квантования и дискретизации, которые рассматриваются как выходные отсчеты MIMO канала и описываются матричным уравнением:
Y=HX+η
где Y - M-мерный вектор выходных отсчетов MIMO канала, Х - N-мерный вектор входных отсчетов MIMO канала, η - M-мерный вектор отсчетов шума, H - канальная матрица MIMO канала размера М×N, при этом полагают, что частота дискретизации равна частоте следования QAM символов, поэтому вектор входных отсчетов MIMO канала Х является также вектором переданных QAM символов,
- по полученным выходным отсчетам MIMO канала оценивают канальную матрицу MIMO канала H, дисперсию шума $${\sigma }_{\eta }^2$$ и среднюю мощность сигнала $${\sigma }_s^2$$ ;
- детектируют QAM символы на основе выходных отсчетов MIMO канала, производя оценку апостериорной вероятности кодированных бит;
- декодируют информационные биты на основе оценки апостериорной вероятности кодированных бит, полученной в детекторе QAM символов, и в соответствии с параметрами канального кода производят оценку апостериорной вероятности информационных бит и уточняют оценки апостериорной вероятности кодированных бит;
- на основе информации о вероятности кодированных бит на выходе канального декодера формируют оценки априорной вероятности о принятых QAM символах и направляют их на дополнительный вход MIMO детектора в качестве априорной информации для формирования новых оценок апостериорных вероятностей кодированных и информационных бит на следующем цикле итеративного процесса;
- после завершения всех итераций на основе оценки апостериорной вероятности информационных бит на основном выходе канального декодера формируют окончательную последовательность жестких оценок информационных бит, восстанавливая оригинальные данные;
отличающийся тем, что
- детектирование QAM символов осуществляют в два этапа, на первом этапе формируют MMSE (Minimum Mean Square Error) оценку переданного символа с помощью линейного MMSE фильтра, описываемого уравнениями:
$${\tilde{X}}_{MMSE}={\overline{X}}_{pr}+G_{MMSE}\left(Y-H{\overline{X}}_{pr}\right)$$
$$G_{MMSE}=V_{pr}H\text{'}{\left(H{V}_{pr}H\text{'}+\frac{{\sigma }_{\eta }^2}{{\sigma }_s^2}I\right)}^{-1}$$
где $${\overline{X}}_{pr}$$ - N-мерный вектор априорных математических ожиданий переданных QAM символов, V_pr - диагональная матрица, характеризующая априорную дисперсию переданных QAM символов, нормированную средней мощностью сигнала $${\sigma }_s^2$$ ,
для каждого компонента $${\tilde{x}}_{MMSE,n}$$ (n=1, 2,…N) вектора MMSE оценок $${\tilde{X}}_{MMSE}$$ , вычисляют дисперсию ошибки оценивания $${\sigma }_{{}_{MMSE,n}}^2$$ , которая совпадает с соответствующим диагональным элементом корреляционной матрицы V_MMSE, вычисляемой по формуле
V_MMSE=V_pr-G_MMSEHV_pr,
на втором этапе путем линейного преобразование каждого n-ого (где n=1, 2,…N) компонента $${\tilde{x}}_{MMSE,n}$$ вектора MMSE оценок $${\tilde{X}}_{MMSE}$$ и соответствующего ему априорного математического ожидания $${\overline{x}}_{pr,n}$$ осуществляют коррекцию и находят новую линейную оценку $${\tilde{x}}_{LMS,n}$$ переданного QAM символа и дисперсию ошибки для этой оценки $${\tilde{\sigma }}_{{}_{LMS,n}}^2$$
$${\tilde{x}}_{LMS,n}=-\frac{{\sigma }_{{}_{MMSE,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}{\overline{x}}_{pr,n}+\frac{{\sigma }_{{}_{pr,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}{\tilde{x}}_{MMSE,n}$$
$${\tilde{\sigma }}_{{}_{LMS,n}}^2=\frac{{\sigma }_{{}_{pr,n}}^2{\sigma }_{{}_{MMSE,n}}^2}{{\sigma }_{{}_{pr,n}}^2-{\sigma }_{{}_{MMSE,n}}^2}$$
где $${\sigma }_{{}_{pr,n}}^2$$ - априорная дисперсия n-ого переданного QAM символа,
- осуществляют демодуляцию оценок QAM символов, полученных в детекторе, и вычисляют апостериорные вероятности переданных бит в соответствии с выражением:
$${\lambda }_{n,k}=\ln \left(\frac{{\displaystyle \sum _{x:f\left(b_{n,k}=1\right)}\exp }\left(-\frac{1}{{\tilde{\sigma }}_{LMS,n}^2}{\left|{\tilde{x}}_{LMS,n}-x\left(b_{n\mathrm{,1}},b_{n\mathrm{,2}},\dots ,b_{n,K}\right)\right|}^2\right){\displaystyle \prod _{t=1}^K\mathrm{Pr}\left(b_{n,t}\right)}}{{\displaystyle \sum _{x:f\left(b_{n,k}=-1\right)}\exp }\left(-\frac{1}{{\tilde{\sigma }}_{LMS,n}^2}{\left|{\tilde{x}}_{LMS,n}-x\left(b_{n\mathrm{,1}},b_{n\mathrm{,2}},\dots ,b_{n,K}\right)\right|}^2\right){\displaystyle \prod _{t=1}^K\mathrm{Pr}\left(b_{n,t}\right)}}\right)$$
где x(b₁,…b_K) - табличная функция с числом состояний 2^K, описывающая сигнальное созвездие переданного QAM символа в зависимости от комбинации кодированных бит b_n,k∈{-1;1}, k=1, 2,…K, n=1, 2,…N, Pr(b_n,k) - априорная вероятность k-го кодированного бита в n-м QAM символе, полученная на предыдущей итерации, λ(b_n,k) - логарифмическое отношение апостериорных вероятностей для k-го бита в n-м QAM символе, K - число бит в одном QAM символе;
- удаляют из вычисленной в демодуляторе апостериорной вероятности кодированных бит априорную информацию для каждого кодированного бита P_pr(b_n,k) путем вычитания логарифмического отношения априорных вероятностей $${\stackrel{⌢}{\lambda }}_{n,k}={\lambda }_{n,k}-{\lambda }_{n,k,pr}$$ ,
где $${\lambda }_{n,k,pr}=\ln \left(\frac{\mathrm{Pr}\left(b_{n,k}\right)=+1}{\mathrm{Pr}\left(b_{n,k}\right)=-1}\right)$$ ;
- на основании этой информации осуществляют декодирование информационных бит в соответствии с параметрами канального кода;
- выделяют из уточненной апостериорной вероятности кодированных бит, полученной после декодирования в виде уточненного логарифмического отношения вероятностей $${\tilde{\lambda }}_{n,k}$$ , внешнюю информацию, удаляя из нее информацию, поступающую на вход декодера в виде мягких оценок кодированных бит $${\stackrel{⌢}{\lambda }}_{n,k}$$ , и на основании этой информации формируют оценку априорной вероятности, поступающую в QAM демодулятор на следующей итерации
$$\mathrm{Pr}\left(b_{n,k}\right)=\frac{e^{b_{n,k}\frac{{\lambda }_{n,k,pr}}{2}}}{e^{\frac{{\lambda }_{n,k,pr}}{2}}+e^{-\frac{{\lambda }_{n,k,pr}}{2}}}$$ , $${\lambda }_{n,k,pr}={\tilde{\lambda }}_{n,k}-{\stackrel{⌢}{\lambda }}_{n,k}$$
- на основе уточненных оценок апостериорной вероятности кодированных бит на выходе канального декодера формируют оценки переданных QAM символов путем весового сложения значений сигнального QAM созвездия, где веса определяются уточненными апостериорными вероятностями кодированных бит в соответствии с выражением:
где
и вычисляют дисперсии этих оценок
,
- сравнивают полученные дисперсии $${\sigma }_n^2$$ с оценками дисперсий после детектирования QAM символов $${\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , полученными на предыдущей итерации и, если новая оценка дисперсии оказывается меньше дисперсии детектирования, т.е. $${\sigma }_n^2<{\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , то осуществляют коррекцию априорного математического ожидания и априорной дисперсии
$${\tilde{x}}_{pr,n}=-\frac{{\sigma }_n^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}{\tilde{x}}_{LMS,n}+\frac{{\tilde{\sigma }}_{{}_{LMS,n}}^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}{\widehat{x}}_n$$
$${\tilde{\sigma }}_{{}_{pr,n}}^2=\frac{{\sigma }_n^2{\tilde{\sigma }}_{{}_{LMS,n}}^2}{{\tilde{\sigma }}_{{}_{LMS,n}}^2-{\sigma }_n^2}$$
если дисперсия $${\sigma }_n^2$$ больше, либо равна дисперсии $${\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , т.е. $${\sigma }_n^2\ge {\tilde{\sigma }}_{{}_{LMS,n}}^2$$ , то априорные математические ожидания и дисперсию оставляют без изменений
$$\begin{array}{c}{\tilde{x}}_{pr,n}={\overline{x}}_{pr,n}\\ {\tilde{\sigma }}_{{}_{pr,n}}^2={\sigma }_{{}_{pr,n}}^2\end{array}$$
- априорные математические ожидания $${\tilde{x}}_{pr,n}$$ и дисперсии $${\tilde{\sigma }}_{{}_{pr,n}}^2$$ , для всех n=1, 2,…N, объединенные в вектор новых априорных математических ожиданий $${\tilde{X}}_{pr}$$ и в новую априорную диагональную корреляционную матрицу $${\tilde{V}}_{pr}$$ , являются входными параметрами MIMO детектора, используемыми на следующей итерации.

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

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

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

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