基于并行Cholesky矩阵分解的MIMO系统信号检测系统及方法

基于并行cholesky矩阵分解的mimo系统信号检测系统及方法
技术领域
1.本发明涉及通信技术领域，更进一步涉及多天线技术领域中的一种基于并行cholesky矩阵分解的多输入多输出mimo(multiple-input multiple-output)系统信号检测系统及方法。本发明能以较低的复杂度与较低的时延完成基于并行cholesky矩阵分解的mimo系统的信号检测。


背景技术：

2.mimo系统利用空间复用技术，通过将原本串行的数据，并行的由多根天线发送，在不增加天线发射功率以及消耗更多频谱资源的前提下，使系统吞吐率成倍提升，被应用于各类无线通信场景之中。近年来，mimo信号检测方法的研究主要围绕k-best信号检测方法展开，该方法可以通过设置保留路径数量的方式，在性能和硬件资源消耗上做出选择，具有极高的灵活性，但存在矩阵分解过程中实现复杂度高、信号搜索过程中展开子节点数量多等问题，导致该方法完成信号检测的任务需要消耗大量的时间和计算资源，这与新一代通信系统低时延的发展趋势不符。
3.乐鑫信息科技(上海)股份有限公司在其申请的专利文献“带有检测中信道矩阵预处理的mimo-ofdm无线信号检测方法和系统”(申请号：201910152084.3，申请公布号：cn 111628952 a)公开了一种基于k-best的无线信号检测方法及其系统。该专利文献公开的方法当接收到的第一个ofdm符号后，对每个子载波的信道矩阵进行预处理，以生成全局动态k值表，确定了每个子载波进行信号检测时保留的路径数量。同时，为了进一步降低方法整体的复杂度，该方法从实数域进行信号检测，避免了复杂的复数运算。虽然该方法有效降低了信号检测过程中展开的子节点数量以及算法的复杂度，但是，该方法仍然存在的不足之处是，在实数域进行信号检测时，需对信道矩阵进行实数化，这会使信道矩阵大小翻倍，导致该方法在信号检测过程中计算复杂度过高，此外，大小翻倍后的信道矩阵在进行信道矩阵分解以及k-best搜索时会带来处理时延翻倍的问题。该专利文献公开的系统由信道矩阵预处理子模块、存储子模块和k-best搜索子模块组成。首先，信道矩阵预处理子模块对第一个ofdm符号的每个子载波的信号矩阵执行信号矩阵预处理，以生成全局动态k值表；接着，存储子模块存储全局动态k值表以及每个子载波的搜索参数；最后，k-best搜索子模块对mimo-ofdm数据包中每个后续ofdm符号进行mimo检测。为了降低系统复杂度，该系统利用实数域替代复数域计算，但是，该系统存在的不足之处是，在实数域进行信号检测时，需要对矩阵实数化，使得信号矩阵大小翻倍，增加了系统在信号检测时计算复杂度高，同时k-bset搜索也会导致系统处理时延翻倍。


技术实现要素：

4.本发明的目的在于针对上述现有技术的不足，提出一种基于并行cholesky矩阵分解的mimo系统信号检测系统及方法，用于解决现有mimo系统信号检测时复杂度高，以及在实数域进行信号检测时，检测时延翻倍的问题。
5.为了实现上述目的，本发明的思路是，本发明的方法通过对实数化后的信道矩阵采用cholesky矩阵模型进行并行矩阵分解，重新设计的信道矩阵分解过程，能以完成对一行元素计算的代价实现对相邻两行元素的并行分解，降低了对实数化后信道矩阵分解时所需计算次数，有效提高了信道矩阵的分解速度，相较于传统信道矩阵分解方法，克服了现有矩阵分解技术中，信道矩阵大小翻倍后带来的信道矩阵分解复杂度过高、分解翻倍时间不足的缺陷，并且，对于分解得到的上三角矩阵，其奇数行对角线元素右侧第一个元素数值固定为0，充分利用了分解得到的矩阵特性。本发明的系统在搜索过程采用k-best算法对相邻两层同时进行子节点展开，并在子节点展开时，将各保留路径依据累积欧氏距离的排序进行划分，每类路径按照se枚举规则顺序展开不同数量的子节点，在对相邻两层展开的子节点进行组合、排序后，最终完成k个最佳路径的选择，相较于传统搜索系统，本发明改进后的系统仅对部分节点进行展开，有效减少了展开的子节点数量，克服了现有方法中子节点展开数量过多的不足。为了进一步降低时延，本发明的系统在进行信号搜索时，在完成相邻两层子节点的并行展开后，通过对子节点进行二次展开，将相邻两层的搜索过程合并为一次，仅需传统方法一半的信号搜索次数，便可完成对全部信号的搜索，有效解决了信道矩阵大小翻倍带来的信号搜索时间翻倍问题。
6.实现本发明目的的技术方案的步骤如下：
7.本发明的系统，包括信道矩阵转换模块、检测结果判决模块、信道矩阵分解模块、k-best并行搜索模块；其中：
8.所述信道矩阵转换模块，用于将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵；
9.所述信道矩阵分解模块，用于利用cholesky矩阵分解方法并行分解迭代矩阵，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数；将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发射天线的数量；利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵；判断当前迭代矩阵是否还存在未被移除的元素，若是，则从当前迭代矩阵中选取前两行元素作为待分解元素，否则，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q；
10.所述k-best并行搜索模块，用于利用k-best方法并行搜索相邻两层的路径，从mimo系统的所有发射天线中，选取一根未检测的发射天线；将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定；将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径；通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径；
11.所述检测结果判决模块，用于判断是否选完mimo系统的所有发射天线，若是，则确定发射信号的检测值，否则，利用k-best方法并行搜索相邻两层的路径，从mimo系统的所有发射天线中，选取一根未检测的发射天线，进行并行搜索；
12.所述检测结果判决模块，用于确定发射信号的检测值。
13.选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。
14.本发明的信号检测方法的步骤包括如下：
15.步骤1，信道矩阵转换模块将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵；
16.步骤2，信道矩阵分解模块利用用cholesky矩阵分解方法并行分解迭代矩阵：
17.步骤2.1，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数；
18.步骤2.2，将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发射天线的数量；
19.步骤2.3，利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵；
20.步骤2.4，判断当前迭代矩阵是否还存在未被移除的元素，若是，则执行步骤2.1，否则，执行步骤2.5；
21.步骤2.5，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q；
22.步骤3，k-best并行搜索模块利用k-best方法并行搜索相邻两层的路径：
23.步骤3.1，从mimo系统的所有发射天线中，选取一根未检测的发射天线；
24.步骤3.2，将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定；
25.步骤3.3，将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径；
26.步骤3.4，通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径；
27.步骤4，检测结果判决模块判断是否选完mimo系统的所有发射天线，若是，则执行步骤5，否则，执行步骤3；
28.步骤5，检测结果判决模块确定发射信号的检测值：
29.选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。
30.本发明与现有技术相比具有以下优点：
31.第一，本发明的方法通过对实数化后的信道矩阵采用cholesky矩阵模型进行并行分解迭代，与其它实数域信号检测方法相比，以完成对一行元素计算的代价实现对相邻两
行元素的并行分解，克服了现有技术在进行实数域信号检测时，信道矩阵大小翻倍后带来的信道矩阵分解复杂度过高的缺陷，有效提高了信道矩阵的分解速度。使得本发明能以传统实数域信道矩阵分解方法一半的计算复杂度，完成同样的信道矩阵分解任务，具有低复杂度、高分解速度的优势。
32.第二，本发明的系统在k-best模块中进行子节点展开时，将各保留路径依据累积欧氏距离的排序进行划分，每类路径按照se枚举规则顺序展开不同数量的子节点，在对相邻两层展开的子节点进行组合、排序后，最终完成k个最佳路径的选择。通过对部分子节点进行节点展开，有效降低了展开子节点的数量，降低了进行子节点展开时的计算复杂度。克服了现有系统中子节点展开数量过多的不足，使得本发明信号检测的复杂度进一步降低，减少了mimo系统中使用的资源数量。
33.第三，本发明的系统为了更进一步降低时延，本系统在进行信号搜索时，在完成了相邻子节点的并行展开后，通过对子节点进行二次展开，将相邻两层的搜索过程合并为一次，有效的解决了信道矩阵大小翻倍后带来的信号搜索时间翻倍问题，使得本发明能更快的完成搜索任务。
附图说明
34.图1是本发明系统框图；
35.图2是本发明算法实现步骤流程图；
36.图3是本发明k-best信号搜索流程图；
37.图4是本发明se枚举示意图；
38.图5是本发明二次展开方案示意图。
具体实施方式
39.下面结合附图和实施例，对本发明做进一步的详细描述。
40.参照图1，对本发明系统做进一步的详细描述。
41.本发明的系统，包括信道矩阵转换模块、检测结果判决模块、信道矩阵分解模块、k-best并行搜索模块；其中：
42.所述信道矩阵转换模块，用于将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵。
43.所述信道矩阵分解模块，用于利用cholesky矩阵分解方法并行分解迭代矩阵，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数；将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发射天线的数量；利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵；判断当前迭代矩阵是否还存在未被移除的元素，若是，则从当前迭代矩阵中选取前两行元素作为待分解元素，否则，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q。
44.所述k-best并行搜索模块，用于利用k-best方法并行搜索相邻两层的路径，从
mimo系统的所有发射天线中，选取一根未检测的发射天线；将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定；将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径；通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径。
45.所述检测结果判决模块，用于判断是否选完mimo系统的所有发射天线，若是，则确定发射信号的检测值，否则，利用k-best方法并行搜索相邻两层的路径，从mimo系统的所有发射天线中，选取一根未检测的发射天线，进行并行搜索。
46.所述检测结果判决模块，用于确定发射信号的检测值。
47.选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。
48.参照图2，对本发明方法的实现步骤做进一步的详细描述。
49.步骤1，信道矩阵转换模块将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵。
50.所述的将mimo系统的信道复矩阵转化为信道实矩阵是由下式实现的：
[0051][0052]
其中，h
′
表示信道实矩阵，re(
·
)表示实部符号，h1表示信道复矩阵的第1个列向量，hn表示信道复矩阵的第n个列向量，n的取值等于mimo系统中发射天线的总数，im(
·
)表示虚部符号。
[0053]
所述的通过信道实矩阵生成迭代矩阵是指，通过信道实矩阵得到迭代矩阵中奇数行元素的数值，选择迭代矩阵奇数行元素中处于奇数列位置的元素，将该元素与其右方、下方以及右下方的元素划分为一组，在组内按照主对角线方向元素数值相同，副对角线方向元素数值互为相反数的方式，逐个得到偶数行元素的数值，其中，迭代矩阵中奇数行元素的数值是由下式得到的：
[0054][0055]
其中，a
i,j
表示迭代矩阵中第i行，第j列的元素，i为在区间[1,2n]内选取的一个奇数值，n的取值等于mimo系统中发射天线的总数，j为在区间[i,2n]内选取的一个实数值，hi′
表示信道实矩阵第i列元素组成的向量，(
·
)
t
表示对矩阵的转置运算。
[0056]
本发明实施例中信道实矩阵的行的总数为4，列的总数为4，通过信道实矩阵得到的迭代矩阵a形式如下：
[0057]
[0058]
步骤2，信道矩阵分解模块利用cholesky矩阵分解方法并行分解迭代矩阵：
[0059]
步骤2.1，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数。
[0060]
所述的通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素是指，对待分解元素中第1行，第1列的元素进行开方运算，得到上三角矩阵r中第2o-1行，第2o-1列元素的数值，通过待分解元素得到上三角矩阵r第2o-1行剩余元素的数值后，使第2o-1行处于奇数列位置的元素，取值等于右上方元素数值的相反数，处于奇数列位置的元素，取值等于左上方元素的数值，其中，上三角矩阵r第2o-1行剩余元素的数值是由下式得到的：
[0061][0062]
其中，r
2o-1,2o+l
表示上三角矩阵r中第2o-1行，第2o+l列的元素，o的取值等于迭代次数，1≤o≤n，n的取值等于mimo系统中发射天线的总数，l为在区间[2,2n-2o]内选取的一个实数值，a
1,l
表示待分解元素中第1行，第l列的元素，r
2o-1,2o-1
表示上三角矩阵r中第2o-1行，第2o-1列的元素。
[0063]
步骤2.2，将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发射天线的数量。
[0064]
步骤2.3，利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵。
[0065]
所述的利用中间矩阵相关元素更新移除分解元素后的矩阵是由下式实现的：
[0066]a″
2x-1,2y-1
＝a2′
x-1,2y-1-(t
2x-1,2y-1
+t
2x,2y
)
[0067]a″
2x-1,2y
＝a
′
2x-1,2y-(t
2x-1,2y-t
2x,2y-1
)
[0068]
其中，a
″
2x-1,2y-1
表示更新后的迭代矩阵第2x-1行，第2y-1列的元素，x为在区间[1,n-o]内选取的一个实数值，n的取值等于mimo系统中发射天线的总数，y为在区间[x,n-o]上选取的一个实数值，a
′
2x-1,2y-1
表示移除分解元素后的矩阵第2x-1行，第2y-1列的元素，t
2x-1,2y-1
表示中间矩阵第2x-1行，第2y-1列的元素，t
2x,2y
表示中间矩阵第2x行，第2y列的元素。
[0069]
步骤2.4，判断当前迭代矩阵是否还存在未被移除的元素，若是，则执行步骤2.1，否则，执行步骤2.5。
[0070]
步骤2.5，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q。
[0071]
所述的通过上三角矩阵r和信道实矩阵生成等效酉矩阵q是由下式实现的：
[0072]
q＝(r-1
)
th′
t
[0073]
其中，q表示等效酉矩阵，r-1
表示上三角矩阵r的逆矩阵，(
·
)
t
表示对矩阵或向量的转置运算，h
′
表示信道实矩阵。
[0074]
步骤3，k-best并行搜索模块利用k-best方法并行搜索相邻两层的路径。
[0075]
参照图3，对本发明采用k-best方法并行搜索相邻两层的路径的步骤做进一步的详细描述。
[0076]
步骤3.1，从mimo系统的所有发射天线中，选取一根未检测的发射天线。
[0077]
步骤3.2，将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定。
[0078]
本发明的实施例中，对信号调制阶数为64qam，保留路径数量为8的k-best搜索算法，设置高置信度路径数量为3条，每条路径展开4个子节点，中置信度路径为3条，每条路径展开2个子节点，低置信度路径为2条，每条展开1个子节点。
[0079]
se枚举规则如图4所示，图4中，-5、-3、-1、1、3、5表示星座图映射后的实部值，表示当前路径下计算出来的信号值，枚举时，对星座映射值到的欧式距离以由小到大的顺序进行排序，并根据排序结果从码本中选择将要被展开的子节点。
[0080]
对信号估计值的计算是由下式实现的：
[0081][0082]
其中，表示当前路径第z个节点的信号估计值，qz表示等效酉矩阵第z行元素组成的向量，表示实数化接收信号向量，∑
·
表示求和操作，m的取值等于mimo系统中接收天线的总数，k为在区间[z+1,2m]上依次选取的实数值，r
z,k
表示上三角矩阵r第z行，第k列的元素。
[0083]
步骤3.3，将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径。
[0084]
本发明的实施例中，并行展开相邻两层子节点后，对子节点进行二次展开的方式如图5所示，首先，对天线i的2i+1层的同一个父节点进行展开操作，得到2i层展开的n个子节点和2i-1层展开的n个子节点；然后，对2i层和2i-1层的n个子节点进行二次展开，并将展开后相邻两层的子节点合并得到等效路径。
[0085]
步骤3.4，通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径。
[0086]
所述的计算每一条二次扩展后路径的累积欧氏距离是由下式实现的：
[0087][0088]
其中，d
m-2
表示待搜索节点数量为m-2的二次扩展后路径的累积欧氏距离，m的取值等于待搜索节点数量，1≤m≤2m，qm表示等效酉矩阵第m行元素组成的向量，表示实数化接收信号向量，表示当前路径第m个节点对应的信号估计值，其取值范围由当前mimo系统使用的信号调制方式决定，r
m,m
表示上三角矩阵r第m行，第m列的元素，|
·
|表示绝对值运算。
[0089]
步骤4，检测结果判决模块判断是否选完mimo系统的所有发射天线，若是，则执行步骤5，否则，执行步骤3。
[0090]
步骤5，检测结果判决模块确定发射信号的检测值：
[0091]
选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。

技术特征：
1.一种基于并行cholesky矩阵分解的mimo系统信号检测系统，包括信道矩阵转换模块、检测结果判决模块；其特征在于，还包括信道矩阵分解模块、k-best并行搜索模块；其中：所述信道矩阵转换模块，用于将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵；所述信道矩阵分解模块，用于利用cholesky矩阵分解方法并行分解迭代矩阵，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数；将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发射天线的数量；利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵；判断当前迭代矩阵是否还存在未被移除的元素，若是，则从当前迭代矩阵中选取前两行元素作为待分解元素，否则，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q；所述k-best并行搜索模块，用于利用k-best方法并行搜索相邻两层的路径，从mimo系统的所有发射天线中，选取一根未检测的发射天线；将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定；将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径；通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径；所述检测结果判决模块，用于判断是否选完mimo系统的所有发射天线，若是，则确定发射信号的检测值，否则，利用k-best方法并行搜索相邻两层的路径，从mimo系统的所有发射天线中，选取一根未检测的发射天线，进行并行搜索；所述检测结果判决模块，用于确定发射信号的检测值。选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。2.根据权利要求1所述检测系统的一种基于并行cholesky矩阵分解的mimo系统信号检测方法，其特征在于，采用cholesky矩阵分解方法并行分解迭代矩阵，采用k-best方法并行搜索相邻两层的路径，该方法的步骤包括如下：步骤1，信道矩阵转换模块将mimo系统的信道复矩阵转化为信道实矩阵，通过信道实矩阵生成迭代矩阵；步骤2，信道矩阵分解模块利用用cholesky矩阵分解方法并行分解迭代矩阵：步骤2.1，从当前迭代矩阵中选取前两行元素作为待分解元素，通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素，同时将待分解元素从当前迭代矩阵中移除，得到移除分解元素后的矩阵，其中，o的取值等于迭代次数；步骤2.2，将上三角矩阵r中第2o-1行，第2o+1列至第2n列元素组成的行向量转置后与自身行向量相乘，组成中间矩阵；其中，o的取值等于迭代次数，n的取值等于mimo系统中发
射天线的数量；步骤2.3，利用中间矩阵相关元素更新移除分解元素后的矩阵，得到更新后的迭代矩阵；步骤2.4，判断当前迭代矩阵是否还存在未被移除的元素，若是，则执行步骤2.1，否则，执行步骤2.5；步骤2.5，通过上三角矩阵r和信道实矩阵生成等效酉矩阵q；步骤3，k-best并行搜索模块利用k-best方法并行搜索相邻两层的路径：步骤3.1，从mimo系统的所有发射天线中，选取一根未检测的发射天线；步骤3.2，将保留的搜索路径按照累积欧式距离的排序划分为3类路径，按照se枚举规则同时展开相邻两层子节点，展开的子节点数量由路径类别以及k值决定；步骤3.3，将展开的下层节点添加至对应保留路径的末端，得到扩展后路径，再将展开的上层节点复制后，添加至对应扩展后路径的末端，对路径做第二次扩展，得到二次扩展后路径；步骤3.4，通过上三角矩阵r和等效酉矩阵q计算出每一条二次扩展后路径的累积欧氏距离，选择其中累积欧氏距离最小的k条路径作为保留路径；步骤4，检测结果判决模块判断是否选完mimo系统的所有发射天线，若是，则执行步骤5，否则，执行步骤3；步骤5，检测结果判决模块确定发射信号的检测值：选中保留路径中累积欧氏距离最小的路径，从根节点开始，每次选择两个值，将所选两个值分别作为实部与虚部，构成一个发射调制信号复数值的检测值，直至到达叶子结点，由所有发射调制信号复数值的检测值组成的复数集合即为所有发射天线的信号检测值。3.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤1中所述的将mimo系统的信道复矩阵转化为信道实矩阵是由下式实现的：其中，h
′
表示信道实矩阵，re(
·
)表示实部符号，h1表示信道复矩阵的第1个列向量，h
n
表示信道复矩阵的第n个列向量，n的取值等于mimo系统中发射天线的总数，im(
·
)表示虚部符号。4.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤1中所述的通过信道实矩阵生成迭代矩阵是指，通过信道实矩阵得到迭代矩阵中奇数行元素的数值，选择迭代矩阵奇数行元素中处于奇数列位置的元素，将该元素与其右方、下方以及右下方的元素划分为一组，在组内按照主对角线方向元素数值相同，副对角线方向元素数值互为相反数的方式，逐个得到偶数行元素的数值，其中，迭代矩阵中奇数行元素的数值是由下式得到的：其中，a
i,j
表示迭代矩阵中第i行，第j列的元素，i为在区间[1,2n]内选取的一个奇数值，n的取值等于mimo系统中发射天线的总数，j为在区间[i,2n]内选取的一个实数值，h
i
′
表示信道实矩阵第i列元素组成的向量，(
·
)
t
表示对矩阵的转置运算。
5.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤2.1中所述的通过待分解元素得到上三角矩阵r中第2o-1行和第2o行元素是指，对待分解元素中第1行，第1列的元素进行开方运算，得到上三角矩阵r中第2o-1行，第2o-1列元素的数值，通过待分解元素得到上三角矩阵r第2o-1行剩余元素的数值后，使第2o-1行处于奇数列位置的元素，取值等于右上方元素数值的相反数，处于奇数列位置的元素，取值等于左上方元素的数值，其中，上三角矩阵r第2o-1行剩余元素的数值是由下式得到的：其中，r
2o-1,2o+l
表示上三角矩阵r中第2o-1行，第2o+l列的元素，o的取值等于迭代次数，1≤o≤n，n的取值等于mimo系统中发射天线的总数，l为在区间[2,2n-2o]内选取的一个实数值，a
1,l
表示待分解元素中第1行，第l列的元素，r
2o-1,2o-1
表示上三角矩阵r中第2o-1行，第2o-1列的元素。6.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤2.3中利用中间矩阵相关元素更新移除分解元素后的矩阵是由下式实现的：a
″
2x-1,2y-1
＝a
′
2x-1,2y-1-(t
2x-1,2y-1
+t
2x,2y
)a
″
2x-1,2y
＝a
′
2x-1,2y-(t
2x-1,2y-t
2x,2y-1
)其中，a
″
2x-1,2y-1
表示更新后的迭代矩阵第2x-1行，第2y-1列的元素，x为在区间[1,n-o]内选取的一个实数值，n的取值等于mimo系统中发射天线的总数，y为在区间[x,n-o]上选取的一个实数值，a
′
2x-1,2y-1
表示移除分解元素后的矩阵第2x-1行，第2y-1列的元素，t
2x-1,2y-1
表示中间矩阵第2x-1行，第2y-1列的元素，t
2x,2y
表示中间矩阵第2x行，第2y列的元素。7.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤2.5中通过上三角矩阵r和信道实矩阵生成等效酉矩阵q是由下式实现的：q＝(r-1
)
t
h
′
t
其中，q表示等效酉矩阵，r-1
表示上三角矩阵r的逆矩阵，(
·
)
t
表示对矩阵或向量的转置运算，h
′
表示信道实矩阵。8.根据权利要求1所述基于并行cholesky矩阵分解的mimo系统信号检测系统，其特征在于，步骤3.4中计算每一条二次扩展后路径的累积欧氏距离是由下式实现的：其中，d
m-2
表示待搜索节点数量为m-2的二次扩展后路径的累积欧氏距离，m的取值等于待搜索节点数量，1≤m≤2m，q
m
表示等效酉矩阵第m行元素组成的向量，表示实数化接收信号向量，表示当前路径第m个节点对应的信号估计值，其取值范围由当前mimo系统使用的信号调制方式决定，r
m,m
表示上三角矩阵r第m行，第m列的元素，|
·
|表示绝对值运算。

技术总结
本发明公开了一种基于并行Cholesky矩阵分解的MIMO系统信号检测系统及方法，主要解决现有MIMO系统信号检测时复杂度高，以及在实数域进行信号检测时，检测时延翻倍的问题。其实现步骤是：信道矩阵转换模块将信道复矩阵转化为信道实矩阵并生成迭代矩阵；信道矩阵分解模块利用Cholesky矩阵分解方法并行分解迭代矩阵；K-Best并行搜索模块利用K-Best并行搜索方法并行搜索相邻两层的路径；检测结果判决模块完成所以天线信号的检测并确定发射信号检测值。本发明对信号检测过程的并行化设计，有效减少了检测的复杂度，降低了检测时延，能以更少的资源数量实现对信号的检测。少的资源数量实现对信号的检测。少的资源数量实现对信号的检测。


技术研发人员：张顺 王祥光 程帅林 段明明 马建鹏 何先灯
受保护的技术使用者：西安电子科技大学
技术研发日：2023.05.25
技术公布日：2023/8/9