基于MIMO系统加权能效最大化的分布式波束赋形优化方法技术方案

本发明专利技术涉及一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法，包括：步骤S1：所有基站初始化发射波束赋形向量、接收波束赋形向量；步骤S2：每个基站向其各自所服务的用户发送该用户的当前发射波束赋形向量；步骤S3：各用户根据自身的发射波束赋形向量更新自身的接收波束赋形向量，并将其反馈给所有基站；步骤S4：各基站接收到的所有用户的接收波束赋形向量，计算自己到所有用户的有效信道信息；步骤S5：根据上述计算结果通过ADMM过程、SCA过程以及MMSE过程三层迭代得到最优发射波束赋形向量和接收波束赋形向量。与现有技术相比，本发明专利技术具有提升通信质量等优点。

Distributed Beam Shaping Optimization Method Based on Maximizing Weighted Energy Efficiency of MIMO System

【技术实现步骤摘要】
基于MIMO系统加权能效最大化的分布式波束赋形优化方法
本专利技术涉及一种通信方法，尤其是涉及一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法。
技术介绍
移动数据通信的持续增长对无线网络的容量提出进一步要求，移动网络的性能主要受到来自相邻小区同频率无线电波的干扰，因此设计先进的干扰协调技术对提升无线网络性能至关重要；其中适用于多小区、多用户MIMO通信的COMP(协作多点传输)技术通过多个基站波束赋形的协同设计，可有效实现各小区间的干扰协调。在优化各基站波束赋形向量时，通常以最小化基站功耗或最大化用户吞吐量作为目标函数，而很少涉及最大化能效(功耗/吞吐量)的问题；此外使用加权能效和作为目标进行优化可以更好的满足异构网络不同小区对能效的不同需求。上述提到的协作波束赋形优化问题通常需要存在一个中心节点(基站)汇总所有基站和用户的信息，或基站之间通过回程链路彼此交换信道状态矩阵，从而集中式的计算最优波束赋形向量；然而实际系统回程链路容量有限，且当基站、用户以及收发双方天线数目非常大时，现有硬件系统难以满足集中式运算的高复杂度要求；此外大多数现有文献只考虑了MISO系统的发射波束赋形优化问题，且在对基站功耗建模时常将电路功耗视为常数，而没有考虑到实际模型功耗也取决于系统速率。
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法。本专利技术的目的可以通过以下技术方案来实现：一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法，包括：步骤S1：所有基站初始化发射波束赋形向量、接收波束赋形向量；步骤S2：每个基站向其各自所服务的用户发送该用户的当前发射波束赋形向量；步骤S3：各用户根据自身的发射波束赋形向量更新自身的接收波束赋形向量，并将其反馈给所有基站；步骤S4：各基站接收到的所有用户的接收波束赋形向量，计算自己到所有用户的有效信道信息；步骤S5：根据上述计算结果通过ADMM过程、SCA过程以及MMSE过程三层迭代得到最优发射波束赋形向量和接收波束赋形向量。所述步骤S1中还初始化松弛变量和引入变量ADMM全局变量v(0)，以及拉格朗日乘子；所述步骤S5具体包括：步骤S51：设ADMM过程、SCA过程以及MMSE过程的迭代指示初始值分别为l＝0，n＝0，m＝0；步骤S52：每个基站执行最里层ADMM迭代过程，更新局部变量为全局变量为v(l+1)，以及拉格朗日乘子为直至更新后的变量与更新前的变量残差小于第一设定阈值，得到ADMM迭代优化后的最优解，并进入中间层SCA迭代过程；步骤S53：各基站根据ADMM迭代的最优解，进行SCA迭代更新，直至SCA迭代更新后的变量与更新前之差小于第二设定阈值，并得到SCA迭代优化后的最优解，其中令：{v(0)，ξ(0)}＝{v(l)，ξ(l)}其中：为用户bk的经过SCA迭代更新后的发射波束赋形向量，为用户bk的经过SCA迭代更新后的引入变量，为经过SCA迭代更新后的松弛变量，v(0)和ξ(0)分别为经SCA迭代更新的ADMM迭代过程中全局变量和拉格朗日乘子的初始值，v(l)，ξ(l)和为ADMM迭代优化的最优解；步骤S54：各基站根据SCA迭代得到的最优解，更新接收波束赋形向量，并判断更新后的接收波束赋形向量与更新前接收波束赋形向量之差是否小于第三设定阈值，若为是，则优化结束，若为否，则返回步骤S2所述步骤S52中，局部变量、全局变量和拉格朗日乘子中任一者更新后与自身当前值之间的残差小于第一设定阈值，即进入中间层SCA迭代过程。所述步骤S3中，用户更新的接收波束赋形向量为：其中：为用户bk更新后的接收波束赋形向量，为基站i到用户bk的信道矩阵，为基站i所服务的用户ij的更新前发射波束赋形向量，(·)H为矩阵的共轭转置，N0为噪声功率谱密度，I为单位矩阵，为基站b到用户bk的信道矩阵，为用户bk的当前的发射波束赋形向量，B为基站的集合，Ki为基站i所服务的用户的集合。所述步骤S4中，基站b到用户ij的有效信道信息具体为：其中：为用户ij更新后的接收波束赋形向量，为基站b到用户ij的信道矩阵。所述步骤S52具体包括：步骤S521：每个基站计算局部变量其中：为更新后的基站b的局部变量，为局部变量sb满足的约束集合，ωb为基站b的能效权重，ηb为代表每个基站能效的松弛变量，为当前的拉格朗日乘子，θb为ADMM优化算法引入的一致性变量，(·)T表示转制，ρ为罚系数，为当前的全局变量，表示欧几里得范数的平方。步骤S522：各基站更新全局变量其中：v(l+1)为更新后的全局变量，步骤S523：更新拉格朗日乘子：其中：为更新基站b的拉格朗日乘子，为更新前基站b的拉格朗日乘子；步骤S524：判断更新后的任何优化变量与更新前优化变量之差是否小于第一设定阈值，若为是，则执行步骤S53，若为否，则l＝l+1，返回步骤S521。所述第一设定阈值为10-5。所述步骤S53具体包括：S531：更新当前SCA优化系数和ADMM参数初始值；S532：判断更新后的发射波束赋形向量与更新前发射波束赋形向量之差是否小于第二设定阈值，若为是，则执行步骤S54，若为否，则n＝n+1，返回步骤S531。所述步骤S54具体包括：步骤S541：更新接收波束赋形向量；步骤S542：判断更新后的接收波束赋形向量与更新前的接收波束赋形向量之差是否小于第三设定阈值，若为是，则优化结束，若为否，则m＝m+1，返回步骤S4。所述局部变量sb满足的约束集合具体为：其中：为经SCA第n次迭代更新的的一阶泰勒近似，为经SCA第n次迭代更新的的一阶泰勒近似，a为给定常数，ε为功率放大效率，P0为基站的电路功耗，Pb为基站b额定发射功率，Γ0为接收用户信噪比阈值，为基站b对用户ij(i≠b)的干扰ICI变量，为用户bk受到的来自其他基站的干扰ICI变量。与现有技术相比，本专利技术具有以下有益效果：基于实际MIMO系统中各小区加权能效最大化的分布式波束赋形算法，结合近似凸优化算法SCA、机器学习中的ADMM算法、以及接收端的最小平方差MSE准则，在MIMO系统的收发双方交替优化发射波束赋形和接收波束赋形向量，并在基站优化发射波束赋形向量时，通过分离局部变量使基站能够独立计算而不需要彼此交换信道状态信息。。附图说明图1为本专利技术方法的主要步骤流程示意图。具体实施方式下面结合附图和具体实施例对本专利技术进行详细说明。本实施例以本专利技术技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于下述的实施例。本专利技术考虑包含B小区的多用户MIMO下行链路系统：其中每个基站b∈B＝{1，...，B}配置NT发射天线、服务Kb用户，共K用户每个用户配置NR接收天线，且每个用户bk仅由单个基站b服务，即设基站b到用户k的信道矩阵为满足平坦衰落信道模型；基站采用线性编码方式，是用户bk的发射数据符号，用户bk的接收信号为表示加性白高斯噪声；设是用户bk的组合编码(接收波束赋形)向量，为用户bk的预编码(发射波束赋形)向量，则用户bk的接收MSE(最小平方误差)为：基于上述定义，基站b的能效表达式EE可由下式给出：其中表示用户bk的接收SINR，ε～(0，1)表示功本文档来自技高网...

【技术保护点】
1.一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法，其特征在于，包括：步骤S1：所有基站初始化发射波束赋形向量、接收波束赋形向量；步骤S2：每个基站向其各自所服务的用户发送该用户的当前发射波束赋形向量；步骤S3：各用户根据自身的发射波束赋形向量更新自身的接收波束赋形向量，并将其反馈给所有基站；步骤S4：各基站接收到的所有用户的接收波束赋形向量，计算自己到所有用户的有效信道信息；步骤S5：根据上述计算结果通过ADMM过程、SCA过程以及MMSE过程三层迭代得到最优发射波束赋形向量和接收波束赋形向量。

【技术特征摘要】
1.一种基于MIMO系统加权能效最大化的分布式波束赋形优化方法，其特征在于，包括：步骤S1：所有基站初始化发射波束赋形向量、接收波束赋形向量；步骤S2：每个基站向其各自所服务的用户发送该用户的当前发射波束赋形向量；步骤S3：各用户根据自身的发射波束赋形向量更新自身的接收波束赋形向量，并将其反馈给所有基站；步骤S4：各基站接收到的所有用户的接收波束赋形向量，计算自己到所有用户的有效信道信息；步骤S5：根据上述计算结果通过ADMM过程、SCA过程以及MMSE过程三层迭代得到最优发射波束赋形向量和接收波束赋形向量。2.根据权利要求1所述的一种基于MIMO系统加权能效最大化的分布式波束赋形方法，其特征在于，所述步骤S1中还初始化松弛变量和引入变量ADMM全局变量v(0)，以及拉格朗日乘子；所述步骤S5具体包括：步骤S51：设ADMM过程、SCA过程以及MMSE过程的迭代指示初始值分别为l＝0,n＝0,m＝0；步骤S52：每个基站执行最里层ADMM迭代过程，更新局部变量为全局变量为v(l+1),以及拉格朗日乘子为直至更新后的变量与更新前的变量残差小于第一设定阈值，得到ADMM迭代优化后的最优解，并进入中间层SCA迭代过程；步骤S53：各基站根据ADMM迭代的最优解，进行SCA迭代更新，直至SCA迭代更新后的变量与更新前之差小于第二设定阈值，并得到SCA迭代优化后的最优解，其中令：{v(0),ξ(0)}＝{v(l),ξ(l)}其中：为用户bk的经过SCA迭代更新后的发射波束赋形向量，为用户bk的经过SCA迭代更新后的引入变量，为经过SCA迭代更新后的松弛变量，v(0)和ξ(0)分别为经SCA迭代更新的ADMM迭代过程中全局变量和拉格朗日乘子的初始值，v(l),ξ(l)和为ADMM迭代优化的最优解；步骤S54：各基站根据SCA迭代得到的最优解，更新接收波束赋形向量，并判断更新后的接收波束赋形向量与更新前接收波束赋形向量之差是否小于第三设定阈值，若为是，则优化结束，若为否，则返回步骤S2。3.根据权利要求2所述的一种基于MIMO系统加权能效最大化的分布式波束赋形方法，其特征在于，所述步骤S52中，局部变量、全局变量和拉格朗日乘子中任一者更新后与自身当前值之间的残差小于第一设定阈值，即进入中间层SCA迭代过程。4.根据权利要求2所述的一种基于MIMO系统加权能效最大化的分布式波束赋形方法，其特征在于，所述步骤S3中，用户更新的接收波束赋形向量为：其中：为用户bk更新后的接收波束赋形向量，为基站i到用户bk的信道矩阵，为基站i所服务的用户ij的更新前发射波束赋形向量，(·)H为矩阵的共轭转置，N0为噪声功率谱密度，I...

【专利技术属性】
技术研发人员：赵生捷，韩丰夏，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海,31