基于分组ML检测和并行迭代干扰抵消的上行MIMO检测方法

基于分组ml检测和并行迭代干扰抵消的上行mimo检测方法
技术领域
1.本发明涉及一种基于分组ml(最大似然)检测和并行迭代干扰抵消的上行mimo(多天线传输)检测方法，属于无线移动通信技术领域。


背景技术：

2.大规模mimo(mutiple-input multiple-output)技术首次在第三代移动无线通信网络研究中被提出，该技术旨在利用发送接收多天线来提高频谱效率以及链路可靠度。因此，大规模mimo技术也成为用于满足用户对第5代移动通信系统更高服务质量(qos)诉求的关键技术。然而天线与天线之间往往存在干扰，所以为进一步提升通信质量往往需要在基站中额外部署一个mimo检测算法以抵消天线之间的干扰。
3.近年来，无线移动网络用户数量巨幅增加，基站与用户之间的交互在一些小城市也已达到艾字节量级。自2015年起至2021年，全球基站与用户的总数据交互量呈指数级增加，用户数的增加对无线移动网络服务者提出了更高频谱效率，更高能量效率，更高传输速率以及更好的移动性的需求。为应对这些需求，一种解决方案便是使用大规模mimo技术，但是随着天线规模的增大和允许接入的上行链路的增加，mimo检测器的复杂度也随之快速提高，为解决此问题需要对应部署高性能低复杂度的检测算法
4.实际上自mimo技术问世近50年来，mimo检测器的研究一直在学界具有极高的关注度。一方面随着芯片工艺的提升，计算机算力的增强使得复杂度较高的高性能mimo检测算法的部署成为可能。另一方面，随着学界对mimo技术研究的深入，大量高性能可物理实现的算法被提出。mimo检测算法可大致分为线性检测算法和非线性检测算法两类，其中主流且常用的算法为基于线性最优推导得到的mmse检测算法，以及基于近似最大似然的固定复杂度的球形译码(fcsd)算法。然而与非线性检测算法性能相比mmse算法检测性能较差，而具有高检测性能的fcsd算法实现复杂度较高。并且，随着发送数据流的增加，两种算法的复杂度都有显著提升。从检测性能上来说，fcsd检测算法的性能与数据流排序准确度有很大关系，但高性能的排序算法诸如snr准则排序方法算法复杂度也较高，相比较而言mmse算法虽然复杂度略低但性能有较明显的衰减并且上述两种算法的算法复杂度都随发送数据流数的增加有显著的增加。
5.综上所述，目前对于大规模mimo检测算法的研究一方面需要解决高性能与低复杂度之间的矛盾，另一方面需要解决发送数据流大幅增加以后带来的复杂度的提高问题。


技术实现要素：

6.本发明所要解决的技术问题是如何提供一种复杂度适中的高性能上行mimo检测方案。
7.本发明为了解决上述技术问题采用以下技术方案：本发明设计了一种基于分组ml检测和并行迭代干扰抵消的上行mimo检测方法，用于针对上行集中式单用户或者多用户mimo接收信号进行实时检测，针对实时所获得的待检测mimo信号，执行如下步骤：
8.步骤a.接收待检测mimo信号，并针对待检测mimo信号进行分组ml检测获得并行迭代干扰抵消的初始值，设置好最大迭代次数t，然后进入步骤b，其中待检测mimo接收信号y可以被表示为：其中，hi表示第i个用户的信道，xi表示第i个用户的发送数据流，n表示加性白噪声，k表示总用户数，在集中式接收检测的基站中，可以按分块矩阵形式简洁的表达成y＝hx+n的形式
9.步骤b.执行并行干扰抵消，首先将其他组数据流对待检测组数据流的干扰减去，然后对待检测组中数据流用分组ml检测方法进行检测，上述过程以组为单位并行执行并视作一次迭代。本次迭代检测的输出作为下一次迭代的并行干扰抵消初始值。
10.步骤c.跳转执行步骤b，并将迭代次数加1，直至迭代次数达到预设值，此时并行迭代干扰抵消检测器输出最终检测结果。
11.作为本发明的一种关键技术方案，所述步骤a中的分组ml检测包括如下步骤：
12.步骤a1：对发送数据流及信道矩阵进行分组，组内数据流数l应满足2≤l≤4。设发送天线数为n
t
，则总共可以分为：p＝n
t
/l组，分组为连续均匀分组。
13.步骤a2：根据接收数据按步骤a1的分组结果进行分组排序，并按排序结果调整对应信道矩阵的列顺序为：
14.步骤a3：对调整过列序的矩阵做qr分解得到根据qr分解结果得到接收向量的有效部分y
eff
＝qhy，对矩阵r按步骤a1的分组结果进行分组则第p个分组可以表示为：r
p
＝[r
(p-1)*l+1
,r
(p-1)*l+2
,
…
,r
(p-1)*l+l
]，ri为矩阵r的第i列向量。
[0015]
步骤a4：组内使用最大似然检测，第t次迭代第p组的检测过程为c
l
表示属于调制空间c的调制符号在l个数据流上并行发送时所生成的所有可能的发送符号向量所构成的空间。检测完成后将第p组对其他数据流的干扰抵消，抵消过程可表示为：
[0016]
步骤a5：按p＝p,p-1,
…
,1的顺序进行检测，重复步骤a4直到所有数据流都被检测，输出最终检测结果。
[0017]
作为本发明的一种优选技术方案，所述步骤a2中的分组排序方法包括如下步骤：
[0018]
步骤a2.1：用h
origin
＝h存储得到的信道矩阵后，对信道矩阵h求逆得到n
t
为发送天线数，为信道逆矩阵h
+
第i行向量。
[0019]
步骤a2.2：根据上述的分组结果分别计算待排序组内数据流的总放大系数score，计算方法如下：当前应当被检测的组为
[0020]
步骤a2.3：删除第p
cur
组数据流在信道矩阵中对应的列向量，该过程可以表示为：并对进行过删除操作的矩阵进行求逆得到新的h
+
。本发明中定义矩阵的减操作(-{})为删除{}中对应的列。
[0021]
步骤a2.4：重复步骤a2.2，a2.3直至所有组都被排序。
[0022]
步骤a2.5：按照排序结果调整原信道矩阵h
origin
的列顺序使得排序靠前的组首先被检测得到
[0023]
作为本发明的一种关键技术方案，所述步骤b中的迭代干扰抵消方法包括如下步骤：
[0024]
步骤b1：以上述的分组ml算法作为检测的初始值上标表示迭代次数。
[0025]
步骤b2：并行的以组为单位进行干扰抵消，对于第p组数据的干扰抵消过程可以表示为：其中表示第t次迭代，t＝0,
…
,t-1时得到得第i组数据流的估计值，具体可以表示为
[0026]
步骤b3：并行的以组为单位进行ml检测，第p组数据的检测结果可以表示为：
[0027]
步骤b4：以b2和b3所构成的干扰抵消与检测为一次迭代，t＝t+1重复执行直至达到预设的迭代次数t＝t通常t≤5。
[0028]
本发明所述一种基于分组ml检测和并行迭代干扰抵消的上行mimo检测方法采用以上技术方案与现有技术相比，具有以下技术效果：
[0029]
通过分组ml检测以及简化的分组排序方法相比于现有的高性能近似mld译码检测器(如基于排序的固定复杂度球形检测方法(fcsd))降低了检测的复杂度，通过采取组内ml检测以及并行迭代干扰抵消的方法确保了检测的准确度相比于最优线性检测方法(mmse)有更好的性能。
附图说明
[0030]
图1是本发明设计方法的检测流程图；
[0031]
图2是本发明设计方法中2数据流分组时，分组ml检测的树形搜索示意图；
[0032]
图3是仿真实施例一：独立复高斯信道下16发送数据流32接收天线所对应的检测性能图；
[0033]
图4是仿真实施例二：相关信道下16发送数据流32接收天线所对应的检测性能图。
具体实施方式
[0034]
下面结合说明书附图对本发明的具体实施方式作详细的说明。
[0035]
本发明提供了一种基于分组ml检测和并行干扰抵消的上行mimo检测方法，该方法主要由分组ml检测算法和并行迭代干扰抵消两个主要部分组成，其中分组ml算法的检测结果会作为并行干扰抵消方法的初始值，并行干扰抵消完成以后对组内元素采用ml检测方法。作为分组ml检测方法排序部分的一种优选方案，本发明提出了一种基于snr准则以组为
单位进行排序的方案，该方案相比于snr准则排序方案在排序过程中求逆次数较少所以复杂度更低，相比于不排序的分组ml检测方案又具有更好的性能。综上，本发明给出了一种新型上行mimo检测方案，可以兼顾计算复杂度和检测性能，相比于线性检测器有更好的检测性能，相比于高性能mld检测器有更低的复杂度，尤其是在具有一定分集度的mimo系统中具有接近最优ml检测器的性能。接下来针对本发明所包含的上述重要技术做进一步详细说明。
[0036]
首先给出本发明仿真过程中所依据的mimo系统模型。在集中式上行多用户mimo系统中若假设多个用户同时向基站发送信号则基站侧的接收信号可以表示为：
[0037][0038]
其中hi表示第i个用户的信道，xi表示第i个用户的发送数据流，n表示加性白噪声，k表示总用户数，在集中式接收检测的基站中，可以按分块矩阵形式简洁的表达成y＝hx+n的形式。关于每个用户的信道矩阵hi，本发明假设为相关或独立的复高斯的平坦瑞利衰落信道。复高斯平坦瑞利衰落信道即信道矩阵中的每一个元素满足零均值高斯分布即：这里表示克罗内克积，在本发明研究的mimo系统中假设每个用户都满流发送，用户发送天线数与发送数据流数相等用lk表示。相关信道即增加了相关性处理的复高斯信道，数学上可以表示为：
[0039][0040]
其中r
t
和rr为发送端和接收端的相关矩阵，两者元素的取值方法类似，相关矩阵r
t
(rr)中的第i行第j列元素的取值方式如下：
[0041][0042]
不同之处在于相关矩阵rr的维度为nr×
nr，相关矩阵r
t
的维度为l
×
l。式(4)中的ρ被称作相关系数满足0＜ρ＜1，相关系数设置的越大，生成的相关信道矩阵的相关性越强。
[0043]
对于上行多用户mimo系统来说，来自同一个用户的数据流之间的相关性往往较高，而不同用户间的数据流的相关性低。依据此结论本发明中所述的分组算法尽量将同一用户的数据流分入一组中。即每一组中的数据流l满足：
[0044][0045]
对应的总分组数p则满足：
[0046][0047]nt
＝k
×
lk表示总发送流数，在本系统模型中假设发送流数与发送天线数相等因此也等于总发送天线数，式(2)中的l
max
为人为设置参数，设置时需要综合考虑性能以及运算复杂度，经本发明研究发现，一般设置为4时可以兼顾上述两方面考虑。对于检测算法来说，数据流间的相关性越高，检测准确的难度越大。例如对于线性mmse检测器来说，其检测性能(检测器输出数据流与真实数据流之间的差别)会随数据流间相关度的提高而大幅增加。因此，本发明对于同组内的高相关数据流采用性能最佳的ml检测方案即
[0048]
本发明提出的分组ml算法在组间使用了串行干扰抵消方法，串行干扰抵消方法以及球形译码算法如fcsd算法对于检测顺序较为敏感，合理的检测顺序可以提高该方法的检测准确度。但是根据最优检测准则(如snr准则)得到检测顺序的排序算法，在确定检测顺序时就会消耗大量运算资源。由于数值计算的原因在确定检测顺序的过程中，每次输出当前的待检测数据流，就需要删除该数据流在信道矩阵对应的列以后对信道矩阵h进行重新求逆，运算量巨大，尤其是在大规模mimo系统中其算法复杂度θ
∝nt4
，但是本分组ml算法中采用的分组排序算法的复杂度θ
∝
p
×nt3
，在n
t
较大时可以显著降低算法复杂度。
[0049]
由于采用简化的排序算法会降低检测的准确度，因此在本发明中还增设了一级并行迭代干扰抵消。并行迭代干扰抵消方法以上述分组ml算法的检测结果作为初始值进行干扰抵消，并且利用了分组排序的结果。作为一种优选方法在迭代干扰抵消中继续对组内采用了ml检测方法，在不显著提高算法复杂度的情况下提升了检测性能。之前提到的排序复杂度可在一定程度上转化为并行迭代干扰抵消产生的运算即在排序时使用更简单的排序算法但是在迭代干扰抵消时增加迭代次数，本发明给出的分组排序加适量次数的并行迭代干扰抵消是经过仿真得出的一种兼顾复杂度和性能的优选方案，本发明在实际部署使用时也可以视情况进行调整。
[0050]
综上，本发明设计了一种基于分组ml检测和并行干扰抵消的上行mimo检测方法，用于对上行大规模mimo场景下的待检测mimo信号进行实时检测，具体实际应用当中，针对实时所获得的待检测上行mimo信号，执行如下步骤：
[0051]
步骤a.接收待检测mimo信号，并针对待检测mimo信号进行分组ml检测获得并行迭代干扰抵消的初始值，设置好最大迭代次数t，然后进入步骤b，其中待检测mimo接收信号y可以被表示为：
[0052]
其中，hi表示第i个用户的信道，xi表示第i个用户的发送数据流，n表示加性白噪声，k表示总用户数，在集中式接收检测的基站中，可以按分块矩阵形式简
[0053]
y＝h1x1+h2x2+
…
+hkxk+n
[0054]
洁的表达成y＝hx+n的形式
[0055]
步骤b.执行并行干扰抵消，首先将其他组数据流对待检测组数据流的干扰减去，然后对待检测组中数据流用分组ml检测方法进行检测，上述过程以组为单位并行执行并视作一次迭代。本次迭代检测的输出作为下一次迭代的并行干扰抵消初始值。
[0056]
步骤c.跳转执行步骤b，并将迭代次数加1，直至迭代次数达到预设值，此时并行迭代干扰抵消检测器输出最终检测结果。
[0057]
作为本发明的一种关键技术方案，所述步骤a中的分组ml检测包括如下步骤：
[0058]
步骤a1：对发送数据流及信道矩阵进行分组，组内数据流数l应满足2≤l≤4。设发送天线数为n
t
，则总共可以分为：p＝n
t
/l组，分组为连续均匀分组。
[0059]
步骤a2：根据接收数据按步骤a1的分组结果进行分组排序，并按排序结果调整对应信道矩阵的列顺序为：
[0060]
步骤a3：对调整过列序的矩阵做qr分解得到根据qr分解结果得到接收向量的有效部分y
eff
＝qhy，对矩阵r按步骤a1的分组结果进行分组则第p个分组可以表示为：r
p
＝[r
(p-1)*l+1
,r
(p-1)*l+2
,
…
,r
(p-1)*l+l
]，ri为矩阵r的第i列向量。
[0061]
步骤a4：组内使用最大似然检测，第t次迭代第p组的检测过程为c
l
表示属于调制空间c的调制符号在l个数据流上并行发送时所生成的所有可能的发送符号向量所构成的空间。检测完成后将第p组对其他数据流的干扰抵消，抵消过程可表示为：
[0062]
步骤a5：按p＝p,p-1,
…
,1的顺序进行检测，重复步骤a4直到所有数据流都被检测，输出最终检测结果。
[0063]
作为本发明的一种优选技术方案，所述步骤a2中的分组排序方法包括如下步骤：
[0064]
步骤a2.1：用h
origin
＝h存储得到的信道矩阵后，对信道矩阵h求逆得到n
t
为发送天线数，为信道逆矩阵h
+
第i行向量。
[0065]
步骤a2.2：根据上述的分组结果分别计算代排序组内数据流的总放大系数score，计算方法如下：当前应当被检测的组为
[0066]
步骤a2.3：删除第p
cur
组数据流在信道矩阵中对应的列向量，该过程可以表示为：并对进行过删除操作的矩阵进行求逆得到新的h
+
。本发明中定义矩阵的减操作(-{})为删除{}中对应的列。
[0067]
步骤a2.4：重复步骤a2.2，a2.3直至所有组都被排序。
[0068]
步骤a2.5：按照排序结果调整原信道矩阵h
origin
的列顺序使得排序靠前的组首先被检测得到
[0069]
作为本发明的一种关键技术方案，所述步骤b中的迭代干扰抵消方法包括如下步骤：
[0070]
步骤b1：以上述的分组ml算法作为检测的初始值上标表示迭代次数。
[0071]
步骤b2：并行的以组为单位进行干扰抵消，对于第p组数据的干扰抵消过程可以表示为：其中表示第t次迭代，t＝0,
…
,t时得到得第i组数据流的估计值，具体可以表示为
[0072]
步骤b3：并行的以组为单位进行ml检测，第p组数据的检测结果可以表示为：
[0073]
步骤b4：以b2和b3所构成的干扰抵消与检测为一次迭代，t＝t+1重复执行直至达到预设的迭代次数t＝t通常t≤5。
[0074]
接下来针对本发明所设计的基于分组ml和并行干扰抵消得上行mimo检测方法，结合仿真进行说明如下：
[0075]
仿真示例一：仿真条件如下表1所示：
[0076]
表1仿真示例一仿真条件
[0077]
用户发送数据流数lk4接收天线数nr32用户数4调制方式qpsk相关系数ρ0迭代次数t5组内数据流数l4分组个数p4
[0078]
图3所示给出了基于以上条件下，本发明提出的检测算法与常用线性检测算法mmse检测算法以及非线性检测算法fcsd算法的性能比较图。可以看出本发明提出的算法的性能在上述场景下性能均好于经典算法。
[0079]
仿真示例二：仿真条件如下表2所示：
[0080]
表2仿真示例二仿真条件
[0081]
用户发送数据流数lk4接收天线数nr32用户数4调制方式qpsk相关系数ρ0.6迭代次数t5组内数据流数l4分组个数p4
[0082]
图3所示给出了基于以上条件下，本发明提出的检测算法与常用线性检测算法mmse检测算法以及非线性检测算法fcsd算法的性能比较图。虽然相关信道使得所有检测算法的性能都有所下降，但本发明提出的算法的性能在上述场景下性能均好于经典算法。
[0083]
下面结合附图对本发明的实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。