发明名称 移动通信网络下设备间通信的最优收发联合处理方法
摘要 本发明公开了一种在移动通信网络下设备间通信的最优收发联合处理方法,按如下步骤进行:首先,时分双工通信系统中的基站通过上行信道估计可以等价地获得它至UE3的下行信道信息以及到D2D接收用户的下行干扰信道信息;UE1采用类似的信道估计方法可以获得它至UE2的传输信道信息和它对UE3的干扰信道信息,并将它们反馈给基站。接着,基站根据获得的所有四类信道信息按照最小化均方误差准则或者最大化信道信道容量准则计算出基站到UE3和UE1到UE2的预编码矩阵和接收矩阵。最后,基站将计算得到的发送预编码矩阵、接收矩阵通过控制信道通知给UE1、UE2和UE3。
申请公布号 CN102882570B 申请公布日期 2015.04.01
申请号 CN201210370443.0 申请日期 2012.09.28
申请人 东南大学;日电(中国)有限公司 发明人 许威;朱道华;李朝峰;雷鸣;赵春明
分类号 H04B7/04(2006.01)I;H04L25/02(2006.01)I 主分类号 H04B7/04(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种移动通信网络下设备间通信的最优收发联合处理方法,其特征在于:该方法按以下步骤进行: 1)D2D是指设备间,D2D链路是指D2D设备对UE1与UE2之间的通信链路,无线移动网络链路是指基站到UE3的下行通信链路,为简便计,采用链路1和链路2分别替代D2D链路和无线移动网络链路,所有相关参数的下标分别采用数字“1”、“2”来区分这两种类型的传输链路,D2D发送用户UE1和D2D接收用户UE2分别拥有N<sub>1</sub>根发送天线和M<sub>1</sub>根接收天线,基站拥有N<sub>2</sub>根发送天线,网络用户UE3拥有M<sub>2</sub>根接收天线,D2D发送用户和基站的发送功率分别为P<sub>1</sub>、P<sub>2</sub>,它们分别发送多流数据x<sub>1</sub>和x<sub>2</sub>给各自的接收机,它们的维度分别为S<sub>1</sub>×1和S<sub>2</sub>×1,S<sub>1</sub>和S<sub>2</sub>分别表示D2D链路和无线移动网络链路传输的数据流数; 2)链路1与链路2采用非正交方式共享时频资源,即D2D设备对之间通信的物理信道与某个无线移动网络用户的物理资源重叠;在TDD系统中,为了使基站获得完整的信道信息来实现最优收发设计,基站可以通过上行信道估计获得基站至UE3的下行信道H<sub>2</sub>以及基站至UE2的下行干扰信道H<sub>2,1</sub>;通过上述相同的信道估计方式,UE1同样可以获得它至UE2的传输信道H<sub>1</sub>和它与UE3之间的干扰信道H<sub>1,2</sub>,最后UE1将H<sub>1</sub>和H<sub>1,2</sub>反馈给基站; 3)基站根据获得的完整的四类信道信息通过最小化接收信号均方误差降低误比特率或者最大化信道信道容量准则优化基站、UE1的预编码矩阵和UE2、UE3的接收矩阵;所述最小化接收信号均方误差迭代算法称作算法一,该算法的迭代计算方法如下:步骤1:令迭代次数变量n=0,按照传输信道H<sub>1</sub>和H<sub>2</sub>奇异值分解结果初始化发送预编码矩阵<img file="FDA0000593819740000011.GIF" wi="260" he="70" />下标u是标识链路的变量,其中u=1和u=2分别表示D2D链路和移动网络下行链路;步骤2:n=n+1,然后根据下面式(2)计算第n+1次迭代过程中链路u的接收矩阵<img file="FDA0000593819740000012.GIF" wi="136" he="74" />步骤3:将步骤2计算得到的<img file="FDA0000593819740000013.GIF" wi="108" he="74" />代入功率约束条件下面式(3)求得参数<img file="FDA0000593819740000014.GIF" wi="144" he="77" />其中<img file="FDA0000593819740000015.GIF" wi="117" he="82" />是保证功率约束条件成立而引入的拉格朗日乘数;步骤4:将步骤2和步骤3中计算得到的<img file="FDA0000593819740000016.GIF" wi="115" he="74" />和<img file="FDA0000593819740000017.GIF" wi="114" he="75" />代入下面式(1)可以计算得到第n+1次迭代过程中的发送预编码矩阵<img file="FDA0000593819740000018.GIF" wi="160" he="77" />步骤5:重复步骤2到步骤4直到W<sub>u</sub>(n+<sup>1</sup>)和F<sub>u</sub>(n+<sup>1</sup>)收敛; <img file="FDA0000593819740000021.GIF" wi="1655" he="125" /><img file="FDA0000593819740000022.GIF" wi="1583" he="123" /><img file="FDA0000593819740000023.GIF" wi="1216" he="239" />式(1)至式(3)中下标u和m都是标识链路的变量且u,m∈{1,2},u≠m,而n表示迭代次数变量,上标H表示矩阵的共轭转置操作,I表示单位阵;<img file="FDA0000593819740000024.GIF" wi="160" he="98" />和<img file="FDA0000593819740000025.GIF" wi="144" he="98" />分别代表了矩阵<img file="FDA0000593819740000026.GIF" wi="140" he="85" />和<img file="FDA0000593819740000027.GIF" wi="141" he="82" />主对角线上的第t个元素;其中,<img file="FDA0000593819740000028.GIF" wi="941" he="82" /><img file="FDA0000593819740000029.GIF" wi="136" he="81" />和<img file="FDA00005938197400000210.GIF" wi="136" he="84" />分别是<img file="FDA00005938197400000211.GIF" wi="525" he="148" />特征分解的对角阵和酉矩阵;所述最大化信道容量算法称作算法二,该算法的迭代计算方法如下: 步骤1:令迭代次数变量n=0,选择足够小的正数τ以及门限,按照传输信道H<sub>1</sub>和H<sub>2</sub>奇异值分解结果初始化发送预编码矩阵<img file="FDA00005938197400000212.GIF" wi="88" he="69" />和接收矩阵<img file="FDA00005938197400000213.GIF" wi="98" he="68" />并将它们代入下面式(4)得到初始化信道容量并记为cap(0);步骤2:n=n+1,然后根据下面式(5)二次规划计算出第n+1次迭代中更新向量z;z的表达式为 <img file="FDA00005938197400000214.GIF" wi="1313" he="242" />其中,<img file="FDA00005938197400000215.GIF" wi="241" he="92" />分别表示取复数的实部和虚部,<img file="FDA00005938197400000216.GIF" wi="114" he="82" />和<img file="FDA00005938197400000217.GIF" wi="110" he="84" />分别表示链路u预编码矩阵和接收矩阵在地n+1次迭代过程中的更新值;u=1、2分别表示D2D链路和移动网络下行链路;vec(·)表示拉直运算,T表示矩阵转置;步骤3:利用步骤2中获得的更新向量并引入线性搜索参数β,β在区间[0,1]中线性变化,同时记<img file="FDA00005938197400000218.GIF" wi="472" he="88" />和<img file="FDA00005938197400000219.GIF" wi="428" he="92" />并将它们代入下面式(4)计算 得到含参的cap(n+1,β),最终求得最大化cap(n+1,β)的β<sup>*</sup>;步骤4:更新<img file="FDA0000593819740000031.GIF" wi="440" he="85" />和<img file="FDA0000593819740000032.GIF" wi="432" he="86" />并得到在n+1迭代过程中最优的信道容量记为cap(n+1)=cap(n+1,β<sup>*</sup>);步骤5:重复步骤2到步骤4直到cap(n+1)‑cap(n)≤门限; 步骤6:最后,根据功率约束条件对W<sub>u</sub>线性加权; 最大化信道容量准则的计算公式如下面式(4)所示,下面式(5)是下面式(4)的近似式且为二次规划问题; <img file="FDA0000593819740000033.GIF" wi="1878" he="380" />其中,式(4)中下标u和m都是标识链路的变量且u,m∈{1,2},u≠m,上标H表示矩阵的共轭转置操作,I表示单位阵;w<sub>u,k</sub>和f<sub>u,k</sub>分别是发送预编码矩阵W<sub>u</sub>和接收矩阵F<sub>u</sub>的第k列向量;ρ<sub>u</sub>是链路u的发送功率约束因子,保证功率恒为P<sub>u</sub>;<img file="FDA0000593819740000034.GIF" wi="77" he="79" />为链路u的噪声功率;<img file="FDA0000593819740000035.GIF" wi="107" he="92" />为链路u的第k个数据流的接收信干噪比,将它泰勒展开后用一阶多项式近似忽略高阶项后可以将(4)式的最优化问题等价的转化为如下的凸优化问题:<img file="FDA0000593819740000036.GIF" wi="1574" he="294" />式(5)中τ1代表了一个足够小的正数τ乘以向量1,1向量中元素全为1;式(5)已经将某些项合并为如下所列中间变量: <img file="FDA0000593819740000037.GIF" wi="1257" he="413" />p<sub>u</sub><sub>,</sub><sub>k</sub>表达式中的<img file="FDA0000593819740000041.GIF" wi="255" he="94" />g<sub>x,u,k</sub>和g<sub>y,u,k</sub>如下所示,<img file="FDA0000593819740000042.GIF" wi="556" he="78" /><img file="FDA0000593819740000043.GIF" wi="1441" he="142" /><img file="FDA0000593819740000044.GIF" wi="1732" he="884" />式(8)中0代表零向量,下标表示向量维度且完全由收、发天线数M<sub>u</sub>、N<sub>u</sub>,数据流个数S<sub>u</sub>决定;式(7)中<img file="FDA0000593819740000045.GIF" wi="82" he="94" />和<img file="FDA0000593819740000046.GIF" wi="88" he="91" />分别为<img file="FDA0000593819740000047.GIF" wi="85" he="93" />的分子和分母,变量a<sub>u,i</sub>,b<sub>u,i</sub>,c<sub>m</sub>和d<sub>u,k</sub>定义如下<img file="FDA0000593819740000048.GIF" wi="1488" he="505" />式(6)中矩阵Q中的行向量为 <img file="FDA0000593819740000049.GIF" wi="1640" he="253" />式(10)中0代表零向量,下标表示向量维度且完全由收、发天线数M<sub>u</sub>、N<sub>u</sub>,数据流个数S<sub>u</sub>决定; 4)基站根据最小化接收信号的均方误差或者最大化信道容量准则计算得到链路1、链路2的发送预编码矩阵和接收矩阵分别为W<sub>1</sub>、F<sub>1</sub>和W<sub>2</sub>、F<sub>2</sub>,基站将W<sub>1</sub>、F<sub>1</sub>通过控制信道传递给D2D通信设备对; 5)最后,UE1利用发送预编码矩阵W<sub>1</sub>加权数据流x<sub>1</sub>发送给UE2,而发送预编码矩阵W<sub>2</sub>则被基站用来加权下行信号x<sub>2</sub>并发送给UE3。 
地址 210096 江苏省南京市四牌楼2号