具体实施方式
下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
LTE系统的物理层实现是基于MIMO-OFDM(Multiple Input MultipleOutput-Orthogonal Frequency Division Multiplexing,多发多收-正交频分复用)技术的,MIMO-OFDM技术作为无线环境下的高速传输技术,主要是将一定频谱宽度内的信道划分为多个正交子信道,将频率选择性衰落信道改良为平坦衰落,MIMO-OFDM技术则可以有效地实现系统容量的大幅扩充。
在包括两个发射天线和两个接收天线的系统中,在具有丰富散射路径的条件下,定义x=(x1,x2)T是发送信号向量,通常情况下该向量是复数形式的(BPSK(Binary Phase Shift Keying,移相键控)调制的信号除外),y=(y1,y2)T为复数形式的接收向量,n=(n1,n2)T是复数形式的噪声向量,则发送信号和接收信号之间的关系可以表示为:y=Hx+n,即每一根天线上接收的都是发送端的两路信号经过相应的信道衰减后与噪声相叠加的信号集合。
基于上述情况,本发明实施例中,通过假定两个发送符号中的一个,并利用V-BLAST(垂直分层空时结构)算法求得另一个,然后进行反向检测以验证假设的正确性,由此约束最大似然搜索的范围,从而在保证解调性能接近最大似然检测算法的情况下,降低检测的复杂度与所用时间,并显著提高检测性能,且提高实际系统中的可实现性,解决现有算法检测效率低的问题。
实施例一
本发明实施例一提出一种信号检测方法,应用于包括两个发射天线和两个接收天线的系统(如LTE的2发2收系统)中,并用于空间复用的场景下,如图1所示,该信号检测方法可以包括以下步骤:
步骤101,选择调制星座图中的各星座点分别作为发射天线1的发送符号x1,并利用V-BLAST算法计算发射天线2上x1对应的发送符号x2。
本发明实施例中,选择调制星座图中的各星座点分别作为发射天线1的发送符号x1,并利用V-BLAST算法计算发射天线2上x1对应的发送符号x2,具体包括但不限于:
步骤1、根据系统调制方式(如LTE的2发2收系统的调制方式)获取调制星座图,该调制星座图中包含|Ω|个星座点,各星座点为复值符号,且|Ω|对应于系统调制阶数(如LTE的2发2收系统的调制阶数)。
例如,在16QAM(Quadrature Amplitude Modulation,正交幅度调制)调制的系统中,调制星座图中有16个调制星座点,其分别为:
需要说明的是,上述QAM调制是一种向量调制,将输入的数据比特先映射(一般采用格雷码)到一个复平面(星座)上,形成复值调制符号,然后将符号的I、Q分量(对应复平面的实部和虚部,也就是水平和垂直方向)分别对应调制在相互正交的两个载波上;例如,16QAM的调制信号具有16个样点,每个样点表示一个复值符号,即为上述16个调制星座点。
步骤2、选择各星座点分别作为发射天线1的发送符号x1,并利用接收信号向量y和信道估计得到的信道矩阵H=(h1,h2)进行V-BLAST计算,得到|Ω|个x1对应的x2;其中,h1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
本发明实施例中,通过信道估计算法可以得到信道矩阵,且将信道矩阵分为两个列向量h1,h2,即信道矩阵H=(h1,h2);其中,信道估计算法包括但不限于现有的最小二乘法、最小均方误差、变换域等算法,对此不作限定。
本发明实施例中,利用接收信号向量y和信道估计得到的信道矩阵H=(h1,h2)进行V-BLAST计算,得到|Ω|个x1对应的x2,具体包括但不限于:
利用如下公式计算得到|Ω|个x1对应的x2:
x
2=Q(h
2 +(y-h
1x
1))=Q(h
2 Hy-h
2 Hh
1x
1)/||h
2||
2);其中,Q为限制函数,且
a为调制星座图中包含的星座点。
具体的V-BLAST算法为:初始化:i=1,G
1=H
+,
循环部分:
i=i+1;其中,(·)
+表示矩阵的M-P逆,<G
i>
j表示矩阵G
i的第j行,同理[H
i]
j表示矩阵H
i的第j列,Q(·)表示限制函数,
重复执行循环部分直到i超过矩阵列数。
本实施例为2发2收的LTE系统,信道估计得到的信道矩阵H=(h1,h2),由于选取各星座点分别作为发射天线1的发送符号x1,即将x1作为已知,使用V-BLAST算法计算x2过程为:x2=Q(h2 +(y-h1x1))=Q(h2 Hy-h2 Hh1x1)/||h2|2),从而可以基于上述公式得到|Ω|个x1及其对应的x2。
步骤102,利用V-BLAST算法计算x2对应的发送符号x1’。
本发明实施例中,利用V-BLAST算法计算x2对应的发送符号x1’,具体包括但不限于:利用如下公式计算得到|Ω|个x2对应的x1’:
x1′=Q(h1 +(y-h2x2))=Q(h1 Hy-h1 Hh2x2)/||h1||2);
其中,调制星座图中包含|Ω|个星座点,且|Ω|对应于系统调制阶数;Q为限制函数,且
a为调制星座图中包含的星座点;h
1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h
2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
步骤103,如果x1与x1’相等,则将(x1,x2)T作为候选发送向量。
在具体实现过程中,通过比较x1与x1’,如果x1与x1’相等,则将(x1,x2)T作为候选发送向量;如果x1与x1’不等,则舍弃相应信息。
本发明实施例中,将(x1,x2)T作为候选发送向量的过程中:通过比较|Ω|组x1及其对应的x1’,如果x1与x1’相等,则说明发射天线1发送该组的复值符号x1时,发射天线2上最有可能发送该组的复值符号x2,且反之亦然,因此,该组两发送符号组成的发送向量(x1,x2)T经过信道后的响应与接收向量的欧氏距离较近,可以将(x1,x2)T作为候选发送向量进行最大似然搜索;如果x1与x1’不等,则说明若发射天线2发送x2,发射天线1发送的并不是相应的x1,假设不成立,说明发送向量(x1,x2)T经过信道后的响应与接收向量的欧氏距离较远,故舍去。
步骤104,利用各星座点计算得到的候选发送向量确定候选集合S,并在候选集合S中进行最大似然信号检测。
本发明实施例中,在候选集合S中进行最大似然信号检测具体包括但不限于:将候选集合S中的各候选发送向量分别带入最大似然检测算法公式
中的x,以分别得到各候选发送向量经过信道矩阵H后的响应与接收信号向量y的欧式距离,并将各欧式距离中最小欧式距离对应的信号向量确定为最大似然信号向量,且最大似然信号向量为最可能的发送信号向量;即在候选集合S中搜索满足条件
的信号向量x
ML,实现最大似然检测;其中,H=(h
1,h
2),且h
1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h
2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
综上所述,本发明实施例中,可以在接收端进行并行检测,在保证解调性能接近最大似然检测算法的情况下,降低检测的复杂度与所用时间,并能够显著提高接收机的检测性能,且提高在实际系统中的可实现性。
如图2所示,为本发明实施例在2发2收的LTE系统下的仿真结果,图中的横坐标表示SNR(信噪比),纵坐标表示BER(误比特率);仿真条件为:信号调制方式为16QAM、信道模型为TS 36.211所定义的信道,信道估计为理想信道估计,涉及的检测算法包括ML检测和本实施例中改进ML检测。
由图2中可以看出,改进算法的检测性能与ML检测相差不大;在复杂度方面,采用最大似然检测算法需要遍历|Ω|2=162=256次,而仿真结果同时表明,使用本发明实施例的信号检测方法,最大似然检测算法平均只需要遍历2.5次,大大减小了搜索的范围,减少了遍历次数;从算法需要的乘加运算次数来看,ML检测至少需要256×4=1024次乘加运算,而本发明实施例提供的算法只需要50次乘加运算,可见本发明实施例大大提高了信号检测的效率,有效地解决了现有的信号检测方法效率低的问题。
实施例二
基于与上述方法同样的发明构思,本发明实施例还提出了一种信号检测设备,应用于包括两个发射天线和两个接收天线的系统中,如图3所示,该信号检测设备包括:
计算模块11,用于选择调制星座图中的各星座点分别作为发射天线1的发送符号x1,并利用垂直分层空时结构V-BLAST算法计算发射天线2上x1对应的发送符号x2;以及,利用V-BLAST算法计算x2对应的发送符号x1’;
确定模块12,用于当所述x1与x1’相等时,将(x1,x2)T作为候选发送向量;
检测模块13,用于利用所述各星座点计算得到的候选发送向量确定候选集合S,并在所述候选集合S中进行最大似然信号检测。
所述计算模块11,具体用于根据系统调制方式获取调制星座图,所述调制星座图中包含|Ω|个星座点,各星座点为复值符号,且|Ω|对应于系统调制阶数;
选择各星座点分别作为发射天线1的发送符号x1,并利用接收信号向量y和信道估计得到的信道矩阵H=(h1,h2)进行V-BLAST计算,得到|Ω|个x1对应的x2;
其中,h1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
所述计算模块11,进一步用于利用如下公式计算得到|Ω|个x1对应的x2:
x2=Q(h2 +(y-h1x1))=Q(h2 Hy-h2 Hh1x1)/||h2||2);其中,Q为限制函数,且a为所述调制星座图中包含的星座点。
所述计算模块11,具体用于利用如下公式计算得到|Ω|个x2对应的x1’:
x1′=Q(h1 +(y-h2x2))=Q(h1 Hy-h1 Hh2x2)/||h1||2);
其中,所述调制星座图中包含|Ω|个星座点,且|Ω|对应于系统调制阶数;Q为限制函数,且 a为所述调制星座图中包含的星座点;h1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
所述确定模块12,具体用于通过比较所述x1与所述x1’,如果所述x1与所述x1’相等,则确定发射天线1发送该组的复值符号x1时,发射天线2上最有可能发送该组的复值符号x2,且发送向量(x1,x2)T经过信道后的响应与接收向量的欧氏距离较近,并将(x1,x2)T作为候选发送向量。
所述检测模块13,具体用于将所述候选集合S中的各候选发送向量分别带入最大似然检测算法公式
中的x,以分别得到各候选发送向量经过信道矩阵H后的响应与接收信号向量y的欧式距离,并将各欧式距离中最小欧式距离对应的信号向量确定为最大似然信号向量,且所述最大似然信号向量为最可能的发送信号向量;
其中,H=(h1,h2),且h1为发射天线1与接收天线1间的信道冲击响应和发射天线1与接收天线2间的信道冲击响应组成的二维列向量,h2为发射天线2与接收天线1间的信道冲击响应和发射天线2与接收天线2间的信道冲击响应组成的二维列向量。
其中,本发明装置的各个模块可以集成于一体,也可以分离部署。上述模块可以合并为一个模块,也可以进一步拆分成多个子模块。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明可以通过硬件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
上述本发明序号仅仅为了描述,不代表实施例的优劣。
以上公开的仅为本发明的几个具体实施例,但是,本发明并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围。