三维大规模MIMO下基于用户划分的信道估计方法

摘要

本发明提出了一种三维大规模MIMO下基于用户划分的信道估计方法，解决了现有技术中因导频污染导致的信道估计精度较低的问题。实施步骤包括：首先构造三维大规模MIMO信道模型；根据给定划分角度将用户划分为中心用户或边缘用户；对中心用户和边缘用户采用软导频复用方法，并将中心用户的信道矩阵投影到边缘区域及相邻小区区域的零空间中，完成上行信道估计。本发明将所有用户划分为中心用户和边缘用户，并在软导频复用的基础上对中心用户进行正交投影。在计算复杂度较低的前提下，有效提升了信道估计的准确性，可用于三维大规模MIMO多小区多用户系统的导频去污染及小区间干扰抑制。

权利要求

1.一种三维大规模MIMO系统下基于用户划分的信道估计方法，其特征在于，包括有以下步骤：(1)信道模型构造：对大规模MIMO多小区多用户系统中的用户，分别构造水平维度和垂直维度的信道矩阵，进而通过矩阵直积，构造三维大规模MIMO信道模型，在三维大规模MIMO信道模型有新增的垂直维度的自由度；

(2)用户划分：将小区中所有用户划分为中心用户和边缘用户，具体划分方法是将三维大规模MIMO信道中用户垂直维度信道协方差矩阵秩为1时的角度定为划分角度，根据每个用户的下倾角角度与垂直方向的划分角度进行对比，对用户进行划分；

(3)导频分配：基站进行基于导频序列的上行信道估计，在导频分配时采用软导频复用方法，所有小区的中心用户使用相同的导频资源；边缘用户使用正交导频资源；

(4)对中心用户进行正交投影：对中心用户进行正交投影，中心用户将其垂直维度信道矩阵投影到边缘区域及相邻小区区域的零空间中；

(5)信道估计：边缘用户根据导频分配的结果进行信道估计，中心用户根据导频分配的结果和正交投影的结果进行信道估计，分别输出边缘用户和中心用户的信道估计结果。

2.根据权利要求1所述的三维大规模MIMO下基于用户划分的信道估计方法，其特征在于，步骤2中所述的用户划分，包括有如下步骤：(2a)确定划分角度：针对大规模MIMO中的one-ring信道模型，根据小区b中用户k的下倾角角度，构造用户垂直维度信道协方差矩阵 并将协方差矩阵的秩 与理想情况下大规模天线信道协方差矩阵的秩rank(RE)进行对比，确定协方差矩阵的秩为1时的用户垂直方向的角度为划分角度θ°；

(2b)划分用户：根据每个用户的下倾角角度θbk，与垂直方向的划分角度θ°进行对比，若

90°≤θbk＜θ°则将用户判定为边缘用户，若θ°≤θbk＜180°则将用户判定为中心用户，完成中心用户和边缘用户的划分。

3.根据权利要求1所述的三维大规模MIMO下基于用户划分的信道估计方法，其特征在于，步骤3中所述的导频分配，包括有如下步骤：(3a)上行信道估计：多小区系统中共有L个小区，c是各个小区的序号，满足c＝1,…,L，Nh是基站端水平天线数，Nv是基站端垂直天线数，Nt是基站端天线总数，满足Nt＝Nh×Nv，以小区b作为参考小区，各个小区用户同时向其发送长度为τ的上行导频序列，此时，小区b基站侧的接收信号 表示为：其中ρP是导频的传输功率，Hbc＝[hbc1,hbc2,…,hbcK]是小区c中所有用户到小区b基站的信道矩阵，K是小区用户数，导频序列 满足ΦΦH＝IK，(·)T表示对矩阵进行转置， 是上行信道加性高斯白噪声矩阵；

小区b的基站对本小区进行基于导频序列的信道估计，得到信道估计矩阵 为：

(3b)软导频复用：软导频复用方法对每个小区中的所有小区中心用户都分配相同的导频序列，而为边缘用户分配与目标小区边缘用户正交的导频序列，通过软导频复用，将小区b基站侧的接收信号 表示为各小区中心用户上行传输信号的和及各小区边缘用户上行传输信号的和：其中，等号右边第一项表示各小区中心用户上行传输信号的和，等号右边第二项表示各小区边缘用户上行传输信号的和， 是小区c中所有中心用户到小区b基站的信道矩阵， 是中心用户的导频， 是小区c中所有边缘用户到小区b基站的信道矩阵，是边缘用户的导频，KC是小区中心用户数，KE是小区边缘用户数，满足K＝KC+KE及Φ＝[ΦCT,ΦET]T， 是小区b接收端对应的加性高斯白噪声，所有小区中心用户的导频资源相同，边缘用户的导频资源相互正交。

4.根据权利要求1所述的三维大规模MIMO下基于用户划分的信道估计方法，其特征在于，步骤4中所述的对中心用户进行正交投影，包括有如下步骤：(4a)选取干扰空间：对小区b的中心用户来说，其干扰空间是垂直维度角度在(π/2,θ°)之间的边缘区域及相邻小区区域对应的垂直维度信道空间；

(4b)中心用户正交投影：对中心用户在干扰空间信道协方差矩阵非零特征值对应的特征向量进行正交投影，使得中心用户的信道矩阵在所有干扰空间非零特征值对应的特征向量的零空间中。

5.根据权利要求4所述的三维大规模MIMO下基于用户划分的信道估计方法，其特征在于，步骤(4b)中所述的中心用户正交投影，包括有如下步骤：

4b1)奇异值分解：对干扰空间的垂直维度信道协方差矩阵 进行奇异值分

解：

其中 是Nv×Nv阶酉矩阵， 是半正定Nv×Nv阶对角矩阵， 是Nv×Nv阶酉矩阵，(·)H表示对矩阵进行共轭转置；

4b2)干扰空间获取：半正定对角矩阵 对角线上的元素 即为信道协方差矩阵的奇异值，取 中非零元素所在列的列数，对应于 中的列，即干扰空间为信道协方差矩阵非零特征值对应的特征向量

4b3)中心用户在干扰空间正交投影：对小区中心用户的信道矩阵，对其在干扰空间信道协方差矩阵非零特征值对应的特征向量 进行正交投影，使得中心用户的信道矩阵在所有非零特征值对应的特征向量的零空间中。

6.根据权利要求1所述的三维大规模MIMO下基于用户划分的信道估计方法，其特征在于，步骤5中所述的信道估计，包括有如下步骤：(5a)中心用户信道估计：通过软导频复用得到的小区b中心用户的信道估计矩阵：令 则最终的中心用户信道估计矩阵为：

其中， 是小区b中所有中心用户到小区b基站的信道矩阵， 是将Nt×1维的向量表示成Nh×Nv维的矩阵，vec(·)表示向量矩阵化， 是Nv×Nv维的单位矩阵， 是对进行垂直维度正交投影后得到的小区b中心用户的信道估计矩阵；

(5b)边缘用户信道估计：由于软导频复用对边缘用户采用正交导频，因此可得到的小区b中边缘用户的信道估计矩阵为：其中， 是小区b中所有边缘用户到小区b基站的信道矩阵，根据小区中心用户和边缘用户的信道估计矩阵，完成对小区中心用户和边缘用户的信道估计。

说明书

三维大规模MIMO下基于用户划分的信道估计方法

技术领域

[0001] 本发明属于无线通信技术领域，特别涉及大规模MIMO信道估计，具体是一种三维大规模MIMO下基于用户划分的信道估计方法，用于降低三维大规模MIMO多小区多用户系统导频污染的影响。

背景技术

[0002] 作为下一代5G通信系统的关键技术之一，大规模MIMO技术引起了广泛的关注和研究。大规模MIMO技术可以充分利用空间自由度，显著地提升频谱效率和数据链路的可靠性。但大规模MIMO多小区系统中存在严重的导频污染问题，对小区中用户特别是小区边缘用户的信道估计产生严重的影响，是制约大规模MIMO系统容量性能的重要因素。

[0003] 为体现大规模MIMO系统的高频谱效率和高能量效率等优点并对抗导频污染带来的影响，首先需要获取准确的信道状态信息。信道估计的结果决定了后续接收端操作的准确性，进而决定了信号传输质量和系统整体性能。常见的信道估计算法根据先验信息的不同，可分为基于导频序列的信道估计、盲信道估计和半盲信道估计三种。盲信道估计和半盲信道估计需要较长的观测数据，计算复杂度较高且可能会出现相位模糊等问题。考虑到大规模MIMO系统的高复杂度和实际应用中的可行性，在大规模MIMO中主要采用传统的最小二乘信道估计方法等基于导频序列的信道估计方法。

[0004] 现有的大规模MIMO多小区系统中，小区中心用户和边缘用户受干扰情况不同，为对受导频污染严重的边缘用户的干扰进行适当的管理和抑制，有采用用户划分的方式，将用户分为小区中心用户和边缘用户，并在此基础上对小区中心用户和边缘用户分别进行信道估计。用户划分方法主要包括门限固定的基于用户水平位置或垂直角度的划分方法，和动态门限的用户划分方法等。门限固定的基于用户水平位置的划分方法是根据水平方向上用户与基站之间距离的某个经验值作为划分门限，将用户和基站之间的距离与门限值进行比较，距离小于门限值的判定为中心用户，大于门限值的判定为边缘用户。该方法计算复杂度较低，但不能很好的利用三维信道空间的用户信息；门限固定的基于用户垂直角度的划分方法是根据垂直方向下倾角角度的某个经验值作为划分门限，可以较好的利用三维空间信道信息，但其性能依赖于经验值的准确性；动态门限的用户划分方法其门限受系统中某种变化因素的影响而动态改变，可以根据系统实际情况较为准确的进行用户划分，但由于需要不断更新门限其计算复杂度和反馈开销较高。

[0005] 上述用户划分方法中，门限固定的用户划分方法其划分结果不能适应系统天线数等参数的变化。动态门限的用户划分方法其门限取值与用户变化情况有关，计算复杂度和反馈开销较高。

发明内容

[0006] 本发明的目的在于针对上述现有技术的不足，提出一种低复杂度且划分方法适应天线数变化的三维大规模MIMO下基于用户划分的信道估计方法。

[0007] 本发明是一种三维大规模MIMO系统下基于用户划分的信道估计方法，其特征在于，包括有以下步骤：

[0008] (1)信道模型构造：对大规模MIMO多小区多用户系统中的用户，分别构造水平维度和垂直维度的信道矩阵，进而通过矩阵直积，构造三维大规模MIMO信道模型，在三维大规模MIMO信道模型有新增的垂直维度的自由度；

[0009] (2)用户划分：将小区中所有用户划分为中心用户和边缘用户，具体划分方法是将三维大规模MIMO信道中用户垂直维度信道协方差矩阵秩为1时的角度定为划分角度，根据每个用户的下倾角角度与垂直方向的划分角度进行对比，对用户进行划分；

[0010] (3)导频分配：基站进行基于导频序列的上行信道估计，在导频分配时采用软导频复用方法，所有小区的中心用户使用相同的导频资源；边缘用户使用正交导频资源，保证边缘用户不受干扰；

[0011] (4)对中心用户进行正交投影：对中心用户进行正交投影，中心用户将其垂直维度信道矩阵投影到边缘区域及相邻小区区域的零空间中，完成中心用户的正交投影，保证中心用户不受干扰；

[0012] (5)信道估计：边缘用户根据导频分配的结果进行信道估计，中心用户根据导频分配的结果和正交投影的结果进行信道估计，分别输出边缘用户和中心用户的信道估计结果。

[0013] 本发明给出了一种基于用户划分的信道估计方案，可用于降低三维大规模MIMO多小区多用户系统导频污染的影响。

[0014] 本发明与现有技术相比具有如下优点：

[0015] 适应性高：对于门限固定的用户划分方法，门限固定的基于用户水平位置的划分方法划分用户时仅利用了水平维度的空间自由度，而本发明利用了三维空间中新增的垂直维度的自由度；门限固定的基于用户垂直位置的划分方法其划分结果不能适应系统天线数等参数的变化，而本发明能够适应垂直天线数不断增加的三维大规模MIMO系统。

[0016] 复杂度低：对于动态门限的用户划分方法，在大规模MIMO多小区多用户系统中，一方面存在用户移动性较强的特点，与用户变化情况有关的动态门限需瞬时更新，计算复杂度较高，另一方面存在天线数较多、用户数较多的特点，与天线数和用户数有关的反馈信道矩阵维度较高，反馈开销较大，因此动态门限的用户划分方法计算复杂度和反馈开销均较高。本发明在用户划分时门限取值与基站端垂直天线数有关，一方面基站端天线变化情况不会瞬时更新，计算复杂度较低，另一方面用户划分门限取值与基站端垂直天线数有关而与用户变化情况无关，且无需用户端向基站端反馈，降低了反馈开销，因此本发明计算复杂度和反馈开销均较低。

[0017] 消除导频污染：对于不进行用户划分的信道估计方法，各个小区边缘用户使用相同的导频资源，在上行传输时基站端无法区分不同小区的用户，造成导频污染和小区间干扰。本发明采用用户划分方法得到小区中心用户和边缘用户，且在分配导频时，对小区边缘用户使用正交的导频资源，保证边缘用户不受干扰，对小区中心用户使用相同的导频资源，且中心用户将其垂直维度信道矩阵投影到边缘区域及相邻小区区域的零空间中，保证中心用户不受干扰，消除了导频污染的影响，有效提升了信道估计的准确性。

附图说明

[0018] 图1为本发明信道估计的实现流程图；

[0019] 图2为本发明的系统模型示意图；

[0020] 图3为本发明的软导频复用示意图；

[0021] 图4为本发明的垂直方向的秩随下倾角角度的变化仿真图；

[0022] 图5为本发明与现有技术的中心用户信道估计性能的仿真对比图；

[0023] 图6为本发明与现有技术的边缘用户信道估计性能的仿真对比图。

具体实施方式

[0024] 下面结合附图对本发明实施方式和效果详细说明：

[0025] 实施例1

[0026] 5G由于其高带宽、低时延、广连接的特点，己经成为目前业界研究的焦点。大规模MIMO技术通过天线数的增加带来了数据传输速率和能量效率的提升，成为下一代5G通信系统的关键技术之一。导频污染是影响大规模MIMO多小区多用户系统性能的主要因素，为降低导频污染影响，有采用用户划分的方式，将用户分为小区中心用户和边缘用户，并在此基础上对小区中心用户和边缘用户分别进行信道估计。常用的用户划分方法中，门限固定的用户划分方法其划分结果不能适应系统天线数等参数的变化，动态门限的用户划分方法其门限取值与用户变化情况有关，计算复杂度和反馈开销较高。针对现有的用户划分方法无法适用于垂直天线数不断增加的三维大规模MIMO场景的问题，本发明提出一种三维大规模MIMO下基于用户划分的信道估计方法。

[0027] 本发明是一种三维大规模MIMO系统下基于用户划分的信道估计方法，估计过程中涉及到基站、用户、小区，参见图1，包括有以下步骤：

[0028] (1)信道模型构造：对大规模MIMO多小区多用户系统中的用户，分别构造水平维度和垂直维度的信道矩阵，进而通过矩阵直积，构造三维大规模MIMO信道模型，在三维大规模MIMO信道模型有新增的垂直维度的自由度。对大规模MIMO多小区多用户系统中的用户，本例中小区数为2，基站端水平天线数为16，垂直天线数为4，即各基站总发送天线数为64，各用户接收天线数为1。将有关参数代入大规模MIMO中的one-ring信道模型，得到水平维度和垂直维度的信道矩阵，进而通过矩阵直积得到三维大规模MIMO信道模型，这种模型能够适应大规模MIMO多小区多用户系统。

[0029] (2)用户划分：将小区中所有用户划分为中心用户和边缘用户，具体划分方法是将三维大规模MIMO信道中用户垂直维度信道协方差矩阵秩为1时的角度定为划分角度，根据每个用户的下倾角角度与垂直维度的划分角度进行对比，若用户的下倾角角度小于垂直维度的划分角度则将用户判定为边缘用户，若用户的下倾角角度大于垂直维度的划分角度则将用户判定为中心用户，完成中心用户和边缘用户的划分。本例中将天线数等有关参数代入大规模MIMO中的one-ring信道模型，得到的划分角度为120°，对各用户若其下倾角角度小于120°则将用户判定为边缘用户，若用户的下倾角角度大于120°则将用户判定为中心用户，完成中心用户和边缘用户的划分。本发明中天线数给定时得到的划分角度固定，随着天线数的变化得到的划分角度也不断变化。

[0030] (3)导频分配：基站进行基于导频序列的上行信道估计，在导频分配时采用软导频复用方法，所有小区的中心用户均使用相同的导频资源。所有边缘用户均使用正交导频资源，保证边缘用户不受干扰。本例中将小区数等有关参数代入软导频复用方法，得到本例中的导频分配方法为将导频集划分为三个正交子集即P＝{P0,P1,P2}，将导频集P0分配给各个小区的中心用户，将导频集P1、P2分别分配给各个小区的边缘用户。各个小区边缘用户的导频序列相互正交，保证边缘用户不受干扰。

[0031] (4)对中心用户进行正交投影：本发明专门对中心用户进行正交投影，中心用户将其垂直维度信道矩阵投影到边缘区域及相邻小区区域的零空间中，完成中心用户的正交投影，保证中心用户不受干扰。本例中针对大规模MIMO中的one-ring信道模型，根据垂直维度下倾角角度范围(π/2,120°)，构造边缘区域及相邻小区区域的垂直维度信道矩阵，对矩阵进行奇异值分解，干扰空间为信道协方差矩阵非零特征值对应的特征向量，中心用户在干扰空间正交投影，完成中心用户的正交投影，保证中心用户不受干扰。

[0032] (5)信道估计：边缘用户根据导频分配的结果进行信道估计；中心用户根据导频分配的结果和正交投影的结果进行信道估计，分别输出边缘用户和中心用户的上行信道估计结果。在下行处理中，时分双工系统具有信道互易性，直接根据上行信道估计结果进行下行预编码以消除小区间干扰。频分双工系统不具有信道互易性，需采用本例中的信道估计方法进行信道估计，并根据信道估计的结果进行下行预编码以消除小区间干扰，提升系统和速率。

[0033] 本发明的技术思路是：首先对大规模MIMO多小区多用户系统中的用户，分别构造水平维度和垂直维度的信道矩阵，进而通过矩阵直积，构造三维大规模MIMO信道模型，这种模型能够适应大规模MIMO多小区多用户系统；然后根据用户下倾角信息将其划分为中心用户或边缘用户；通过软导频复用，所有小区的中心用户使用相同的导频资源，并在信道估计时，将估计到的信道矩阵投影到边缘区域及相邻小区区域的零空间中，避免对边缘区域的干扰并提升估计精度；所有边缘用户使用正交导频资源，以使其信道估计精度不受导频污染的影响。

[0034] 本发明给出了一个整体的信道估计技术方案，用户划分标准与垂直方向上的角度有关，且该角度的取值随垂直方向上的天线数变化而动态改变，解决了现有用户划分方法中计算复杂度高且划分结果不能适应系统天线数等参数的变化的技术问题。在用户划分的基础上，通过对中心用户分配相同的导频资源并在干扰空间进行正交投影，对边缘用户分配相互正交的导频资源，并对中心用户和边缘用户进行基于导频的信道估计，能够有效提升信道估计的准确性，解决了现有技术中因导频污染导致的信道估计精度较低的技术问题，有效降低导频污染的影响。

[0035] 实施例2

[0036] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1，步骤2中所述的用户划分，包括有如下步骤：

[0037] (2a)确定划分角度：针对大规模MIMO中的one-ring信道模型，根据小区b中用户k的下倾角角度，构造用户垂直维度信道协方差矩阵 并将协方差矩阵的秩与理想情况下大规模天线信道协方差矩阵的秩rank(RE)进行对比，确定协方差矩阵的秩为
1时的用户垂直方向的角度为划分角度θ0。本发明通过实验数据也证明取用户垂直维度信道协方差矩阵秩为1时的下倾角角度作为用户划分角度是可行的。

[0038] (2b)划分用户：根据每个用户的下倾角角度θbk，与垂直方向的划分角度θ0进行对比，若90°≤θbk＜θ0则将用户判定为边缘用户，若θ0≤θbk＜180则将用户判定为中心用户，完成中心用户和边缘用户的划分。本发明选择这样的用户划分方式，是因为边缘区域垂直方向上的下倾角范围本身就很小。三维大规模MIMO系统相较于二维大规模MIMO系统多出一个垂直方向上的维度，在带来空间自由度扩张的同时也增加了信道估计复杂度和反馈开销。选取垂直方向信道秩为1时的下倾角角度作为区域划分角度，在时分双工系统下可以降低信道维度进而降低预编码运算的复杂度，在频分双工系统下边缘区域用户只需要估计水平方向这一个维度的信道，降低了信道估计的计算复杂度和反馈开销。本发明通过实验数据也证明取用户垂直维度信道协方差矩阵秩为1时的下倾角角度作为用户划分角度是可行的。

[0039] 实施例3

[0040] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-2，步骤3中所述的导频分配，包括有如下步骤：

[0041] (3a)上行信道估计：参见图2，多小区系统中共有L个小区，c是各个小区的序号，满足c＝1,…,L，Nh是基站端水平天线数，Nv是基站端垂直天线数，Nt是基站端天线总数，满足Nt＝Nh×Nv，以小区b作为参考小区，各个小区用户同时向其发送长度为τ的上行导频序列，此时，小区b基站侧的接收信号 表示为：

[0042]

[0043] 其中ρP是导频的传输功率，Hbc＝[hbc1,hbc2,…,hbcK]是小区c中所有用户到小区b基站的信道矩阵，K是小区用户数，导频序列 满足ΦΦH＝IK，(·)T表示对矩阵进行转置， 是上行信道加性高斯白噪声矩阵。参见图2，图2中小区数L＝3，每个小区中各有一个基站，θ表示用户下倾角，基站侧不仅接收到来自本小区用户的上行传输信号，也接收到来自其他小区用户的上行传输信号。

[0044] 小区b的基站对本小区进行基于导频序列的信道估计，得到信道估计矩阵 为：

[0045]

[0046] (3b)软导频复用：软导频复用方法对每个小区中的所有小区中心用户都分配相同的导频序列，而为边缘用户分配与目标小区边缘用户正交的导频序列，通过软导频复用，将小区b基站侧的接收信号 表示为各小区中心用户上行传输信号的和及各小区边缘用户上行传输信号的和：

[0047]

[0048] 其中，等号右边第一项表示各小区中心用户上行传输信号的和，等号右边第二项表示各小区边缘用户上行传输信号的和， 是小区c中所有中心用户到小区b基站的信道矩阵， 是中心用户的导频， 是小区c中所有边缘用户到小区b基站的信道矩阵， 是边缘用户的导频，KC是小区中心用户数，KE是小区边缘用户数，满足
是小区b接收端对应的加性高斯白噪声，所有小区中
心用户的导频资源相同，边缘用户的导频资源相互正交。本发明选择这样的导频分配方式，是因为使用相同导频资源的各边缘用户在上行传输时基站端无法区分不同小区的用户，导致导频污染问题，降低系统和速率。而各个小区中心用户所在区域距离较远，考虑到大尺度衰落等路径损耗因素，各个小区中心用户之间的干扰较小，且为了降低导频开销，因此使用相同的导频资源。本发明通过实验数据也证明本发明所提方法能够降低边缘用户信道估计误差。

[0049] 实施例4

[0050] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-3，步骤4中所述的对中心用户进行正交投影，包括有如下步骤：

[0051] (4a)选取干扰空间：对小区b的中心用户来说，其干扰空间是垂直维度角度在(π/2,θ0)之间的边缘区域及相邻小区区域对应的垂直维度信道空间。

[0052] (4b)中心用户正交投影：对中心用户在干扰空间信道协方差矩阵非零特征值对应的特征向量进行正交投影，使得中心用户的信道矩阵在所有干扰空间非零特征值对应的特征向量的零空间中。本发明通过软导频复用对不同小区的边缘用户分配相互正交的导频，保证边缘用户不受干扰。而此时对于中心用户，一方面由于划分用户时，大规模MIMO中one-ring信道模型的角度扩展导致下倾角角度在划分角度附近的中心用户和边缘用户之间产生干扰，因此中心用户会受到边缘用户的干扰。另一方面由于各个小区中心用户使用相同的导频资源，因此中心用户会受到其他小区中心用户的干扰。从而有必要对于边缘区域及相邻小区区域对应的干扰空间需要通过正交投影，保证中心用户不受干扰。本发明通过实验数据也证明本发明所提方法能够降低中心用户信道估计误差。

[0053] 实施例5

[0054] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-4，步骤(4b)中所述的中心用户正交投影，包括有如下步骤：

[0055] 4b1)奇异值分解：对干扰空间的垂直维度信道协方差矩阵 进行奇异值分解：

[0056]

[0057] 其中 是Nv×Nv阶酉矩阵， 是半正定Nv×Nv阶对角矩阵， 是Nv×Nv阶酉矩阵，(·)H表示对矩阵进行共轭转置。

[0058] 4b2)干扰空间获取：半正定对角矩阵 角线上的元素 即为信道协方差矩阵的奇异值，取 中非零元素所在列的列数，对应于 中的列，即干扰空间为信道协方差矩阵非零特征值对应的特征向量

[0059] 4b3)中心用户在干扰空间正交投影：对小区中心用户的信道矩阵，对其在干扰空间信道协方差矩阵非零特征值对应的特征向量 进行正交投影，使得中心用户的信道矩阵在所有非零特征值对应的特征向量的零空间中。本发明通过中心用户在干扰空间正交投影，保证中心用户不受干扰，提升中心用户信道估计的准确性。

[0060] 实施例6

[0061] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-5，步骤5中所述的信道估计，包括有如下步骤：

[0062] (5a)中心用户信道估计：通过软导频复用得到的小区b中心用户的信道估计矩阵：

[0063]

[0064] 令 则最终的中心用户信道估计矩阵为：

[0065]

[0066] 其中， 是小区b中所有中心用户到小区b基站的信道矩阵， 是将Nt×1维的向量 表示成Nh×Nv维的矩阵，vec(·)表示向量矩阵化， 是Nv×Nv维的单位矩阵， 是对 进行垂直维度正交投影后得到的小区b中心用户的信道估计矩阵；

[0067] (5b)边缘用户信道估计：由于软导频复用对边缘用户采用正交导频，因此可得到的小区b中边缘用户的信道估计矩阵为：

[0068]

[0069] 其中， 是小区b中所有边缘用户到小区b基站的信道矩阵，根据小区中心用户和边缘用户的信道估计矩阵，完成对小区中心用户和边缘用户的信道估计。本发明得到不受干扰的小区中心用户的信道估计和不受干扰的边缘用户的信道估计，降低信道估计的误差，保证后续预编码处理的准确性，提升系统和速率。

[0070] 下面给出一个更加详细的例子，对本发明进一步说明

[0071] 实施例7

[0072] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-6，

[0073] 下面结合附图对本发明实施方式和效果作进一步说明：

[0074] 参见图2，本发明的系统模型包含多个六边形小区，基站配置二维有源天线阵列，每个小区多个用户，本发明研究基站对用户的信道估计。

[0075] 参见图1，实现步骤如下：

[0076] 步骤1，信道模型构造：对大规模MIMO多小区多用户系统中的用户，分别构造水平维度和垂直维度的信道矩阵，进而通过矩阵直积，构造三维大规模MIMO信道模型，在三维大规模MIMO信道模型有新增的垂直维度的自由度。

[0077] 本步骤的具体实现如下：

[0078] (1a)构造系统模型：参见图2，假设该多小区多用户系统有L个六边形蜂窝小区，每个小区用户数为K。每个小区基站水平方向天线数为Nh，垂直方向天线数为Nv，即总的基站天线数为Nt＝Nh×Nv，每个用户接收天线数为1。本例中小区数L＝3，用户数K＝20，每个小区基站水平方向天线数为16，垂直方向天线数为8，即总的基站天线数为128，每个用户接收天线数为1，每个小区中各有一个基站，θ表示用户下倾角，基站侧不仅接收到来自本小区用户的上行传输信号，也接收到来自其他小区用户的上行传输信号。

[0079] 采用适用于大规模MIMO场景的信道模型如one-ring信道等，分别作为水平和垂直方向上的相关矩阵。水平和垂直方向上任意两天线(下标为n1,n2)间的相关性为：

[0080]

[0081]

[0082] 对其进行克罗内克积操作，可生成小区b与小区c中的用户k之间的信道相关矩阵Rbck如下：

[0083]

[0084] 其中，φbck,max和φbck,min分别表示最大方位角和最小方位角，θbck,max和θbck,min分别表示最大下倾角和最小下倾角，λ是信号波长， 是天线间距， 表示克罗内克积。

[0085] (1b)构造大规模MIMO三维信道：对相关矩阵进行奇异值分解，得到水平方向相关矩阵 及垂直方向相关矩阵

[0086]

[0087]

[0088] 由Karhunen-Loeve变换得到小区b与小区c中的用户k之间的信道矩阵的表达式：

[0089]

[0090]

[0091] 其中wbck独立同分布的服从均值为0，方差为1的高斯分布，βbck是大尺度衰落系数。

[0092] 步骤2，用户划分：本发明将小区中所有用户划分为中心用户和边缘用户，具体划分方法是将三维大规模MIMO信道中用户垂直维度信道协方差矩阵秩为1时的角度定为划分角度，根据每个用户的下倾角角度与垂直方向的划分角度进行对比，对用户进行划分。

[0093] 本步骤的具体实现如下：

[0094] (2a)确定划分角度：

[0095] 以one-ring信道模型为例，理想情况下大规模天线信道相关矩阵的秩rank(RE)的公式如下：

[0096]

[0097] 其中θmin是下倾角的最小值，θmax是下倾角的最大值，Nv是垂直方向天线数，d是天线间距，λ是信号波长。

[0098] 小区b中的用户k垂直维度信道协方差矩阵 的公式如下：

[0099]

[0100] 其中θbbk,min是用户k下倾角的最小值，θbbk,max是用户k下倾角的最大值， 是天线间距。

[0101] 根据垂直维度信道协方差矩阵求得对应的的秩 与理想情况下大规模天线信道相关矩阵的秩rank(RE)进行对比，确定协方差矩阵的秩为1时的垂直方向的角度
0
为划分角度θ。本发明通过实验数据也证明取用户垂直维度信道协方差矩阵秩为1时的下倾角角度作为用户划分角度是可行的。

[0102] (2b)划分用户：根据每个用户的下倾角角度θbk，与垂直方向的划分角度θ0进行对比，若90°≤θbk＜θ0则将用户判定为边缘用户，若θ0≤θbk＜180°则将用户判定为中心用户，完成中心用户和边缘用户的划分。

[0103] 步骤3，导频分配：基站进行基于导频序列的上行信道估计，在导频分配时采用软导频复用方法，所有小区的中心用户使用相同的导频资源；边缘用户使用正交导频资源，保证边缘用户不受干扰。

[0104] 本步骤的具体实现如下：

[0105] (3a)上行信道估计：多小区系统中共有L个小区，c是各个小区的序号，满足c＝1,…,L，Nh是基站端水平天线数，Nv是基站端垂直天线数，Nt是基站端天线总数，满足Nt＝Nh×Nv，以小区b作为参考小区，各个小区用户同时向其发送长度为τ的上行导频序列，此时，小区b基站侧的接收信号 表示为：

[0106]

[0107] 其中ρP是导频的传输功率，Hbc＝[hbc1,hbc2,…,hbcK]是小区c中所有用户到小区b基站的信道矩阵，K是小区用户数，导频序列 满足ΦΦH＝IK，(·)T表示对矩阵进行转置， 是上行信道加性高斯白噪声矩阵。

[0108] 小区b的基站对本小区进行基于导频序列的信道估计，得到信道估计矩阵 为：

[0109]

[0110] 可以看出，小区b的信道估计结果会受到其他小区中使用相同导频的用户信道的干扰，即存在导频污染问题。针对小区中心用户受导频污染较轻，而小区边缘用户受严重的导频污染的情况，采用软导频复用策略，对小区中心用户及小区边缘用户采取不同的导频分配方法。

[0111] (3b)软导频复用：

[0112] 软导频复用策略对每个小区中的所有小区中心用户都复用相同的导频序列，而为边缘用户分配与目标小区边缘用户正交的导频序列。参见图3，以三小区软导频复用为例，对所有小区中心区域导频复用因子设为1，对所有边缘区域导频复用因子设为1/3。将导频集划分为四个正交子集即P＝{P0,P1,P2,P3}，将导频集P0分配到三个小区的中心位置附近，将其余三个导频集分别分配到三个小区的边缘位置附近。对于相邻小区的边缘交界处，不同用户的导频始终保持正交性，这将大大降低小区之间的导频污染程度。

[0113] 通过软导频复用，将小区b基站侧的接收信号 表示为各小区中心用户上行传输信号的和及各小区边缘用户上行传输信号的和：

[0114]

[0115] 其中，等号右边第一项表示各小区中心用户上行传输信号的和，等号右边第二项表示各小区边缘用户上行传输信号的和， 是小区c中所有中心用户到小区b基站的信道矩阵， 是中心用户的导频， 是小区c中所有边缘用户到小区b基站的信道矩阵， 是边缘用户的导频，KC是小区中心用户数，KE是小区边缘用户数，满足K＝KC+KE及Φ＝[ΦCT,ΦET]T， 是小区b接收端对应的加性高斯白噪声，所有小区中心用户的导频资源相同，边缘用户的导频资源相互正交。本发明通过实验数据也证明本发明所提方法能够降低边缘用户信道估计误差。

[0116] 步骤4，对中心用户进行正交投影：对中心用户进行正交投影，中心用户将其垂直维度信道矩阵投影到边缘区域及相邻小区区域的零空间中，完成中心用户的正交投影。

[0117] 本步骤的具体实现如下：

[0118] (4a)选取干扰空间：对小区b的中心用户来说，其干扰空间是垂直维度角度在(π/2,θ0)之间的边缘区域及相邻小区区域对应的垂直维度信道空间。

[0119] (4b)中心用户正交投影：对中心用户在干扰空间信道协方差矩阵非零特征值对应的特征向量进行正交投影，使得中心用户的信道矩阵在所有干扰空间非零特征值对应的特征向量的零空间中。

[0120] 4b1)奇异值分解：对干扰空间的垂直维度信道协方差矩阵 进行奇异值分解：

[0121]

[0122] 其中 是Nv×Nv阶酉矩阵， 是半正定Nv×Nv阶对角矩阵， 是Nv×Nv阶酉矩阵，(·)H表示对矩阵进行共轭转置。

[0123] 4b2)干扰空间获取：半正定对角矩阵 对角线上的元素 即为信道协方差矩阵 的奇异值，取 中非零元素所在列的列数，对应于 中的列，即干扰空间为
信道协方差矩阵非零特征值对应的特征向量

[0124] 4b3)中心用户在干扰空间正交投影：对小区中心用户的信道矩阵，对其在干扰空间信道协方差矩阵非零特征值对应的特征向量 进行正交投影，使得中心用户的信道矩阵在所有非零特征值对应的特征向量的零空间中。本发明通过实验数据也证明本发明所提方法能够降低中心用户信道估计误差。

[0125] 步骤5，信道估计：边缘用户根据导频分配的结果进行信道估计，中心用户根据导频分配的结果和正交投影的结果进行信道估计，分别输出边缘用户和中心用户的信道估计结果。

[0126] 本步骤的具体实现如下：

[0127] (5a)中心用户信道估计：通过软导频复用得到的小区b中心用户的信道估计矩阵：

[0128]

[0129] 令 则最终的中心用户信道估计矩阵为：

[0130]

[0131] 其中， 是小区b中所有中心用户到小区b基站的信道矩阵， 是将Nt×1维的向量 表示成Nh×Nv维的矩阵，vec(·)表示向量矩阵化， 是Nv×Nv维的单位矩阵， 是对进行垂直维度正交投影后得到的小区b中心用户的信道估计矩阵；

[0132] (5b)边缘用户信道估计：由于软导频复用对边缘用户采用正交导频，因此可得到的小区b中边缘用户的信道估计矩阵为：

[0133]

[0134] 其中， 是小区b中所有边缘用户到小区b基站的信道矩阵，根据小区中心用户和边缘用户的信道估计矩阵，完成对小区中心用户和边缘用户的信道估计。

[0135] 下面通过仿真及其分析，对本发明的技术效果再做说明

[0136] 实施例8

[0137] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-7，

[0138] 本发明的效果可通过仿真进一步说明：

[0139] 仿真条件

[0140] 仿真使用MATLAB R2015b软件实现，设置发送水平天线数为16，垂直天线数为8，即各基站总发送天线数为128，各用户接收天线数为1，小区数为3，小区半径为500m，基站高度为35m，每小区用户数为20等参数。

[0141] 仿真内容和结果

[0142] 仿真1.在不同垂直天线数下，对本发明的垂直方向的秩随下倾角角度的变化情况进行仿真，结果如图4所示，图4为本发明的垂直方向的秩随下倾角角度的变化仿真图。

[0143] 图4横坐标为垂直方向下倾角角度，纵坐标为垂直方向信道矩阵的秩。带叉虚线表示理论垂直方向8天线下大规模MIMO信道矩阵秩的曲线，虚线表示理论垂直16方向天线下大规模MIMO信道矩阵秩的曲线，点划线表示理论垂直方向24天线下大规模MIMO信道矩阵秩的曲线，带圈虚线表示理论垂直32方向天线下大规模MIMO信道矩阵秩的曲线，带叉实线表示实际垂直方向8天线下信道矩阵的秩的曲线，带圈实线表示实际垂直方向8天线下信道矩阵的秩的曲线。本发明不仅做了方案的设计和理论的分析，而且还进行了实验与探讨。从图4可见，理论垂直方向8天线下大规模MIMO信道矩阵秩的曲线和实际垂直方向8天线下信道矩阵的秩的曲线在纵坐标为1时存在交点，从而该交点处对应的横坐标值即为秩为1时对应的下倾角角度。理论垂直方向32天线下大规模MIMO信道矩阵秩的曲线和实际垂直方向32天线下信道矩阵的秩的曲线在纵坐标为1时存在交点，从而该交点处对应的横坐标值即为秩为1时对应的下倾角角度。由此得出结论，随着下倾角度的增大，不同垂直天线数下实际垂直方向信道矩阵的秩与理论秩的曲线均存在交点，从而取各天线下秩为1时对应的下倾角角度作为用户划分角度是可行的。

[0144] 本发明不仅做了方案的设计和理论的分析，而且还进行了实验与探讨。

[0145] 实施例9

[0146] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-7，仿真数据和内容同实施例8

[0147] 仿真2.在上述仿真条件下，对本发明与最小二乘信道估计方法和软导频复用信道估计方法的中心用户信道估计误差进行仿真，结果如图5所示，图5为本发明与现有技术的中心用户信道估计性能的仿真对比图。

[0148] 图5横坐标为信噪比，纵坐标为中心用户信道估计归一化均方误差。带圈实线表示中心用户采用本发明所提方法的信道估计误差曲线，菱形实线表示中心用户采用最小二乘信道估计方法的曲线，十字实线表示中心用户采用软导频复用信道估计方法的曲线。本发明不仅做了方案的设计和理论的分析，而且还进行了实验与探讨。从图5可见，中心用户采用本发明所提方法的信道估计误差可达10-4，而中心用户采用最小二乘信道估计方法和中心用户采用软导频复用信道估计方法的信道估计误差仅能达到10-3，信道估计误差越小，对信道估计的准确性越高，误差减少一个数量级，可保证后续接收端操作的精确性大大提升。由此得出结论，随着信噪比的增大，本发明和最小二乘信道估计方法的归一化均方误差都在逐渐减小，由于本发明中的中心用户进行了正交投影，所提方法的中心用户不受边缘区域和相邻小区的干扰，因此比起最小二乘信道估计方法和软导频复用信道估计方法的信道估计误差低。

[0149] 实施例10

[0150] 三维大规模MIMO下基于用户划分的信道估计方法同实施例1-7，仿真数据和内容同实施例8

[0151] 仿真3.在上述仿真条件下，对本发明与最小二乘信道估计方法和软导频复用信道估计方法的的边缘用户信道估计误差进行仿真，结果如图6所示，图6为本发明与现有技术的边缘用户信道估计性能的仿真对比图。

[0152] 图6横坐标为信噪比，纵坐标为边缘用户信道估计归一化均方误差。带圈实线表示边缘用户采用本发明所提方法的信道估计误差曲线，菱形实线表示边缘用户采用最小二乘信道估计方法的曲线，十字实线表示边缘用户采用软导频复用信道估计方法的曲线。本发明不仅做了方案的设计和理论的分析，而且还进行了实验与探讨。从图6可见，边缘用户采用本发明所提方法的信道估计误差和边缘用户采用软导频复用信道估计方法的信道估计误差可达10-3，而边缘用户采用最小二乘信道估计方法的信道估计误差仅能达到10-2，信道估计误差越小，对信道估计的准确性越高，误差减少一个数量级，可保证后续接收端操作的精确性大大提升。由此得出结论，从图6可见，随着信噪比的增大，本发明和最小二乘信道估计方法的归一化均方误差都在逐渐减小，由于本发明中的边缘用户垂直维度秩为1，因此三维信道等效于二维信道，信道估计误差和软导频复用方法性能相同；由于进行软导频复用消除了导频污染，比起受相邻小区干扰的最小二乘信道估计方法相比起估计误差更低。

[0153] 本发明提出了一种三维大规模MIMO下基于用户划分的信道估计方法，解决了现有技术中因导频污染导致的信道估计精度较低的技术问题。实施步骤包括：对大规模MIMO多小区多用户系统中的用户，构造三维大规模MIMO信道模型；然后根据给定划分角度将用户划分为中心用户或边缘用户；对中心用户和边缘用户采用软导频复用方法，并将中心用户的信道矩阵投影到边缘区域及相邻小区区域的零空间中，最后完成上行信道估计。本发明首先对用户进行了划分，得到中心用户和边缘用户，并在软导频复用的基础上对中心用户进行正交投影。在计算复杂度较低的前提下，有效提升了信道估计的准确性，可用于三维大规模MIMO多小区多用户系统的导频去污染及小区间干扰抑制。

[0154] 以上描述仅是本发明的一个具体实例，并不构成对本发明的任何限制。显然对于本领域的专业人员来说，在了解了本发明内容和原理后，都可能在不背离本发明原理、结构的情况下，进行形式和细节上的各种修正和改变，但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。