发明名称 一种基于禁忌搜索的小基站容量与覆盖优化方法
摘要 本发明涉及一种基于禁忌搜索的容量与覆盖优化方法,属于无线通信技术领域,该方法包括:根据小基站初始功率,计算容量与覆盖目标优化函数初始值,采用禁忌搜索算法进行迭代优化,包括:计算当前解的邻域中使目标优化函数最优的解向量;若该解向量使目标优化函数值为目前为止最优,或非目前为止最优但处于非禁忌状态,则记为当前解并加入禁忌表,否则将该解向量从当前解的邻域中排除后返回计算邻域中最优的解向量;完成迭代优化后,更新小基站发射功率,输出此时小基站系统的容量与覆盖目标优化函数值。可实现小基站系统容量与覆盖的自优化,降低优化迭代的计算复杂度,实现网络整体性能的优化提升,且能达到全局最优。
申请公布号 CN105764068A 申请公布日期 2016.07.13
申请号 CN201610204356.6 申请日期 2016.04.01
申请人 清华大学 发明人 曾捷;粟欣;林小枫;朱晓鹏;肖驰洋;张黎;肖立民;许希斌
分类号 H04W16/18(2009.01)I;H04W52/24(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 廖元秋
主权项 一种基于禁忌搜索的容量与覆盖优化方法,其特征在于,该方法包括以下步骤:1)根据小基站初始发射功率,计算容量与覆盖目标优化函数的初始值,具体计算步骤如下:设该小基站系统包含K个宏基站l、M个小基站j和N个随机分布的用户终端UE<sub>i</sub>,其中,l=1,2,…,l,…,K,j=1,2,…j,…,M,i=1,2,…,i,…,N;小基站j的发射功率为p<sub>j</sub>,小基站j与用户终端UE<sub>i</sub>之间的传输信道增益为g<sub>ij</sub>,则从小基站j传输到用户终端UE<sub>i</sub>的参考信号接收功率P<sub>rx</sub>(i,j)用式(1)表达:P<sub>rx</sub>(i,j)=p<sub>j</sub>g<sub>ij</sub>  (1)设系统噪声为σ<sup>2</sup>,则小基站j经由下行链路至用户终端UE<sub>i</sub>的信干噪比SINR<sub>i</sub>表达式如式(2):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>SINR</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>p</mi><mi>j</mi></msub></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>g</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>p</mi><mi>k</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000957064120000011.GIF" wi="590" he="159" /></maths>其中g<sub>ik</sub>为小基站k与用户终端UE<sub>i</sub>之间的传输信道增益,p<sub>k</sub>为小基站k的发射功率,k≠j;假定每个导频功率都占据了总发射功率的均等固定的部分,并且每个专用信道也占据总功率的同等比例部分,则用户终端UE<sub>i</sub>的归一化吞吐量t<sub>i</sub>表示为式(3):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>t</mi><mi>i</mi></msub><mo>=</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>SINR</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>p</mi><mi>j</mi></msub></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>g</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>p</mi><mi>k</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000957064120000013.GIF" wi="1083" he="231" /></maths>则小基站j的归一化总吞吐量T<sub>j</sub>为表达式(4):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>j</mi></msub><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>UE</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>U</mi><mi>j</mi></msub></mrow></munder><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>p</mi><mi>j</mi></msub></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>g</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>p</mi><mi>k</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mo>(</mo><mn>4</mn><mo>)</mo><mo>,</mo></mrow>]]></math><img file="FDA0000957064120000012.GIF" wi="822" he="231" /></maths>其中,U<sub>j</sub>代表小基站j服务的所有用户终端,设容量的评估指标为整个网络中用户终端的平均吞吐量,其公式表达如式(5):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munder><mi>&Sigma;</mi><mrow><msub><mi>UE</mi><mi>i</mi></msub><mo>&Element;</mo><msub><mi>U</mi><mi>j</mi></msub></mrow></munder><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>p</mi><mi>j</mi></msub></mrow><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>+</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></munder><msub><mi>g</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>p</mi><mi>k</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000957064120000021.GIF" wi="1078" he="230" /></maths>以每个宏基站覆盖小区中用户吞吐量的累积分布的最低的p%来表示覆盖,用T<sub>l,p%</sub>表示宏基站l覆盖小区中用户吞吐量累积分布的最低的p%,引入折中系数γ来平衡覆盖和容量两个性能指标,0<γ<1,则M个小基站发射功率为p<sub>j</sub>'时,p<sub>j</sub>'={p<sub>1</sub>,p<sub>2</sub>,...,p<sub>j</sub>,...,p<sub>M</sub>},网络系统容量与覆盖的目标优化函数F(p<sub>j</sub>')定义如式(6):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>F</mi><mo>(</mo><mrow><msup><msub><mi>p</mi><mi>j</mi></msub><mo>&prime;</mo></msup></mrow><mo>)</mo><mo>=</mo><mi>&gamma;</mi><mfrac><mn>1</mn><mrow><mfrac><mi>K</mi><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>j</mi></msub></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&gamma;</mi><mo>)</mo></mrow><mfrac><mn>1</mn><mrow><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>T</mi><mrow><mi>l</mi><mo>,</mo><mi>p</mi><mi>%</mi></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000957064120000022.GIF" wi="878" he="207" /></maths>2)采用禁忌搜索算法生成新的小基站发射功率p<sub>j</sub>',并计算目标优化函数新的取值F(p<sub>j</sub>'),具体步骤如下:2‑a)初始化禁忌搜索的参数:设定禁忌表为空,设定禁忌长度(记为P,典型值为2,5,10等,为正整数)、邻域宽度、小基站功率的取值范围(即所有可取的功率值)、搜索的最大迭代次数m_iter;令m=1,m=1,2,...,m,...,m_iter,第m次迭代开始时,令l=1,设宏基站l覆盖小区中小基站的发射功率为X<sup>(m,l)</sup>,记为当前解;2‑b)计算在当前解X<sup>(m,l)</sup>的邻域范围内,使得目标优化函数F(p<sub>j</sub>')最优的解向量;并判断该解向量是否使得目标优化函数值F(p<sub>j</sub>')为目前为止最优,若最优则执行步骤2‑e),否则执行步骤2‑c);2‑c)判断该解向量是否处于禁忌状态,若是,则将该解向量从当前解的邻域中排除,并执行步骤2‑d),否则执行步骤2‑e);2‑d)判断当前解X<sup>(m,l)</sup>的邻域中是否尚有处于非禁忌状态的解向量,若尚有,则返回执行步骤2‑b);否则,将当前解的邻域中最早被标记为禁忌状态的解向量改置为非禁忌状态,并返回执行步骤2‑b);2‑e)将该解向量记为当前解,并将p<sub>j</sub>'中对应表示宏基站l覆盖小区中小基站的元素值更新为X<sup>(m,l)</sup>中对应值;更新禁忌表各对应元素的禁忌状态,即使得目标优化函数F(p<sub>j</sub>')最优的解向量在禁忌表的对应元素值置为P,其余已处于禁忌状态的解向量在禁忌表的对应元素值减一;判断是否遍历所有K个宏基站,若已遍历则执行步骤2‑g),未完成遍历则执行步骤2‑f);2‑f)令l=l+1,将宏基站覆盖小区内小基站的发射功率记为X<sup>(m,l)</sup>,同时记为当前解,返回执行2‑b);2‑g)判断m是否小于m_iter,若小于,则令m=m+1,并循环执行步骤2);否则结束循环,执行步骤3);3)完成m_iter次迭代优化后,小基站发射功率已更新为p<sub>j</sub>',输出此时小基站系统的容量与覆盖目标优化函数值F(p<sub>j</sub>');结束本次小基站系统容量与覆盖优化。
地址 100084 北京市海淀区清华园1号