发明名称 一种基于改进粒子群的超密集异构网络最优功率协调方法
摘要 本发明提供了一种基于改进粒子群的超密集异构网络最优功率协调方法。该方法联合调整微微站的发送功率,避免密集部署的微微站以最大功率发送对宏站、邻居微微站构成的严重干扰,最大化系统吞吐量。方法考虑用户服务小区在功率调整过程中的变化,并通过邻域局部搜索和多次初始化过程改进粒子群,加快收敛速度,提高功率解的质量。提出的改进通过邻域局部搜索过程保证功率解的局部最优性,并通过多次初始化过程进一步保证得到的功率解的全局最优性,能够以更少的迭代次数获得更高的系统吞吐量,在不影响子载波分配自由度、不损失资源利用率的基础上协调同层和跨层干扰,从功率调整的角度最大程度地提升系统频谱效率和总体吞吐量。
申请公布号 CN104066096A 申请公布日期 2014.09.24
申请号 CN201410315609.8 申请日期 2014.07.03
申请人 东南大学 发明人 潘志文;蒋慧琳;刘楠;尤肖虎
分类号 H04W16/14(2009.01)I;H04W52/04(2009.01)I 主分类号 H04W16/14(2009.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 一种基于改进粒子群的超密集异构网络最优功率协调方法,包括如下步骤:第一步,采集网络构成信息,初始化改进粒子群的参数。采集网络中的宏站个数M、微微站个数I及用户个数U。将站点集合记为C={C<sub>m</sub>,C<sub>p</sub>},其中宏站集合C<sub>m</sub>={m<sub>1</sub>,m<sub>2</sub>,...,m<sub>M</sub>},微微站集合C<sub>p</sub>={p<sub>1</sub>,p<sub>2</sub>,...,p<sub>I</sub>}。初始化改进粒子群的参数:迭代次数记为t,重新初始化次数记为s,初始化最大迭代次数和最大重新初始化次数分别为t<sub>iter</sub>、t<sub>res</sub>,局部搜索半径调整因子ξ&gt;1,当前迭代次数t=0,当前重新初始化次数s=0,当前局部搜索半径r<sub>Q</sub>(t)=1,搜索成功次数N<sub>s</sub>(t)=1,搜索失败次数N<sub>f</sub>(t)=1,成功失败比值门限η<sub>th</sub>>1。第二步:设置N种微微站候选发送功率集合<img file="FDA0000532169770000011.GIF" wi="427" he="76" />及功率调整尺度集合<img file="FDA0000532169770000012.GIF" wi="438" he="75" />每个候选发送功率集合和功率调整尺度集合包含I个微微站的候选发送功率和功率调整尺度,即<img file="FDA0000532169770000013.GIF" wi="698" he="79" /><img file="FDA0000532169770000014.GIF" wi="694" he="84" />微微站的候选发送功率和功率调整尺度需满足<img file="FDA0000532169770000015.GIF" wi="785" he="105" />其中<img file="FDA0000532169770000016.GIF" wi="89" he="85" />为微微站能容忍的最大发送功率。第三步:计算各候选功率集合对应的用户可达速率<img file="FDA0000532169770000017.GIF" wi="103" he="78" />及系统吞吐量T<sup>n</sup>(t)。针对当前每个候选功率集合<img file="FDA0000532169770000018.GIF" wi="154" he="75" />根据公式(1)为用户u(u∈U)计算其接收到的来自各个站点的参考信号接收功率RSRP<sub>u,c</sub>(u∈U,c∈C),RSRP<sub>u,c</sub>=P<sub>c</sub>G<sub>u,c</sub>   (1)其中G<sub>u,c</sub>为用户u与站点c间的信道增益。将RSRP最大的小区作为用户u的服务小区,记为C<sub>u</sub>,并根据公式(2)、(3)计算各候选功率集合下各用户的可达速率<img file="FDA0000532169770000019.GIF" wi="109" he="74" />及系统吞吐量<img file="FDA00005321697700000110.GIF" wi="236" he="96" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>r</mi><mi>u</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>W</mi><mi>u</mi><mi>n</mi></msubsup><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msubsup><mi>P</mi><msub><mi>C</mi><mi>u</mi></msub><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>G</mi><mrow><mi>u</mi><mo>,</mo><msub><mi>C</mi><mi>u</mi></msub></mrow></msub></mrow><mrow><msub><mi>&Sigma;</mi><mrow><msub><mi>C</mi><mi>I</mi></msub><mo>&Element;</mo><mi>C</mi><mo>,</mo><msub><mi>C</mi><mi>I</mi></msub><mo>&NotEqual;</mo><msub><mi>C</mi><mi>u</mi></msub></mrow></msub><msubsup><mi>P</mi><msub><mi>C</mi><mi>I</mi></msub><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>G</mi><mrow><mi>u</mi><mo>,</mo><msub><mi>C</mi><mi>I</mi></msub></mrow></msub><mo>+</mo><msub><mi>N</mi><mn>0</mn></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005321697700000111.GIF" wi="1222" he="191" /></maths><img file="FDA0000532169770000021.GIF" wi="1064" he="94" />其中<img file="FDA00005321697700000221.GIF" wi="74" he="68" />为用户u在功率集合<img file="FDA0000532169770000022.GIF" wi="72" he="60" />下在其服务小区分得的带宽,N<sub>0</sub>为噪声功率。第四步:采集当前自身最优发送功率集合<img file="FDA00005321697700000220.GIF" wi="119" he="78" />和全局最优发送功率集合<img file="FDA0000532169770000023.GIF" wi="141" he="77" />根据公式(4)选择各候选发送功率集合n目前为止的自身最优功率集合,<img file="FDA0000532169770000024.GIF" wi="1180" he="113" />并根据公式(5)选择目前为止所有候选功率集合经历的全局最优发送功率集合。<img file="FDA0000532169770000025.GIF" wi="1127" he="116" />第五步:采集当前全局最优功率集合在N个候选功率集合中对应的标号g(t)。<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>arg</mi><mi>max</mi></mrow><mi>n</mi></munder><msubsup><mi>T</mi><mi>s</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532169770000026.GIF" wi="981" he="98" /></maths>第六步:更新除全局最优集合外,其他候选发送功率的调整尺度<img file="FDA0000532169770000027.GIF" wi="252" he="77" />和功率集合<img file="FDA0000532169770000028.GIF" wi="292" he="79" />根据公式(7)、(8)更新非当前全局最优候选发送功率集合<img file="FDA0000532169770000029.GIF" wi="254" he="83" />和功率调整尺度集合<img file="FDA00005321697700000210.GIF" wi="280" he="84" /><img file="FDA00005321697700000211.GIF" wi="1396" he="85" /><img file="FDA00005321697700000212.GIF" wi="1075" he="78" />其中r<sub>1</sub><sup>n</sup>(t)、<img file="FDA00005321697700000222.GIF" wi="103" he="64" />为[0,1]区间内的随机数,为保证寻优过程的收敛性,限制惯性权重ω及加速系数c<sub>1</sub>、c<sub>2</sub>的取值范围:<img file="FDA00005321697700000213.GIF" wi="946" he="236" />第七步:根据局部随机搜索项,更新当前全局最优功率集合的功率调整尺度<img file="FDA00005321697700000214.GIF" wi="219" he="85" />和发送功率集合<img file="FDA00005321697700000215.GIF" wi="253" he="83" /><img file="FDA00005321697700000216.GIF" wi="1092" he="82" /><img file="FDA00005321697700000217.GIF" wi="1190" he="77" />其中<img file="FDA00005321697700000218.GIF" wi="461" he="84" />为邻域局部随机搜索项,<img file="FDA00005321697700000219.GIF" wi="447" he="84" />第八步:用户根据更新的功率集合重新选择服务小区,计算各候选集合获得的吞吐量。用户根据更新了的功率集合及公式(1)重新选择服务小区,并根据公式(2)、(3)计算各候选功率集合对应的系统吞吐量。第九步:更新自身最优发送功率集合<img file="FDA0000532169770000031.GIF" wi="220" he="77" />全局最优发送功率集合<img file="FDA0000532169770000032.GIF" wi="174" he="77" />和全局最优功率集合标号g(t+1)。根据公式(4)、(5)、(6)更新当前各候选功率集合的自身最优发送功率、全局最优发送功率及全局最优功率集合标号。第十步:更新搜索到更优值的成功次数N<sub>s</sub>和失败次数N<sub>f</sub>。方法为:判断全局最优功率集合的标号是否发生改变,若g(t+1)=g(t)且T<sup>g(t+1)</sup>(t+1)>T<sup>g(t)</sup>(t),即全局最优标号不变且搜索到更优值,则N<sub>s</sub>=N<sub>s</sub>+1;若g(t+1)=g(t)且T<sup>g</sup>(<sup>t+1)</sup>(t+1)≤T<sup>g(t)</sup>(t),即全局最优标号不变且未搜索到更优值,则N<sub>f</sub>=N<sub>f</sub>+1;若g(t+1)≠g(t),即全局最优标号发生改变,则重置成功和失败次数,N<sub>s</sub>=N<sub>f</sub>=1。第十一步:更新局部搜索半径r<sub>Q</sub>(t+1)。根据(12)式更新局部搜索半径,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>r</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>&delta;</mi><msub><mi>r</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>if</mi><msub><mi>N</mi><mi>s</mi></msub><mo>/</mo><msub><mi>N</mi><mi>f</mi></msub><mo>></mo><msub><mi>&eta;</mi><mi>th</mi></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>/</mo><mi>&delta;</mi></mtd><mtd><mi>if</mi><msub><mi>N</mi><mi>s</mi></msub><mo>/</mo><msub><mi>N</mi><mi>f</mi></msub><mo>&lt;</mo><mn>1</mn><mo>/</mo><msub><mi>&eta;</mi><mi>th</mi></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd><mtd><mi>otherwise</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532169770000033.GIF" wi="1168" he="245" /></maths>也就是说,如果成功与失败比值超过一定门限,则增大在全局最优粒子附近搜索的搜索半径以尝试搜索更多可行功率集合;若果比值低于门限,说明在当前搜索半径内搜索到更优解的几率小,应减小全局最优粒子附近搜索的搜索半径以在收敛的小范围内详尽搜索更优功率集合;若成功失败比值等于门限,说明可以在当前搜索半径下继续搜索。第十二步:判断是否结束本次初始化的迭代。若t+1<t<sub>iter</sub>,且所有候选功率集合与全局最优功率集合间的距离之和大于门限值ε,则更新迭代次数,即t=t+1,并回到第六步,计算各更新的候选功率集合对应的系统吞吐量并更新自身及全局最优发送功率集合;否则,进行第十三步。第十三步:判断重新初始化结束条件。若s<t<sub>res</sub>,设置重新初始化次数s=s+1,t=0,重新初始化N‑1种候选发送功率集合,与当前全局最优功率集合一起作为新的N种候选发送功率集合,并重新初始化N种功率调整尺度,回到第二步,计算各新的候选集合对应的系统吞吐量并更新自身及全局最优发送功率集合;否则,进行第十四步。第十四步:停止,按照得到的全局最优发送功率集合设置各微微站的发送功率。
地址 210096 江苏省南京市四牌楼2号