技术领域
[0001] 本
发明涉及移动通信系统技术领域，具体为一种基于MacQueen聚类的 MIMO预编码码本构造方法。
背景技术
[0002] 多输入多输出(Multiple Input Multiple Output，MIMO)系统架构作为4G-LTE中的关键技术得到了广泛的应用，最新出炉的5G标准，仍然继续采用MIMO技术作为其物理层架构。在采用MIMO技术时，因为基站端部署了多根天线，且小区内也存在着大量的用户，所以MIMO系统存在着严重的流间干扰和用户间干扰，预编码作为解决这一问题的有效手段成为了MIMO系统架构的核心功能模
块并得到了广泛的应用，同时也受到了学术界和各大通信企业的深入研究。
[0003] 预编码目前主要有线性预编码和非线性预编码两种，非线性预编码复杂度高实现困难，而线性预编码复杂度低、原理简单，普遍应用在实际的通信系统中。在已有线性预编码方案中又可以分为非码本的和基于码本的预编码，而实际通信系统中受反馈链路带宽的限制，往往无法实现理想的信道状态信息(CSI)反馈，对于部分CSI的反馈，3GPP规定采用基于码本的预编码方案。基于码本的预编码方案主要分为两个部分：码本的构造与码字选取，码本构造的方法有很多，例如DFT码本、Householder码本、层结构码本、Lloyd矢量量化码本等。其中Lloyd矢量量化码本是采用Lloyd聚类
算法来构造码本，但该方法生成的码本误码率性能好坏受限于初始码本的选取和误差
门限的设置，因此导致了码本误码率性能的不稳定。
发明内容
[0004] 本发明的目的在于提供一种基于MacQueen聚类的MIMO预编码码本构造方法，本发明通过采用MacQueen
聚类算法来替代Lloyd矢量量化码本中的 Lloyd聚类算法，使得在计算过程中并不涉及类似于Lloyd
迭代中的误差门限设置，避免了误差门限值的选取对码本误码率性能的影响，同时在每个信道矩阵H分类之后都更新了质心，为下一次信道矩阵的分类提供了更加准确的分类依据，因此本方法使MIMO系统的误码率性能得到进一步的提高，提高通信系统的可靠性。
[0005] 为实现上述目的，本发明提供如下技术方案：一种基于MacQueen聚类的 MIMO预编码码本构造方法，包括如下步骤：
[0006] 步骤1：随机选取Q个信道样本H；
[0007] 步骤2：对Q个信道样本中的前N(Q>>N)个信道矩阵H进行SVD分解，得到相应的N个预编码矩阵f，并将这N个预编码矩阵作为初始码本 F0＝{f1,f2,…,fN}，每个码字fi(i＝1,2,…,N)即是N个类的初始质心；
[0008] 步骤3：选取下一个信道矩阵H，计算其与每个码字矩阵fi相乘的 Frobenius范数的平方，并将该信道矩阵H分配到使Frobenius范数的平方最大的码字矩阵fi所对应的类中；
[0009] 步骤4：重新计算步骤3中更新了的类的质心fi；
[0010] 步骤5：重复步骤3和步骤4，直到Q个信道样本全部分配完毕；
[0011] 步骤6：将最终N个类所对应的质心作为预编码码本F。
[0012] 优选的，对Q个信道样本中的前N(Q>>N)个信道矩阵H进行SVD分解，得到相应的N个预编码矩阵f，其中预编码矩阵的获得按照以下公式进行：
[0013] H＝U∑V
[0014] f＝V(j)
[0015] 其中U、∑、V矩阵分别为信道矩阵SVD分解得到的左奇异矩阵，奇异值组成的对
角矩阵，右奇异矩阵，V(j)为矩阵V的前j列组成的矩阵。
[0016] 优选的，把信道矩阵按顺序逐个进行分类，且在每分配完一个信道矩阵H 之后，都要重新计算增加了新信道矩阵的类的质心f，其中质心f的计算按照以下公式进行：
[0017]
[0018] 其中 表示 中每个信道矩阵H进行SVD分解得到的预编码矩阵f组成的类。
[0019] 优选的，所述步骤3中选取下一个信道矩阵H，是从第N+1个信道矩阵H 开始，依次为第N+2,N+3,…,Q个信道矩阵H。
[0020] 优选的，所述步骤3中计算H与每个码字矩阵fi相乘的Frobenius范数的平方的计算表达式为：
[0021]
[0022] 其中 表示质心为fk的信道矩阵H组成的类。
[0023] 本发明提供了一种基于MacQueen聚类的MIMO预编码码本构造方法，具备以下有益效果：
[0024] 本发明通过采用MacQueen聚类算法来替代Lloyd矢量量化码本中的 Lloyd聚类算法，使得在计算过程中并不涉及类似于Lloyd迭代中的误差门限设置，避免了误差门限值的选取对码本误码率性能的影响，同时在每个信道矩阵H分类之后都更新了质心，为下一次信道矩阵的分类提供了更加准确的分类依据，因此本方法使MIMO系统的误码率性能得到进一步的提高，提高通信系统的可靠性。
附图说明
[0025] 图1为MIMO-OSTBC系统
框图；
[0027] 图3为本发明提供的基于MacQueen聚类的MIMO预编码码本构造方法流程图；
[0028] 图4为预编码码本大小L＝64时，本发明提供的基于MacQueen聚类的MIMO 预编码码本构造方法与Lloyd矢量量化码本，DFT码本的误码率(BER)性能曲线对比示意图；
[0029] 图5为预编码码本大小L＝16时，本发明提供的基于MacQueen聚类的MIMO 预编码码本构造方法与Lloyd矢量量化码本的误码率(BER)性能曲线对比示意图。
具体实施方式
[0030] 下面将结合本发明
实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
[0031] 在本实例中仿真采用的通信系统为MIMO-OSTBC系统，采用4-QAM调制， Alamouti空时编码，发射天线数Tx为4，接收天线数Rx为2，系统框图如图 1所示，其中预编码部分采用的是基于码本的预编码方案，其具体工作流程如图2所示。
[0032] 实施例1
[0033] 一种基于MacQueen聚类的MIMO预编码码本构造方法，包括如下步骤：
[0034] 步骤1：随机选取Q＝10000个信道矩阵H；
[0035] 步骤2：将前N＝64个信道矩阵进行SVD分解H＝UΣV，并根据f＝V(2)，即选取右奇异矩阵V的前两列作为预编码矩阵，得到相应的64个预编码矩阵 f，并将这64个预编码矩阵作为初始码本F0＝{f1,f2,…,f64}，码本大小L＝64，每个码字fi(i＝1,2,…,64)即是64个类的初始质心， 表示质心为fi的信道矩阵H组成的类；
[0036] 步骤3：选取第65个信道矩阵H65，计算该矩阵与每个码字矩阵fi相乘的 Frobenius范数的平方，即
[0037]
[0038] 然后将信道矩阵H65分配到使value值最大的码字矩阵f对应的类中,为了表述的更清楚，假设当i＝m时，value值最大，则H65被分配到 类中；
[0039] 步骤4：按照下述公式重新计算步骤3中更新了的 类的质心fm：
[0040]
[0041] 其中公式表示计算 类中每个f与 类中所有的信道矩阵H计算的值，并求它们的期望，然后选出使期望值最小的f作为该 类的新质心， 表示 中每个信道矩阵H进行SVD分解得到的预编码矩阵f 组成的类；
[0042] 步骤5：分别选取第66,67,…,10000个信道矩阵重复步骤3和步骤4；
[0043] 步骤6：将最终分配完成后，每个 类的质心fi作为码字组成我们所需要的预编码码本F。
[0044] 实施例2
[0045] 一种基于MacQueen聚类的MIMO预编码码本构造方法，包括如下步骤：
[0046] 步骤1：随机选取Q＝10000个信道矩阵H；
[0047] 步骤2：将前N＝16个信道矩阵进行SVD分解H＝UΣV，并根据f＝V(2)，即选取右奇异矩阵V的前两列作为预编码矩阵，得到相应的64个预编码矩阵 f，并将这64个预编码矩阵作为初始码本F0＝{f1,f2,…,f16}，码本大小L＝16，每个码字fi(i＝1,2,…,16)即是16个类的初始质心， 表示质心为fi的信道矩阵H组成的类；
[0048] 步骤3：选取第17个信道矩阵H17，计算该矩阵与每个码字矩阵fi相乘的 Frobenius范数的平方，即
[0049]
[0050] 然后将信道矩阵H17分配到使value值最大的码字矩阵f对应的类中,为了表述的更清楚，假设当i＝m时，value值最大，则H17被分配到 类中；
[0051] 步骤4：按照下述公式重新计算步骤3中更新了的 类的质心fm：
[0052]
[0053] 其中公式表示计算 类中每个f与 类中所有的信道矩阵H计算的值，并求它们的期望，然后选出使期望值最小的f作为该 类的新质心， 表示 中每个信道矩阵H进行SVD分解得到的预编码矩阵f 组成的类；
[0054] 步骤5：分别选取第18,19,…,10000个信道矩阵重复步骤3和步骤4；
[0055] 步骤6：将最终分配完成后，每个 类的质心fi作为码字组成我们所需要的预编码码本F。
[0056] 对本发明实施例1与已有预编码码本在上述通信系统中的误码率性能进行仿真对比，得到图4，对本发明实施例1与已有预编码码本在上述通信系统中的误码率性能进行仿真对比，得到图5。
[0057] 结论：
[0058] (1)由从图4和图5可以看出，在相同条件下，与已有的预编码码本方案相比，本发明可以在不增加反馈开销，也就是相同码本大小的情况下，进一步降低系统的误码率，并且在
信噪比较高的条件下，本发明展现出了更优异的误码率性能。
[0059] (2)从本文中所描述的具体实例可以看出，本发明方法在计算过程中并不涉及类似于Lloyd迭代中的误差门限设置，避免了误差门限值的选取对码本误码率性能的影响，同时在每个信道矩阵H分类之后都更新了质心，使得分类更加准确，使MIMO系统的误码率性能得到进一步的提高，从而增强通信系统的可靠性。
[0060] 尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、
修改、替换和变型，本发明的范围由所附
权利要求及其等同物限定。
