CN102811117B - Mimo系统的译码方法及装置 - Google Patents

Abstract

本发明涉及MIMO系统的译码方法及装置，上述方法根据系统收到的调制符号r，计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影对的实部和虚部进行量化估计，得到的实部和虚部的估计值s_2,I、s_2,Q，对s_2,I、s_2,Q进行星座映射，得到s₂在星座图中的位置，根据s₂的位置，计算s₁的估计值与实际值之间的误差D，根据D，更新s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，根据min_dist_0[k]和min_dist_1[k]，计算s₁的每个比特的对数似然比，得到译码结果；上述装置包括投影模块、估计模块、映射模块、误差计算模块、最小误差更新模块以及对数似然比计算模块。本发明译码复杂度低。

Description

MIMO系统的译码方法及装置

技术领域

本发明涉及通信技术领域，尤其涉及一种MIMO(Multiple-Input Multiple-Out-put，多输入多输出)系统的译码方法及装置。

背景技术

多天线传输技术对提高无线通信系统的容量和可靠性具有至关重要的作用，分布式接入方式也是目前高速分组传输研究的热点。传输信号由多个天线同时发送和接收，发送端和接收端之间的无线信道由传统的SISO(Single-InputSingle-Out-put，单输入单输出)系统变成MIMO系统，并使通信系统具有了除传统时间、频率、码道资源以外的空间资源。理论研究表明，MIMO信道是多个SISO信道的叠加，其容量和min(N_T,N_R)成正比，其中N_T和N_R分别是发送和接收天线的数目。这说明当增加发送和接收天线的数目时，可有效地提高信道容量。所以，MIMO系统为提高无线网络的信息吞吐量、扩大覆盖区域和提高传输质量提供了一条具有巨大的潜力的道路。MIMO技术能够在空间中产生独立的并行信道同时传输多路数据流，这样就有效地提高了系统的传输速率，即在不增加系统带宽的情况下成倍地提高通信系统的容量和频谱利用率。MIMO-OFDM(Orthogonal Frequency Division Multiplexing，正交频分复用)系统结合OFDM技术和MIMO技术的优势，在提高无线链路的传输速率和可靠性方面具有巨大潜力。

MIMO-OFDM系统中检测和译码技术也是现在研究的热点。由于采用多天线同时发射，存在着共信道干扰，接收机的检测技术的优劣以及复杂度直接影响系统的性能和应用前景。整个接收部分联合(ML/MAP)接收机能取得最优性能，但复杂度太高，目前的硬件处理能力尚不能满足运算要求，往往只能在天线数目和调制阶数都很小时才能应用。线性接收方法(ZF、MMSE)的复杂度低但性能差，性能介于ML和线性接收之间的有干扰抵消算法和球形译码算法。干扰抵消算法需要从接收信号中减去先检测出的数据部分，因此存在错误传播的现象，性能受干扰抵消次序影响。球形译码算法是极大似然算法的简化，通过动态改变搜索的圆心与半径减少搜索的次数，在高信噪比时性能逼近ML的性能，复杂度较低，但在低信噪比时搜索的时间比较长复杂度仍然很高。如何减少最优检测算法的复杂度并避免传统检测算法应用的局限性对系统的实现很有意义。

发明内容

本发明的目的是，提供一种MIMO系统的译码方法及装置，以优化现有的译码方法复杂度太高的问题。

本发明提供了一种MIMO系统的译码方法，

根据系统接收到的调制符号r，计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

对上述的实部和虚部进行量化估计，得到上述的实部和虚部的估计值s_2,I、s_2,Q；

对上述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

根据上述估计得到的调制符号s₂在星座图中的位置，计算上述调制符号s₁的估计值与实际值之间的误差D；

根据上述误差D，更新上述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；

根据上述最小误差min_dist_0[k]和min_dist_1[k]，计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

优选地，上述计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影步骤具体包括：

计算上述另一发射终端的信道响应h₂的模方||h₂||²；

计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H；

计算上述h₂ ^H与上述系统接收到的调制符号r的乘积；

计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积；

根据上述另一发射终端的信道响应h₂的模方||h₂||²，上述h₂ ^H与上述系统接收到的调制符号r的乘积以及上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积，计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

优选地，通过如下公式计算上述当前发射终端上发射的调制符号s₁在另一发射终端上发射的调制符号s₂的投影

其中，h₂表示上述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置，||h₂||²为h₂的模方。

优选地，上述计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影步骤具体包括：

计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H；

计算上述h₂ ^H与上述系统接收到的调制符号r的乘积；

计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积；

根据上述h₂ ^H与上述系统接收到的调制符号r的乘积以及上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积，计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

优选地，通过如下公式计算上述当前发射终端上发射的调制符号s₁在另一发射终端上发射的调制符号s₂的投影

其中，h₂表示上述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置。

优选地，通过如下公式对上述的实部和虚部进行估计：

其中，表示的实部，表示的虚部，表示向下取整。

优选地，上述对的实部和虚部进行量化估计，得到实部和虚部的估计值s_2,I、s_2,Q步骤具体包括以下步骤：

计算上述另一发射终端的信道响应h₂的模方||h₂||²；

根据如下公式计算上述的实部的中间估计值

根据如下公式计算上述的虚部的中间估计值

判断上述是否小于-2||h₂||²，若是，则s_2,I＝0；否则，

判断上述是否小于0，若是，则s_2,I＝1；否则，

判断上述是否小于2||h₂||²，若是，则s_2,I＝2；否则，s_2,I＝3；

判断上述是否小于-2||h₂||²，若是，则s_2,Q＝0；否则，

判断上述是否小于0，若是，则s_2,Q＝1；否则，

判断上述是否小于2||h₂||²，若是，则s_2,Q＝2；否则，s_2,Q＝3。

优选地，通过如下公式计算上述调制符号s₁的估计值与实际值之间的误差D：

D＝||r-s₁h₁-s₂h₂||²

其中，s₁表示上述当前发射终端发射的调制符号，h₁表示上述当前发射终端的信道响应，s₂表示上述另一发射终端发射的调制符号，h₂表示上述另一发射终端的信道响应。

优选地，上述方法还包括：

上述系统初始化时，将上述最小误差min_dist_0[k]和min_dist_1[k]中各元素的值初始化为正无穷大。

优选地，上述根据上述误差D，更新上述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]步骤具体包括以下步骤：

判断上述调制符号s₁的第k个比特的值，

若

则判断上述误差D是否小于min_dist_0[k]，若是，则令min_dist_0[k]＝D；否则，不更新min_dist_0[k]；

若

则判断上述误差D是否小于min_dist_1[k]，若是，则令min_dist_1[k]＝D；否则，不更新min_dist_1[k]。

优选地，通过如下公式计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比：

其中，LLR₁ ^k表示上述当前发射终端发射的调制符号s₁的第k个比特的对数似然比。

本发明还提供了一种MIMO系统的译码装置，包括投影模块、估计模块、映射模块、误差计算模块、最小误差更新模块以及对数似然比计算模块，上述投影模块，用于根据系统接收到的调制符号r，计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

上述估计模块，用于对上述的实部和虚部进行量化估计，得到上述的实部和虚部的估计值s_2,I、s_2,Q；

上述映射模块，用于对上述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

上述误差计算模块，用于根据上述估计得到的调制符号s₂在星座图中的位置，计算上述调制符号s₁的估计值与实际值之间的误差D；

上述最小误差更新模块，用于根据上述误差D，更新上述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；

上述对数似然比计算模块，用于根据上述最小误差min_dist_0[k]和min_dist_1[k]，计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

优选地，上述装置还包括初始化模块，用于在上述MIMO系统初始化时，将上述最小误差min_dist_0[k]和min_dist_1[k]中各元素的值初始化为正无穷大。

优选地，上述投影模块，还用于计算上述另一发射终端的信道响应h₂的模方||h₂||²，计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H，计算上述h₂ ^H与上述系统接收到的调制符号r的乘积以及计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积；

上述估计模块，还用于计算上述另一发射终端的信道响应h₂的模方||h₂||²；计算上述的实部的中间估计值计算上述的虚部的中间估计值以及判断上述是否小于-2||h₂||²，并在上述小于-2||h₂||²时，令上述的估计值s_2,I等于0，在上述大于等于-2||h₂||²时，判断上述是否小于0，并在上述小于0时，令上述的估计值s_2，I等于1，在上述大于等于0时，判断上述是否小于2||h₂||²，在上述小于2||h₂||²时，令上述的估计值s_2,I等于2，在上述大于等于2||h₂||²时，令上述的估计值s_2,I等于3；判断上述是否小于-2||h₂||²，并在上述小于-2||h₂||²时，令上述的估计值s_2，Q等于0，在上述大于等于-2||h₂||²时，判断上述是否小于0，并在上述小于0时，令上述的估计值s_2,Q等于1，在上述大于等于0时，判断上述是否小于2||h₂||²，在上述小于2||h₂||²时，令上述的估计值s_2,Q等于2，在上述大于等于2||h₂||²时，令上述的估计值s_2,Q等于3；

上述最小误差更新模块，还用于判断上述调制符号s₁的第k个比特的值，并在上述时，判断上述误差D是否小于min_dist_0[k]，在上述误差D小于min_dist_0[k]时，令min_dist_0[k]的值等于上述误差D，在上述误差D大于等于min_dist_0[k]时，不更新min_dist_0[k]；或者在上述判断上述误差D是否小于min_dist_1[k]，在上述误差D小于min_dist_1[k]时，令min_dist_1[k]的值等于上述误差D，在上述误差D大于等于min_dist_1[k]时，不更新min_dist_1[k]。

本发明优化了传统MAP算法复杂度太高的问题，本发明计算出的译码结果相当于软解调的软信息，可直接用于译码，无需解调，对于多天线、高阶和低阶调制都有很好的译码效果。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是本发明MIMO系统的译码方法优选实施例的流程图；

图2是16QAM星座示意图；

图3是估计得到的调制符号s₂的星座示意图；

图4是本发明MIMO系统的译码装置优选实施例的原理框图。

具体实施方式

为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚、明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

如图1所示，是本发明MIMO系统的译码方法优选实施例的流程图；本实施例以16QAM为例，其星座示意图如图2所示，其中，b₁b₀,b₃b₂表示16QAM星座图中每个调制符号对应的4个bit，c为归一化因子，16个调制符号分别标记为C₁,C₂,…,C₁₆，因对于本发明而言，调制符号标记和调制符号如何对应并不重要，故图中未示出16个调制符号的具体位置。每个调制符号有4个比特。假设MIMO系统有N_RX根接收天线，系统模型如下：

其中，表示MIMO系统的各接收天线接收到的调制符号；s₁表示MIMO系统的一个发射终端发射的调制符号，表示该发射终端在接收天线N_RX的信道响应；s₂表示MIMO系统的另一个发射终端发射的调制符号，表示该发射终端在接收天线N_RX的信道响应；表示MIMO系统的各接收天线上的噪声。

令：

则上述系统模型变为：

r＝s₁h₁+s₂h₂+n

其中，r表示MIMO系统接收到的调制符号；s₁表示MIMO系统的一个发射终端发射的调制符号，s₁∈{C₁,C₂,…,C₁₆}，h₁表示该发射终端的信道响应；s₂表示MIMO系统的另一个发射终端发射的调制符号，s₂∈{C₁,C₂,…,C₁₆}，h₂表示该发射终端的信道响应；n表示噪声。

基于以上模型，本实施例包括以下步骤：

步骤S001：初始化当前发射终端发射的调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，将其中各元素的值设置为正无穷大，即令：

min_dist_0[k]＝+∝

min_dist_1[k]＝+∝

其中，k表示调制符号的第k个比特位，在16QAM中，k＝0,1,2,3。

步骤S002：根据系统接收到的调制符号r，计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

在实施例一中，本步骤具体包括：

步骤S10：计算上述另一发射终端的信道响应h₂的模方||h₂||²并保存；

步骤S11：计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H并保存；

步骤S12：计算上述h₂ ^H与上述系统接收到的调制符号r的乘积A并保存，即

A＝h₂ ^Hr

步骤S13：计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积B并保存；即

B＝h₂ ^Hh₁

步骤S14：根据上述另一发射终端的信道响应h₂的模方||h₂||²，上述h₂ ^H，以及上述A和B，通过如下公式计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

其中，h₂表示上述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置，||h₂||²为h₂的模方。

实施例一中，对于不同的调制符号s₁，每次只要作一次复数乘法、一次复数减法和一次复数除实数的除法即可得到相应的投影

由于实施例一中，需要做一次除法运算，开销比较大，还可以对实施例进行进一步的优化。

在实施例二中，本步骤具体包括：

步骤S10’：计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H并保存；

步骤S11’：计算上述h₂ ^H与上述系统接收到的调制符号r的乘积A并保存；

步骤S12’：计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积B并保存；

步骤S13’：根据上述h₂ ^H与A、B，通过如下公式计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

其中，h₂表示上述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置。

实施例二与实施例一相比，对于不同的调制符号s₁，仅需要作一次复数乘法、一次复数减法和一次复数与实数的乘法即可得到相应的投影

步骤S003：对上述的实部和虚部进行量化估计，得到实部和虚部的估计值s_2,I、s_2,Q；

本步骤中，

若步骤S002中采用实施例一计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影则本步骤通过如下步骤对上述的实部和虚部进行估计：

步骤S20：根据如下公式计算上述的实部的中间估计值

步骤S21：根据如下公式计算上述的虚部的中间估计值

步骤S22：根据上述通过如下公式计算上述的实部的估计值s_2I为和虚部的估计值s_2,Q

其中，表示的实部，表示的虚部，表示向下取整。

若步骤S002中采用实施例二计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影则本步骤通过如下步骤对上述的实部和虚部进行估计：

步骤S20’：计算上述另一发射终端的信道响应h₂的模方||h₂||²并保存；

步骤S21’：根据如下公式计算上述的实部的中间估计值

步骤S22’：根据如下公式计算上述的虚部的中间估计值

步骤S23’：判断上述是否小于-2||h₂||²，若是，则执行步骤S24’，否则，执行步骤S25’；

步骤S24’：s_2,I＝0，执行步骤S30’；

步骤S25’：判断上述是否小于0，若是，则执行步骤S26’，否则，执行步骤S27’；

步骤S26’：s_2,I＝1，执行步骤S30’；

步骤S27’：判断上述是否小于2||h₂||²，若是，则执行步骤S28’，否则，执行步骤S29’；

步骤S28’：s_2,I＝2，执行步骤S30’；

步骤S29’：s_2,I＝3，执行步骤S30’；

步骤S30’：判断上述是否小于-2||h₂||²，若是，则执行步骤S31’，否则，执行步骤S32’；

步骤S31’：s_2,Q＝0，估值结束；

步骤S32’：分别判断上述是否小于0，若是，则执行步骤S33’，否则，执行步骤S34’；

步骤S33’：s_2,Q＝1，估值结束；

步骤S34’：判断上述是否小于2||h₂||²，若是，则执行步骤S35’，否则，执行步骤S36’；

步骤S35’：s_2,Q＝2，估值结束；

步骤S36’：s_2,Q＝3，估值结束。

步骤S004：对上述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

如图3所示，是估计得到的调制符号s₂的星座示意图，I路即横轴通过s_2,I映射，s_I表示调制符号的实部，Q路即纵轴通过s_2,Q映射，s_Q表示调制符号s的虚部。比如，当s_2,I＝1，s_2,Q＝3时，则估计得到的调制符号s₂被映射到图中的第二象限中左上角的点O。

步骤S005：根据上述估计得到的调制符号s₂在星座图中的位置，计算上述调制符号s₁的估计值与实际值之间的误差D；

本步骤通过如下公式计算上述调制符号s₁的估计值与实际值之间的误差D：

D＝||r-s₁h₁-s₂h₂||²

其中，s₁表示上述当前发射终端发射的调制符号，h₁表示上述当前发射终端的信道响应，s₂表示上述另一发射终端发射的调制符号，h₂表示上述另一发射终端的信道响应。

步骤S006：根据上述距离D，更新上述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；

本步骤具体为：

步骤S0061：判断上述调制符号s₁的第k个比特的值，若则执行步骤S0062；若执行步骤S0065；

步骤S0062：判断上述误差D是否小于min_dist_0[k]，若是，执行步骤S0063，否则，执行步骤S0064；

步骤S0063：令min_dist_0[k]＝D；

步骤S0064：不更新min_dist_0[k]；

步骤S0065：判断所述误差D是否小于min_dist_1[k]，若是，则执行步骤S0066，否则，执行步骤S0067；

步骤S0066：令min_dist_1[k]＝D；

步骤S0067：不更新min_dist_1[k]。

步骤S007：根据上述最小误差min_dist_0[k]和min_dist_1[k]，计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

本步骤通过如下公式计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比：

其中，LLR₁ ^k表示上述当前发射终端发射的调制符号s₁的第k个比特的对数似然比。

如图4所示，是本发明MIMO系统的译码装置的优选实施例原理框图，本实施中，译码装置包括投影模块02、估计模块03、映射模块04、误差计算模块05、最小误差更新模块06以及对数似然比计算模块07，其中，

初始化模块01，用于在上述MIMO系统初始化时，将上述最小误差min_dist_0[k]和min_dist_1[k]中各元素的值初始化为正无穷大；

投影模块02，用于根据系统接收到的调制符号r，计算计算上述另一发射终端的信道响应h₂的模方||h₂||²并保存，计算上述另一发射终端的信道响应h₂的共轭转置h₂ ^H并保存，计算上述h₂ ^H与上述系统接收到的调制符号r的乘积并保存以及计算上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积并保存，以及根据上述另一发射终端的信道响应h₂的模方||h₂||²，上述h₂ ^H与上述系统接收到的调制符号r的乘积以及上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积，计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影或者根据上述h₂ ^H与上述系统接收到的调制符号r的乘积以及上述h₂ ^H与上述当前发射终端的信道响应h₁的乘积，计算上述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

估计模块03，用于对上述的实部和虚部进行量化估计，得到上述的实部和虚部的估计值s_2，I、s_2，Q；具体为：用于计算上述另一发射终端的信道响应h₂的模方||h₂||²并保存；计算上述的实部的中间估计值计算上述的虚部的中间估计值以及判断上述是否小于-2||h₂||²，并在上述小于-2||h₂||²时，令上述的估计值s_2，I等于0，在上述大于等于-2||h₂||²时，判断上述是否小于0，并在上述小于0时，令上述的估计值s_2,I等于1，在上述大于等于0时，判断上述是否小于2||h₂||²，在上述小于2||h₂||²时，令上述的估计值s_2,I等于2，在上述大于等于2||h₂||²时，令上述的估计值s_2,I等于3；判断上述是否小于-2||h₂||²，并在上述小于-2||h₂||²时，令上述的估计值s_2,Q等于0，在上述大于等于-2||h₂||²时，判断上述是否小于0，并在上述小于0时，令上述的估计值s_2,Q等于1，在上述大于等于0时，判断上述是否小于2||h₂||²，在上述小于2||h₂||²时，令上述的估计值s_2,Q等于2，在上述大于等于2||h₂||²时，令上述的估计值s_2,Q等于3；

映射模块04，用于对上述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

误差计算模块05，用于根据上述估计得到的调制符号s₂在星座图中的位置，计算上述调制符号s₁的估计值与实际值之间的误差D；

最小误差更新模块06，用于更新上述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；具体为：判断上述调制符号s₁的第k个比特的值，并在上述时，判断上述误差D是否小于min_dist_0[k]，在上述误差D小于min_dist_0[k]时，令min_dist_0[k]的值等于上述误差D，在上述误差D大于等于min_dist_0[k]时，不更新min_dist_0[k]；或者在上述判断上述误差D是否小于min_dist_1[k]，在上述误差D小于min_dist_1[k]时，令min_dist_1[k]的值等于上述误差D，在上述误差D大于等于min_dist_1[k]时，不更新min_dist_1[k]。

对数似然比计算模块07，用于根据上述最小误差min_dist_0[k]和min_dist_1[k]，计算上述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

上述说明示出并描述了本发明的优选实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (14)

1.一种多输入多输出MIMO系统的译码方法，其特征在于，

根据系统接收到的调制符号r，计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

对所述的实部和虚部进行量化估计，得到所述的实部和虚部的估计值s_2,I、s_2,Q；

对所述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

根据所述估计得到的调制符号s₂在星座图中的位置，计算所述调制符号s₁的估计值与实际值之间的误差D；

根据所述误差D，更新所述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；

根据所述最小误差min_dist_0[k]和min_dist_1[k]，计算所述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

2.根据权利要求1所述的方法，其特征在于，所述计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影步骤具体包括：

计算所述另一发射终端的信道响应h₂的模方||h₂||²；

计算所述另一发射终端的信道响应h₂的共轭转置h₂ ^H；

计算所述h₂ ^H与所述系统接收到的调制符号r的乘积；

计算所述h₂ ^H与所述当前发射终端的信道响应h₁的乘积；

根据所述另一发射终端的信道响应h₂的模方||h₂||²，所述h₂ ^H与所述系统接收到的调制符号r的乘积以及所述h₂ ^H与所述当前发射终端的信道响应h₁的乘积，计算所述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

3.根据权利要求1或2所述的方法，其特征在于，通过如下公式计算所述当前发射终端上发射的调制符号s₁在另一发射终端上发射的调制符号s₂的投影

${\tilde{s}}_2={h_2}^H(r-s_1{h}_1)/\left|\right|h_2|{|}^2$

其中，h₂表示所述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置，||h₂||²为h₂的模方。

4.根据权利要求1所述的方法，其特征在于，所述计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影步骤具体包括：

计算所述另一发射终端的信道响应h₂的共轭转置h₂ ^H；

计算所述h₂ ^H与所述系统接收到的调制符号r的乘积；

计算所述h₂ ^H与所述当前发射终端的信道响应h₁的乘积；

根据所述h₂ ^H与所述系统接收到的调制符号r的乘积以及所述h₂ ^H与所述当前发射终端的信道响应h₁的乘积，计算所述当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

5.根据权利要求1或4所述的方法，其特征在于，通过如下公式计算所述当前发射终端上发射的调制符号s₁在另一发射终端上发射的调制符号s₂的投影

${\tilde{s}}_2=h_2^H(r-s_1{h}_1)\sqrt{10}$

其中，h₂表示所述另一发射终端的信道响应，h₂ ^H表示h₂的共轭转置。

6.根据权利要求3所述的方法，其特征在于，通过如下公式对所述的实部和虚部进行估计：

其中，表示的实部，表示的虚部，表示向下取整。

7.根据权利要求5所述的方法，其特征在于，所述对的实部和虚部进行量化估计，得到实部和虚部的估计值s_2,I、s_2,Q步骤具体包括以下步骤：

计算所述另一发射终端的信道响应h₂的模方||h₂||²；

根据如下公式计算所述的实部的中间估计值

根据如下公式计算所述的虚部的中间估计值

判断所述是否小于-2||h₂||²，若是，则s_2,I＝0；否则，

判断所述是否小于0，若是，则s_2,I＝1；否则，

判断所述是否小于2||h₂||²，若是，则s_2,I＝2；否则，s_2,I＝3；

判断所述是否小于-2||h₂||²，若是，则s_2,Q＝0；否则，

判断所述是否小于0，若是，则s_2,Q＝1；否则，

判断所述是否小于2||h₂||²，若是，则s_2,Q＝2；否则，s_2,Q＝3。

8.根据权利要求1所述的方法，其特征在于，通过如下公式计算所述调制符号s₁的估计值与实际值之间的误差D：

D＝||r-s₁h₁-s₂h₂||²

其中，s₁表示所述当前发射终端发射的调制符号，h₁表示所述当前发射终端的信道响应，s₂表示所述另一发射终端发射的调制符号，h₂表示所述另一发射终端的信道响应。

9.根据权利要求1所述的方法，其特征在于，所述方法还包括：

所述系统初始化时，将所述最小误差min_dist_0[k]和min_dist_1[k]中各元素的值初始化为正无穷大。

10.根据权利要求1或9所述的方法，其特征在于，所述根据所述误差D，更新所述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]步骤具体包括以下步骤：

判断所述调制符号s₁的第k个比特的值，

若

则判断所述误差D是否小于min_dist_0[k]，若是，则令min_dist_0[k]＝D；否则，不更新min_dist_0[k]；

若

则判断所述误差D是否小于min_dist_1[k]，若是，则令min_dist_1[k]＝D；否则，不更新min_dist_1[k]。

11.根据权利要求1所述的方法，其特征在于，通过如下公式计算所述当前发射终端发射的调制符号s₁的每个比特的对数似然比：

${\mathrm{LLR}}_1^k=\min \_dist\_1\[k\]-\min \_dist\_0\[k\]$

其中，LLR₁ ^k表示所述当前发射终端发射的调制符号s₁的第k个比特的对数似然比。

12.一种MIMO系统的译码装置，其特征在于，包括投影模块、估计模块、映射模块、误差计算模块、最小误差更新模块以及对数似然比计算模块，所述投影模块，用于根据系统接收到的调制符号r，计算当前发射终端发射的调制符号s₁在另一发射终端发射的调制符号s₂上的投影

所述估计模块，用于对所述的实部和虚部进行量化估计，得到所述的实部和虚部的估计值s_2,I、s_2,Q；

所述映射模块，用于对所述s_2,I、s_2,Q进行星座映射，得到调制符号s₂在星座图中的位置；

所述误差计算模块，用于根据所述估计得到的调制符号s₂在星座图中的位置，计算所述调制符号s₁的估计值与实际值之间的误差D；

所述最小误差更新模块，用于根据所述误差D，更新所述调制符号s₁的估计值与实际值之间的最小误差min_dist_0[k]和min_dist_1[k]，其中，k表示调制符号s₂的第k个比特；

所述对数似然比计算模块，用于根据所述最小误差min_dist_0[k]和min_dist_1[k]，计算所述当前发射终端发射的调制符号s₁的每个比特的对数似然比，得到译码结果。

13.根据权利要求12所述的装置，其特征在于，所述装置还包括初始化模块，用于在所述MIMO系统初始化时，将所述最小误差min_dist_0[k]和min_dist_1[k]中各元素的值初始化为正无穷大。

14.根据权利要求12或13所述的装置，其特征在于，

所述投影模块，还用于计算所述另一发射终端的信道响应h₂的共轭转置h₂ ^H，计算所述h₂ ^H与所述系统接收到的调制符号r的乘积以及计算所述h₂ ^H与所述当前发射终端的信道响应h₁的乘积；

所述估计模块，还用于计算所述另一发射终端的信道响应h₂的模方||h₂||²；计算所述的实部的中间估计值计算所述的虚部的中间估计值以及判断所述是否小于-2||h₂||²，并在所述小于-2||h₂||²时，令所述的估计值s_2,I等于0，在所述大于等于-2||h₂||²时，判断所述是否小于0，并在所述小于0时，令所述的估计值s_2,I等于1，在所述大于等于0时，判断所述是否小于2||h₂||²，在所述小于2||h₂||²时，令所述的估计值s_2,I等于2，在所述大于等于2||h₂||²时，令所述的估计值s_2,I等于3；判断所述是否小于-2||h₂||²，并在所述小于-2||h₂||²时，令所述的估计值s_2,Q等于0，在所述大于等于-2||h₂||²时，判断所述是否小于0，并在所述小于0时，令所述的估计值s_2,Q等于1，在所述大于等于0时，判断所述是否小于2||h₂||²，在所述小于2||h₂||²时，令所述的估计值s_2,Q等于2，在所述大于等于2||h₂||²时，令所述的估计值s_2,Q等于3；

所述最小误差更新模块，还用于判断所述调制符号s₁的第k个比特的值，并在所述时，判断所述误差D是否小于min_dist_0[k]，在所述误差D小于min_dist_0[k]时，令min_dist_0[k]的值等于所述误差D，在所述误差D大于等于min_dist_0[k]时，不更新min_dist_0[k]；或者在所述判断所述误差D是否小于min_dist_1[k]，在所述误差D小于min_dist_1[k]时，令min_dist_1[k]的值等于所述误差D，在所述误差D大于等于min_dist_1[k]时，不更新min_dist_1[k]。