CN103973602A - 信号检测方法和装置 - Google Patents

信号检测方法和装置 Download PDF

Info

Publication number
CN103973602A
CN103973602A CN201310031384.9A CN201310031384A CN103973602A CN 103973602 A CN103973602 A CN 103973602A CN 201310031384 A CN201310031384 A CN 201310031384A CN 103973602 A CN103973602 A CN 103973602A
Authority
CN
China
Prior art keywords
measurement value
branched measurement
increment
judgement amount
constellation point
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.)
Granted
Application number
CN201310031384.9A
Other languages
English (en)
Other versions
CN103973602B (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.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
China Academy of Telecommunications Technology CATT
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 China Academy of Telecommunications Technology CATT filed Critical China Academy of Telecommunications Technology CATT
Priority to CN201310031384.9A priority Critical patent/CN103973602B/zh
Priority to PCT/CN2014/071178 priority patent/WO2014114246A1/zh
Publication of CN103973602A publication Critical patent/CN103973602A/zh
Application granted granted Critical
Publication of CN103973602B publication Critical patent/CN103973602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明实施例公开了一种信号检测方法和装置,涉及无线通信领域,用于降低信号检测的复杂度,提高检测效率。本发明中,在采用QRDM-List算法对接收信号进行信号检测的过程中,在对除最先搜索层之外的其他传输层进行符号搜索时,只需对部分判决量通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量,对于其他判决量,仅需采用星座符号的硬判决方法确定与其欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量,与现有技术中需要对所有判决量都通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量相比,本发明的信号检测复杂度较低,进而提高了信号检测效率。

Description

信号检测方法和装置
技术领域
本发明涉及无线通信领域,尤其涉及一种信号检测方法和装置。
背景技术
在空间复用传输的多天线系统中,在实现中可采用的一种近似最大似然的检测算法为基于查找表的QR分解和M算法的(List based QR decompositionwith M algorithm,QRDM-List)检测算法,该算法采用逐层选取最有可能多个的星座符号保留,并选取有可能的星座符号的组合(该星座符号的组合是一个矢量,在本检测算法中称为分支),作为输出的检测结果。
假设空间复用传输的MIMO系统的传输模型如下公式1:
r=Hx+n…………………………………….(1)
其中,r为NR×1的接收信号矢量,H为NR×NL的信道矩阵,x表示NL×1的发送信号向量,n表示NR×1的噪声向量。
对于基于QR分解的接收端检测算法,如公式1转化为下式2:
y = Q H r
= Q H ( Hx + n )
= Rx + n ~ . . . ( 2 )
上式2中,Q矩阵是酉矩阵,R矩阵是上三角矩阵,如下式3所示:
对式1两端同时左乘QH,可以得到下式4:
y = Q H ( Hx + n )
= Rx + n ~ . . . ( 4 )
其中,y=QHr,
QRDM-List检测算法是一种基于存储表的简化最大似然检测算法。该算法在每一层搜索的过程中根据判决量的位置,通过查找预先存储的星座图上的各区域和最近的星座点的映射表格的方法,确定NC2(NC2<Q)个与判决量位置较近的星座点,对这些星座点计算对应的分支度量值。假设每一层保留NC1个分支,这样每一层需要计算NC1×NC2个分支度量值,并在这NC1×NC2个分支度量值中选取最小的NC1个度量值和分支保留。即在QRDM-List检测算法中,在保留分支下并没有遍历所有星座点计算分支度量值,而是通过查表的方法,只计算一部分星座点的分支度量值。
QRDM-List检测算法和QR分解和M算法的最大似然的(QR detection withM algorithm Maximum Likelihood Detection,QRD-M)检测算法的区别仅在树形搜索部分,之前的信道矩阵的QR分解和之后的软比特计算方法都是相同的。该算法的执行步骤如下:
步骤0:预先针对各种不同的调制方式,初始化查找表,通过获得的判决量确定,该判决量所述的区间,每个区间对应于NC2个保留星座点。
步骤1:对信道矩阵H进行QR分解,输出Q矩阵和R矩阵。其中Q矩阵是酉矩阵,R矩阵是上三角矩阵。初始化Xleft为NL×NC1的空矩阵,bm为1×NC1的向量。
步骤2:对接收信号向量y左乘QH,得到等效接收向量
步骤3:对NL层传输的符号从第NL层到第1层进行逐层搜索。
步骤3.1:对于第NL层,即第一个被搜索的符号,将Q个星座点带入下式进行分支度量值的计算:
d 1 , q = | y N L - r N L , N L x ^ q | 2 ;
其中,表示当前层传输的符号,来自星座集合。在Q个分支度量值中保留NC1个度量值最小的存入bm,并将对应的星座点存入Xleft的第NL行。
对于第i,i∈{2,…,NL}个被搜索的符号,即第NL-i+1层的符号,先根据已搜索层的保留分支计算在对应保留分支的情况下的判决量:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
上式中,nc1∈{1,…,NC1}。根据和调制方式,通过查表的方法确定和欧氏距离最近的NC2个星座点这样在该层需要进行NC1×NC2个分支度量值的增量计算
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
上式中,nc2∈{1,…,NC2}。累积分支度量值采用分支度量值的增量和保存的累积分支度量值bm中对应的元素进行相加:
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
每一层在NC1×NC2个分支中保留最小的NC1个累积分支度量值,将对应的累积分支度量值存入bm,将对应的已搜索层保留分支和当前层保留的星座点写入Xleft
最后,所有符号搜索完成之后输出的bm为保留的最小的NC1个分支对应的分支度量值,保留分支Xleft包含对应的bm最小的NC1个分支。
QRDM-List的树形搜索示意图如图1所示,该图示中假设NC1=16,NC2=4。
在未编码系统中,将Xleft中BM最小值对应的矢量转化为比特输出即为检测结果。然而在编码的系统中,需要生成每个比特的对数似然比(也称软比特)输入到译码器。该算法通过树形搜索输出保留分支并计算对应的分支度量值,采用保留分支中该比特为0的最小累积分支度量值和保留分支中该比特为1的最小累积分支度量值的差值,作为该软比特的输出。如下式所示:
b ~ ( n l ) ( n b ) &ap; min x &Element; X left n l , n b , 0 ( BM ) - min x &Element; X left n l , n b , 1 ( BM ) ;
上式中,表示保留分支中第nl层第nb比特为0的分支集合,表示留分支中第nl层第nb比特为1的分支集合,nl∈{1,…,NL},nb∈{1,…,log2Q}。NL表示传输的层数,log2Q表示调制阶数。
综上,现有的QRDM-List检测算法存在两方面的缺点,第一,该算法对所有的保留分支都进行基于查表法的进一步搜索,而各条保留分支是正确检测结果的概率则是有差别的,分支度量值越小对应的分支正确的概率越大,这样对各分支采用相同的处理,显然会计算了更多的分支度量值,从而引入了不必要的复杂度;第二,由于引入了一些不必要的分支度量值的计算,从而增加了排序部分的输入长度,这样也会增加排序部分复杂度和处理时延。
发明内容
本发明实施例提供一种信号检测方法和装置,用于降低信号检测的复杂度,提高检测效率。
一种信号检测方法,该方法包括:
在采用基于QR分解和M算法的查找表最大似然检测QRDM-List算法对接收信号进行信号检测的过程中,对多个传输层逐层进行符号搜索,并且,按照如下方法对除最先搜索层之外的其他传输层进行符号搜索:
确定已搜索传输层的每个保留星座符号组合对应的判决量;
选取确定的各判决量中的部分判决量,对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;
对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
一种信号检测装置,该装置包括:
符号搜索单元,用于在采用基于QR分解和M算法的查找表最大似然检测QRDM-List算法对接收信号进行信号检测的过程中,对多个传输层逐层进行符号搜索;
所述符号搜索单元包括判决量确定单元、判决量选取单元、第一度量值增量确定单元、第二度量值增量确定单元和结果确定单元;
所述判决量确定单元,用于对除最先搜索层之外的其他传输层进行符号搜索时,确定已搜索传输层的每个保留星座符号组合分别对应的判决量;
所述判决量选取单元,用于选取确定的各判决量中的部分判决量;
所述第一度量值增量确定单元,用于对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;
所述第二度量值增量确定单元,用于对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
所述结果确定单元,用于根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
本发明实施例提供的方案中,在采用QRDM-List算法对接收信号进行信号检测的过程中,在对除最先搜索层之外的其他传输层进行符号搜索时,只需对部分判决量通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量,对于其他判决量,仅需采用星座符号的硬判决方法确定与其欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量,与现有技术中需要对所有判决量都通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量相比,本发明的信号检测复杂度较低,进而提高了信号检测效率。
附图说明
图1为现有技术中QRM-List算法的树形搜索过程示意图;
图2为本发明实施例提供的方法流程示意图;
图3为本发明实施例中的QRDM-List&HD检测算法的执行示意图;
图4为本发明实施例中的QRDM-List&HD检测算法树形搜索的执行示意图;
图5为本发明实施例中的QRM-List&HD检测算法的树形搜索过程示意图;
图6为本发明实施例中的16QAM的星座区域划分示意图;
图7为本发明实施例提供的装置结构示意图。
具体实施方式
为了降低信号检测的复杂度,提高检测效率,本发明实施例提供一种信号检测方法。
参见图2,本发明实施例提供的信号检测方法中,采用基于QR分解和M算法的查找表最大似然检测(QRDM-List)算法对接收信号进行信号检测,在信号检测过程中,对多个传输层逐层进行符号搜索,其中对除最先搜索层之外的其他传输层进行符号搜索的方法包括如下步骤:
步骤20:确定已搜索传输层的每个保留星座符号组合对应的判决量;
步骤21:选取确定的各判决量中的部分判决量;
步骤22:对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
步骤23:根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
步骤20中,确定已搜索传输层的每个保留星座符号组合对应的判决量,具体实现可以如下:
按照如下公式确定已搜索传输层的第nc1个保留星座符号组合对应的判决量nc1∈{1,…,NC1},NC1为保留星座符号组合的个数:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;ri,i为所述R矩阵中第i行第i列的元素。
步骤21中,选取确定的各判决量中的部分判决量,具体实现可以如下:
选取确定的各判决量中累积分支度量值最小的设定数目个判决量,该设定数目小于判决量的总个数。这里,累积分支度量值最小的设定数目个判决量是指在上一次符号搜索过程中得到的累积分支度量值最小的设定数目个保留星座符号组合对应的判决量。
步骤23中,根据确定的每个分支度量值的增量得到本传输层的符号搜索结果,具体实现可以如下:
对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值,该累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留星座符号组合的分支度量值之和;
从得到的各累积分支度量值中选取NC1个最小的累积分支度量值,将选取的NC1个最小的累积分支度量值对应的星座点组合作为本传输层的保留星座符号组合并保存到Xleft中,将选取的NC1个最小的累积分支度量值作为各保留星座符号组合分别对应的分支度量值并保存到bm中,其中Xleft为NL行NC1列的矩阵;bm为长度为NC1的向量;NL为传输层的总层数,NC1为保留星座符号组合的个数。
具体的,在确定星座点对应的分支度量值的增量时,若该星座点是通过查表确定的,则确定该星座点对应的分支度量值的增量的方法如下:
按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量nc1∈{1,…,NC3},1≤NC3≤NC1,NC1为保留星座符号组合的个数,nc2∈{1,…,NC2},NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数:
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的第nc2个星座点;ri,i为所述R矩阵中第i行第i列的元素。
相应的,在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是通过查表确定的星座点对应的分支度量值的增量,则按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp((nc1-1)×NC2+nc2):
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
具体的,在确定星座点对应的分支度量值的增量时,若该星座点是通过星座符号硬判决方法确定的,则确定该星座点对应的分支度量值的增量的方法如下:
按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量nc1∈{NC3+1,…,NC1},1≤NC3≤NC1,NC1为保留星座符号组合的个数,NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数,NC3为选取的部分判决量的个数:
d i , NC 3 &times; NC 2 + nc 1 - NC 3 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,..,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的一个星座点;ri,i为所述R矩阵中第i行第i列的元素。
相应的,在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是通过星座符号硬判决方法确定的星座点对应的分支度量值的增量,则按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp(NC3×NC2+nc1-NC3):
BM tmp ( NC 3 &times; NC 2 + nc 1 - NC 3 ) = bm ( nc 1 ) + d i , NC 3 &times; NC 2 + nc 1 - NC 3 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
本发明中,每个传输层的保留星座符号组合的个数为16,通过查表确定的星座点的个数为4,选取的部分判决量的个数为4。
下面结合具体实施例对本发明进行说明:
本实施例给出一种结合存储表和硬判的简化最大似然检测算法,QRDM-List&HD(List&Hard Decision based QR decomposition with Malgorithm)算法。该算法根据保留分支的度量值大小确定对该保留分支下一步进行基于查表法的保留还是进行硬判仅保留一条路径,对于分支度量值较小的保留分支进行查表法的扩展获得当前传输层用于计算分支度量值的多条保留分支,对于分支度量值较大的分支进行硬判获得当前层用于计算分支度量值的一条保留分支。即通过进一步删减分支度量值较大的分支度量值的扩展分支,减少了需要计算的分支度量值,从而也减少了排序输入的序列长度。
在执行QRDM-List&HD算法前,需要预先针对各种不同的调制方式,初始化查找表,该查找表中保存有星座图的区域标识与对应区域中保留星座点的映射关系;其中区域标识对应的保留星座点的个数为NC2
如图3所示,信号检测流程如下:
步骤1:对信道矩阵H进行QR分解,输出Q矩阵和R矩阵。其中Q矩阵是酉矩阵,R矩阵是上三角矩阵。初始化Xleft为NL×NC1的空矩阵,bm为1×NC1的向量;
对接收信号向量y左乘QH,得到等效接收向量
步骤2:从第NL个传输层到第1个传输层逐层进行基于查表和硬判决的树形符号搜索。如图4所示,本步骤的执行如下:
A、对于第NL个传输层,即第1个被搜索的符号,将星座集合中的Q个星座点分别带入下式进行分支度量值的计算:
d 1 , q = | y N L - r N L , N L x ^ q | 2 ;
其中,表示当前的星座点。表示等效接收向量中的第NL个元素,表示R矩阵中第NL行第NL列的元素。
从计算出的Q个分支度量值中选取出最小的NC1个分支度量值,将该NC1个分支度量值按从小到大的顺序存入bm,并将该NC1个分支度量值分别对应的星座点存入Xleft的第NL行。
B、对于第NL-i+1个传输层,即第i(i∈{2,…,NL})个被搜索的符号,先计算已搜索传输层的每个保留分支即每个保留星座符号组合(一个保留星座符号组合是Xleft中的一列星座点)对应的判决量:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
上式中,nc1∈{1,…,NC1}。对于这NC1个判决量分为两种不同的情况进行处理。
情况1:选取NC1个判决量中累积分支度量值最小的NC3(1≤NC3≤NC1)个判决量,对于选取的每个判决量,通过查表的方法确定与该判决量(nc1∈{1,…,NC3)的欧氏距离最近的NC2个星座点并按照下式计算该NC2个星座点中每个星座点对应的分支度量值的增量:
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
这样在当前传输层的前NC3个保留分支下,共需要进行NC3×NC2个分支度量值的增量计算。
通过查表的方法确定与判决量的欧氏距离最近的NC2个星座点的方法为:在预先生成的查找表中查找数值与该判决量相同的区域标识,将该区域标识对应的保留星座点作为与判决量的欧氏距离最近的NC2个星座点。
上式中,nc1∈{1,…,NC3},nc2∈{1,…,NC2}。
对于得到的每个分支度量值的增量,按照下式确定该分支度量值的增量对应的累积分支度量值:
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留分支的分支度量值之和,即累积分支度量值是将分支度量值的增量和保存的bm中对应的元素进行相加得到的。
情况2:对于NC1个判决量中未被选取的后NC1-NC3个判决量(nc1∈{NC3+1,…,NC1}),通过星座符号硬判决的方法确定与(nc1∈{NC3+1,…,NC1})的欧式距离最近的1个星座点并按照下式计算该星座点对应的分支度量值的增量;
d i , NC 3 &times; NC 2 + nc 1 - NC 3 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
这样在当前传输层的后NC1-NC3个保留分支下,共需要进行NC1-NC3个分支度量值的增量计算。
上式中,nc1∈{NC3+1,…,NC1};
对于得到的每个分支度量值的增量,按照下式确定该分支度量值的增量对应的累积分支度量值:
BM tmp ( NC 3 &times; NC 2 + nc 1 - NC 3 ) = bm ( nc 1 ) + d i , NC 3 &times; NC 2 + nc 1 - NC 3 ;
累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留分支的分支度量值之和,即累积分支度量值是将分支度量值的增量和保存的bm中对应的元素进行相加得到的。
两种情况下,共计算了NC3×NC2+NC1-NC3个分支度量值的增量和累积分支度量值。每一传输层在NC3×NC2+NC1-NC3个累积分支度量值中保留最小的NC1个累积分支度量值,将保留的NC1个累积分支度量值按从小到大的顺序存入bm,将保留的NC1个累积分支度量值对应的星座点组合保存到Xleft
最后,所有传输层的符号搜索完成之后输出bm和Xleft,Xleft中的每一列表示NL各传输层可能传输的星座点组合,bm的每个元素表示Xleft中对应列的星座点组合的分支度量值。
步骤3:根据步骤2中的树形符号搜索结果进行软解调即软比特的计算,输出得到的软比特,得到信号检测结果。
一个较优实施例如下:
对于LTE系统采用4层空间复用传输的模式下,以16QAM调制方式下的检测为例,假设NC1=16,NC2=4,NC3=4存储的映射表格和星座划分如图5和表1所示。表1中的候选列表一栏中存储的是图5中星座点的编号,表1中的区域一栏存储的是图5中的区域标识。
区域 候选列表
1 13,14,15,16
2 9,11,13,15
3 9,10,11,12
4 5,6,13,14
5 1,5,9,13
6 1,2,9,10
7 5,6,7,8
8 1,3,5,7
9 1,2,3,4
表116QAM下的查找表
信号检测流程执行过程如下:
步骤0:根据图5和表1的方式划分星座区间和建立各区域标识与保留星座点的映射表。
步骤1:对信道矩阵H进行QR分解,输出Q矩阵和R矩阵。其中Q矩阵是酉矩阵,R矩阵是上三角矩阵。初始化Xleft为4×16的空矩阵,bm为1×16的向量。
步骤2:对接收信号向量y左乘QH,得到等效接收向量y=[y1,…y4]T
步骤3:对4层传输的符号从第4层到第1层进行逐层搜索。
步骤3.1:对于第4层,即第一个被搜索的符号,将16个星座点带入下式进行分支度量值的计算:
d 1 , q = | y 4 - r 4,4 x ^ q | 2
其中,表示当前层传输的符号,来自星座集合。计算全部16个分支度量值按从小到大的顺序存入bm,并将对应的星座点存入Xleft的第4行。(因为星座点数目和保留星座点个数相等,所以不必进行排序,全部保留)
对于第i(i∈{2,…,4})个被搜索的符号,即第4-i+1层的符号,先根据已搜索层的保留分支计算在对应保留分支的情况下的判决量:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 4 r i , j X j , nc 1 left ) / r i , i ;
上式中,nc1∈{1,…,16}。对于这16个判决量分为两种不同的情况进行处理。
情况1:对于其中的16个判决量中的前4个判决量通过查表的方法确定和nc1∈{1,…,4}欧氏距离最近的4个星座点这样在该层前4个保留分支下,需要进行4×4=16个分支度量值的增量计算。
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2
上式中,nc1∈{1,…,4},nc2∈{1,…,4}。累积分支度量值采用分支度量值的增量和保存的累积分支度量值bm中对应的元素进行相加,即:
BM tmp ( ( nc 1 - 1 ) &times; 4 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; 4 + nc 2 ;
情况2:对于其中的16个判决量中的后12个判决量通过硬判决的方法确定和nc1∈{5,…,16}欧式距离最近的1个星座点这样在该层后12个保留分支下,需要进行12个分支度量值的增量计算。
d i , 12 + nc 1 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
上式中,nc1∈{5,…,16},累积分支度量值采用分支度量值的增量和保存的累积分支度量值bm中对应的元素进行相加,即:
BM tmp ( 12 + nc 1 ) = bm ( nc 1 ) + d i , 12 + nc 1 ;
两种情况下,共计算了28个分支度量值。每一层在28个分支中保留最小的16个累积分支度量值,将对应的累积分支度量值按从小到大的顺序存入bm,将对应的已搜索层保留分支和当前层保留的星座点写入Xleft
如图6所示,为已搜索传输层的16个保留分支,c1到c4是通过查表得到的与保留分支的判决量的欧氏距离最近的4个星座点,…,c13到c16是通过查表得到的与保留分支的判决量的欧氏距离最近的4个星座点;c17是通过星座符号的应判决方法得到的与保留分支的判决量的欧氏距离最近的一个星座点,…,c28是通过星座符号的应判决方法得到的与保留分支的判决量的欧氏距离最近的一个星座点;是对当前传输层进行符号搜索后得到的保留分支。
最后,所有符号搜索完成之后输出的bm为保留的最小的16个分支对应的分支度量值,保留分支Xleft包含对应的bm最小的16个分支。
假设,采用QRDM-List算法的情况下NC1=16,NC2=4,采用QRDM-List检测算法情况下NC1=16,NC2=4,NC3=4,在64QAM的调制方式下,给出两种算法的复杂度对比。下面主要在分支度量值的计算和排序两方面进行分析。
分支度量值计算量的比较:
对于非第一层的搜索,QRDM-List算法需要进行4×16=64个分支度量值的计算,QRDM-List&HD算法需要进行4×4+12=28个分支度量值的计算。后者的运算量缩减了超过50%。
排序的比较:
对于非第一层的搜索,QRDM-List算法需要进行64选16的排序。QRDM-List&HD算法需要进行28选16的排序。排序实现复杂度和时延得到了缩减。
参见图7,本发明实施例提供一种信号检测装置,该装置包括:
符号搜索单元70,用于在采用基于QR分解和M算法的查找表最大似然检测QRDM-List算法对接收信号进行信号检测的过程中,对多个传输层逐层进行符号搜索;
所述符号搜索单元70包括判决量确定单元71、判决量选取单元72、第一度量值增量确定单元73、第二度量值增量确定单元74和结果确定单元75;
所述判决量确定单元71,用于对除最先搜索层之外的其他传输层进行符号搜索时,确定已搜索传输层的每个保留星座符号组合分别对应的判决量;
所述判决量选取单元72,用于选取确定的各判决量中的部分判决量;
所述第一度量值增量确定单元73,用于对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;
所述第二度量值增量确定单元74,用于对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
所述结果确定单元75,用于根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
进一步的,所述判决量选取单元72用于:
选取确定的各判决量中累积分支度量值最小的设定数目个判决量,该设定数目小于判决量的总个数。
进一步的,所述判决量确定单元71用于:
按照如下公式确定已搜索传输层的第nc1个保留星座符号组合对应的判决量nc1∈{1,…,NC1},NC1为保留星座符号组合的个数:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;ri,i为所述R矩阵中第i行第i列的元素。
进一步的,所述结果确定单元75包括:
累积度量值确定单元76,用于对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值,该累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留星座符号组合的分支度量值之和;
结果保存单元77,用于得到的各累积分支度量值中选取NC1个最小的累积分支度量值,将选取的NC1个最小的累积分支度量值对应的星座点组合作为本传输层的保留星座符号组合并保存到Xleft中,将选取的NC1个最小的累积分支度量值作为各保留星座符号组合分别对应的分支度量值并保存到bm中,其中Xleft为NL行NC1列的矩阵;bm为长度为NC1的向量;NL为传输层的总层数,NC1为保留星座符号组合的个数;将当前的Xleft和bm作为本传输层的符号搜索结果。
进一步的,所述第一度量值增量确定单元73,用于按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量nc1∈{1,…,NC3},1≤NC3≤NC1,NC1为保留星座符号组合的个数,nc2∈{1,…,NC2},NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数:
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,..,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的第nc2个星座点;ri,i为所述R矩阵中第i行第i列的元素。
进一步的,所述累积度量值确定单元76用于:
在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是所述第一度量值增量确定单元确定的,则按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp((nc1-1)×NC2+nc2):
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
进一步的,所述第二增量确定单元74用于:
按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量nc1∈{NC3+1,…,NC1},1≤NC3≤NC1,NC1为保留星座符号组合的个数,NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数,NC3为选取的部分判决量的个数:
d i , NC 3 &times; NC 2 + nc 1 - NC 3 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,..,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的一个星座点;ri,i为所述R矩阵中第i行第i列的元素。
进一步的,所述累积度量值确定单元76用于:
在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是所述第二度量值增量确定单元确定的,则按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp(NC3×NC2+nc1-NC3):
BM tmp ( NC 3 &times; NC 2 + nc 1 - NC 3 ) = bm ( nc 1 ) + d i , NC 3 &times; NC 2 + nc 1 - NC 3 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
进一步的,每个传输层的保留星座符号组合的个数为16,通过查表确定的星座点的个数为4,选取的部分判决量的个数为4。
综上,本发明的有益效果包括:
本发明实施例提供的方案中,在采用QRDM-List算法对接收信号进行信号检测的过程中,在对除最先搜索层之外的其他传输层进行符号搜索时,只需对部分判决量通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量,对于其他判决量,仅需采用星座符号的硬判决方法确定与其欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量,与现有技术中需要对所有判决量都通过查表确定与其欧氏距离最近的多个星座点,并确定每个星座点对应的分支度量值的增量相比,本发明的信号检测复杂度较低,进而提高了信号检测效率。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (18)

1.一种信号检测方法,其特征在于,该方法包括:
在采用基于QR分解和M算法的查找表最大似然检测QRDM-List算法对接收信号进行信号检测的过程中,对多个传输层逐层进行符号搜索,并且,按照如下方法对除最先搜索层之外的其他传输层进行符号搜索:
确定已搜索传输层的每个保留星座符号组合对应的判决量;
选取确定的各判决量中的部分判决量,对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;
对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
2.如权利要求1所述的方法,其特征在于,所述选取确定的各判决量中的部分判决量,具体包括:
选取确定的各判决量中累积分支度量值最小的设定数目个判决量,该设定数目小于判决量的总个数。
3.如权利要求1所述的方法,其特征在于,所述确定已搜索传输层的每个保留星座符号组合对应的判决量,具体包括:
按照如下公式确定已搜索传输层的第nc1个保留星座符号组合对应的判决量nc1∈{1,…,NC1},NC1为保留星座符号组合的个数:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;ri,i为所述R矩阵中第i行第i列的元素。
4.如权利要求1所述的方法,其特征在于,所述根据确定的每个分支度量值的增量得到本传输层的符号搜索结果,具体包括:
对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值,该累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留星座符号组合的分支度量值之和;
从得到的各累积分支度量值中选取NC1个最小的累积分支度量值,将选取的NC1个最小的累积分支度量值对应的星座点组合作为本传输层的保留星座符号组合并保存到Xleft中,将选取的NC1个最小的累积分支度量值作为各保留星座符号组合分别对应的分支度量值并保存到bm中,其中Xleft为NL行NC1列的矩阵;bm为长度为NC1的向量;NL为传输层的总层数,NC1为保留星座符号组合的个数;将当前的Xleft和bm作为本传输层的符号搜索结果。
5.如权利要求4所述的方法,其特征在于,在确定星座点对应的分支度量值的增量时,若该星座点是通过查表确定的,则确定该星座点对应的分支度量值的增量的方法包括:
按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量nc1∈{1,…,NC3},1≤NC3≤NC1,NC1为保留星座符号组合的个数,nc2∈{1,…,NC2},NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数:
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的第nc2个星座点;ri,i为所述R矩阵中第i行第i列的元素。
6.如权利要求5所述的方法,其特征在于,在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是通过查表确定的星座点对应的分支度量值的增量,则按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp((nc1-1)×NC2+nc2):
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
7.如权利要求4所述的方法,其特征在于,在确定星座点对应的分支度量值的增量时,若该星座点是通过星座符号硬判决方法确定的,则确定该星座点对应的分支度量值的增量的方法包括:
按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量nc1∈{NC3+1,…,NC1},1≤NC3≤NC1,NC1为保留星座符号组合的个数,NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数,NC3为选取的部分判决量的个数:
d i , NC 3 &times; NC 2 + nc 1 - NC 3 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的一个星座点;ri,i为所述R矩阵中第i行第i列的元素。
8.如权利要求7所述的方法,其特征在于,在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是通过星座符号硬判决方法确定的星座点对应的分支度量值的增量,则按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp(NC3×NC2+nc1-NC3):
BM tmp ( NC 3 &times; NC 2 + nc 1 - NC 3 ) = bm ( nc 1 ) + d i , NC 3 &times; NC 2 + nc 1 - NC 3 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
9.如权利要1-8中任一所述的方法,其特征在于,每个传输层的保留星座符号组合的个数为16,通过查表确定的星座点的个数为4,选取的部分判决量的个数为4。
10.一种信号检测装置,其特征在于,该装置包括:
符号搜索单元,用于在采用基于QR分解和M算法的查找表最大似然检测QRDM-List算法对接收信号进行信号检测的过程中,对多个传输层逐层进行符号搜索;
所述符号搜索单元包括判决量确定单元、判决量选取单元、第一度量值增量确定单元、第二度量值增量确定单元和结果确定单元;
所述判决量确定单元,用于对除最先搜索层之外的其他传输层进行符号搜索时,确定已搜索传输层的每个保留星座符号组合分别对应的判决量;
所述判决量选取单元,用于选取确定的各判决量中的部分判决量;
所述第一度量值增量确定单元,用于对于选取的每个判决量,通过查表确定与该判决量的欧氏距离最近的多个星座点,并确定该多个星座点中每个星座点对应的分支度量值的增量;
所述第二度量值增量确定单元,用于对于未被选取的每个判决量,采用星座符号的硬判决方法确定与该判决量的欧氏距离最近的一个星座点,并确定该星座点对应的分支度量值的增量;
所述结果确定单元,用于根据确定的每个分支度量值的增量得到本传输层的符号搜索结果。
11.如权利要求10所述的装置,其特征在于,所述判决量选取单元用于:
选取确定的各判决量中累积分支度量值最小的设定数目个判决量,该设定数目小于判决量的总个数。
12.如权利要求10所述的装置,其特征在于,所述判决量确定单元用于:
按照如下公式确定已搜索传输层的第nc1个保留星座符号组合对应的判决量nc1∈{1,…,NC1},NC1为保留星座符号组合的个数:
x ~ i , nc 1 = ( y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left ) / r i , i ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;ri,i为所述R矩阵中第i行第i列的元素。
13.如权利要求10所述的装置,其特征在于,所述结果确定单元包括:
累积度量值确定单元,用于对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值,该累积分支度量值是该分支度量值的增量与该分支度量值的增量对应的已搜索传输层的保留星座符号组合的分支度量值之和;
结果保存单元,用于得到的各累积分支度量值中选取NC1个最小的累积分支度量值,将选取的NC1个最小的累积分支度量值对应的星座点组合作为本传输层的保留星座符号组合并保存到Xleft中,将选取的NC1个最小的累积分支度量值作为各保留星座符号组合分别对应的分支度量值并保存到bm中,其中Xleft为NL行NC1列的矩阵;bm为长度为NC1的向量;NL为传输层的总层数,NC1为保留星座符号组合的个数,将当前的Xleft和bm作为本传输层的符号搜索结果。
14.如权利要求13所述的装置,其特征在于,所述第一度量值增量确定单元,用于按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量nc1∈{1,…,NC3},1≤NC3≤NC1,NC1为保留星座符号组合的个数,nc2∈{1,…,NC2},NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数:
d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , nc 2 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的第nc2个星座点;ri,i为所述R矩阵中第i行第i列的元素。
15.如权利要求14所述的装置,其特征在于,所述累积度量值确定单元用于:
在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是所述第一度量值增量确定单元确定的,则按照如下公式确定由第nc1个判决量得到的第nc2个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp((nc1-1)×NC2+nc2):
BM tmp ( ( nc 1 - 1 ) &times; NC 2 + nc 2 ) = bm ( nc 1 ) + d i , ( nc 1 - 1 ) &times; NC 2 + nc 2 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
16.如权利要求13所述的装置,其特征在于,所述第二增量确定单元用于:
按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量nc1∈{NC3+1,…,NC1},1≤NC3≤NC1,NC1为保留星座符号组合的个数,NC2为通过查表确定的与判决量的欧氏距离最近的星座点的个数,NC3为选取的部分判决量的个数:
d i , NC 3 &times; NC 2 + nc 1 - NC 3 = | y i - &Sigma; j = i + 1 N L r i , j X j , nc 1 left - r i , i x ^ i , nc 1 , 1 | 2 ;
其中,yi为等效接收信号向量中的第i个元素,从第NL-1层到第1层,i的取值分别为2,...,NL;NL为传输层的总层数;ri,j为预先对接收信号向量进行QR分解得到的R矩阵中第i行第j列的元素;Xleft为NL行NC1列的矩阵,Xleft初始为空矩阵,并在每次对传输层进行符号搜索后将得到的保留星座符号组合保存在Xleft中;为Xleft中第j行第nc1列的元素;为由第nc1个判决量得到的一个星座点;ri,i为所述R矩阵中第i行第i列的元素。
17.如权利要求16所述的装置,其特征在于,所述累积度量值确定单元用于:
在对于得到的每个分支度量值的增量,确定该分支度量值的增量对应的累积分支度量值时,若该分支度量值的增量是所述第二度量值增量确定单元确定的,则按照如下公式确定由第nc1个判决量得到的一个星座点对应的分支度量值的增量所对应的累积分支度量值BMtmp(NC3×NC2+nc1-NC3):
BM tmp ( NC 3 &times; NC 2 + nc 1 - NC 3 ) = bm ( nc 1 ) + d i , NC 3 &times; NC 2 + nc 1 - NC 3 ;
其中,bm(nc1)为向量bm中的第nc1个元素。
18.如权利要10-17中任一所述的装置,其特征在于,每个传输层的保留星座符号组合的个数为16,通过查表确定的星座点的个数为4,选取的部分判决量的个数为4。
CN201310031384.9A 2013-01-28 2013-01-28 信号检测方法和装置 Active CN103973602B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310031384.9A CN103973602B (zh) 2013-01-28 2013-01-28 信号检测方法和装置
PCT/CN2014/071178 WO2014114246A1 (zh) 2013-01-28 2014-01-23 信号检测方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310031384.9A CN103973602B (zh) 2013-01-28 2013-01-28 信号检测方法和装置

Publications (2)

Publication Number Publication Date
CN103973602A true CN103973602A (zh) 2014-08-06
CN103973602B CN103973602B (zh) 2017-07-21

Family

ID=51226937

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310031384.9A Active CN103973602B (zh) 2013-01-28 2013-01-28 信号检测方法和装置

Country Status (2)

Country Link
CN (1) CN103973602B (zh)
WO (1) WO2014114246A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106561008A (zh) * 2015-09-30 2017-04-12 北京展讯高科通信技术有限公司 一种降低信号接收复杂度的方法及接收机以及移动终端
CN107027132A (zh) * 2016-02-01 2017-08-08 电信科学技术研究院 一种信号检测方法及装置
WO2018086405A1 (zh) * 2016-11-10 2018-05-17 电信科学技术研究院 一种调制方式检测方法和装置

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2548876B (en) * 2016-03-31 2019-12-25 Toshiba Res Europe Limited Low-complexity LLR computation for nonuniform QAM constellations
TWI726347B (zh) * 2019-07-03 2021-05-01 國立陽明交通大學 無線通訊裝置以及訊號偵測方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080279299A1 (en) * 2007-05-10 2008-11-13 Comsys Communication & Signal Processing Ltd. Multiple-input multiple-output (mimo) detector incorporating efficient signal point search and soft information refinement
CN101542994A (zh) * 2006-11-24 2009-09-23 Nxp股份有限公司 在多天线系统的接收机中用来产生软比特信息的方法和装置
CN101582750A (zh) * 2009-06-02 2009-11-18 北京天碁科技有限公司 一种基于宽度优先搜索的球形译码检测方法
US20120106684A1 (en) * 2010-10-29 2012-05-03 Intel Mobile Communications Technology Dresden GmbH Method and arrangement for soft-decision sphere decoding
CN102723975A (zh) * 2012-06-11 2012-10-10 北京大学 Mimo系统的信号检测方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101188478A (zh) * 2006-11-17 2008-05-28 中兴通讯股份有限公司 分层空时码系统的信号检测方法
CN101888287B (zh) * 2009-05-14 2014-12-03 株式会社Ntt都科摩 一种多输入多输出接收机信号检测方法及装置
US8488721B2 (en) * 2009-08-20 2013-07-16 Electronics And Telecommunications Research Institute Adaptive QRD-M algorithm based signal detecting method by using constellation set grouping in spatial multiplexing multiple-input multiple-output system
CN101662341B (zh) * 2009-09-25 2012-08-01 意法·爱立信半导体(北京)有限公司 一种多输入多输出信号检测方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101542994A (zh) * 2006-11-24 2009-09-23 Nxp股份有限公司 在多天线系统的接收机中用来产生软比特信息的方法和装置
US20080279299A1 (en) * 2007-05-10 2008-11-13 Comsys Communication & Signal Processing Ltd. Multiple-input multiple-output (mimo) detector incorporating efficient signal point search and soft information refinement
CN101582750A (zh) * 2009-06-02 2009-11-18 北京天碁科技有限公司 一种基于宽度优先搜索的球形译码检测方法
US20120106684A1 (en) * 2010-10-29 2012-05-03 Intel Mobile Communications Technology Dresden GmbH Method and arrangement for soft-decision sphere decoding
CN102723975A (zh) * 2012-06-11 2012-10-10 北京大学 Mimo系统的信号检测方法及装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106561008A (zh) * 2015-09-30 2017-04-12 北京展讯高科通信技术有限公司 一种降低信号接收复杂度的方法及接收机以及移动终端
CN106561008B (zh) * 2015-09-30 2020-04-24 北京展讯高科通信技术有限公司 一种降低信号接收复杂度的方法及接收机以及移动终端
CN107027132A (zh) * 2016-02-01 2017-08-08 电信科学技术研究院 一种信号检测方法及装置
WO2017133408A1 (zh) * 2016-02-01 2017-08-10 电信科学技术研究院 一种信号检测方法及装置
CN107027132B (zh) * 2016-02-01 2019-05-21 电信科学技术研究院 一种信号检测方法及装置
WO2018086405A1 (zh) * 2016-11-10 2018-05-17 电信科学技术研究院 一种调制方式检测方法和装置
TWI650984B (zh) * 2016-11-10 2019-02-11 電信科學技術研究院 一種調製方式檢測方法和裝置

Also Published As

Publication number Publication date
WO2014114246A1 (zh) 2014-07-31
CN103973602B (zh) 2017-07-21

Similar Documents

Publication Publication Date Title
JP4429945B2 (ja) Mimo多重通信装置および信号分離方法
EP2520055B1 (en) Branch processing of search tree in a sphere decoder
CN103973602A (zh) 信号检测方法和装置
Roger et al. Fully parallel GPU implementation of a fixed-complexity soft-output MIMO detector
US20170163400A1 (en) Data processing method and apparatus for maximum likelihood ml receiver
CN106357312B (zh) 格约减辅助广度优先树搜索mimo检测方法
EP3304836B1 (en) Parallel processing of sphere decoders and other vector finding approaches using tree search
CN104901910A (zh) 一种mimo系统的检测方法及装置
CN106656406A (zh) 一种非正交多址接入中信号检测方法及装置
CN101562464B (zh) 一种基于深度优先搜索的球形译码检测方法
TWI650984B (zh) 一種調製方式檢測方法和裝置
CN101557269B (zh) 一种基于超大规模集成电路的球形译码检测方法
CN103532609A (zh) 基于K-Best算法的MIMO信号检测方法及检测器
EP2882153B1 (en) Mimo maximum likelihood detector (mld) accelerator
CN103384180B (zh) 用于多输入多输出系统的空时译码方法和装置
CN103684565B (zh) 一种确定软比特信息的方法及装置
CN103427958A (zh) 信号检测方法和装置
CN109660473B (zh) 一种球形译码检测方法及装置、计算机可读存储介质
CN104735003B (zh) 欧式距离计算方法、模块和多输入多输出译码装置
CN103188003B (zh) 基于迭代检测的低复杂度并行干扰消除方法及系统
CN104639495A (zh) 星座软解映射的方法和装置
CN103873203B (zh) 一种信号检测的方法及装置
CN107027132B (zh) 一种信号检测方法及装置
CN102882815B (zh) 一种多输入多输出数据检测方法及装置
CN107104714A (zh) 一种无需qr分解的mimo检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210601

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY