CN107171704B - 一种大规模mimo***的上行功率控制方法及装置 - Google Patents
一种大规模mimo***的上行功率控制方法及装置 Download PDFInfo
- Publication number
- CN107171704B CN107171704B CN201710434668.0A CN201710434668A CN107171704B CN 107171704 B CN107171704 B CN 107171704B CN 201710434668 A CN201710434668 A CN 201710434668A CN 107171704 B CN107171704 B CN 107171704B
- Authority
- CN
- China
- Prior art keywords
- optimization function
- user equipment
- cell
- convex optimization
- convex
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005457 optimization Methods 0.000 claims abstract description 149
- 230000006870 function Effects 0.000 claims description 172
- 238000005562 fading Methods 0.000 claims description 38
- 239000011159 matrix material Substances 0.000 claims description 37
- 230000005540 biological transmission Effects 0.000 claims description 27
- 230000009466 transformation Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 238000012549 training Methods 0.000 description 5
- 238000010295 mobile communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000010432 diamond Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0426—Power distribution
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W52/00—Power management, e.g. TPC [Transmission Power Control], power saving or power classes
- H04W52/04—TPC
- H04W52/06—TPC algorithms
- H04W52/14—Separate analysis of uplink or downlink
- H04W52/146—Uplink power control
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W52/00—Power management, e.g. TPC [Transmission Power Control], power saving or power classes
- H04W52/04—TPC
- H04W52/18—TPC being performed according to specific parameters
- H04W52/26—TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service]
- H04W52/267—TPC being performed according to specific parameters using transmission rate or quality of service QoS [Quality of Service] taking into account the information rate
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Quality & Reliability (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供一种大规模MIMO***的上行功率控制方法及装置,该方法包括:S1,根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;S2,根据所述凸性优化函数,获取小区内用户设备的功率分配结果。这样在满足本小区内用户设备速率需求的同时,使得大规模MIMO***内所有用户设备的速率之和达到最大,充分发挥了大规模MIMO***的潜力。
Description
技术领域
本发明涉及无线通信技术领域,更具体地,涉及一种大规模MIMO***的上行功率控制方法及装置。
背景技术
多入多出(Multiple-Input Multiple-Output,MIMO)技术是指在发射端和接收端分别使用多个发射天线和接收天线,使信号通过发射端与接收端的多个天线传送和接收,从而改善通信质量。它能充分利用空间资源,通过多个天线实现多发多收,在不增加频谱资源和天线发射功率的情况下,可以成倍的提高***信道容量,显示出明显的优势,被视为下一代移动通信的核心技术。
大规模MIMO(Massive MIMO)时分双工***作为一种新颖的蜂窝网络结构成为当前的研究热点,在大规模MIMO***中,基站端有数量巨大的低功率小天线,天线数目远远超过同时调度的单天线用户设备数量,基站和用户设备之间通过时分双工进行通信。大规模天线阵列提供了很高的阵列增益和很高的空间分辨率。在较好的传播条件下,线性大规模MIMO***可以有效地消除小区内的干扰和噪声。
目前,现有的大规模MIMO***在功率分配时,将总功率平均分配给所有用户设备,这样不能充分挖掘大规模MIMO***在单用户功率控制方面的潜力;并且大规模MIMO***内所有用户设备的和速率也未达到最大,没有发挥大规模MIMO***的优势。
发明内容
针对上述的技术问题,本发明提供一种大规模MIMO***的上行功率控制方法及装置。
第一方面,本发明实施例提供一种大规模MIMO***的上行功率控制方法,包括:S1,根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;S2,根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
其中,在所述S1之前还包括:S0,根据获取的小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和大尺度衰落系数,获取***和速率上界。
其中,所述S0包括:S01,根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数,获取本小区基站的接收信号;S02,根据所述基站接收的信号获取信道状态矩阵,并根据所述信道状态矩阵获取所述小区内各用户设备的信噪比;S03,根据所述信噪比及詹森不等式,获得所述***和速率上界。
其中,所述S1中的将所述非凸性优化函数转化为凸性优化函数包括:S11,根据对数下界不等式和变量指数变换,将所述非凸性优化函数转化为半非凸性优化函数;S12,通过一阶泰勒级数展开,将所述半非凸性优化函数转化为所述凸性优化函数;其中,所述半非凸性优化函数为目标函数为凸、约束条件为非凸的函数。
其中,所述S2包括:S21,获取所述凸性优化函数的解,并根据所述解获取所述非凸性优化函数的解;S22,根据所述非凸性优化函数的解更新对数下界不等式的参数,然后再次基于获取的所述凸性优化函数的解,获取所述非凸性优化函数的解;S23,直到获取的所述非凸性优化函数第t次与第t-1次的解的差的范数平方小于预设阈值时,得到小区内用户设备的功率分配结果,其中t≥2。
其中,所述信噪比为:
其中,γlk为第l小区内的第k个用户设备的信噪比;plk为第l小区内的第k个用户设备的发射功率;为信道估计误差的自相关矩阵;为小区间干扰的自相关矩阵;为噪声的自相关矩阵;为信道状态矩阵; 为第i小区的第k个用户设备到第l个小区基站的路损和阴影衰落效应,即大尺度衰落系数;pik为第i小区内的第k个用户设备的发射功率。
相应地,所述和速率上界为:
其中,所述非凸性优化函数为:
所述凸性优化函数为:
第二方面,本发明提供一种大规模MIMO***的上行功率控制装置,包括:获取模块,用于根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数,获取***和速率上界;转化模块,用于根据所述和速率上界,建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;分配模块,用于根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
第三方面,本发明提供一种大规模MIMO***的上行功率控制装置,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至7任一所述的方法。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至7任一所述的方法。
本发明提供的一种大规模MIMO***的上行功率控制方法及装置,通过根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;然后根据所述凸性优化函数,获取小区内用户设备的功率分配结果。这样在满足本小区内用户设备速率需求的同时,使得大规模MIMO***内所有用户设备的速率之和达到最大,充分发挥了大规模MIMO***的潜力。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的大规模MIMO***的上行功率控制方法的流程图;
图2为图1所述的大规模MIMO***的上行功率控制方法中的大规模MIMO***的结构示意图;
图3为图1所述的大规模MIMO***的上行功率控制方法的根据凸性优化函数获取小区内用户设备的功率分配结果的流程图;
图4为本发明另一实施例提供的大规模MIMO***的上行功率控制方法中***和速率随基站天线数的变化图;
图5为本发明另一实施例提供的大规模MIMO***的上行功率控制方法中优化后的速率增益随基站天线数的变化图;
图6为本发明另一实施例提供的大规模MIMO***的上行功率控制方法中本小区内用户设备的总发射功率随基站天线数的变化图
图7为本发明实施例提供的大规模MIMO***的上行功率控制装置的结构框图;
图8为本发明另一实施例提供的大规模MIMO***的上行功率控制装置的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的大规模MIMO***的上行功率控制方法的流程图,如图1所示,该方法包括:S1,根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;S2,根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
其中,MIMO(Multiple-Input Multiple-Output)技术指在发射端和接收端分别使用多个发射天线和接收天线,使信号通过发射端与接收端的多个天线传送和接收,从而改善通信质量。它能充分利用空间资源,通过多个天线实现多发多收,在不增加频谱资源和天线发射功率的情况下,可以成倍的提高***信道容量,显示出明显的优势、被视为下一代移动通信的核心技术。
在本发明实施例中,以一个具有多个小区的大规模MIMO***为例,***中由六边形形状的L个小区组成,每个小区的中心有一个基站,基站有M根天线,且假定每个小区内有K个单天线用户,满足大规模MIMO***条件,即M>>K。另外,假定每个小区使用相同的频率资源,即频率再使用因子为1,整个***操作在时分双工(Time Division Duplexing,TDD)模式。图2为一个有3个小区的大规模MIMO***,假定在上行接收,每个基站预先不知道对用户的信道状态信息(CSI),并采用破零(ZF)滤波处理。CSI的获取通过导频训练获得,即在每个信道相干间隔(T)的开始,每个用户都向所在的基站发射一个长度为τ的导频序列,这里显然满足T>τ。小区内的用户使用正交的导频序列集,且正交导频序列集在所有相邻的小区中重复使用。
其中,凸函数就是一个定义在某个向量空间的凸子集C(区间)上的实值函数。
具体地,在控制大规模MIMO***的上行功率时,根据MIMO***的和速率上界建立非凸性优化函数,为了较方便的获得非凸性优化函数的解,将该非凸性优化函数转化为凸性优化函数。对凸性优化函数进行求解,根据求得的凸性优化函数的解获取到非凸性优化函数的解,即可以获取到小区内用户设备的功率分配结果。
本发明实施例与现有技术相比,通过根据***和速率上界建立非凸性优化函数,并将非凸性优化函数转化为凸性优化函数,然后通过求解凸性优化函数获取到非凸性优化函数的解,进而得到小区内用户设备的功率分配结果,这样在满足本小区内用户设备速率需求的同时,使得大规模MIMO***内所有用户设备的速率之和达到最大,充分发挥了大规模MIMO***的潜力。
在上述实施例的基础上,在所述S1之前还包括:S0,根据获取的小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和大尺度衰落系数,获取***和速率上界。
其中,大尺度衰落系数是指***内所有用户设备到本小区基站的路损和阴影衰落效应对应的参数。
具体地,在根据***和速率上界建立非凸性优化函数之前,需要先获取到***的和速率上界。在每个小区的用户设备同时向各自所在对的小区基站发射导频训练序列后,各个小区的基站都会获取到小区内用户设备的发射功率,并根据发射功率计算出大尺度衰落系数,然后***内其他小区基站会将小区内用户设备的发射功率和大尺度衰落系数反馈至本小区。然后根据本小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和所有用户设备到本小区基站的大尺度衰落系数,可以获取到***的和速率上界。
在本发明实施例中,通过根据本小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和所有用户设备到本小区基站的大尺度衰落系数,获取到***的和速率上界,这样保证了在***的速率之和达到最大时进行小区内用户设备的功率分配,充分发挥的大规模MIMO***的潜能。
在上述实施例的基础上,所述S0包括:S01,根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数,获取本小区基站的接收信号;S02,根据所述基站接收的信号获取信道状态矩阵,并根据所述信道状态矩阵获取所述小区内各用户设备的信噪比;S03,根据所述信噪比及詹森不等式,获得所述***和速率上界。
其中,信道状态信息(Channel State Information,CSI)在无线通信领域就是通信链路的信道属性,它描述了信号在每条传输路径上的衰弱因子,即信道增益矩阵H中每个元素的值,如信号散射(Scattering),环境衰弱(fading,multipath fading or shadowingfading),距离衰减(power decay of distance)等信息。CSI可以使通信***适应当前的信道条件,在多天线***中为高可靠性高速率的通信提供了保障。
其中,信噪比英文名称叫做SNR或S/N(SIGNAL-NOISE RATIO),又称为讯噪比,是指一个电子设备或者电子***中信号与噪声的比例。这里面的信号指的是来自设备外部需要通过这台设备进行处理的电子信号,噪声是指经过该设备后产生的原信号中并不存在的无规则的额外信号(或信息),并且该种信号并不随原信号的变化而变化。
其中,詹森不等式是关于凸性(convexity)的不等式,它给出积分的凸函数值和凸函数的积分值间的关系。其中,凸性是非常好的性质,在最优化问题里面,线性和非线性不是本质的区别,只有凸性才是。如果最优化的函数是凸的,那么局部最优就意味着全局最优,否则无法推得全局最优。
具体地,在根据获取的小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和***内所有用户设备到本小区基站的大尺度衰落系数,获取***和速率上界时,先根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和所有用户设备到本小区基站的大尺度衰落系数,获取本小区基站的接收信号。然后对本小区基站的接收信号进行处理得到所有用户设备到本小区基站的信道矩阵,并根据基站接收的信号和该信道矩阵可以获得本小区内各用户设备的信噪比。最后通过应用詹森不等式和根据该信噪比,以获得***的和速率上界。
在本发明实施例中,通过根据获得的信噪比和应用詹森不等式获得***的和速率上界,这样在满足小区内用户设备的功率需求的同时,可以使***的速率之和达到最大,充分体现了大规模MIMO***的优势。
在上述实施例的基础上,所述S1中的将所述非凸性优化函数转化为凸性优化函数包括:S11,根据对数下界不等式和变量指数变换,将所述非凸性优化函数转化为半非凸性优化函数;S12,通过一阶泰勒级数展开,将所述半非凸性优化函数转化为所述凸性优化函数;其中,所述半非凸性优化函数为目标函数为凸、约束条件为非凸的函数。
其中,对数下界不等式是指有最小值的对数不等式。
其中,指数变换是指将变量用其对应的指数替换。
其中,泰勒级数用无限项连加式——级数来表示一个函数,这些相加的项由函数在某一点的导数求得。泰勒级数是以于1715年发表了泰勒公式的英国数学家布鲁克·泰勒(Sir Brook Taylor)的名字来命名的。通过函数在自变量零点的导数求得的泰勒级数又叫做迈克劳林级数,以苏格兰数学家科林·麦克劳林的名字命名。泰勒级数在近似计算中有重要作用。
具体地,在根据获取的和速率上界建立好非凸性欧化函数之后,需要将建立的非凸性优化函数转化为凸性优化函数,首先根据对数下界不等式和变量的指数变换,将非凸性优化函数转化为半非凸性优化函数,即目标函数为凸、约束条件为非凸的函数。然后再通过一阶泰勒级数展开,将半非凸性优化函数转化为凸性优化函数,这样就可以通过使用离线优化工具包求解凸性优化函数,例如,用CVX获取凸性优化函数的解,然后根据获取的凸性优化函数的解来得到非凸性优化函数的解,进而获得本小区内用户设备的功率分配结果。
在本发明实施例中,通过根据对数不等式、变量的指数变换以及一阶泰勒级数展开,将非凸性优化函数转化为凸性优化函数,这样可以较简单和快速的得到凸性优化函数的解,进而得到非凸性优化函数的解,达到简化求解过程的目的,并且使得获得本小区内用户设备的功率分配更准确。
在上述实施例的基础上,结合图3,所述S2包括:S21,获取所述凸性优化函数的解,并根据所述解获取所述非凸性优化函数的解;S22,根据所述非凸性优化函数的解更新对数下界不等式的参数,然后再次基于获取的所述凸性优化函数的解,获取所述非凸性优化函数的解;S23,直到获取的所述非凸性优化函数第t次与第t-1次的解的差的范数平方小于预设阈值时,得到小区内用户设备的功率分配结果,其中t≥2。
其中,范数是数学中的一种基本概念,在泛函分析中,范数是一种定义在赋范线性空间中函数,满足相应条件后的函数都可以被称为范数。
具体地,在根据***和速率上界建立好非凸性优化函数,并将非凸性优化函数转化为凸性优化函数之后,先求解转化得到的凸性优化函数,并根据获得的凸性优化函数的解得到之前非凸性优化函数的解。然后根据获取的非凸性优化函数的解更新对数下界不等式中的参数,然后根据更新后的对数下界不等式再次对凸性优化函数进行求解,并根据求得的解再次对凸性优化函数进行求解。这样一直循环直到获取的非凸性优化函数第t次与第t-1次的解的差的范数平方小区预设阈值,其中t≥2,例如,非凸性优化函数第t次与第t-1次的解的差的范数平方小于10-6,则可以结束循环,将该第t次获得的非凸性优化函数的解作为最终解,即可以得到小区内用户设备的功率分配结果。
在本发明实施例中,通过不断的更新对数不等式的参数,并且对每次获得的非凸性优化函数的解进行检测,直到相近两次非凸性优化函数的解满足条件时,才停止更新对数不等式,并将最后一次非凸性优化函数的解作为本小区内用户设备的功率分配结果,这样可以提高本小区内用户设备功率分配的准确性,使得功率分配更合理。
在上述实施例的基础上,所述信噪比为:
其中,γlk为第l小区内的第k个用户设备的信噪比;plk为第l小区内的第k个用户设备的发射功率;为信道估计误差的自相关矩阵;为小区间干扰的自相关矩阵;为噪声的自相关矩阵;为估计的信道矩阵; 为第i小区的第k个用户设备到第l个小区基站的路损和阴影衰落效应,即大尺度衰落系数;pik为第i小区内的第k个用户设备的发射功率;相应地,所述和速率上界为:
具体地,在获取到小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数之后,则小区内的基站接收的数据为:
其中,Pl是对角矩阵,其对角元素[Pl]k,k=plk表示第l个小区内的第k个用户的数据发射功率,si是第i个小区内K个用户发射的数据向量,满足nl表示第l个小区基站M根天线上接收到的独立加性高斯噪声向量,假定噪声的功率已被归一化,即
利用估计的信道矩阵,基站执行破零(zero forcing,ZF)检测后输出的信号可写成,
其中,I1,I2和I3分别代表由信道估计误差带来的干扰,小区间干扰和噪声。它们的均值都为零,自相关矩阵计算如下,
根据基站输出的信号可以获得第l小区内的第k个用户设备的信噪比,用γlk表示,即
根据香农公式,于是第l个小区遍历的上行可达和速率可表示为:
其中,M-K+1是形状参数影响着分布的形状不仅包括位置移动及伸展/收缩,δllk是速率参数,该值越大概率密度函数的形状拉伸的越厉害,该值越小概率密度的形状函数越集中。
利用f(x)=log2(1+x)是x的凹函数的结论,并对和速率表达式应用詹森不等式,可以得到,
由(11)和(12)可得:
根据(10)式容易计算出E(x)=(M-K+1)δllk。将这个结果带入(13),可得到上行速率的上界:
其中,在信道矩阵估计阶段,每个小区的用户设备同时都向各自所在的小区基站发射长度为τ的导频训练序列,用符号表示每个小区内K个用户的正交导频序列集,满足ΦHΦ=IK。用户设备发射完导频序列后,第l小区的基站接收的训练信号可以表示成,
其中,Yp,l的下标p表示导频(pilot)信号,下标l指第l个小区;表示第i小区内的K个用户到第l个小区基站的信道矩阵;表示用户到基站的瑞利衰落效应,其中的元素都是均值为零、方差为1的独立同分布(i.i.d.)复高斯随机变量,即对角阵表示第i小区内的K个用户到第l个小区基站的路损和阴影衰落效应,并假定是已知的;对角阵Q的对角元素[Qi]k,k=qik表示第i小区内的第k个用户发射的导频功率;矩阵Nl表示第l个小区基站的M根天线上接收的噪声,满足
为了估计信道状态信息(Channel State Information,CSI),基站首先对接收到的信号乘以正交的导频序列集Φ并进行预处理可得,
由于用户设备位置独立,所以Gll的各列向量gllk,k=1,2,…,K也是独立的,因而可以取出Gll中的某个列向量单独进行估计,于是
ek表示第k个元素为1,其余元素为0的K×1维的列向量。注意:由于正交变换不影响元素的分布特性,故上式中省略了Nl后的正交矩阵Φ。
利用MMSE估计理论,可得gllk的MMSE估计结果为:
在上述实施例的基础上,所述非凸性优化函数为:
所述凸性优化函数为:
具体地,在根据***和速率上界建立好非凸性优化函数时,对于其他小区固定而已知的发射功率,通过下行控制优化用户的发射功率以获得最大小区吞吐量,则非凸性优化函数可建模为,
由于上面的目标函数
将下面的对数下界不等式,
log2(1+z)≥ξlog2z+μ (24)
当近似常数作如下选择时,上面的不等式在z=z0时将非常的紧,
其中,
给定ξk和μk,由于(29)式中第一项是线性函数和凹函数之和(注意log-sum-exp是凸函数),因此,问题的目标按说关于变量是凹函数。很明显,(28)式中的第一个约束条件也是凸函数,但第二个约束条件还是非凸的,此时得到的是半非凸性优化函数。则对第二个约束条件的左边用一阶泰勒级数展开,通过迭代逼近该约束条件,于是可转换为下面的凸性优化函数
通过不断地进行迭代,直到非凸性优化函数的解收敛,例如,非凸性优化函数第t次与第t-1次的解的差的范数平方小于预设阈值,例如,预设阈值为10-6。这样就可以得到最终的非凸性优化函数的解,即可以得到本小区用户设备的功率分配结果。
以下对本发明实施例进行举例说明,但不限制本发明的保护范围。假设有3个小区,每个小区有5个用户设备,以第1个小区为例,设第1个小区的用户设备到本基站的大尺度衰落系数为β11=[1,0.8,0.4,0.1,0.05],第2个小区的用户设备到第1个小区的基站的大尺度衰落系数为β12=10-2×[0.17,0.15,0.14,0.105,0.08],第3个小区的用户设备到第1个小区的基站的大尺度衰落系数为β13=10-2×[0.18,0.16,0.13,0.07,0.06]。第1个小区内的每个用户设备发射功率都为30dB,第2个和第3个小区的5个用户设备的发射功率取相同的一组值,即5dB,10dB,15dB,20dB,30dB。每个用户的最大发射功率为50dB,最小发射功率为5dB。可以得到,如图4所示的***和速率随基站天线数的变化情况;如图5所示的优化功率控制获得的速率增益随基站天线数的变化情况;如图6所示的优化功率控制后的用户设备的发射功率之和随基站天线数的变化情况。
图4中带星号的虚线是每个用户设备以30dB的功率发射时获得的仿真和速率;菱形的虚线是每个用户设备以30dB的功率发射时得到的可达和速率的上界;带星号的实线是通过本发明实施例提供的方法进行优化功率控制以后的得到的仿真和速率,带菱形的实线是通过本发明实施例提供的方法进行优化功率控制以后的得到的可达和速率上界。从图4可以看出,无论是否进行功率控制,增加基站天线数都可以提高***的可达和速率;且得到的闭式的可达和速率上界表达式在整个基站天线范围内都能精确地匹配仿真得到的速率;基站对用户设备的发射功率进行优化控制以后,在考虑的基站天线数范围内都可以使得***获得的可达和速率得到明显地提高。
图5给出了速率增益ΔR=Ropt-R随基站天线数的变化情况,其中,Ropt表示基站对用户设备进行优化功率控制后***获得的可达和速率,R表示基站对用户设备没有进行优化功率控制,即基站都按照30dB发射功率时***获得的可达和速率。从图5可以明显地看出,当基站的天线数增加时,功率控制的优势也越来越明显,即功率控制后获得的速率增益随着基站天线数的增加也在不断增加。
图6给出了优化功率控制后第l个小区所有用户设备总的发射功率随基站天线数的变化情况。从图6可以看出,随着基站天线数的增加,所有用户设备总的发射功率在不断减小,特别是基站天线数比较小时,增加天线数会使总的发射功率得到明显降低。应该注意的是,即使是总的发射功率降低了,从图5可知,***获得的可达和速率是在增大的,所以说,使用本发明实施例提供的功率控制方法,不仅可以提高***的可达和速率,还可以减少用户设备的辐射功率。
图7为本发明实施例提供的大规模MIMO***的上行功率控制装置的结构框图,如图7所示,该功率控制装置包括:获取模块701、转化模块702及分配模块703。
具体地,在控制大规模MIMO***的上行功率时,每个小区的用户设备同时向各自所在对的小区基站发射导频训练序列,然后获取模块701会获取到小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率以及所有用户设备到本小区基站的大尺度衰落系数。获取模块701则根据本小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和大尺度衰落系数,可以获取到***的和速率上界。然后转化模块702根据和速率上界建立非凸性优化函数,且将非凸性优化函数转化为凸性优化函数。最后分配模块703对凸性优化函数进行求解,进而获取到非凸性优化函数的解,即得到小区内用户设备的功率分配结果。
本发明实施例与现有技术相比,通过转化模块根据***和速率上界建立非凸性优化函数,并将非凸性优化函数转化为凸性优化函数,然后通过分配模块求解凸性优化函数获取到非凸性优化函数的解,进而得到小区内用户设备的功率分配结果,这样在满足本小区内用户设备速率需求的同时,使得大规模MIMO***内所有用户设备的速率之和达到最大,充分发挥了大规模MIMO***的潜力。
图8为本发明另一实施例提供的大规模MIMO***的上行功率控制装置的结构框图,如图8所示,该识别装置包括:处理器801、存储器802和总线803;
其中,处理器801和存储器802通过总线83完成相互间的通信;所述处理器801用于调用所述存储器802中的程序指令,以执行上述各方法实施例所提供的方法,例如包括:根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
在另一实施例中,提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行上述各方法实施例所提供的方法,例如包括:根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种大规模MIMO***的上行功率控制方法,其特征在于,包括:
S1,根据***和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;
S2,根据所述凸性优化函数,获取小区内用户设备的功率分配结果;
在所述S1之前还包括:S0,根据获取的小区内用户设备的发射功率、***内其他小区基站反馈的用户设备的发射功率和大尺度衰落系数,获取***和速率上界。
2.根据权利要求1所述的方法,其特征在于,所述S0包括:
S01,根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数,获取本小区基站的接收信号;
S02,根据所述基站接收的信号获取信道状态矩阵,并根据所述信道状态矩阵获取所述小区内各用户设备的信噪比;
S03,根据所述信噪比及詹森不等式,获得所述***和速率上界。
3.根据权利要求1所述的方法,其特征在于,所述S1中的将所述非凸性优化函数转化为凸性优化函数包括:
S11,根据对数下界不等式和变量指数变换,将所述非凸性优化函数转化为半非凸性优化函数;
S12,通过一阶泰勒级数展开,将所述半非凸性优化函数转化为所述凸性优化函数;
其中,所述半非凸性优化函数为目标函数为凸、约束条件为非凸的函数。
4.根据权利要求1所述的方法,其特征在于,所述S2包括:
S21,获取所述凸性优化函数的解,并根据所述解获取所述非凸性优化函数的解;
S22,根据所述非凸性优化函数的解更新对数下界不等式的参数,然后再次基于获取的所述凸性优化函数的解,获取所述非凸性优化函数的解;
S23,直到获取的所述非凸性优化函数第t次与第t-1次的解的差的范数平方小于预设阈值时,得到小区内用户设备的功率分配结果,其中t≥2。
7.一种大规模MIMO***的上行功率控制装置,其特征在于,包括:
获取模块,用于根据小区内用户设备的发射功率、***内其他小区基站反馈的用户设备发射功率和大尺度衰落系数,获取***和速率上界;
转化模块,用于根据所述和速率上界建立非凸性优化函数,并将所述非凸性优化函数转化为凸性优化函数;
分配模块,用于根据所述凸性优化函数,获取小区内用户设备的功率分配结果。
8.一种大规模MIMO***的上行功率控制装置,其特征在于,包括:至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中,
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至6任一所述的方法。
9.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至6任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710434668.0A CN107171704B (zh) | 2017-06-09 | 2017-06-09 | 一种大规模mimo***的上行功率控制方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710434668.0A CN107171704B (zh) | 2017-06-09 | 2017-06-09 | 一种大规模mimo***的上行功率控制方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107171704A CN107171704A (zh) | 2017-09-15 |
CN107171704B true CN107171704B (zh) | 2020-04-14 |
Family
ID=59825570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710434668.0A Active CN107171704B (zh) | 2017-06-09 | 2017-06-09 | 一种大规模mimo***的上行功率控制方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107171704B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108880734B (zh) * | 2018-04-28 | 2020-05-15 | 哈尔滨工程大学 | 量子回溯搜索优化的CCFD-Massive MIMO***功率分配方法 |
CN110149128B (zh) * | 2019-05-09 | 2021-01-15 | 中国科学院计算技术研究所 | 一种mimo-noma***中功率调整方法及装置 |
CN110166088B (zh) * | 2019-05-15 | 2021-09-28 | 南京邮电大学 | 以用户为中心的无小区mimo***的功率控制算法 |
CN111148208B (zh) * | 2020-01-06 | 2022-04-19 | 武汉虹信科技发展有限责任公司 | 一种功率平衡的方法及装置 |
CN114285444B (zh) * | 2021-09-30 | 2023-03-28 | 南京邮电大学 | 一种去蜂窝大规模mimo***的功率优化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104393956A (zh) * | 2014-11-26 | 2015-03-04 | 北京邮电大学 | 一种用于无线携能通信***最大化和速率预编码方法 |
CN105392161A (zh) * | 2015-10-15 | 2016-03-09 | 北京航空航天大学 | 协作异构网络中的用户接入与功率控制联合优化方法 |
CN105430728A (zh) * | 2015-12-08 | 2016-03-23 | 同济大学 | 用于微小区的基于最优化理论的上行链路功率控制方法 |
CN105467365A (zh) * | 2015-12-08 | 2016-04-06 | 中国人民解放军信息工程大学 | 一种改善mimo雷达doa估计性能的低旁瓣发射方向图设计方法 |
US9369925B2 (en) * | 2012-10-30 | 2016-06-14 | Futurewei Technologies, Inc. | System and methods to achieve optimum efficiency-jain fairness in wireless systems |
CN106603133A (zh) * | 2016-12-28 | 2017-04-26 | 北京邮电大学 | 基于迫零预编码的非匹配信道功率分配方法及*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8218422B2 (en) * | 2008-06-03 | 2012-07-10 | Nec Laboratories America, Inc. | Coordinated linear beamforming in downlink multi-cell wireless networks |
US8934557B2 (en) * | 2010-06-30 | 2015-01-13 | Telefonaktiebolaget L M Ericsson (Publ) | Statistical joint precoding in multi-cell, multi-user MIMO |
-
2017
- 2017-06-09 CN CN201710434668.0A patent/CN107171704B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9369925B2 (en) * | 2012-10-30 | 2016-06-14 | Futurewei Technologies, Inc. | System and methods to achieve optimum efficiency-jain fairness in wireless systems |
CN104393956A (zh) * | 2014-11-26 | 2015-03-04 | 北京邮电大学 | 一种用于无线携能通信***最大化和速率预编码方法 |
CN105392161A (zh) * | 2015-10-15 | 2016-03-09 | 北京航空航天大学 | 协作异构网络中的用户接入与功率控制联合优化方法 |
CN105430728A (zh) * | 2015-12-08 | 2016-03-23 | 同济大学 | 用于微小区的基于最优化理论的上行链路功率控制方法 |
CN105467365A (zh) * | 2015-12-08 | 2016-04-06 | 中国人民解放军信息工程大学 | 一种改善mimo雷达doa估计性能的低旁瓣发射方向图设计方法 |
CN106603133A (zh) * | 2016-12-28 | 2017-04-26 | 北京邮电大学 | 基于迫零预编码的非匹配信道功率分配方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN107171704A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107171704B (zh) | 一种大规模mimo***的上行功率控制方法及装置 | |
Ashikhmin et al. | Interference reduction in multi-cell massive MIMO systems I: Large-scale fading precoding and decoding | |
KR101424697B1 (ko) | 다중 사용자 다중 입출력 무선통신 시스템에서 간섭 제거를위한 신호 처리 장치 및 방법 | |
US9559759B2 (en) | System and method for massive MIMO communication | |
JP5666581B2 (ja) | Mu−mimo通信システムの送信機のためのプリコーディング方法 | |
EP2248276B1 (en) | Precoding for multiple antennas | |
US9450787B2 (en) | System and method for early termination in iterative null-space directed singular value decomposition for MIMO | |
KR102334011B1 (ko) | 무선 통신 시스템에서 머신 러닝 기반 제한된 피드백 방법 및 장치 | |
US9923614B2 (en) | Multiple-antenna data transmission method, base station, user equipment, and system | |
EP3185434B1 (en) | Method and device for beamforming | |
US9918317B2 (en) | Apparatus configured to approximate a power coefficient in a cell-free massive MIMO wireless system and method of performing same | |
EP2795807B1 (en) | Downlink transmission in a mu-mimo system | |
US8559545B2 (en) | Channel information feedback method and apparatus thereof | |
WO2017132984A1 (en) | Method and apparatus of topological pilot decontamination for massive mimo systems | |
US9072116B2 (en) | Systems and methods for reducing complexity in modulation coding scheme (MCS) adaptation | |
Müller et al. | Analysis of blind pilot decontamination | |
Abdelfatah et al. | A Study on the Basics Processes of Massive MIMO. | |
EP4304237A1 (en) | Information reporting method and apparatus, first device, and second device | |
US10063396B2 (en) | Method and apparatus of topological pilot decontamination for massive MIMO systems | |
US8995401B2 (en) | Multiple antenna method and apparatus for reducing inter-cell interference in multi-user wireless systems | |
KR101857671B1 (ko) | Mimo 송신기에서 re 그룹을 형성하는 방법 | |
US9680550B2 (en) | Method and apparatus for transmitting feedback about channel having spatial correlation | |
KR101857670B1 (ko) | Mimo 수신기에서 re 그룹을 형성하여 수신 신호를 처리하는 방법 | |
CN104253639A (zh) | 获取信道质量指示的方法及装置 | |
Wu et al. | Enhanced time-shifted pilot based channel estimation in massive MIMO systems with finite number of antennas |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |