技术领域
[0001] 本
申请涉及无线通信技术领域信道状态信息反馈方法，特别涉及一种基于PCA演进的大规模MIMO信道状态信息压缩反馈方法。
背景技术
[0002] 多输入多输出(Multiple-Input Multiple-Output，MIMO)技术是指在发射端和接收端分别使用多个发射天线和接收天线，使
信号通过发射端与接收端的多个天线传送和接收，从而改善通信
质量。它能充分利用空间资源，通过多个天线实现多发多收，在不增加
频谱资源和天线发射功率的情况下，可以成倍的提高系统信道容量，被视为下一代移动通信的核心技术。
正交频分复用(Orthogonal Frequency Division Multiplexing，OFDM)是一种有效的调制方案，它将信道分成若干正交子信道，将高速数据信号转换成并行的低速子数据流，调制到在每个子信道上进行传输。接收端采用相关技术来将正交信号分开，减少了子信道之间的相互干扰。每个子信道上的信号带宽小于信道的相关带宽，因此每个子信道上的可以看成平坦性衰落，从而可以消除符号间干扰。
[0003] 移动通信大规模MIMO-OFDM系统中，为实现使用预编码技术将
能量集中到目标移动终端上，一般由接收端将估计出来的信道信息通过反馈信道反馈给发射端，让发送端实时的获得准确的下行信道的信道信息。但随着大规模MIMO技术推进，天线数量会大幅增加，如果将下行链路信道信息直接完全反馈给发射端，会极大增加容量损失，所以对信道信息的压缩反馈方式成为研究的热点。
[0004] 主成分分析(Principal Component Analysis，PCA)，也称主分量分析或矩阵数据分析，一种数据分析的技术，主要思想是将高维数据投影到较低维空间，提取多元事物的主要因素，揭示其本质特征。它可以高效地找出数据中的主要部分，将原有的复杂数据
降维，去除整个数据中的噪音和冗余，将把相关的变量变为若干不相关的综合指标变量，是一种最小均方意义上的最优变换，目的是去除输入随机向量之间的相关性，突出原始数据中的隐含特性。主成分分析方法的优势在于
数据压缩以及对多维数据进行降维，操作简单，且没有参数限制，可以方便的应用于各个场合。因此，将PCA用于大规模MIMO系统的信道状态信息压缩中，可以减少反馈开销。但是，现提出的使用PCA对信道信息压缩的方法大多只考虑了空间域或频域一个域内的相关性，且目前对于主要成分提取的选择矩阵需要动态更新，这无疑从另一方面增加了反馈量并提高了计算复杂度。
[0005] 综上，现在已经出现的大规模MIMO系统中信道状态信息反馈方法，针对它们计算方法复杂度高、反馈量大、反馈
精度较低等问题，现提出基于主成分分析演进的大规模MIMO-OFDM信道信息反馈方法，考虑信道信息在
空域频域两个维度的相关性，结合分簇思想，将大量信道状态信息压缩反馈到发射端并进行重构操作。
发明内容
[0006] 本发明旨在至少解决
现有技术中存在的技术问题，特别创新地提出了一种基于PCA演进的大规模MIMO-OFDM系统信道反馈方法。
[0007] 发明目的：为了实现在大规模MIMO-OFDM系统中，对信道状态信息进行压缩后的反馈，本发明提出了一种基于PCA演进的信道信息压缩方法来反馈高度相关的信道状态信息。该方法首先考虑空域和频域两个域内的信道特性，通过获取下行信道相关特性参数建立信道模型；然后将矢量化后的信道状态信息即信道向量进行分簇，提高信道状态信息的相关性和压缩性，以便对每一簇内的高维信道状态信息降低到低维数据；其次，在收发两端计算特性参数信道模型每一簇内信道状态信息的协方差矩阵以及稀疏矩阵，在接收端，通过选择矩阵提取稀疏向量主要成分，使
选定的稀疏信道向量满足各元素之和最大原则；最后将选定的信道信息及选择矩阵的码字索引反馈到发射端进行信道状态信息的重构。
[0008] 为了实现本发明的上述目的，本发明提供了一种基于PCA演进的大规模MIMO-OFDM系统信道反馈方法，其特征包括：
[0009] S1，获取下行信道相关特性参数，建立信道模型。
[0010] S2，将信道向量hk进行分簇，根据信道特性参数计算出每一簇信道状态信息的协方差矩阵，并通过协方差矩阵的特征值分解获得稀疏矩阵。
[0011] S3，预定义一个收发两端都已知的二进制选择矩阵码本，根据选定的稀疏信道向量各元素之和最大原则选取最优的选择矩阵。
[0012] 所述的基于PCA演进的大规模MIMO-OFDM系统信道反馈方法，其特征在于，所述S1包括：
[0013] 本发明考虑一个大规模MIMO-OFDM系统，发射端为配置Nt根均匀天线的线性阵列，接收端为单天线用户，考虑频域有Nc个
子载波；在大规模MIMO-OFDM系统中，发射端与第k个用户间在空域和频域的信道向量hk表示为：
[0014]
[0015] 其中N＝Nt×Nc，vec(A)表示将一个a×b的矩阵A矢量化，变成ab×1的列向量，hk(n)为第k个用户第n个子载波的空域信道向量，建立模型为：
[0016]
[0017] 在大规模MIMO系统中，由于发射端配置天线数较多时，天线间距较小，天线间呈现较强的空间相关性，因此，发射端与第k个单天线用户间的信道可建模为：
[0018]
[0019] 其中， 表示发送相关矩阵，n表示第n个子载波，n∈{1,…,Nc}；表示均值为0，方差为1的独立分布复高斯随机矢量，Hiid,k(n)的第(u,v)个
元素代表发射天线u和接收天线v之间由路径损耗和小尺度衰落组成的信道增益。
[0020] 由于发射端配置Nt根均匀线性阵列天线，空间相关矩阵可以通过Jakes模型得到，第p根天线与第q根天线之间的相关系数表示为：
[0021]
[0022] 其中J0(·)表示第一类零阶贝塞尔函数，dpq表示第p根天线与第q根天线之间的距离，λ表示载波
波长。
[0023] 所述的基于PCA演进的大规模MIMO-OFDM系统信道反馈方法，其特征在于，所述S2包括：
[0024] S2-1，将信道向量hk进行分簇。
[0025] 接收端将连续的信道状态信息分为G个簇，即hkT＝[hk,1T…hk,GT]，其中G≤N；hk,g表示第g个簇的信道状态信息，用G个簇来代替N个信道状态信息元素，则每一簇包含 个信道状态信息元素。
[0026] S2-2，根据特性参数计算出信道每一簇信道状态信息的协方差矩阵，并通过协方差矩阵的特征值分解获得稀疏矩阵；第k个用户第g簇的信道协方差矩阵表示为：
[0027]
[0028] 对第g簇信道状态信息的协方差矩阵进行特征值分解，获得稀疏矩阵：
[0029]
[0030] 其中， 表示hk,g的协方差矩阵， 表示对协方差矩阵 进行特征值分解，即Ck,g＝Ψk,gΛΨk,gH；Λ是对
角矩阵，其对角元素是协方差矩阵 的特征值，Ψk,g的列向量是协方差矩阵 的
特征向量；通过卡洛南变换，Ψk,g同时表示稀疏矩阵。
[0031] 所述的基于PCA演进的大规模MIMO-OFDM系统信道反馈方法，其特征在于，所述S3包括：
[0032] S3-1，预定义一个收发两端都已知的二进制选择矩阵码本CP：
[0033] CP＝{P1…PL}
[0034] 其中Pi表示选择矩阵码本的第i个码字，即选择矩阵，L表示码本CP包含L个选择矩阵。
[0035] S3-2，选取最优的选择矩阵Pi。
[0036] 通过选择矩阵Pi从稀疏向量sk,g中选择M个最重要的反馈信息，其中
[0037] s′k,g＝Pisk,g
[0038] 全局搜索选择矩阵码本，当选定的稀疏信道向量s′k,g各元素之和最大时，接收端确定使用的选择矩阵Pi。
[0039] i＝argimax|sk,gPi|
[0040] 确定选择矩阵Pi后，通过s′k,g＝Pisk,g从稀疏向量sk,g中选择M个最重要的反馈信息。
[0041] S3-3，采用随机矢量量化码本，量化选定的稀疏信道向量s′k,g，s′k,g码字索引通过以下公式获得：
[0042] j＝argjmax|s′k,gwj|
[0043] 其中wj是随机矢量量化码本的第j个码字；然后将s′k,g码字索引以及选择矩阵Pi的码字索引一起反馈到发射端。
[0044] 综上所述，由于采用了上述技术方案，本发明的有益效果是：
[0045] 根据发射端为大规模MIMO均匀线性阵列，接收端为单天线用户这一场景建立了相关的信道模型，将空域和频域连续的N个信道状态信息进行分簇。首先根据信道特性计算出每一簇信道状态信息的协方差矩阵和稀疏矩阵；然后，在收发两端预定义一个相同的二进制选择矩阵码本CP，根据选定的稀疏信道向量各元素之和最大原则确定选择矩阵，获得选定的稀疏信道向量；最后将选定的稀疏信道向量随机矢量化以及选择矩阵索引反馈到发射端并进行重构操作。本发明将连续的N个信道状态信息进行分簇，提高了信道状态信息相关性及基于PCA的压缩性能；通过信道特性得到信道协方差矩阵和稀疏矩阵以及在收发两端预定义的二进制选择矩阵码本，减少了大量的反馈量和降低了复杂度。
[0046] 本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
[0047] 本发明的上述和/或附加的方面和优点从结合下面附图对
实施例的描述中将变得明显和容易理解，其中：
[0048] 图1是本发明系统结构图；
具体实施方式
[0050] 下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
[0051] 在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“
水平”、“顶”、“底”“内”、“外”等指示的方位或
位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
[0052] 在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
[0053] 本发明提出了一种由PCA演进的大规模MIMO-OFDM系统信道状态信息压缩反馈方法，根据考虑场景建立了相关的信道模型，将连续的信道状态信息进行分簇，提高了信道状态信息的相关性以及基于PCA的压缩性能，通过信道统计特性以及在收发两端预定义的二进制选择矩阵码本，减少了大量的反馈量和降低了复杂度。
[0054] 结合附图1和附图2对本发明进行详细说明，主要包括以下步骤：
[0055] 步骤1：开始。
[0056] 步骤2：获取下行信道相关特性参数，建立相关信道模型。
[0057] 本发明考虑大规模MIMO-OFDM系统，发射端为配置Nt根均匀线性阵列的天线，接收端为单天线用户，有Nc个子载波。在大规模MIMO系统中，由于发射端配置天线数较多时，天线间距较小，天线间呈现较强的空间相关性，因此，发射端与第k单天线用户间的信道可建模为：
[0058]
[0059] 其中， 表示发送相关矩阵，n表示第n个子载波，n∈{1,…,Nc}；表示均值为0，方差为1的独立分布复高斯随机矢量，Hiid,k(n)的第(u,v)个
元素代表发射天线u和接收天线v之间由路径损耗和小尺度衰落组成的信道增益。
[0060] 由于发射端配置Nt根均匀线性阵列天线，空间相关矩阵可以通过Jakes模型得到，第p根天线与第q根天线之间的相关系数表示为：
[0061]
[0062] 其中J0(·)表示第一类零阶贝塞尔函数，dpq表示第p根天线与第q根天线之间的距离，λ表示载波波长。
[0063] 大规模MIMO-OFDM系统中，发射端与第k个用户间在空域和频域的信道向量hk表示为：
[0064]
[0065] 其中N＝Nt×Nc，vec(A)表示将一个a×b的矩阵A矢量化，变成ab×1的列向量，hk(n)为第k个用户第n个子载波的空域信道向量，建立模型为：
[0066]
[0067] 步骤3：对信道向量hk进行分簇。
[0068] 接收端将连续的N个信道状态信息分为G个簇，即 其中G≤N；hk,g表示第g个簇的信道状态信息，用G个簇来代替N个信道状态信息元素，则每一簇包含 个信道状态信息元素。
[0069] 步骤4：根据特性参数和步骤3，首先计算出特性参数信道每一簇信道状态信息的协方差矩阵 然后将每一簇信道状态信息的协方差矩阵进行特征值分解，即Ck,g＝Ψk,gΛΨk,gH，获得稀疏矩阵Ψk,g，并将Ψk,g储存在发射端和接收端。具体实施过程如下：
[0070] (1)根据特性参数计算信道每一簇信道状态信息的协方差矩阵，第k个用户第g簇信道状态信息的协方差矩阵表示为：
[0071]
[0072] (2)对第g簇信道状态信息的协方差矩阵进行特征值分解，获得稀疏矩阵：
[0073]
[0074] 其中， 表示对协方差矩阵 进行特征值分解，即Ck,g＝Ψk,gΛΨk,gH。Λ是对角矩阵，其对角元素是协方差矩阵 的特征值，Ψk,g的列向量是协方差矩阵 的
特征向量，通过卡洛南变换，Ψk,g同时表示稀疏矩阵。
[0075] 由于步骤4中的信道向量hk是通过信道特性参数获得的，则只需确定信道状态信息分簇数G，那么收发两端的协方差矩阵就是已知的，进而收发两端的稀疏矩阵Ψk,g也是已知的。将稀疏矩阵Ψk,g储存在发射端和接收端，则接收端不需再向基站反馈Ψk,g，因此减少了反馈量。
[0076] 步骤5：信道估计。
[0077] 首先估计第k个用户所有个子载波的空域相关信道，即获得1×Nt的Hk(n)，n∈{1,…,Nc}；然后根据公式(4)将Hk(n)矢量化转变为Nt×1的空域信道向量hk(n)；再根据公式(3)构造信道向量hk。本发明假定采用理想的信道估计，不考虑信道估计误差。
[0078] 步骤6：信道的稀疏。
[0079] 首先根据步骤3将信道估计构造的信道向量hk进行分簇；然后接收端第g簇信道状态信息hk,g通过稀疏矩阵Ψk,g得到一个稀疏向量sk,g，sk,g通过提取hk,g的主要成分来代表hk,g，表示为
[0080]
[0081] 步骤7：预定义选择矩阵码本和选取最优选择矩阵。
[0082] 首先预定义一个收发两端都已知的二进制选择矩阵码本CP，然后通过选定的稀疏信道向量s′k,g内各元素之和最大准则选取最优选择矩阵Pi。具体实施过程如下：
[0083] (1)预定义一个收发两端都已知的二进制选择矩阵码本CP：
[0084] CP＝{P1…PL} (8)
[0085] 其中Pi表示选择矩阵码本的第i个码字，即选择矩阵；L表示码本CP包含L个选择矩阵。
[0086] (2)通过选择矩阵Pi从稀疏向量sk,g中选择M个最重要的反馈信息，其中
[0087] s′k,g＝Pisk,g (9)
[0088] 全局搜索选择矩阵码本，当选定的稀疏信道向量s′k,g内各元素之和最大时，接收端确定使用的选择矩阵Pi。
[0089] i＝argimax|sk,gPi| (10)
[0090] 其中Pi表示选择矩阵码本的第i个码字。传统主成分分析方法需要动态更新每簇和不同时刻的选择矩阵Pi并将其反馈到发射端，增加了计算复杂度。本发明预先在收发两端设定一个二进制选择矩阵码本CP，只需反馈选择矩阵的索引到基站，减少了计算复杂度和反馈开销。
[0091] 步骤8：确定选择矩阵Pi后，通过公式(9)从稀疏向量sk,g中提取M个主要成分，即主要反馈信息。步骤8将高维的信道状态信息变成了低维的信道状态信息。
[0092] 步骤9：采用随机矢量量化码本，量化低维的选定的稀疏信道向量s′k,g，然后将s′k,g码字索引以及步骤7选择矩阵Pi的码字索引一起反馈到发射端。s′k,g码字索引通过以下公式获得：
[0093] j＝argjmax|s′k,gwj| (11)
[0094] 其中wj是随机矢量量化码本的第j个码字。与直接量化高维的信道状态信息相比，随机矢量量化码本能够设计得小一些，减少了反馈开销。
[0095] 步骤10：发射端重构信道状态信息。
[0096] 发射端通过反馈获得的s′k,g码字索引j以及选择矩阵Pi码字索引i进行信道状态信息重构，具体实施过程如下：
[0097] (1)发射端从随机矢量量化码本中找到对应第g簇s′k,g的码字，从选择矩阵码本CP＝{P1…PL}中找到对应第g簇选择矩阵Pi的码字，根据以下公式：
[0098]
[0099] 恢复出第g簇的高维信道状态信息向量
[0100] (2)发射端通过选择矩阵Pi和稀疏矩阵Ψk,g作逆变换，根据以下公式：
[0101]
[0102] 恢复第g簇的信道状态信息
[0103] (3)发射端根据分簇序号g依次对 进行重新排列，重构出下行链路信道状态信息hk。
[0104] 步骤11：结束。
[0105] 本发明与传统PCA用于信道信息压缩反馈方法相比，其创新之处在于建立信道模型时考虑了空域和频域两个域的信道特性，并将得到的信道向量进行分簇，再对每一簇信道状态信息基于卡洛南变换得到稀疏向量，最后通过选择矩阵提取稀疏向量的主要成分并反馈，本发明中预定义了静态的选择矩阵并以码本形式反馈，减小反馈量。
[0106] 在本
说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0107] 尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、
修改、替换和变型，本发明的范围由
权利要求及其等同物限定。
