基于天线选择和干扰对齐的多小区大规模mimo干扰消除方法

基于天线选择和干扰对齐的多小区大规模mimo干扰消除方法
【专利摘要】本发明公开了一种基于天线选择和干扰对齐的多小区大规模MIMO干扰消除方法，属于移动无线通信技术领域。本发明包括：基站端获取信道矩阵并分别计算各发射天线所对应的列向量的范数，选择范数最大的前Ns根天线作为活跃天线，并关闭其它天线；基于活跃天线所对应的列向量作为新的信道矩阵并计算每个基站的干扰对齐预编码向量，根据干扰对齐预编码向量和新的信道矩阵计算接收端的接收迫零向量并据此分配各数据流的发送功率，然后再基于发送功率对接收迫零向量进行更新，最后，当前基站基于其干扰对齐预编码向量、当前接收迫零向量、各数据流的发送功率生成接收端对应的发送数据并发送。本发明的实施，能降低计算复杂度，提高传输速率及降低误比特率。
【专利说明】
基于天线选择和干扰对齐的多小区大规模MIMO干扰消除方法
技术领域
[0001] 本发明属于移动无线通信技术领域，特别设及一种基于天线选择和干扰对齐的多 小区大规模MIMO干扰消除方法
【背景技术】
[0002] 大规模MIMO技术，通过在基站端配置大量的天线可W显著提高系统的频谱效率和 能效效率，在基站天线数目趋于无穷时不同用户之间的信道矩阵趋于正交，从而可W消除 用户间干扰和小区间干扰，由此受到了广泛的研究和关注。然而，实际情况中的天线数目不 可能无穷多，不同用户之间的信道矩阵并非是理想正交的，即用户间干扰依旧存在，尤其是 小区边缘用户。
[0003] 干扰对齐技术，是在发送端利用已经获知的信道状态信息设计预编码矩阵，使得 干扰信号在接收端对齐，从而得到了无干扰的期望信号，同时压缩了干扰信号所占的子空 间，增大了期望信号的维度，提高了系统的传输速率。目前干扰对齐算法一般分为线性干扰 对齐和迭代干扰对齐，而前者由于计算复杂度较低在实际应用具有一定的优势。
[0004] 天线选择技术是按照一定的准则，从所有的服务天线中选择其中一部分天线进行 信号发射。天线选择准则可W是最大化信道容量准则，最大信噪比准则和最小误码率则准 等等。在实际系统中，考虑到实际物理限制，实现复杂度和成本等因素，数量庞大的天线可 能会成为工程实现的瓶颈。通过天线选择技术可W有效降低系统信号处理的复杂度和成 本，并且在一定程度上保留了大规模天线阵列带来的绝大部分优良性能，可W实现在大规 模MIMO系统性能和复杂度的有效折衷。
[0005] 功率分配是资源分配的内容之一，根据实际信道情况进行功率分配可W提高通信 系统的性能，比如注水算法可W使得MIMO系统的信道容量最大化。如果不考虑功率分配，即 等功率分配，将会使得所传输的每个数据流的有效信道增益严重失衡，很难使系统的性能 达到最优。文南犬('Interference Alignment Based On Antenna Selection For Massive MIMO System"(Zhiyuan Shi,Xiaopeng Zhu,Yifeng Zhao,Lianfen Huang,Computer SciencefeEducation(ICCSE),2015 10th International Conference on,22-24July 2015，pp.606-610)在大规模MIMO系统中联合采用了天线选择和干扰对齐技术来消除干扰， 但是并没有考虑功率分配的问题，而且天线选择算法的复杂度较高，系统性能仍有较大的 提升空间。

【发明内容】

[0006] 本发明的发明目的在于:针对上述存在的问题，提出了一种基于天线选择，干扰对 齐和动态功率分配的大规模MIMO干扰消除方法，其利用干扰对齐技术大幅度提高了用户的 传输速率和误比特率性能，同时通过天线选择技术降低了信道矩阵维度，从而降低了运算 复杂度，有利于工程实现。
[0007] 本发明的基于天线选择和干扰对齐的多小区大规模MIMO干扰消除方法，包括下列 步骤：
[0008] 步骤I:基站端通过上行导频获取所有用户的信道状态信息CSI，即获取信道矩阵 出jEcWxm(cWxm表示NXM矩阵），其中出徒示用户j至蜡站i的信道矩阵，M表示基站端的天线 数目，N为用户端的天线数目。
[0009] 步骤2:在每个基站，根据信道容量最大化准则和基站各自所服务的用户的信道矩 阵进行天线选择。所述天线选择依据的准则为：
I其中4表示基 站端所选择的天线子集，hi表示发射天线i所对应的列向量。即分别计算信道矩阵出1的所有 列向量的范数，选择范数最大的前化根天线作为活跃天线，Ns为所需选择天线数目，并关闭 剩余M-化根天线。
[0010] 步骤3:基站端根据步骤2天线选择的结果重新确定实际采用的信道矩阵，对于每 一个信道矩阵出j，选择其中的化列列向量作为新的信道矩阵山曰3山即山曰3^=比^山刮， 其中，比1北]表示信道矩阵出J的第k列。
[0011] 步骤4:基站端根据天线选择后的信道矩阵H_asu计算干扰对齐预编码向量。
[0012] 本发明所采用的干扰对齐方法为步骤401-403所述：
[OOU] 401:对联合矩阵进行特征值分解，选取特征值最大 的D(D为发送的数据流数目）列特征向量作为基站1的干扰对齐预编码向量Vi并进行向量归 一化处理；即 vi(:，;r) = vi(: ,rVl |vi(:，!) I I，r = l,. . .，D，符号（：，!)表示取第 r 列运算
[0014] 402:根据公式V';=八;.;| * V;计算基站2的干扰对齐预编码向量V2并进行向量归 一化处理；
[0015] 40 3:根据公式1? =总*巧1。*巧i十算基站3的干扰对齐预编码向量V. 3并进行向量 归一化处理。
[0016] 步骤5:计算接收端的接收迫零向量。
[0017] 所述迫零处理采用的迫零矩阵根据下列公式得到：
[001 引
[0019]
[0020]
[0021] 其中，Wj为DXN的用户j的接收迫零矩阵(D为用户j发送的数据流数目，N为用户j 的天线数目），即Wj由其中对应的d行行向量组成，.曰为等效联合发送信道矩阵，U_&为等效 联合矩阵的伪逆矩阵，（?）H表示矩阵的共辆转置，下同。
[0022] 步骤6:分别计算每个用户的接收迫零矩阵的列向量的模值，并根据模值的大小进 行动态功率分配：
[0023] 先计算用户j的迫零向量的模值向量
'其中符号diag(A) 表示取矩阵A的对角元素。
[0024] 再根据模值进行数据流的发送功率分配
I其中，P表示所有 基站的发射总功率，Pk表示当前基站向用户j发射第r流数据的功率，山r表示用户j的第r个 迫零向量的模值，即矩阵dj的第j列所对应的模值 J
[00巧]步骤7:根据当前功率分配情况更新迫零矩阵
，其中pj = diag (Pj1，PW，...，PW)。步骤8:结合上述天线选择、干扰对齐W及功率分配的计算结果生成发送 数据:Xj = Vi*pj*Sj，其中Xj表示当前基站对用户j的N X 1发送数据向量，Pj为当前基站关于 用户j的功率分配对角矩阵，对角元素为每个数据流分配的功率Wr, Sj为发送给用户j的原 始发送数据。
[0026] 综上所述，由于采用了上述技术方案，本发明的有益效果是：与现有技术相比，本 发明通过天线选择方法降低了信道矩阵维度，从而使得可用相对较低的复杂度实现了大规 模MIMO系统干扰消除，提高了用户的传输速率，同时通过动态功率分配降低了系统的误比 特率性能，提高了系统传输的可靠性。
【附图说明】
[0027] 图1为本发明的实施例场景示意图。
[0028] 图2为本发明实施例的实施结构框图。
[0029] 图3为本发明所采取天线选择算法与性能最佳的渐消算法的可达率对比曲线。
[0030] 图4为本发明实施案例4天线情况下的误比特率性能对比曲线。
[0031] 图5为本发明实施案例8天线情况下的误比特率性能对比曲线。
[0032] 图6为本发明实施案例4天线情况下的可达率对比曲线。
[0033] 图7为本发明实施案例8天线情况下的可达率对比曲线。
【具体实施方式】
[0034] 为使本发明的目的、技术方案和优点更加清楚，下面结合实施方式和附图，对本发 明作进一步地详细描述。
[00对实施例
[0036] 参见图1，本实例的系统场景为基于C-RAN架构的大规模MIMO多小区系统，系统采 用TDD制式。在C-RAN架构下可W轻易的无延时的共享所有小区的信道信息，大大提高了联 合处理的可实施性。在本案例中考虑3个小区，每个小区一个基站，每个基站服务一个小区 边缘用户，基站端的天线数目为128,用户端的天线数目为N，本实例中对N取4和8进行仿真 分析。本实施例中的信道状态信息为零均值，方差为1的复高斯变量。所有基站的发射总功 率均为P。
[0037] 针对上述实施例，本发明提供的基于天线选择和干扰对齐W及功率分配的大规模 MIMO干扰消除方法的具体实现步骤为：
[0038] 步骤SI:基站端通过上行导频获取所有用户的信道状态信息CSI，即获取信道矩阵 出j G CWXi28，i，j = 1，2，3，其中出j表示用户j到基站i的信道矩阵。
[0039] 步骤S2:在每个基站，根据信噪比最大化准则和基站各自所服务的用户的信道矩 阵出i，i = 1，2，3进行天线选择。
[0040] 即根据公J
分别计算信道矩阵出1，此沸曲3 的所有列向量的模值，选择其中模值最大的咐良天线作为活跃天线，关闭剩余128-咐良天线。
[0041] 步骤S3:基站端根据步骤S2天线选择的结果重新确定实际采用的信道矩阵，对于 每一个信道矩阵，选择其中的N列列向量作为新的信道矩阵H_asij。
[0042] 山曰34 =比1化]1^￡*其中，比1化]表示信道矩阵出撕第巧!]。
[0043] 步骤S4:基站端根据天线选择后的信道矩阵H_asu计算干扰对齐预编码向量VI。首 先对联合矩阵W,;| *//,, * W,,1 * W，, */"/,:1 * W。进行特征值分解，选取特征值最大的N/2(本实 施例中，发送的数据流数目D为N/2)列特征向量作为基站1的干扰对齐预编码向量VI，并根 据公式vi(:，r) = vi(:，r)/ M vi(:，r) I I，r= 1，. . .，N/2进行向量归一化处理;然后根据公式 V3==パ:>パl.;*，，l计算基站2的干扰对齐预编码向量V2并进行向量归一化处理;最后根据公 式V'; = /7。| *//|: *巧计算基站3的干扰对齐预编码向量V3并进行向量归一化处理。
[0044] 步骤S5:计算接收端的接收迫零向量。
[0045] 首先根据
计算U_Hj，再根据Wj = U_Hj((j-1 )蝴+1: j*D，：），j = 1，2，3得到用户j的接收迫零矩阵Wj，即Wj 由U_Hj中对应的N/2行行向量组成。
[0046] 步骤S6:分别计算每个用户的接收迫零矩阵的列向量的模值，并根据模值进行数 据流的发送功率分配，目[
庚中Cl^为用户j的迫零向量的模值向量dj
的第r列的模值，T为所有用户迫零向量的模值之和。
[0047] 步骤S7:根据功率分配情况更新迫零矩阵。目H ，其中
[0048] 步骤S8:当前基站端结合上述天线选择、干扰对齐W及功率分配的计算结果生成 关于原始发送数据Sj的发送数据Xj : Xj = Vi*p片Sj。
[0049] 步骤S9:接收端对接收信号yj进行迫零处理即可获取到对应的发送数据，即yj = Wi*：Tj，其中rj为NX 1的接收向量。
[0050] 为了评估本发明的性能，采用蒙特卡洛仿真方法对信道容量和误比特率进行了仿 真。
[0051] 图3为本发明采用的天线选择算法与性能最优的渐消算法的可达率对比曲线，其 中QPSK表示正交频移键控，AS表示天线选择，IA表示干扰对齐，ZF表示迫零预编码，MF表示 匹配滤波预编码。可W看到，在本实例所假设的信道系数为零均值，方差为1的复高斯变量 场景下，本发明所采取的简单的天线选择算法相对于渐消算法的信道容量，在低信噪比时 性能几乎一致，而高信噪比下在4天线时仅有约1-化ps/Hz的性能损失，在8天线时仅有约2- SbpsAlz的性能损失。而本发明所采取的天线选择算法复杂度则要低的多，有利于工程实 现。本发明所采取的天线选择算法仅需计算每根发射天线所对应的传输向量的模值并取其 中模值最大的Ns根，而信道容量最大化渐消算法需要进行N-Ns次迭代，每次迭代需要计算 去除每根天线所产生的信道容量性能损失，复杂度明显要高于本发明所采取的天线选择算 法。
[0052] 图4和图5为不同算法的误比特率性能对比，可W看到本发明发明提出的基于干扰 对齐和动态功率分配的干扰消除方法，在误比特率性能方面几乎接近于理想环境下的QPSK 调制系统，多点协作迫零预编码和未进行功率分配(等功率分配）的干扰对齐性能接近，在 高信噪比时本发明所提出的方法要明显好于多点协作迫零预编码，在相同的误比特率下有 3-4地的性能提升。
[0053] 图6和图7为不同算法的可达率对比曲线，比较对象为本发明所提出的基于天线选 择，干扰对齐和动态功率分配的干扰消除方法，和基于天线选择，干扰对齐的干扰消除方 法，基于天线选择的多点协作迫零预编码算法，W及不进行天线选择的非协作迫零预编码 算法。仿真结果显示，不进行天线选择的非协作迫零预编码系统是干扰受限的，基于天线选 择的协作迫零预编码是最优的，基于等功率分配的干扰对齐方法与之接近。本发明的方法 相对非协作迫零预编码方法有非常明显的优势，而在加入动态分配算法后在可达率性能上 有一定的损失，但损失较小，相对于多点协作迫零预编码方法在4天线时性能损失约为 3bits/sec/Hz，在8天线时性能损失约为加 its/sec/化，运仍属于可接受范围。因为本发明 所提方法的复杂度相对多点协作预编码方法要低得多，本发明所提方法中最为复杂的部分 为矩阵规模为NX N的特征值分解和矩阵求逆运算，复杂度均为0(妒），而多点协作迫零预编 码需要进行全局信道矩阵的求逆运算，计算复杂度为〇(化N)3)，其中K为总用户数。
[0054] W上所述，仅为本发明的【具体实施方式】，本说明书中所公开的任一特征，除非特别 叙述，均可被其他等效或具有类似目的的替代特征加 W替换;所公开的所有特征、或所有方 法或过程中的步骤，除了互相排斥的特征和/或步骤W外，均可W任何方式组合。
【主权项】
1.基于天线选择和干扰对齐的多小区大规模ΜΙΜΟ干扰消除方法，其特征在于，包括下 列步骤： 步骤1:基站端获取信道矩阵出j，其中用户标识符j = 1，2，3，基站标识符i = 1，2，3，且出j 为NXM矩阵，其中N为用户端的天线数目，Μ为基站端的天线数目； 步骤2:基于信道矩阵出J，基站端分别计算各发射天线所对应的列向量的范数，选择范 数最大的前化根天线作为活跃天线，并关闭剩余Μ-化根天线； 步骤3:基站端根据从信道矩阵出J中选择化根活跃天线所对应的列向量作为新的信道 矩阵 H_asij; 步骤4:基站端根据信道矩阵H_asu计算每个基站的干扰对齐预编码向量vi: 401:对联合矩阵i * //,, a //2,1 * //;| * //。| * //。进行特征值分解，选取特征值最大的D列 特征向量作为基站1的干扰对齐预编码向量VI并进行向量归一化处理，其中D为发送的数据 流数目；即 V1(:，;r) = vi(:，;r)/| |vi(:，!) I I，r = l,. . .，D，符号（：，!)表示取第 r 列运算； 402:根据公式V。= *//|; *崎:计算基站2的干扰对齐预编码向量V2并进行向量归一化 处理； 403:根据公式Vg = *//|, *，'i计算基站3的干扰对齐预编码向量V.3并进行向量归一化 处理； 步骤5:基站端根据公式W戸11_&(〇-1州+1。'*0，：）计算用户前接收迫零向量巧^，其 中。_//, =(/7,"/7,) '/7," '且馬二巧 _化、*v'2 巧-化ν;,.*ν'3)，符号（·）H表 示矩阵的共辆转置； 步骤6:基站端基于接收迫零向量W汾配各数据流的发送功率化，即P,,. = 户，Γ = 1，...，D， 其4>表示所有基站的发射总功率，山r表示矩阵山的第r 列的非零元素的值，符号diag( ·)表示求取对角线； 步骤7 :根据步骤6分配的发送功率wr更新接收迫零向量，即令= 1.2,3，其 中pj = diag(pji，pi2,. . .，p扣）； 步骤8:当前基站端i生成发送给用户j发送数据Xj并发送，其中Xj = Vi*Pj*Sj，所述S徒 示发送给用户j的原始发送数据。
【文档编号】H04B7/06GK105978615SQ201610325499
【公开日】2016年9月28日
【申请日】2016年5月17日
【发明人】杨全校, 杨海芬, 林水生, 周亮, 袁强, 尹腾飞
【申请人】电子科技大学