CN107210981B - 用于解码通过通信信道接收到的数据块的方法和接收器 - Google Patents

用于解码通过通信信道接收到的数据块的方法和接收器 Download PDF

Info

Publication number
CN107210981B
CN107210981B CN201680007298.6A CN201680007298A CN107210981B CN 107210981 B CN107210981 B CN 107210981B CN 201680007298 A CN201680007298 A CN 201680007298A CN 107210981 B CN107210981 B CN 107210981B
Authority
CN
China
Prior art keywords
data block
qrd
las
value
size
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201680007298.6A
Other languages
English (en)
Other versions
CN107210981A (zh
Inventor
金耕进
菲利普·奥尔利克
秋浓俊昭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Publication of CN107210981A publication Critical patent/CN107210981A/zh
Application granted granted Critical
Publication of CN107210981B publication Critical patent/CN107210981B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045Arrangements at the receiver end
    • H04L1/0054Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03248Arrangements for operating in conjunction with other apparatus
    • H04L25/03286Arrangements for operating in conjunction with other apparatus with channel-decoding circuitry
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03CMODULATION
    • H03C3/00Angle modulation
    • H03C3/38Angle modulation by converting amplitude modulation to angle modulation
    • H03C3/40Angle modulation by converting amplitude modulation to angle modulation using two signal paths the outputs of which have a predetermined phase difference and at least one output being amplitude-modulated
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0242Channel estimation channel estimation algorithms using matrix methods
    • H04L25/0246Channel estimation channel estimation algorithms using matrix methods with factorisation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03203Trellis search techniques
    • H04L25/03216Trellis search techniques using the M-algorithm
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03891Spatial equalizers

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Artificial Intelligence (AREA)
  • Error Detection And Correction (AREA)
  • Radio Transmission System (AREA)

Abstract

一种方法,该方法解码通过通信信道接收到的数据块。该方法利用QR分解(QRD‑M)法确定对所述块中的比特的初始估计,该QRD‑M方法通过减小所述初始估计与所述块中的所述比特之间的累积距离来顺序地解码所述块中的所述比特。接着,所述方法利用似然上升(LAS)法确定所述块,该LAS方法从所述初始估计开始迭代地更新所述块的所述比特。

Description

用于解码通过通信信道接收到的数据块的方法和接收器
技术领域
本发明总体上涉及数字通信,并且更具体地,涉及解码通过通信信道接收到的数据块。
背景技术
多输入多输出(MIMO)系统使用多个发送天线和接收天线。与利用单一天线的系统相比,MIMO系统可以在不分配附加频率或发送功率的情况下与天线数量成比例地增加信道传送容量。MIMO系统的信道容量主要取决于在接收器中使用的信号检测方法,以便恢复所传送符号的块。重要的是设计MIMO系统的信号检测方法,以使实现高性能以及低的复杂性和检测延迟。对于经受时变衰落和多路径延迟的双选择性信道,这是一个特殊的问题,其中,发生载波间干扰(ICI)干扰和符号间干扰(ISI)两者。
MIMO系统的信号检测方法的示例可以包括:最大似然(ML)检测法、球形解码算法、QR分解与M-算法(QRD-M)算法等。
尽管ML检测法在MIMO系统中提供最优性能,但当发送天线的数量增加并且使用较高阶调制方法时,运算复杂性指数地增加。因此,存在的缺点在于事实上不使用ML检测法。
与ML检测法相比,球形解码算法提供与ML法相似的性能和显著降低的平均运算复杂性。然而,球形解码算法因信道矩阵的条件数和噪声分散而瞬间改变复杂性。结果,在最糟情况下,球形解码算法表现了与ML方法类似的运算复杂性。换句话说,球形解码算法的运算复杂性具有较大标准偏差和随机性。因此,难以将球形解码算法应用于移动基站具有有限功率和低检测等待时间容限的应用。
提供QRD-M法作为性能与复杂性之间的折衷方案。在QRD-M法中,无论信道条件或噪声功率如何,检测信号所需的计算量都是固定的。因此,检测信号的QRD-M法在每个处理中考虑更多的信息,由此使其可以进一步降低运算复杂性。换句话说,当存在条件好的信道环境或低噪声功率时,QRD-M法减少剩余候选符号的数量,由此减少与在每个分支处计算的累积距离有关的运算。然而,QRD-M法的检测性能取决于所选择的候选的数量,并且候选数量越多,运算复杂性就变得越大。
发明内容
在双重选择性信道上发送的单个符号块上的信道变化因ICI和ISI两者而不可忽略。最大似然检测(MLD)在高数据速率MIMO系统中是不可行的。几种检测方法可以实现近MLD性能。例如,QR分解与M-算法,即,QRD-M法,通过在树搜索的每一级处选择具有最小累积度量(metrics)的M个候选来降低复杂性。
然而,为了实现QRD-M算法的近MLD性能,应当使得M较大,由此使需要更高计算复杂性,这是不合需要的。用于低复杂性大MlMO检测的另一种似然上升搜索(LAS)法实现了近似最大似然(ML)性能,但仅适用于具有数百个天线的系统。在具有空间限制的通信终端中难于放置数百个天线。另外,LAS法仅针对相对较大的传送数据块具有好的性能。
本发明的一些实施方式基于这样的认识,即,对于通过多个通信信道传送的数据块的不同尺寸,不同类型的检测器可以比其它检测器更优化,尤其是考虑到接收器的有限计算资源。例如,对于要由具有预定计算能力的接收器检测的一个尺寸的块来说,QRD-M可以比LAS更优化。小块的尺寸和M值由解码器的计算能力决定。随着块尺寸的增加,LAS法可以变得更优化。
本发明的一些实施方式基于这样的认识,即,可以使用QRD-M与LAS的组合来检测不同尺寸的数据块。例如,具有可能的次优M的QRD-M的输出可以用作针对LAS的迭代的输入。这种组合对于通过具有有限计算能力的解码器解码中等尺寸的数据块来说能更好工作。而且,一些实施方式基于另一种认识,即,改变M值和LAS的迭代次数,可以改变QRD-M或LAS方法超过另一个方法的优势。
因此,一些实施方式基于这样的认识,即,用于检测块的方法应当被选择为块尺寸的函数。本发明的各种实施方式描述了针对不同尺寸的传送数据块的检测器的不同组合,并且提供了一种用于基于块的尺寸和/或接收器的计算能力来选择检测器中的一个或组合的方法。
因此,一个实施方式公开了一种用于解码通过通信信道接收到的数据块的方法。该方法包括:利用QR分解(QRD-M)法确定对所述块中的比特的初始估计,其中,所述QRD-M法通过减小所述初始估计与所述块中的所述比特之间的累积距离来顺序地解码所述块中的所述比特;以及利用似然上升(LAS)法确定所述块,其中,所述LAS法从所述初始估计开始迭代地更新所述块的所述比特,其中,所述方法的步骤由接收器的处理器来执行。
另一实施方式公开了一种接收器,该接收器包括:至少一个天线,所述至少一个天线用于通过通信信道接收数据块;以及处理器,该处理器用于解码所述数据块,其中,所述处理器被设置用于利用QR分解(QRD-M)法确定对所述块中的比特的初始估计,并且利用似然上升(LAS)法确定所述块,该LAS方法从所述初始估计开始迭代地更新所述块的所述比特。
附图说明
图1是采用本发明的一些原理的通信系统的示意图。
图2是根据本发明的一个实施方式的解码器的实现的示例。
图3是根据本发明的一些实施方式的自适应地改变解码方法的优势的方法的框图。
图4是根据本发明的一些实施方式的用于确定用于QRD-M法的M值的方法的框图。
图5是根据本发明的一些实施方式的用于选择用于LAS法的迭代次数N的方法的框图。
图6是根据本发明的一些实施方式的示例性无线通信系统的框图。
图7是根据本发明的一些实施方式的QRD-M法的树结构的示意图。
图8是根据本发明的一些实施方式的QRD-M法的框图。
图9是根据本发明的一个实施方式的作为用于LAS法的迭代的函数的对数似然度的图。
图10示出了根据本发明的一个实施方式的LAS法的框图。
具体实施方式
图1示出了采用本发明的一些实施方式的系统100的示意图。系统100包括经由信道150通信的发送器110和接收器120。该系统可以是多输入多输出(MIMO)系统,其中,发送器和/或接收器使用一个或多个天线111、121。发送器和接收器都可以包括用于处理传送数据的处理器和存储器以及连接到天线的至少一个RF链。
与利用单一天线的系统相比,MIMO系统可以在不分配附加频率或发送功率的情况下与天线数量成比例地增加信道150的传送容量。然而,本发明的一些实施方式考虑了在发送器和接收器处具有单个天线的点对点通信系统。
一些实施方式基于这样的一般认识,即,通过通信信道传送的数据块可以通过检测器的组合而不是单个检测器来更好地检测。检测器的组合例如可以基于所发送的数据块的尺寸和/或接收器的计算能力来选择。
例如,接收器120包括第一检测器130和第二检测器140。第一检测器提供对所发送的数据块115的初始估计。该初始估计135可以是次优的,但不是输出该次优检测,该初始估计135用作针对第二检测器140的输入,第二检测器140修正该初始估计以检测145该数据块。按这种方式,检测器甚至可以被设计成次优的,例如,由于接收器的有限计算资源,只要检测器的组合提供可接受的检测精度即可。
图2示出了根据本发明的一个实施方式的接收器120的解码器200的实现的示例。该实施方式基于这样的具体认识,即,QR分解与M(QRD-M)法和似然上升搜索(LAS)法的组合可以一起用于检测不同尺寸的数据块。例如,具有可能次优M的QRD-M的输出可以用作针对LAS的迭代的输入。
该实施方式对于具有有限计算能力的接收器来说是有利的。QRD-M方法通过在树搜索的每一级处选择具有最小累积度量的M个候选来降低复杂性。最优M值随着要检测的数据块的尺寸增加而增加。然而,QRD-M中最大M值由解码器的处理器的计算能力决定。对于具有足够小的尺寸的块115来说,最大可能M值可以足以最优地解码该数据块。在这种情况下,QRD-M法优于LAS。
例如,考虑最优的ML检测器,其具有以下形式
Figure GDA0001356681040000041
其中,长度N的矢量
Figure GDA0001356681040000042
表示对发送块的估计,并且其中,N对应于正被检测的块的尺寸。fy|x,H(·)是接收块y、给定发送块x以及信道H的似然函数。注意,上述变量的以下维度:x是长度N的矢量,其中,每个要素表示根据尺寸为
Figure GDA0001356681040000043
的字母表绘制的调制符号。例如,如果在发送器处采用QPSK调制,则
Figure GDA0001356681040000044
H是复值信道系数的N×N矩阵。
由此,ML检测器测试所有可能的传送块x,并针对给定信道选择最大化
Figure GDA0001356681040000051
的那一个。需要测试的可能发送块的数量取决于字母表的尺寸
Figure GDA0001356681040000052
和块的尺寸N两者,并且可以明确地写成
Figure GDA0001356681040000053
对于N=3并且
Figure GDA0001356681040000054
的示例树图来说,总共有43=64个可能块,因此作为潜在解,总共有64个长度为3的路径。ML检测器测试所有64个路径。
QRD-M法将检测问题处理为树形搜索问题,其中,
Figure GDA0001356681040000055
个可能块中的每一个都由贯穿具有N+1级的树的、长度为N的路径表示。QRD-M检测器通过将树中搜索的路径的数量迭代地修剪至仅M来操作。即,在每一级处,计算
Figure GDA0001356681040000056
个分支度量,其也是由QRD-M检测器测试的长度为N的路径的数量。因此,如在此使用的,用于QRD-M检测器的最优M值是这样的值,即,其允许QRD-M检测器获取到完全的ML检测器性能,即,最优M值使得
Figure GDA0001356681040000057
求解M生成
Figure GDA0001356681040000058
在QRD-M检测器和ML检测器相当的意义上,这是最优M值。
对于小块尺寸和低调制阶数来说,从计算和存储器的角度来看,设定
Figure GDA0001356681040000059
是可行的。然而,QRD-M检测器的益处是其用远小于
Figure GDA00013566810400000510
的M值获取近ML性能的能力。根据经验,设定
Figure GDA00013566810400000511
是足够大的,使得性能非常接近ML检测器。注意,虽然设定
Figure GDA00013566810400000512
生成具有近ML性能的检测器,但在ML检测器提供有关可达到的检测性能的下限的意义上,其仍是次优设定。针对任何M值和调制字母表尺寸
Figure GDA00013566810400000513
QRD-M法最终计算出
Figure GDA00013566810400000514
个分支度量。为了设定
Figure GDA00013566810400000515
分支度量计算的数量变为
Figure GDA00013566810400000516
并且是为获取近ML性能所需的QRD-M分支度量计算的典型数量。
一般来说,为了使QRD-M算法实现近MLD性能,应当使M较大,由此需要更高的计算复杂性,这是不合需要的。因此,本发明的一些实施方式针对QRD-M法使用次优M值。如在此使用的,QRD-M法的次优M值为
Figure GDA00013566810400000517
对于相对较大的传送数据块来说,用于低复杂性大MlMO检测的另一种似然上升搜索(LAS)法实现了近似最大似然(ML)的性能。LAS是一种贪婪搜索算法,其同样通过发现具有较大似然值的新的估计
Figure GDA00013566810400000518
来尝试迭代地改进其对传送块的估计
Figure GDA00013566810400000519
也就是说,对于i=0,...Niter-1来说,
Figure GDA0001356681040000061
其中,i是迭代计数,Niter是LAS算法中的最大迭代数,而
Figure GDA0001356681040000062
是在LAS算法的第i次迭代时对传送块的估计。在LAS算法的每次迭代中,检查/测试
Figure GDA0001356681040000063
中的N个符号中的每一个,以查看使用来自调制字母表
Figure GDA0001356681040000064
中的另一个值来替换该符号是否将增加似然函数。
因为LAS是一种贪婪算法,所以LAS可以收敛至局部最大值而不是全局最优最大值。然而,LAS在维度上的增加降低了一给定起始点(初始估计
Figure GDA0001356681040000065
)收敛至局部最小值的概率。换句话说,在每次迭代时LAS在维度上的增加提高了发现改进该似然度的维度的概率。针对i.i.d.信道(矩阵H的元素是独立且相同分布的),可以看出,随着块尺寸趋于无穷大,LAS的性能收敛至ML性能。虽然无限长的块尺寸在实践中不可行,但通常在200-400范围内的块尺寸足以观察到ML性能。
然而,一些实施方式基于这样的认识,即,如果选择LAS的起始点以使得初始估计
Figure GDA0001356681040000066
与少数位置中的ML解不同,则LAS可以对传送数据块的较小尺寸有效。这确保了LAS算法在全局最优解附近开始。需要初始估计值接近ML解的这种要求,因为该小块尺寸抑制了LAS算法用于从接近局部最小值的初始起始点恢复的能力。此外,本发明的一些实施方式基于这样的认识,即,甚至利用次优M值的QRD-M的输出可以提供这样的初始估计。
因此,在一个实施方式中,一种用于解码通过通信信道接收到的数据块的解码器包括第一检测器,该第一检测器用于利用QRD-M法210来确定对所述块115中的比特的初始估计215,该QRD-M法通过减小初始估计与所述块中的比特之间的累积距离来顺序地解码所述块,并且包括第二检测器,该第二检测器用于利用LAS法220确定所述块,该LAS法从初始估计开始迭代地更新所述块的比特。
在一些实施方式中,要发送的数据块的尺寸预先获知。在那些实施方式中,用于QRD-M检测器的M值和用于LAS检测器的迭代数N针对该尺寸和解码器的计算能力进行了优化。例如,假设该解码器具有足够大以存储q个分支度量的存储器,则可以将M参数可被设定成的最大值是q。该解码器的LAS部分中的最大迭代数Niter接着可以根据解码器的处理速度进行设定。
然而,在某些情况下,数据块的尺寸在收到之前是未知的。例如,在诸如根据IEEE802.11标准设计的可变长度数据传送中,其中,块的尺寸是可变的,并且该尺寸被提供在数据包报头中。
针对那些情况,本发明的一些实施方式针对块的平均尺寸确定M值和Niter。然而,另选实施方式根据块的尺寸来改变M值和LAS的迭代数Niter。那些实施方式基于这样的认识,即,对于较小的块尺寸,最大可能M值可以足够解码该块,所以不需要严重依赖于LAS。与此相反,对于较大的块尺寸,LAS方法不需要来自QRD-M的良好初始估计来检测块,从而可以进一步降低QRD-M法的复杂性。
图3示出了根据一些实施方式的方法的框图,其自适应地选择M的值和/或LAS的迭代数Niter来改变解码器中QRD-M和LAS方法相对于彼此的优势。该方法基于块的尺寸来确定用于QRD-M法的M值315和用于LAS法的迭代数335中的至少一个,其中,M是所述QRD-M法的每个阶段中的幸存路径的数量。接下来,该方法利用具有M值315的QRD-M来解码320该块,以确定该块的初始估计325,并且利用具有从初始估计325开始的Niter335次迭代来解码330块。
图4示出了用于确定用于QRD-M法的M值的方法的框图。该方法根据解码数据的接收器的计算能力405,确定410用于QRD-M法的最大M值,其可以包括存储器406的尺寸和接收器的处理器的计算速度或处理速度407中的至少一个。在某些实现中,这个步骤离线地执行。例如,如果接收器具有足以保持q个分支度量的存储器和至少足够快以在检测间隔中计算bmax个分支度量的处理速度407,则针对给定字母表尺寸和块尺寸,可以由硬件支持的最大M值为
Figure GDA0001356681040000071
本质上,最大M值由存储器和硬件的处理限制决定。
在在线解码期间,接收器的解码器将块尺寸430与阈值440进行比较420,并且如果块尺寸小于451阈值,则选择450用于QRD-M法的最大M值415,否则452选择小于最大M值的M值。在当块尺寸增加时M值减小的意义上,这个实施方式是反直觉的。然而,当这种复杂性不是必需时,该实施方式允许降低QRD-M法的复杂性。在一些实施方式中,M值被选择为一,以将QRD-M法的效果降低到其最小值。
在一些实施方式中,通过考虑q是对最大M值的硬件约束来选择阈值。由此,接收器计算420
Figure GDA0001356681040000072
并检查以查看其是否大于bmax。如果是,则在QRD-M算法中使用的M值被设定成Mmax 415。Mmax大于
Figure GDA0001356681040000073
是可能的。在这种情况下,检测器可以考虑给定的系统参数:块尺寸N、字母表尺寸
Figure GDA0001356681040000074
以及处理能力bmax使得QRD_M单独可以处理所接收的块。在这种情况下,当
Figure GDA0001356681040000075
时,检测器将顺序地执行QRD-M检测和LAS检测,其中QRD-M检测的输出充当用于LAS检测的初始估计
Figure GDA0001356681040000076
在这种情况下,M值取决于块尺寸,并且可以计算
Figure GDA0001356681040000081
其中,
Figure GDA0001356681040000082
指示顶函数(ceiling function)。上面示出的值256根据这样的事实提出,即,针对大块(>200个符号),当M小于最优经验值
Figure GDA0001356681040000083
时,LAS算法比QRD-M算法执行得更好。另外,随着块尺寸的增加,一些实施方式允许LAS算法处理更多的检测处理。随着块尺寸增加而减小M参数值的其它函数关系也是合适的。然而,独立于块尺寸地设定M=Mmax可导致检测器的QRD-M部分320的不必要的复杂性。
另外或者另选地,对于选择用于QRD-M法的M值来说,一些实施方式选择用于LAS方法330的迭代数。因为LAS法发现似然函数的局部最小值,所以一些实施方式不适应该迭代数Niter。本质上,在多次迭代后,LAS法收敛至最小值,并且额外的迭代不再导致改进的性能。因此,一些实施方式测试每个新的估计块
Figure GDA0001356681040000084
以查看任何组成部分是否已经从先前估计
Figure GDA0001356681040000085
改变。如果是,则潜在地需要进一步的迭代。
然而,一些实施方式限制总允许迭代数,以最小化针对每个块执行的计算数量,和或强调ORD-M法相对于LAS法的优势。
图5示出了根据块尺寸和/或M值来选择用于LAS法的最大迭代数N的框图。在该实施方式中,解码器将块尺寸430与阈值540进行比较520,如果块尺寸小于551阈值,则选择550最小N值515(例如,N=1)用于LAS法,否则552选择最大迭代数N,例如,需要不大于5-10次迭代以供收敛。
示例性实施方式
本发明的一个实施方式考虑了在发送器和接收器处具有单个天线的点对点通信系统。该实施方式考虑针对用于K个块符号的准静态信道的频率平坦衰落。
图6示出了包括采用本发明的一些原理的发送器110和接收器120的示例性无线通信系统的大体框图。发送器调制620输入比特610以生成符号625。调制620的示例包括但不限于BPSK(二进制相移键控)、QPSK(正交相移键控)、M进制QAM(正交幅度调制)。符号625可以包括一次采取N个的顺序符号的块。例如,矢量x(n)可以表示包括N个调制符号的第n个块。调制符号可以被假设成具有零均值和单位功率,并且在统计学上是独立的,由此,以下成立:
E{x(n))}=0和E{x(n)x(n)H}=IN
其中,运算符E{}表示期望值,IN是N×N单位矩阵。第n个块中的第j个符号传送由x(n)中的xj(n)表示,即,u(n)也表示针对该信道的N次使用符号的集合。
在一些实施方式中,通过N×N酉预编码矩阵P来对数据块进行预编码630。矩阵P可以是任何N×N酉矩阵,诸如离散傅里叶变换(DFT)或Walsh-Hadamard矩阵。一个实施方式还可以将随机旋转640应用于这样的预编码器输出,即,其采用具有沿其对角线的随机复数条目的N×N对角矩阵形式。该旋转640可以具有以下形式
Figure GDA0001356681040000091
其中,e是表示自然对数的基数的数,j是虚数单位,即,
Figure GDA0001356681040000092
并且θi(i=1,2,...N)是超过{0,2π}的i.i.d均匀随机变量。旋转矩阵具有在较高信噪比下改进误码率性能的效果。
因此,在应用预编码矩阵P和相位旋转矩阵Λ之后,针对信道H 650的输入645是u=PΛx(n)。因为P和Λ是酉(PΛ)(PΛ)H=IN。该信道H 650可以被建模为N×N对角矩阵,H=diag(h1,h2,...hN)其中,N个元素是具有单位方差的复数正规随机变量。将白高斯噪声z(n)添加660至信道,并且包括数据块的接收信号670是
y(n)=HPΛx(n)+z(n)=Heqx(n)+z(n)。
上述等式示出了该系统可以被看作具有N个虚拟空间流的N×N虚拟MIMO系统。因为所接收的矢量信号y(n)借助于块符号x(n)表达,所以根据本发明的各种实施方式,接收器120采用解码器对来自y(n)的块符号x(n)进行逐块检测。
QRD-M方法的示例
图7示出了具有QPSK调制和M=4的3×3MIMO系统的树结构700的示意图。粗体线指示借助于QRD-M算法的幸存路径。在这个示例中,例示了用于发送器和接收器处配备有三个天线的系统的检测的树结构。该阶数与天线数量相同,所以有三个阶段。该示例将QPSK用于传送符号,并且将QRD-M用于M=4,即,在检测处理中有四条幸存路径。在每个阶段710、720及730,针对所有可用路径计算累积的欧几里得距离度量,然后仅选择在所有路径当中具有最小累积欧几里德距离的M个路径(在这种情况下为四个路径)。
在第一阶段710中,因为当使用QPSK调制时,针对四个潜在符号只有四条路径,所以所有路径711被选择为幸存路径。在第二阶段720,仅从幸存路径计算累积欧氏距离的进一步计算。即,存在对应于所有路径的十六个累积欧氏距离,所以在第二阶段仅选择四个幸存路径721、722、723及724。在第三阶段730,QRD-M方法选择具有最小累积欧几里德距离的路径(例如,731)。该路径经由QRD-M给出最终数据决定,因为从根节点向下直至叶节点的、具有最小累积欧几里德距离的长度N路径表示最有可能被发送的符号序列。
图8示出了根据本发明的一个实施方式的QRD-M法的框图。该实施方式针对等效信道矩阵Heq应用QR分解801,以确定[Q,R]=QR(Heq),即,QR=Heq,其中,Q是酉矩阵而R是上三角矩阵,并且将y(n)乘以810Hermitian矩阵QH以生成
Figure GDA0001356681040000101
该实施方式还初始化计数器k=0、一组M个累加器acck(1)=0,...,acck(M)=0),以及一组幸存路径
Figure GDA0001356681040000102
如果820k小于821块尺寸N 820,则针对可能的分支数,
Figure GDA0001356681040000103
该实施方式针对j=1,2,...,ΔM计算840累积欧几里德距离
Figure GDA0001356681040000104
其中,
Figure GDA0001356681040000105
Figure GDA0001356681040000106
的第(N-K)个元素,R(N-K)是R的第(N-K)个行矢量,并且xj是在第j分支上的星座符号。
接下来,本实施方式在所有ΔM个分支当中选择850具有最小的累积欧几里德距离的仅M个分支,如M个分支(j1,j2,...,jM)及它们对应的累积欧几里得距离
Figure GDA0001356681040000107
该实施方式使860计数器k加一,并利用其先前幸存路径Sk={Sk-1;(j1,j2,...,jM)}和累积欧几里德距离
Figure GDA0001356681040000108
来更新幸存路径,并且该迭代返回870至测试820。当k=N 880时,该实施方式追溯在k=N具有最小累积欧几里得距离的幸存路径Sk,以使可以进行数据决定
Figure GDA0001356681040000109
LAS方法的示例
图9示出了根据本发明的一些实施方式的作为用于LAS法的迭代的函数的对数似然度的图。如前所述,LAS算法是一种贪婪搜索算法,其寻求通过在每次迭代中发现具有较大似然度的数据矢量x来改进其当前估计。下面描述如何执行搜索的细节,并且图9旨在示出针对每次迭代,候选解具有比先前值大的似然值。实际上,如果当前候选解无法改进先前解,则算法终止,因为它已经找到了全局最优解或者已经稳定在局部最大值。
图10示出了根据一个实施方式的LAS法的实现的图。本实施方式假设在发送器中使用BPSK调制,这相当于设定
Figure GDA0001356681040000111
根据针对所接收的信号y(n)的表达式,利用等效信道矩阵将矩阵G定义1010为G=2Re(Heq HHeq),并且形成202符号块的初始估计
Figure GDA0001356681040000112
初始估计的一个示例是
Figure GDA0001356681040000113
在k=0时,其中,k指示迭代数,形成1020块符号的初始估计。对于给定迭代,k1030小于1031Niter,该实施方式计算N乘1矢量g(k),其第j个元素由gj k表示1040。该计算使用等效信道Heq和先前数据矢量估计
Figure GDA0001356681040000114
来确定
Figure GDA0001356681040000115
接下来,该实施方式考虑先前迭代中的所估计的符号,将gj k与矩阵G的第k个对角元素的量值的进行比较,更新1050
Figure GDA0001356681040000116
的第j个符号的估计,并且进行1070至下一个符号。当所需的迭代数Niter达到1032时,获取1035数据检测。
本发明的上述实施方式可以按许多方式中的任一种来实现。例如,这些实施方式可以利用硬件、软件或其组合来实现。当按软件来实现时,软件代码可以在任何合适处理器或处理器集合上执行,而不管设置在单一计算机中还是在多个计算机当中分布。这种处理器可以被实现为集成电路,在集成电路组件中具有一个或更多个处理器。然而,处理器可以利用采用任何合适格式的电路来实现。
而且,本发明的实施方式可以被具体实施为已经提供了其一示例的方法。作为该方法的一部分执行的动作可以按任何合适方式来排序。因此,即使在示例性实施方式中被示出为顺序动作,也可以构造按与所例示相比不同的次序来执行动作的实施方式,其可以包括同时执行一些动作。
在权利要求书中使用诸如“第一”、“第二”的序数词来修改权利要求部件本身并不意味着一个权利要求要素相对于另一个要素的任何优先级、优先权或次序或执行方法的动作的时间次序。而是仅仅被用作用于区分具有特定名称的一个权利要求要素与具有相同名称(但供序数词使用)的另一要素的标记,以区分这些权利要求要素。

Claims (9)

1.一种用于解码通过通信信道发送的数据块的方法,该方法包括:
接收通过通信信道发送的数据块;
确定所述数据块的尺寸;
基于所述数据块的尺寸来选择用于QRD-M法的M值和用于似然上升(LAS)法的迭代数中的至少一个,其中,M是该QRD-M法的每个阶段中的幸存路径的数量;
利用所述QRD-M法确定对所述数据块中的比特的初始估计,其中,所述QRD-M法通过减小所述初始估计与所述数据块中的所述比特之间的累积距离来顺序地解码所述数据块中的所述比特;以及
利用该似然上升(LAS)法确定所述数据块,其中,该似然上升(LAS)法从所述初始估计开始迭代地更新所述数据块的所述比特,其中,所述方法的步骤由接收器的处理器执行。
2.根据权利要求1所述的方法,所述方法还包括:
根据所述接收器的存储器的大小和所述处理器的计算速度,选择用于所述QRD-M法的最大M值;
如果所述数据块的尺寸小于阈值,则选择用于所述QRD-M法的所述最大M值;以及
如果所述数据块的尺寸不小于阈值,则选择小于所述最大M值的M值。
3.根据权利要求2所述的方法,其中,所述最大M值为
Figure FDA0002539570340000011
其中,N是所述似然上升(LAS)法的迭代数,而
Figure FDA0002539570340000012
是针对所述数据块中的所述比特的值的字母表的尺寸。
4.根据权利要求2所述的方法,其中,所述最大M值是所述似然上升(LAS)法的、能够存储在所述接收器的存储器中的最大分支数。
5.根据权利要求2所述的方法,所述方法还包括:
根据利用所述数据块传送的数据包报头来确定所述数据块的尺寸。
6.根据权利要求1所述的方法,所述方法还包括:
如果所述数据块的尺寸小于阈值,则最小化所述似然上升(LAS)法的迭代数。
7.一种接收器,该接收器包括:
至少一个天线,所述至少一个天线用于通过通信信道接收数据块;以及
处理器,该处理器用于解码所述数据块,其中,所述处理器被设置用于:利用QRD-M法确定对所述数据块中的比特的初始估计,该QRD-M法通过减小所述初始估计与所述数据块中的所述比特之间的累积距离来顺序地解码所述数据块;以及利用似然上升(LAS)法解码所述数据块,该似然上升(LAS)法从所述初始估计开始迭代地更新所述数据块的所述比特,
其中,所述处理器确定所述数据块的尺寸,并且基于所述数据块的尺寸来选择所述QRD-M法的每个阶段中的幸存路径的数量M值和用于所述似然上升(LAS)法的迭代数中的至少一个。
8.根据权利要求7所述的接收器,该接收器还包括:
存储器,所述存储器存储基于所述接收器的所述存储器的大小和所述处理器的处理速度确定的最大M值,其中,如果所述数据块的尺寸小于阈值,则所述处理器从所述存储器选择用于所述QRD-M法的所述最大M值,否则确定小于所述最大M值的M值。
9.根据权利要求8所述的接收器,其中,所述处理器确定所述最大M值为
Figure FDA0002539570340000021
其中,N是所述似然上升(LAS)法的迭代数,而
Figure FDA0002539570340000022
是针对所述数据块中的所述比特的值的字母表的尺寸。
CN201680007298.6A 2015-01-26 2016-01-15 用于解码通过通信信道接收到的数据块的方法和接收器 Active CN107210981B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US14/604,871 US9525514B2 (en) 2015-01-26 2015-01-26 System and method for decoding block of data received over communication channel
US14/604871 2015-01-26
PCT/JP2016/051760 WO2016121625A1 (en) 2015-01-26 2016-01-15 Method for decoding block of data received over communication channel and receiver

Publications (2)

Publication Number Publication Date
CN107210981A CN107210981A (zh) 2017-09-26
CN107210981B true CN107210981B (zh) 2020-08-11

Family

ID=55315677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680007298.6A Active CN107210981B (zh) 2015-01-26 2016-01-15 用于解码通过通信信道接收到的数据块的方法和接收器

Country Status (5)

Country Link
US (1) US9525514B2 (zh)
EP (1) EP3251312B1 (zh)
JP (1) JP6272574B2 (zh)
CN (1) CN107210981B (zh)
WO (1) WO2016121625A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10355762B2 (en) * 2017-07-07 2019-07-16 Mitsubishi Electric Research Laboratories, Inc System and method for adaptive beamforming communication
CN109347772A (zh) * 2018-11-20 2019-02-15 重庆邮电大学 广义空间调制系统中改进的qrd-m检测算法
WO2022061788A1 (en) * 2020-09-25 2022-03-31 Intel Corporation Versatile systolic array for maximum likelihood mimo detectors
CN116938394B (zh) * 2023-09-18 2024-01-30 北京智芯微电子科技有限公司 Mimo检测方法、接收方法及其系统

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE513813C2 (sv) * 1993-02-05 2000-11-06 Ericsson Ge Mobile Communicat Anordning med SAW-filter för duplexöverföring
US6201802B1 (en) * 1997-08-29 2001-03-13 Qualcomm Inc. Method and apparatus for analyzing base station timing
JP2002050980A (ja) * 2000-08-04 2002-02-15 Matsushita Electric Ind Co Ltd 高周波スイッチおよびそれを用いた無線通信機
US7359684B2 (en) * 2001-11-06 2008-04-15 Texas Instruments Incorporated Wireless communication system with variable intermediate frequency transmitter
JP2003152815A (ja) * 2001-11-14 2003-05-23 Hitachi Ltd 通信用半導体集積回路
JP4154949B2 (ja) * 2002-08-06 2008-09-24 松下電器産業株式会社 Sawフィルタ
US6903665B2 (en) * 2002-10-30 2005-06-07 Spacebridge Semiconductor Corporation Method and apparatus for error control coding in communication systems using an outer interleaver
GB2416437B (en) * 2004-07-19 2007-12-27 Renesas Tech Corp A communication semiconductor integrated circuit device and a wireless communication system
US20070171312A1 (en) * 2005-12-20 2007-07-26 Samsung Electronics Co., Ltd. Frequency converter
US8116411B2 (en) * 2007-08-06 2012-02-14 Indian Institute Of Science Method to detect data transmitted from multiple antennas and system thereof
US8245104B2 (en) * 2008-05-02 2012-08-14 Lsi Corporation Systems and methods for queue based data detection and decoding
US8438463B2 (en) * 2010-12-08 2013-05-07 Research In Motion Limited Decoding tail-biting convolutional codes
US8503584B2 (en) * 2010-12-21 2013-08-06 Lsi Corporation Efficient implementation of M-algorithm based on QR decomposition for higher-order constellations
US8737540B1 (en) * 2011-07-15 2014-05-27 Qualcomm Atheros, Inc. System and method for providing reduced complexity maximum likelihood MIMO detection
CN103095511A (zh) * 2011-10-28 2013-05-08 华为技术有限公司 一种在IPsec机制下的网络测试方法,装置及系统
JP2013197751A (ja) * 2012-03-16 2013-09-30 Fujitsu Ltd 無線装置、無線装置制御方法、無線装置制御プログラム

Also Published As

Publication number Publication date
JP6272574B2 (ja) 2018-01-31
US20160218827A1 (en) 2016-07-28
US9525514B2 (en) 2016-12-20
CN107210981A (zh) 2017-09-26
JP2017535213A (ja) 2017-11-24
WO2016121625A1 (en) 2016-08-04
EP3251312B1 (en) 2021-10-06
EP3251312A1 (en) 2017-12-06

Similar Documents

Publication Publication Date Title
KR101124863B1 (ko) 다중 소스로부터의 통신신호를 처리하는 장치 및 방법
JP4429945B2 (ja) Mimo多重通信装置および信号分離方法
EP3104564B1 (en) Anticipated termination for sequential decoders
TWI495305B (zh) 利用依序檢索的二階等化之方法及接收機
US8488721B2 (en) Adaptive QRD-M algorithm based signal detecting method by using constellation set grouping in spatial multiplexing multiple-input multiple-output system
CN107210981B (zh) 用于解码通过通信信道接收到的数据块的方法和接收器
KR20080034050A (ko) 다중 안테나를 사용하는 통신 시스템에서 신호 검출 장치 및 방법
CN107094063B (zh) 半穷举迭代块解码方法和设备
US9054845B2 (en) Radio communication apparatus and signal processing method
US8391421B2 (en) EDA-based detection of communication signals
JP6180333B2 (ja) 無線周波数受信機において信号を復号化する方法
US20150180682A1 (en) Turbo Equalisation
EP3104565A1 (en) Sequential decoding with stack reordering
JP5531089B2 (ja) Mimo受信信号を検出するための深さ優先探索による木探索方法
CN109167648B (zh) 候选星座点集合生成方法及mimo空间复用检测方法
WO2021068007A2 (en) Apparatus and method of multiple-input multiple-output detection with successive transmission layer detection and soft interference cancellation
WO2008151164A1 (en) N-candidate depth-first decoding
EP3836499B1 (en) Device and method for training a model
CN109039539B (zh) 候选星座点集合生成方法及mimo空间复用检测方法
Lai et al. Channel-aware local search (CA-LS) for iterative MIMO detection
TWI401905B (zh) Multiple Input Multiple Output Detection Method and Detector
US7844008B2 (en) Delayed branch decision in quadrature decomposition with M-searching
Fan et al. A new dynamic K-best SD algorithm for MIMO detection
WO2013137760A1 (en) Device and method for detecting transmitted signals in a mimo communication system
KR100989888B1 (ko) 다중 입력 다중 출력 시스템의 신호 검출 방법

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant