发明名称 一种基于带内自回传的干扰抑制方法
摘要 本发明公开了一种基于带内自回传的干扰抑制方法,特征是在采用带内自回传技术的高密度小基站系统中,通过对自回传资源分配与接入链路预编码的联合设计,实现在接入链路以用户为中心构建小基站协作簇,簇内进行部分协作抑制小基站间的同频干扰,自适应的根据自回传链路速率调整协作簇的大小,在联合设计过程中考虑小基站的回传负载,减小不同小基站自回传资源分配的差异。采用本发明提出的方法能够均衡小基站的回传负载,降低协作开销,提高系统的下行接入链路吞吐量。
申请公布号 CN105072647A 申请公布日期 2015.11.18
申请号 CN201510470746.3 申请日期 2015.07.31
申请人 中国科学技术大学 发明人 张静雷;邱玲
分类号 H04W28/08(2009.01)I;H04B7/04(2006.01)I 主分类号 H04W28/08(2009.01)I
代理机构 安徽省合肥新安专利代理有限责任公司 34101 代理人 汪祥虬
主权项 一种基于带内自回传的干扰抑制方法,在采用带内自回传技术的高密度小基站系统(S8)中,系统包含一个宏基站(S1)和多个无线带内自回传小基站(S2)以及若干用户(S3),所述带内自回传小基站的无线传输分为接入链路(S5)与自回传链路(S4)两个部分;宏基站作为一个中心处理单元通过专用的光纤回传链路(S6)从核心网网关(S7)获取用户的数据并通过自回传链路下发给小基站,小基站通过接入链路将数据发送给用户;其特征在于:宏基站收集其到所有小基站的自回传链路信道以及小基站到用户的接入链路信道的信道状态信息,利用虚拟完全协作法,确定自回传资源分配因子的初始搜索区间上界,进而利用区间消去法确定三个自回传资源分配因子搜索值;在确定的资源分配因子搜索值条件下,基于分块坐标下降法,根据接收矢量更新公式、辅助变量更新公式以及利用优化工具求解的考虑小基站回传负载的组稀疏优化问题,通过交替迭代得到组稀疏化的预编码矩阵;利用区间消去法确定最优的自回传资源分配因子和相应的组稀疏化预编码矩阵;宏基站根据上述自回传资源分配与接入链路预编码矩阵的联合设计结果,设置自回传资源分配,并通过自回传链路下发用户数据和组稀疏化的预编码矩阵给小基站;在接入链路,小基站根据稀疏化的预编码构建以用户为中心可控大小的小基站协作簇进行部分协作;具体操作步骤如下:第一步,宏基站收集其到所有小基站的回传信道以及小基站到用户的接入信道的信道状态信息,利用虚拟完全协作法,将小基站全部无线时频资源用于接入链路资源传输,所有小基站为系统内的所用用户进行完全协作,利用考虑单基站功率限制的完全协作方法求得对应最大的协作增益的系统所有用户接入链路和速率R<sup>Full‑JP</sup>,利用资源分配比例因子搜索区间上界公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>&beta;</mi><mo>^</mo></mover><mo>=</mo><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mfrac><msup><mi>R</mi><mrow><mi>F</mi><mi>u</mi><mi>l</mi><mi>l</mi><mo>-</mo><mi>J</mi><mi>P</mi></mrow></msup><mrow><msubsup><mi>R</mi><mi>l</mi><mrow><mi>B</mi><mi>L</mi></mrow></msubsup><mo>+</mo><msup><mi>R</mi><mrow><mi>F</mi><mi>u</mi><mi>l</mi><mi>l</mi><mo>-</mo><mi>J</mi><mi>P</mi></mrow></msup></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000772636260000011.GIF" wi="1255" he="159" /></maths>确定自回传资源分配因子β的初始搜索区间的上界<img file="FDA0000772636260000012.GIF" wi="70" he="82" />自回传资源分配因子β指无线自回传小基站自回传链路占用的时频资源与小基站全部时频资源的比例,<img file="FDA0000772636260000013.GIF" wi="80" he="77" />是通过小基站反馈得到的小站l的回传链路速率;利用区间消去法,设置搜索区间的下界a=0,搜索区间的上界<img file="FDA0000772636260000021.GIF" wi="163" he="81" />区间长度为<img file="FDA0000772636260000022.GIF" wi="321" he="81" />根据区间消去法初始搜索值设置公式:β<sub>1</sub>=a+L/4β<sub>2</sub>=b‑L/4                (2)β<sub>m</sub>=(a+b)/2设置资源分配因子的三个搜索值β<sub>1</sub>,β<sub>2</sub>,β<sub>m</sub>;初始化相关参数,随机生成满足小基站的单基站功率约束的初始化预编码矩阵;利用用户速率计算公式<img file="FDA0000772636260000023.GIF" wi="1654" he="198" />求得用户速率R<sub>k</sub>,将其作为迭代算法中用户速率的初始化迭代值<img file="FDA0000772636260000024.GIF" wi="85" he="83" />其中小基站的发射天线数M,用户接受天线数为N,<img file="FDA0000772636260000025.GIF" wi="210" he="79" />表示基站l到用户k的预编码矢量,<img file="FDA0000772636260000026.GIF" wi="484" he="95" />为所有小基站到用户k的预编码矢量,<img file="FDA0000772636260000027.GIF" wi="351" he="82" />为所有小基站用户的集合,<img file="FDA0000772636260000028.GIF" wi="296" he="82" />为系统内所有小基站的集合;<img file="FDA0000772636260000029.GIF" wi="255" he="78" />为所有小基站到用户k的信道,通过小基站向宏基站的反馈获得,<img file="FDA00007726362600000210.GIF" wi="66" he="79" />为将用户k的噪声建模为加性高斯白噪声的均方差;初始化预编码l<sub>2</sub>范数近似的权值<img file="FDA00007726362600000211.GIF" wi="156" he="79" />初始化辅助变量值w<sub>k</sub>=1;第二步,在自回传资源分配比例因子β分别取β<sub>1</sub>,β<sub>2</sub>,β<sub>m</sub>的条件下,利用干扰信道加权和速率与最小化和加权均方误差的等价关系,采用分块坐标下降法,分别更新用户接收矢量,辅助变量以及预编码矢量:利用最优的用户接收矢量更新公式<img file="FDA00007726362600000212.GIF" wi="1387" he="171" />,更新每个用户的线性接受矢量;利用均方差计算公式<img file="FDA00007726362600000213.GIF" wi="1705" he="133" />计算每个用户的均方差,其中e<sub>k</sub>表示用户k的均方差;将当前辅助变量值赋值保存为w′<sub>k</sub>,利用辅助变量更新公式w<sub>k</sub>=e<sub>k</sub><sup>‑1</sup>                 (6)更新辅助变量w<sub>k</sub>;利用现有的凸优化工具,根据考虑小基站回传负载的组稀疏优化问题公式,<img file="FDA0000772636260000031.GIF" wi="1614" he="436" />更新用户的预编码矢量,其中<img file="FDA0000772636260000032.GIF" wi="161" he="103" />为基站l到用户k的预编码适量<img file="FDA0000772636260000033.GIF" wi="62" he="79" />的l<sub>0</sub>范数<img file="FDA0000772636260000034.GIF" wi="114" he="97" />的凸近似,而<img file="FDA0000772636260000035.GIF" wi="52" he="79" />为凸近似权值;基于组稀疏优化,根据组稀疏预编码范数权值更新公式,更新预编码l<sub>2</sub>范数近似的权值;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mi>k</mi><mi>l</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mo>|</mo><mo>|</mo><msubsup><mi>v</mi><mi>k</mi><mi>l</mi></msubsup><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000772636260000036.GIF" wi="1159" he="165" /></maths>在上述三个更新过程中交替迭代,直到满足迭代终止条件公式<img file="FDA0000772636260000037.GIF" wi="1317" he="156" />其中γ为迭代终止条件门限;利用用户速率计算公式求得用户速率R<sub>k</sub>,将其作为迭代算法中用户速率的初始化迭代值<img file="FDA0000772636260000038.GIF" wi="91" he="83" />计算系统的下行接入链路和速率<img file="FDA0000772636260000039.GIF" wi="303" he="117" />根据预编码组的功率<img file="FDA00007726362600000310.GIF" wi="139" he="105" />将其小于指定门限的预编码组中的元素置为零得到稀疏化的预编码矩阵;第三步,利用区间消去法确定最优的自回传资源分配因子和相应的稀疏化的预编码矩阵;利用区间消去法,更新自回传资源分配比例因子搜索值并计算搜索区间:令系统下行接入链路和速率关于资源分配因子β函数<img file="FDA00007726362600000311.GIF" wi="459" he="117" />比较资源分配因子三个搜索值β<sub>1</sub>,β<sub>2</sub>,β<sub>m</sub>分别对应的系统下行接入链路和速率f(β<sub>1</sub>),f(β<sub>2</sub>),f(β<sub>m</sub>),首先比较资源分配因子搜索值β<sub>1</sub>,β<sub>m</sub>对应的系统下行接入链路和速率f(β<sub>1</sub>),f(β<sub>m</sub>):若f(β<sub>1</sub>)&lt;f(β<sub>m</sub>),令b=β<sub>m</sub>,β<sub>m</sub>=β<sub>1</sub>;否则,比较资源分配因子搜索值β<sub>2</sub>,β<sub>m</sub>对应的系统下行接入链路和速率f(β<sub>2</sub>),f(β<sub>m</sub>):若f(β<sub>2</sub>)&lt;f(β<sub>m</sub>),令a=β<sub>m</sub>,β<sub>m</sub>=β<sub>2</sub>;否则,令a=β<sub>1</sub>,b=β<sub>2</sub>;更新区间长度L=b‑a,根据设定终止的区间值,若更新后的区间长度L不满足目标终止的区间值要求,则根据更新后的资源分配因子搜索值转至第二步;否则得到最优的自回传资源分配因子和相应的预编码矩阵,输出对应最大的接入链路和速率<img file="FDA0000772636260000041.GIF" wi="295" he="117" />的自回传资源分配因子β<sub>m</sub>和相应的稀疏化的预编码矩阵;第四步,宏基站根据上述自回传资源分配与接入链路预编码矩阵的联合设计结果设置自回传资源分配,下发用户数据和预编码矩阵给小基站;对于预编码矩阵中<img file="FDA0000772636260000042.GIF" wi="60" he="78" />被置为零的组,用户k的协作基站簇中将不包含基站l,宏基站不下发用户k的数据给小基站l;进而所有收到用户k的数据的小基站以相干联合处理的方式向用户k联合传输数据;小基站通过自回传链路获得宏基站下发的用户数据和相应的预编码矩阵;根据稀疏化的预编码矩阵,小基站构建以用户为中心的可控大小的协作簇进行部分协作。
地址 230026 安徽省合肥市包河区金寨路96号