发明名称 协作异构网络中的用户接入与功率控制联合优化方法
摘要 本发明公开了一种协作异构网络中的用户接入与功率控制联合优化方法,属于无线通信技术领域。本方法首先获取各基站到所有用户的下行信道信息,计算预编码向量和等效信道能量,然后通过回传链路共享给中心处理器;然后中心处理器收集用户数据率需求、各基站的最大功率、发射天线数,以及基站端通过回传链路共享的等效信道能量,将所涉及的用户接入和功率控制联合优化问题,转化为一个凸优化问题进行求解;中心处理器将求解得到的各基站的功率控制结果和用户接入结果通过回传链路通知每个基站。本发明能够有效地平衡网络中各小区的业务负载,能够大幅降低系统的总发射功率,并能够支持更高的用户数据率需求,用户可任意接入多个基站。
申请公布号 CN105392161B 申请公布日期 2017.05.03
申请号 CN201510671784.5 申请日期 2015.10.15
申请人 北京航空航天大学 发明人 赵景;韩圣千;杨晨阳
分类号 H04W28/08(2009.01)I;H04W52/02(2009.01)I 主分类号 H04W28/08(2009.01)I
代理机构 北京永创新实专利事务所 11121 代理人 祗志洁
主权项 一种协作异构网络中的用户接入与功率控制联合优化方法,其特征在于,包括以下步骤:步骤1:对每个基站,获取它到所有用户的下行信道信息,计算预编码向量和等效信道能量,然后通过回传链路共享给中心处理器;基站包括宏基站和微基站;设某个宏基站的覆盖范围内随机部署B个微基站,共同服务K个用户;宏基站有N<sub>M</sub>个天线,最大发射功率为P<sub>M</sub>;微基站有N<sub>S</sub>个天线,最大发射功率为P<sub>S</sub>,用户端只有一个天线;设宏基站到用户k的信道向量为<img file="FDA0001231676540000011.GIF" wi="227" he="54" />微基站到用户k的信道向量为<img file="FDA0001231676540000012.GIF" wi="226" he="55" />宏基站到用户k的预编码向量为<img file="FDA0001231676540000013.GIF" wi="257" he="55" />微基站到用户k的预编码向量为<img file="FDA0001231676540000014.GIF" wi="227" he="55" />所有基站到用户k的预编码向量为<img file="FDA0001231676540000015.GIF" wi="1083" he="67" />p<sub>kb</sub>为基站b分配给用户k的发射功率,v<sub>kb</sub>为归一化的基站b到用户k的发射预编码向量,b=0,1,2,…,B,当b取0时表示宏基站,当b不为0时表示微基站;<img file="FDA0001231676540000016.GIF" wi="109" he="55" />表示N<sub>M</sub>×1维复向量空间,<img file="FDA0001231676540000017.GIF" wi="101" he="55" />表示N<sub>S</sub>×1维复向量空间,<img file="FDA0001231676540000018.GIF" wi="205" he="55" />表示(N<sub>M</sub>+BN<sub>S</sub>)×1维复向量空间;k=1,2,…,K;每个基站获得它到所有用户之间的信道信息,采用基于信漏噪比最大化准则的发射预编码,则基站b根据已知的信道信息得到v<sub>kb</sub>为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>=</mo><mfrac><mrow><msup><mrow><mo>(</mo><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>K</mi></msubsup><msub><mi>h</mi><mrow><mi>i</mi><mi>b</mi></mrow></msub><msubsup><mi>h</mi><mrow><mi>i</mi><mi>b</mi></mrow><mi>H</mi></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>h</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></mrow><mrow><mo>|</mo><mo>|</mo><msup><mrow><mo>(</mo><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>K</mi></msubsup><msub><mi>h</mi><mrow><mi>i</mi><mi>b</mi></mrow></msub><msubsup><mi>h</mi><mrow><mi>i</mi><mi>b</mi></mrow><mi>H</mi></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>h</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mi>b</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>B</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001231676540000019.GIF" wi="1481" he="231" /></maths>h<sub>ib</sub>为基站b到用户i的信道向量;用户k收到来自基站b的期望信号等效信道能量表示为<img file="FDA00012316765400000110.GIF" wi="298" he="63" />由基站b发送给用户i的干扰信号等效信道能量表示为<img file="FDA00012316765400000111.GIF" wi="281" he="63" />v<sub>ib</sub>为归一化的微基站b到用户i的发射预编码向量;采用功率控制变量p<sub>kb</sub>来表示用户k到基站b的接入情况,当用户k未接入基站b时,p<sub>kb</sub>=0,否则p<sub>kb</sub>>0;定义用户k的接入基站集合为<img file="FDA00012316765400000112.GIF" wi="628" he="63" />步骤2:中心处理器收集用户数据率需求、各基站的最大功率、发射天线数,以及基站端通过回传链路共享的等效信道能量,然后对所有协作宏、微基站的功率控制和所有用户的接入方式进行优化,具体包括:步骤2.1,在满足用户数据率要求、单基站功率约束以及给定每个用户最多接入基站个数的限制下,以最小化系统总发射功率为准则,得到用户接入和功率控制联合优化模型;所述的用户接入和功率控制联合优化模型表示为:<img file="FDA0001231676540000021.GIF" wi="1643" he="454" />其中,(2.1)是用户k的最小信干噪比约束,γ<sub>k</sub>是用户k要求的最小信干噪比,(2.2)和(2.3)分别表示宏基站和微基站的最大发射功率约束,(2.4)限制用户k最多只能接入L<sub>k</sub>个协作基站,<img file="FDA0001231676540000027.GIF" wi="66" he="67" />表示集合<img file="FDA0001231676540000023.GIF" wi="51" he="54" />中的元素个数;步骤2.2,将步骤2.1中的非凸优化问题转化为凸优化问题,再进行求解,获得用户接入结果和各基站的功率控制结果;定义函数<img file="FDA0001231676540000024.GIF" wi="322" he="110" />其中t→0是一个极小值;利用函数f(p<sub>kb</sub>)对约束(2.4)进行连续化表示,表示为:<img file="FDA0001231676540000025.GIF" wi="772" he="125" />进而将公式(2)表示为:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><munder><mi>min</mi><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></munder></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mfrac><mrow><msub><mi>&gamma;</mi><mi>k</mi></msub><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>K</mi></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msub><mi>&alpha;</mi><mrow><mi>k</mi><mi>i</mi><mi>b</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mi>b</mi></mrow></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msub><mi>&alpha;</mi><mrow><mi>k</mi><mi>k</mi><mi>b</mi></mrow></msub><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></mrow></mfrac><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>3.1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>M</mi></msub></mrow></mtd><mtd><mrow><mo>(</mo><mn>3.2</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>S</mi></msub><mo>,</mo><mi>b</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>B</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>3.3</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><mfrac><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mrow><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>+</mo><mi>t</mi></mrow></mfrac><mo>&le;</mo><msub><mi>L</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>3.4</mn><mo>)</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001231676540000026.GIF" wi="1732" he="639" /></maths>利用单项式近似的方法将约束(3.1)和(3.4)中不等式左侧的分母部分由现在的多项式转变为关于变量p<sub>kb</sub>的单项式,将优化问题转变为标准的几何规划问题,同时,引入一组松弛变量β<sub>k</sub>,k=1,...,K,β<sub>k</sub>≥1,将得到的几何规划问题进一步松弛为:<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><munder><mi>min</mi><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></munder></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>+</mo><msubsup><mi>&lambda;&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>&beta;</mi><mi>k</mi></msub></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mfrac><mrow><msub><mi>&gamma;</mi><mi>k</mi></msub><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>k</mi></mrow><mi>K</mi></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msub><mi>&alpha;</mi><mrow><mi>k</mi><mi>i</mi><mi>b</mi></mrow></msub><msub><mi>p</mi><mrow><mi>i</mi><mi>b</mi></mrow></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup></mrow><mrow><mrow><mo>(</mo><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><msub><mi>&alpha;</mi><mrow><mi>k</mi><mi>k</mi><mi>b</mi></mrow></msub><msubsup><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup></mrow><mo>)</mo></mrow><munderover><mi>&Pi;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></munderover><msup><mrow><mo>(</mo><mfrac><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><msubsup><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup></mfrac><mo>)</mo></mrow><msub><mi>a</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></msup></mrow></mfrac><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>4.1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>M</mi></msub></mrow></mtd><mtd><mrow><mo>(</mo><mn>4.2</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></msubsup><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mi>S</mi></msub><mo>,</mo><mi>b</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>B</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>4.3</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>0</mn></mrow><mi>B</mi></msubsup><mfrac><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><mrow><mrow><mo>(</mo><mrow><msubsup><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>+</mo><mi>t</mi></mrow><mo>)</mo></mrow><msup><mrow><mo>(</mo><mfrac><msub><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub><msubsup><mi>p</mi><mrow><mi>k</mi><mi>b</mi></mrow><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup></mfrac><mo>)</mo></mrow><msub><mi>c</mi><mrow><mi>k</mi><mi>b</mi></mrow></msub></msup></mrow></mfrac><mo>&le;</mo><msub><mi>&beta;</mi><mi>k</mi></msub><msub><mi>L</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>4.4</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><msub><mi>&beta;</mi><mi>k</mi></msub><mo>&GreaterEqual;</mo><mn>1</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>K</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>4.5</mn><mo>)</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001231676540000031.GIF" wi="1750" he="880" /></maths>其中,<img file="FDA0001231676540000032.GIF" wi="78" he="69" />表示第l次迭代得到的功率值,<img file="FDA0001231676540000033.GIF" wi="381" he="155" />是利用单项式近似方法得到的中间变量;目标函数里的λ是一个常数;步骤3:中心处理器将各基站的功率控制结果和用户接入结果通过回传链路通知每个基站,每个基站根据用户接入情况从核心网请求用户数据,并根据收到的功率控制结果和步骤1中设计的预编码来发送接入的用户的数据。
地址 100191 北京市海淀区学院路37号