CN102572864B - 最大化吞吐量的多小区联合波束成形设计方法 - Google Patents
最大化吞吐量的多小区联合波束成形设计方法 Download PDFInfo
- Publication number
- CN102572864B CN102572864B CN201110383525.4A CN201110383525A CN102572864B CN 102572864 B CN102572864 B CN 102572864B CN 201110383525 A CN201110383525 A CN 201110383525A CN 102572864 B CN102572864 B CN 102572864B
- Authority
- CN
- China
- Prior art keywords
- overbar
- alpha
- beta
- sigma
- base station
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
一种最大化吞吐量的多小区联合波束成形设计方法,属于无线通信技术领域,基于拉格朗日对偶理论得到对偶上行问题,采用广义特征值分解和内点法迭代求解该对偶上行问题,确定下行单位波束成形向量和最优的上行信干噪比,再由上下行信干噪比对偶相等得到下行功率分配。该方法能够获得较高的系统吞吐量,并且在用户数较多时能够节省计算时间。
Description
技术领域
本发明涉及的是一种无线通信技术领域的方法,具体是一种最大化吞吐量的多小区联合波束成形设计方法。
背景技术
近年来,随着无线通信技术的飞速发展,无线通信大业务量、高速率和高频谱效率的要求日益迫切。在下一代无线通信系统中,频率复用因子为1,存在严重的同信道干扰(Co-Channel Interference,CCI)。作为抑制CCI的关键技术,波束成形技术得到了高速发展和广泛应用。但是,传统的单小区处理将小区间干扰当作噪声,制约了系统性能的大幅提升。多小区联合处理利用基站间信息交互,协作设计传输方案,很大程度上抑制CCI,提升系统性能,因而成为了下一代无线通信系统的关键技术之一。
目前,多小区联合波束成形设计主要有三种准则:1)保证用户QoS的条件下最小化总发射功率;2)在功率约束下最大化系统吞吐量;3)保证用户间公平性,满足功率约束下优化最差用户性能。
现有技术中公开了H Dahrouj and W.Yu的文献“Coordinated beamforming for themulticell multi-antenna wireless system(多小区多天线无线系统的协作波束成形)”,IEEETransactions on Wireless Communications,vol.9,no.5,pp.1748-1759,May.2010,它针对第一种设计准则利用拉格朗日上下行对偶提出一种通过固定点迭代求解上行参数从而确定下行波束成形向量的设计方法。D.Nguyen and T.Le-Ngoc的文献“Multiuserdownlink beamforming in multicell wireless systems:A game theoretical approach(多小区无线系统中多用户下行波束成形设计的博弈论方法)”,IEEE Transactions on SignalProcessing,vol.55,no.7,pp.3326-3338,Jul.2011,对于同样的问题提出了基于博弈论的能够分布式处理的波束成形设计方法。
D.W.H.Cai,T.Q.S.Quek and C.W.Tan的文献“Coordinated max-min SIRoptimization in multicell downlink-duality and algorithm(多小区下行协作最大最小信干比优化的对偶问题和算法)”,in Proceedings of IEEE International Conference onCommunications,Jun.2011,采用拉格朗日对偶和非线性Perron-Frobenius定理给出在多个加权功率约束条件下的最大化最小加权信干比的波束成形设计方法。
目前比较多的研究是针对第一种和第三种准则的,优化系统吞吐量的问题由于很难转化为凸优化形式,因而针对此问题的算法比较少。
L Venturino,N.Prasad and X.Wang的文献“Coordinated linear beamforming indownlink multi-cell wireless networks(多小区无线网络的下行协作线性波束成形)”in42 nd Asilomar Conference on Signals,Systems and Computers,pp.356-360,Oct.2008,在每个基站功率约束下最大化系统加权和速率,通过对优化问题的Karush-Kuhn-Tucker(KKT)条件的分析,确定最优波束成形向量的结构,然后通过迭代的方法求解KKT条件确定相应系数。该算法称为ICBF(Iterative Coordinated Beam-Forming)算法,其计算时间随用户数近似线性增加,不适用于用户数较多的场景。
J.Yang and D.K.Kim的文献“Multi-cell uplink-downlink beamforming throughputduality based on Lagrangian duality with per-base station power constraints(基站功率约束下的基于拉格朗日对偶的多小区上下行吞吐量对偶波束成形)”IEEE CommunicationsLetters,vol.12,no.4,pp.277-279,Apr.2008,推导了在每个基站功率约束下系统下行吞吐量的拉格朗日对偶问题和波束成形向量结构,但是并没有给出求解算法,并且假设每个小区单用户,不存在小区内干扰。
发明内容
本发明针对现有技术存在的上述问题,提供一种最大化吞吐量的多小区联合波束成形设计方法,本发明基于拉格朗日对偶,能够获得较高的系统吞吐量,并且在用户数较多时能够节省大量计算时间。
本发明是通过以下技术方案实现的:
一种最大化吞吐量的多小区联合波束成形设计方法,基于拉格朗日对偶推导出多小区多用户下行吞吐量最大化问题的上行对偶问题,通过求解该上行对偶问题进行波束成形设计;
所述波束成形设计方法具体包括以下步骤:
1)设置系统参数:小区数M,每个小区用户数K,第n个基站到用户(第m个小区的第k个用户)的信道状态信息其中:m,n=1,...,M,k=1,...,K,基站天线数Nt,第m个基站的功率约束Pm,其中:m=1,..,M,零均值复高斯加性噪声的协方差N0;
2)定义和分别为对偶上行问题的用户的发射功率和第m个基站的接收噪声功率,为对应用户的下行单位波束成形向量。构造上行发射功率向量和噪声功率向量以及矩阵T,其中:
3)初始化向量和取一个很小的正实数ε>0,最大迭代次数imax,迭代次数i=0;
4)固定和通过广义特征值分解计算单位波束成形向量
其中:vmax(·)为归一化最大广义特征向量,
5)固定和T(i),采用内点法求解子优化问题得到新的上行发射功率向量内点法初始值为
所述的子优化问题是:
6)固定和T(i),采用内点法求解子优化问题得到新的噪声功率向量内点法初始值为
所述的子优化问题是:
其中:
7)增加迭代次数i=i+1,检验终止条件是否满足,如果不满足则返回步骤4),如果满足则进行步骤8);
所述的终止条件是:
8)如果i≤imax,最优的上行发射功率向量最优的噪声功率向量由广义特征值分解计算下行最优单位波束成形向量如同步骤4),其中:m=1,...,M,k=1,...,K,然后进行下行功率分配。
定义是对应用户的下行发射功率,是对应用户的最优上行信干噪比,构造下行发射功率向量最优上行信干噪比向量 其中:
构造MK×MK维矩阵F为
其中:子矩阵Fmn为
下行发射功率向量p由矩阵F和γ*确定,p=F-1γ*N0。若向量p满足每个基站的功率约束,则算法结束,输出p和T*,若不满足,则采用现有技术中的ICBF算法重新计算波束成形向量并输出;
如果i>imax,采用现有技术中的ICBF算法重新计算波束成形向量并输出。
本发明所提供的这种最大化吞吐量的多小区联合波束成形设计方法,采用广义特征值分解和内点法迭代求解对偶上行问题,确定下行单位波束成形向量和最优的上行信干噪比,再由上下行信干噪比对偶相等得到下行功率分配。该方法基于拉格朗日对偶,能够获得较高的系统吞吐量,并且在用户数较多时能够节省大量计算时间。
附图说明
图1为实施例的场景示意图;
图2为发送信噪比为10dB时分别采用本实施例方法和现有技术中的ICBF算法的吞吐量性能比较图。
图3为发送信噪比为10dB时分别采用本实施例方法和现有技术中的ICBF算法的计算时间性能比较图。
图4为每个小区用户数K=3时分别采用本实施例方法和现有技术中的ICBF算法的吞吐量性能比较图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明所公开的最大化吞吐量的多小区联合波束成形设计方法,是基于拉格朗日对偶推导出多小区多用户下行吞吐量最大化问题的上行对偶问题,通过求解该上行对偶问题进行波束成形设计的,所述的下行优化问题是:
其中:是用户(第m个小区的第k个用户)的单位波束成形向量,是用户的发射功率,是基站n到用户的信道状态信息,N0是零均值复高斯加性噪声的协方差,Pm是基站m的功率约束,Nt是基站天线数。
所述的上行对偶问题是:
其中:和可以看作对偶上行发射功率和接收端噪声功率。
如图1所示的两小区仿真场景,BS表示小区基站,图中示意了两个BS1、BS2,其之间的距离为1000m。U表示每个小区K个位置随机分布的用户,图中每个基站分别示意了三个用户,分别是BS1基站的用户,BS2基站的用户。图2-4中的每个数据点是对1000种节点分布情况的性能平均值。
本发明设计方法包括以下步骤:
第一步、设置系统参数:小区数M,每个小区用户数K,基站天线数Nt,第m个基站的功率约束Pm,其中:m=1,...,M,零均值复高斯加性噪声的协方差N0,第n个基站到用户的信道状态信息为
其中:表示快衰落,每一项为服从均值为0方差为1的复高斯分布的随机变量,10ξ/10表示对数阴影衰落,ξ为均值为0标准差为σξ的实高斯分布随机变量,表示路径损耗,表示基站n到用户的距离,单位为m。
本实施例中,M=2,Nt=4,σξ=8dB,α=-3,β=1.35×107,N0=1。
本实施例中,各个基站功率约束相等,即Pm=PBS,m=1,...,M,发送信噪比为
第二步、定义和分别为对偶上行问题中用户的发射功率和第m个基站的接收噪声功率,为对应用户的下行单位波束成形向量。构造上行发射功率向量和噪声功率向量以及矩阵T,其中:
第三步、初始化向量和取一个很小的正实数ε>0,最大迭代次数imax,迭代次数i=0。
本实施例中, ε=10-5,imax=100;
第四步、固定和通过广义特征值分解计算单位波束成形向量
其中:vmax(·)为归一化最大广义特征向量,
第五步、固定和T(i),采用内点法求解子优化问题得到内点法的初始值为
所述的子优化问题是:
第六步、固定和T(i),采用内点法求解子优化问题得到内点法的初始值为
所述的子优化问题是:
其中:
第七步、增加迭代次数i=i+1,检验终止条件或者i>imax是否满足,如果不满足则返回第四步,如果满足则进行第八步;
第八步、如果i≤imax,最优的上行发射功率向量最优的噪声功率向量由广义特征值分解计算下行最优单位波束成形向量如同第四步,其中:m=1,...,M,k=1,...,K,然后进行下行功率分配。
定义是对应用户的下行发射功率,是对应用户的最优上行信干噪比,构造下行发射功率向量最优上行信干噪比向量 其中:
构造MK×MK维矩阵F为
其中:子矩阵Fmn为
下行发射功率向量p由矩阵F和γ*确定,p=F-1γ*N0。若向量p满足每个基站的功率约束,则算法结束,输出p和T*,若不满足,则采用现有技术中的ICBF算法重新计算波束成形向量并输出;
如果i>imax,采用现有技术中的ICBF算法重新计算波束成形向量并输出。
图2为发送信噪比为10dB,每个小区用户数K从1到10变化时分别采用本实施例方法和现有技术中的ICBF算法的吞吐量性能比较图。
图3为发送信噪比为10dB,每个小区用户数K从1到10变化时分别采用本实施例方法和现有技术中的ICBF算法的计算时间性能比较图。
图4为每个小区用户数K=3,发送信噪比从0dB到25dB变化时分别采用本实施例方法和现有技术中的ICBF算法的吞吐量性能比较图。
由图2-3可见:在发送信噪比为10dB时采用本实施例的最大化吞吐量的多小区联合波束成形设计方法能够获得较高的系统吞吐量,并且在用户数较多时能够节省大量计算时间。由图4可见:在每个小区用户数K=3时采用本实施例的最大化吞吐量的多小区联合波束成形设计方法能够获得较高的系统吞吐量。
Claims (4)
1.一种最大化吞吐量的多小区联合波束成形设计方法,基于拉格朗日对偶推导出多小区多用户下行吞吐量最大化问题的上行对偶问题,通过求解该上行对偶问题进行波束成形设计,其特征在于,包括以下步骤:
1)设置系统参数:小区数M,每个小区用户数K,第n个基站到第m个小区的第k个用户的信道状态信息其中:m,n=1,...,M,k=1,...,K,基站天线数Nt,第m个基站的功率约束Pm,其中:m=1,...,M,零均值复高斯加性噪声的协方差N0;
2)定义和分别为对偶上行问题的用户的发射功率和第m个基站的接收噪声功率,为对应用户的下行单位波束成形向量,构造上行发射功率向量和噪声功率向量以及矩阵T,其中:
3)初始化向量和取一个很小的正实数ε>0,最大迭代次数imax,迭代次数i=0;
4)固定和通过广义特征值分解计算单位波束成形向量
其中:为归一化最大广义特征向量,
5)固定和采用内点法求解子优化问题得到新的上行发射功率向量内点法初始值为
6)固定和T(i),采用内点法求解子优化问题得到新的噪声功率向量内点法初始值为
7)增加迭代次数i=i+1,检验终止条件是否满足,如果不满足则返回步骤4),如果满足则进行步骤8);
8)如果i≤imax,最优的上行发射功率向量最优的噪声功率向量由广义特征值分解计算单位波束成形向量如同步骤4),其中:m=1,...,M,k=1,...,K,然后进行下行功率分配,定义是对应用户的下行发射功率,得到下行发射功率向量若向量p满足每个基站的功率约束,则算法结束,输出p和T*,T*表示取得最优值的T,T表示下行单位波束成形向量构成的矩阵,若不满足,则采用现有技术中的ICBF算法重新计算波束成形向量并输出;
如果i>imax,采用现有技术中的ICBF算法重新计算波束成形向量并输出。
2.根据权利要求1所述的最大化吞吐量的多小区联合波束成形设计方法,其特征是,所述步骤5)中的子优化问题是指:
3.根据权利要求1所述的最大化吞吐量的多小区联合波束成形设计方法,其特征是,所述步骤6)中的子优化问题是指:
其中:
4.根据权利要求1所述的最大化吞吐量的多小区联合波束成形设计方法,其特征是,所述步骤7)中的终止条件是指:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110383525.4A CN102572864B (zh) | 2011-11-25 | 2011-11-25 | 最大化吞吐量的多小区联合波束成形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110383525.4A CN102572864B (zh) | 2011-11-25 | 2011-11-25 | 最大化吞吐量的多小区联合波束成形设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102572864A CN102572864A (zh) | 2012-07-11 |
CN102572864B true CN102572864B (zh) | 2015-04-22 |
Family
ID=46417069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110383525.4A Expired - Fee Related CN102572864B (zh) | 2011-11-25 | 2011-11-25 | 最大化吞吐量的多小区联合波束成形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102572864B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102905320B (zh) * | 2012-09-24 | 2015-04-01 | 上海交通大学 | 最大化多小区下行加权和速率的单调优化方法 |
CN103684560B (zh) * | 2013-12-04 | 2017-02-08 | 上海交通大学 | 多小区多用户系统中基于用户公平性的鲁棒预编码方法 |
CN104486775B (zh) * | 2015-01-08 | 2018-03-27 | 西安电子科技大学 | 基于用户位置的跟踪预测的有源天线阵列波束优化方法 |
CN105933979B (zh) * | 2016-04-12 | 2019-06-25 | 东南大学 | 多小区bdma传输功率分配方法 |
CN106559115A (zh) * | 2016-11-17 | 2017-04-05 | 河海大学 | 一种多用户mimo中继系统中的用户调度方法 |
CN106788637B (zh) * | 2016-12-22 | 2020-04-28 | 电子科技大学 | 一种d2d通信中传输模式和传输波束的联合优化方法 |
CN109451584B (zh) * | 2018-11-28 | 2021-02-09 | 电子科技大学 | 一种多天线数能一体化通信网络的最大化上行吞吐量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101442551A (zh) * | 2008-11-13 | 2009-05-27 | 上海交通大学 | 基于ieee802.15.4协议的传感器节点占空比独立自适应调节方法 |
CN101610107A (zh) * | 2005-12-20 | 2009-12-23 | 高通股份有限公司 | 用于反向链路发射波束成形的方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070189151A1 (en) * | 2006-02-10 | 2007-08-16 | Interdigital Technology Corporation | Method and apparatus for performing uplink transmission in a multiple-input multiple-output single carrier frequency division multiple access system |
-
2011
- 2011-11-25 CN CN201110383525.4A patent/CN102572864B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101610107A (zh) * | 2005-12-20 | 2009-12-23 | 高通股份有限公司 | 用于反向链路发射波束成形的方法和装置 |
CN101442551A (zh) * | 2008-11-13 | 2009-05-27 | 上海交通大学 | 基于ieee802.15.4协议的传感器节点占空比独立自适应调节方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102572864A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102572864B (zh) | 最大化吞吐量的多小区联合波束成形设计方法 | |
Omar et al. | Multiobjective optimization in 5G hybrid networks | |
Jia et al. | Analysis and optimization of an intelligent reflecting surface-assisted system with interference | |
CN108234101B (zh) | 能效最大化导频信号设计方法及大规模多天线系统 | |
Munir et al. | Energy efficient resource allocation in 5G hybrid heterogeneous networks: A game theoretic approach | |
CN104579443A (zh) | 一种基于协作多小区Massive MIMO系统的线性预编码方法 | |
Lin et al. | Dynamic base station formation for solving NLOS problem in 5G millimeter-wave communication | |
CN104104425A (zh) | 一种基于期望与干扰信号关系的多用户mimo自适应接收方法 | |
CN107241799B (zh) | 异构网络中能量有效的资源分配和干扰对齐联合方法 | |
CN113114317A (zh) | 一种irs辅助下行多用户通信系统的相移优化方法 | |
CN104579444B (zh) | 一种异构蜂窝网络中干扰对齐方法 | |
Akhtar et al. | Analysis of decoupled association in HetNets using soft frequency reuse scheme | |
CN101989869B (zh) | 一种多天线蜂窝网络的联合预编码与功率控制方法 | |
Luo et al. | On the study of sustainability and outage of swipt-enabled wireless communications | |
CN104954055A (zh) | 一种低复杂度的多用户信能同传系统能效优化方法 | |
CN115811734B (zh) | 混合卫星-地面无线网络中基于物理层安全的安全速率最大化方法 | |
CN112584403B (zh) | Noma小小区最大速率最小功率的联合优化方法 | |
Ni et al. | Performance analysis for large intelligent surface assisted vehicular networks | |
Heath | Comparing massive MIMO and mmWave MIMO | |
CN105577591B (zh) | 一种异构网中基于全双工通信的跨层串行干扰删除方法 | |
CN103595453B (zh) | 多小区时分复用无线系统波束成型方法 | |
Zhao et al. | Antenna location design for distributed-antenna based secondary systems | |
Rahman et al. | Multi-cell full-duplex wireless communication for dense urban deployment | |
CN104967470B (zh) | 一种单用户认知网络预编码优化设计与空间对齐方法 | |
Zhu et al. | Uplink pilot-to-data power ratio design based on user joint optimization algorithm in multi-cell massive mimo system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150422 Termination date: 20171125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |