一种近似MMSE性能的高效大规模MIMO检测方法及系统技术方案

本发明专利技术公开了一种近似MMSE性能的高效大规模MIMO检测方法及系统，通过引入预条件技术，本发明专利技术能够显著加快传统GS方法的迭代速率，从而使得本发明专利技术提出的大规模MIMO检测算法在恶劣传播环境(如发射/接收天线数相近或者空间相关性较大的信道)中依然能够快速逼近精确MMSE检测算法的性能。数值模拟结果表明，本发明专利技术提出的大规模MIMO检测算法在恶劣传播环境中表现出的误码率性能要优于基于Neumann级数、GS方法、CG方法的传统大规模MIMO检测算法。另一方面，本发明专利技术提供的系统创新性地发掘了GS迭代在元素更新过程中的循环移位特性，从而使得其能够以较低的硬件消耗和延迟进行GS迭代操作。

【技术实现步骤摘要】
一种近似MMSE性能的高效大规模MIMO检测方法及系统
本专利技术属于计算机通信和数字电路领域，涉及一种近似MMSE性能的高效大规模MIMO检测方法及系统。
技术介绍
大规模多输入多输出(MIMO)被公认为是第5代(5G)无线通信系统的重要技术之一。该技术通过在基站和用户端配备大量天线来提供更高的频谱效率，更快的峰值数据速率以及比传统小规模MIMO系统更好的能量效率。然而，伴随着天线数目的大量增加，在大规模MIMO系统中基带算法的复杂度也在急剧增加。其中，上行链路的最佳多用户检测方法，例如最大似然(ML)检测和最大后验(MAP)方法在计算复杂度方面将变得难以承受(由于它们的指数复杂度)。因此，更加可行和高效的检测器设计吸引了大量关注。近年来，研究人员将目光转向了线性检测方法，如传统的迫零(ZF)和最小均方误差(MMSE)，这是因为它们在大规模MIMO系统中有着次优检测性能和低复杂性的特性。值得注意的是，对于大规模MIMO系统中的MMSE检测方法，其主要的计算复杂度在于一个高阶矩阵的求逆。假设M为单天线用户数，若采用精确的矩阵求逆方法，如Cholesky分解法，则计算复杂度为O(M3)。这意味着如果M的数量极大时，精确的MMSE检测将需要巨大的计算量和硬件消耗。近几年，有国内外研究人员相继提出了基于Gauss-Seidel(GS)方法、Neumann级数、共轭梯度(CG)的大规模MIMO检测方法，获得了接近MMSE算法的性能表现。这些方法的共同点在于都是传统的迭代数值计算方法，虽然在一定程度上减小了计算复杂度，但是对于恶劣的传播环境(如发射/接收天线数相近或者空间相关性较大的信道)，它们的性能表现将下降甚至无法收敛。
技术实现思路
专利技术目的：针对现有技术的不足，本专利技术提出了一种近似MMSE性能的高效大规模MIMO检测方法及系统。技术方案：一种近似MMSE性能的高效大规模MIMO检测方法，包括以下步骤：步骤1：预处理；将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中Gram矩阵G＝HHH，No为噪声方差，IM为M维单位矩阵，(.)H为共轭转置操作；步骤2：计算标准化矩阵和标准化向量其中D为W的对角元素矩阵，使得系数矩阵对角线元素为1；步骤3：预条件；构造预条件矩阵P＝S+IM，计算系数矩阵和常数向量其中S为一个和有关的矩阵：步骤4：根据步骤3输出的系数矩阵和常数向量设置迭代初始解为x(0)＝0，并开始迭代操作，输出检测结果；算法伪代码如下：K轮迭代后，x(K)即为待检测信号的估计结果。本专利技术还提供一种近似MMSE性能的高效大规模MIMO检测系统，包括：预条件模块，用于完成预处理，将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中Gram矩阵G＝HHH，No为噪声方差，IM为M维单位矩阵，(.)H为共轭转置操作；然后计算标准化矩阵和标准化向量其中D为W的对角元素矩阵，使得系数矩阵对角线元素为1；最后构造预条件矩阵P＝S+IM，计算系数矩阵和常数向量其中S为一个和有关的矩阵；GS迭代模块，用于完成根据预条件模块输出的系数矩阵和常数向量设置迭代初始解为x(0)＝0，并进行迭代操作，输出检索结果，算法伪代码如下：K轮迭代后，x(K)即为待检测信号的估计结果。进一步的，所述的预条件模块包括6个脉动阵列构成的矩阵乘法器、2个加法器阵列以及1个求倒数单元；其中，用2个脉动阵列计算匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中脉动阵列的处理单元为基本的复数乘法累加器；用另2个脉动阵列计算标准化矩阵和标准化向量其中D-1由求倒数单元计算得到，求倒数单元由查找表生成，脉动阵列的处理单元依然为基本的复数乘法累加器；用剩余的2个脉动阵列计算系数矩阵和常数向量其中预条件矩阵P＝S+IM。进一步的，所述的GS迭代模块包含M-1个复数乘法器、加法器和寄存器，其进行每一轮GS迭代需要M个时钟周期，M为发射天线数。工作原理：考虑到大规模MIMO系统上行链路MMSE检测中滤波矩阵W为Hermitian正定阵且主对角线占优，本专利技术采用的GS迭代方法在多次迭代后一定收敛。在恶劣传播环境(如发射/接收天线数相近或者空间相关性较大的信道)中，本专利技术采用的预条件方法一定能减小迭代矩阵的谱半径，从而达到加速收敛的效果。有益效果：与现有技术相比，本专利技术重点考虑在恶劣传播环境(如发射/接收天线数相近或者空间相关性较大的信道)中如何以较低的计算复杂度达到近似MMSE性能的检测效果。通过采用预条件处理的GS迭代方法，本专利技术能够在同样迭代次数的条件下取得比传统基于Neumann级数、GS方法、CG方法的大规模MIMO检测算法更好的误码率性能，尤其是在恶劣传播环境(如发射/接收天线数相近或者空间相关性较大的信道)下，并且在较少的迭代次数之后得到近似MMSE性能的检测效果。另一方面，本专利技术提供的系统创新性地发掘了GS迭代在元素更新过程中的循环移位特性，从而使得其能够以较低的硬件消耗和延迟进行GS迭代操作。此外，这一特性还使得对应的控制电路的设计变得十分容易。附图说明图1为采用本专利技术信号检测方法和其他传统检测方法的误码率对比图(发射天线数为32，接收天线数为128，相关系数为0时)；图2为采用本专利技术信号检测方法和其他传统检测方法的误码率对比图(发射天线数为32，接收天线数为128，相关系数为0.3时)；图3为采用本专利技术信号检测方法和其他传统检测方法的误码率对比图(发射天线数为16，接收天线数为128，相关系数为0.3时)；图4为采用本专利技术信号检测算法和其他传统检测算法的误码率对比图(发射天线数为8，接收天线数为128，相关系数为0.3时)；图5本专利技术系统示意图；图6为本专利技术系统中GS迭代模块示意图；图7为本专利技术系统中GS迭代模块时序调度示意图(系统发射天线数为4时)。具体实施方式下面将结合附图，对本专利技术的实施案例进行详细的描述；本实施例中建立一个大规模MIMO上行链路系统进行模拟操作。在大规模MIMO上行链路中，一般有N＞＞M(基站天线数N远大于发射天线数，即用户数M)。首先M个不同用户产生的并行传输比特流分别通过信道编码进行编码，然后映射到星座符号，并采取星座图集合能量归一化。让x＝[x1，x2，x3，...，xM]T表示信号向量，x中包含了分别从M个用户产生的传输符号，采用64-QAM方式映射。H表示维度是N×M信道矩阵，故上行链路基站端的接收信号向量y可以表示为y＝Hx+n其中y的维度为N×1，n为N×1维的加性白噪声向量，其元素服从零均值方差为No的高斯分布。上行链路多用户信号检测任务就是从接收机接收向量y＝[y1，y2，y3，...，yN]T估计传输信号符号x。假设H已知，采用最小均方误差(MMSE)线性检测理论，对传输信号向量的估计表示为该估计过程等效为求解线性方程组基于上述模型，本专利技术实施例公开一种近似MMSE性能的高效大规模MIMO检测方法，包括如下步骤：步骤1：预处理，将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中Gram矩阵本文档来自技高网...

【技术保护点】
一种近似MMSE性能的高效大规模MIMO检测方法，其特征在于，包括以下步骤：步骤1：预处理；将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出y

【技术特征摘要】
1.一种近似MMSE性能的高效大规模MIMO检测方法，其特征在于，包括以下步骤：步骤1：预处理；将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中Gram矩阵G＝HHH，No为噪声方差，IM为M维单位矩阵，(.)H为共轭转置操作；步骤2：计算标准化矩阵和标准化向量其中D为W的对角元素矩阵，使得系数矩阵对角线元素为1；步骤3：预条件；构造预条件矩阵P＝S+IM，计算系数矩阵和常数向量其中S为一个和有关的矩阵：步骤4：根据步骤3输出的系数矩阵和常数向量设置迭代初始解为x(0)＝0，并开始迭代操作，输出检测结果；算法伪代码如下：K轮迭代后，x(K)即为待检测信号的估计结果。2.一种近似MMSE性能的高效大规模MIMO检测系统，其特征在于，包括：预条件模块，用于完成预处理，将信道矩阵H和接收信号向量y输入检测器，得到匹配滤波器输出yMF＝HHy和规则化Gram矩阵W＝G+NoIM，其中Gram矩阵G＝HHH，No为噪声方差，IM为M维单位矩阵，(.)H为共轭转置操作；然后计算标准化矩阵和标准化向量其中D为W的对角元素矩阵，使得系数矩阵对...

【专利技术属性】
技术研发人员：张川，吴至榛，尤肖虎，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32