一种mimo全双工中继系统信源中继联合预编码方法
【专利摘要】本发明公开了一种MIMO全双工中继系统信源中继联合预编码方法，包括以下步骤：首先利用导频信道估计方法得到信源?中继、中继?终端、以及中继发射端?中继接收端的信道矩阵，设定系统功率上限及背景噪声值；然后设定中继放大矩阵其中xt,xr分别为发射波束成形和接收波束成形；其次将非凸的多变量系统速率最大化问题转化为仅关于xr的无约束问题，利用梯度上升法求解得到xr，再根据xr计算得到信源和中继的预编码矩阵；最后系统使用计算得到的源端预编码矩阵和中继放大矩阵产生源端和中继处的发送信号，以完成MIMO全双工中继系统的信息传输。本发明通过利用中继放大矩阵和梯度上升法设计MIMO全双工中继系统的预编码方案，从而极大地提升MIMO全双工中继系统的传输速率。
【专利说明】
一种ΜI MO全双工中继系统信源中继联合预编码方法
技术领域
[00011本发明属于数字移动通信系统与技术领域，具体为一种MIMO(Multiple-Input Mu 11 ip 1 e-output)全双工中继系统信源中继联合预编码方法设计。
【背景技术】
[0002] 通信技术的迅猛发展使人们对数据传输速率和通信质量的要求越来越高。在新型 的mMO无线中继通信系统中，一方面，中继技术可以提高系统容量和增加小区覆盖范围因 而成为下一代移动通信的关键技术之一;而另一方面，传统中继都是工作在半双工模式下， 这意味着一个时隙被分成两个甚至更多时隙用于中继的收发而导致系统的频谱效率将会 降低至少50%。而当下，中继的同频全双工模式具有显著的理论研究及应用价值，该模式也 即在同时工作的收发信道上使用相同的频率资源。
[0003] 然而，"Kim D，Lee Η，Hong D · A Surνey 〇f In-band Fu 11-dup 1 eX Transmission : From the Perspective of PHY and MAC Layers[J]. IEEE Communications Surveys&Tutorials，2015,17(4): 1-1 ·"提出由于在从中继输出到接收端 输入过程中信号泄露，全双工操作的主要限制是回环干扰(也称为回环自干扰）。具体而言， 全双工操作的主要缺点是全双工终端生成的回环干扰和来自远处的特定的接收信号之间 存在明显的功率差。大回环路线跨越大部分在接收端的模数转换器的动态范围，因此抑制 该部分的干扰对于全双工中继的操作执行非常关键。
[0004] 在现有的研究中，大部分工作均没有考虑到信息传送在中继的时延误差，因此，在 考虑中继时延、多天线全双工中继系统特有的中继回环干扰约束的条件下，本发明首先利 用在中继的预编码放大矩阵消除回环干扰，再利用中继放大矩阵Q = x;xf和梯度上升法开 展ΜΙΜ0无线中继系统的高效传输机制研究，实现在信源功率、中继功率等多约束下的传输 速率最大化问题。
【发明内容】
[0005] 本发明的目的在于针对现有技术的不足，提供一种ΜΜ0全双工中继系统信源中继 联合预编码方法，包括以下步骤：
[0000] (1)初始化:系统确定源节点的天线数Ns、中继节点接收端的天线数Nr、中继节点发 射端的天线数Ντ、信宿节点的天线数ND、信源的功率上限Ps、中继的功率上限PR;设定中继处 噪声功率σ〗，信宿处噪声功率g ;
[0007] (2)利用导频方法进行信道估计得到信源-中继、中继发射端-中继接收端、中继- 信宿的信道矩阵，分别记为HsR、HRR、HRD ;设定中继放大矩阵.Q = Xfxf，其中X t、Xr分别为一维 发射波束成形向量和一维接收波束成形向量，将ΜΙΜ0全双工中继系统的速率最大化问题等 价为仅关于的无约束优化问题：
[0008]
[0009] 其中Amax(Xr)为矩阵?μΠΗ盂的最大特征值，
[0010] (3)利用梯度上升法得到步骤2所述无约束问题的最优解Xr，并按如下公式计算信 源预编码矩阵V、中继发送波束成形 Xt，中继接收波束成形xr:
[0011]
[0012] ⑷在信源端根据V对传输符号进行预编码，同时根据Xt，Xr计算得到中继放大矩阵 Q = x,:xf，利用Q对中继接收端的接收信号进行放大处理，完成ΜΜ0全双工中继系统的预编 码设计，最终实现信源与信宿之间以全双工的方式进行信息传输。
[0013] 进一步地，所述步骤3中梯度上升方法，具体包括以下步骤：
[0014] (3.1)初始化迭代次数t = 0、迭代次数上限tmax，初始化Xr为矩阵最大特征 值所对应的特征向量；
[0015] (3.2)更新迭代次数:t = t+Ι;令山为Hm ΠHlo最大特征值Amax所对应特征向量， 由下式计算得到矩阵Η? Π 最大特征值的梯度，即：
[0016]
[0017] 继而根据计算得到步骤2所述无约束问题目标函数关于xr的梯度grad_ f;
[0018] (3.3)判断是否满足grad_f彡ε或t彡tmax，其中ε表示判定阈值，其值在〇 . 〇〇 1~ 0.00000 1之间，若满足，则输出最终解Xr，否则根据线搜索方法确定步长β，更新Xr = Xr+β X grad_f，并重复步骤3 · 2~3 · 3 〇
[0019]本发明有益效果:本发明方法首先设定中继放大矩阵Q = X,xf，利用该性质将原非 凸复杂模型简化为单变量无约束问题，从而大大简化设计复杂度，再基于利用梯度上升法 最终设计出一种ΜΜ0全双工信源中继联合预编码方法，该方法可以在达到系统功率要求的 同时提升ΜΙΜΟ全双工中继系统的传输速率。
【附图说明】
[0020]图1是本发明实施例采用该方法的系统模型图。
[0021 ]图2是本发明实施例采用该方法的具体流程图。
[0022]图3是本发明实施例的系统平均传输速率与信噪比的关系图。
[0023]图4是本发明实施例的系统平均传输速率与各节点天线数量的关系图。
【具体实施方式】
[0024]为了使本发明的目的和效果更加清楚，下面对ΜΜ0全双工中继系统及本发明方法 进行详细描述。
[0025]如图1所示，本发明考虑一个三节点的ΜΜ0全双工中继网络，其中信源S将在一个 全双工中继R的帮助下传输信息到终端D。在这个传输网络中，信源和信宿分别配置Ns和Nd根 天线，同时中继的接收端和发送端分别有Nr和Ντ根天线。假设所有的信道都遵循独立瑞利衰 落，令Hsr表示信源和中继之间的信道，令Hrd表示中继和信宿终端之间的信道。值得注意的 是本模型考虑回环自干扰(SI)的影响，因此令Hrr表示从中继发送端到中继自身接收端的自 干扰信道。
[0026] 同时，在中继实施全双工操作需要一定的执行时间，这就造成中继存在一定的时 延，在这里假设这个时延为周期τ。但是，时延τ远远小于数据传输时隙，因此，在中继的时延 对可实现的速率的影响是忽略不计的。本发明将在信源和中继采取线性处理来增强系统传 输性能，在信源端利用预编码矩阵V e Cv〃"处理后发送信息到中继接收端，然后中继接收 端接收信息再经过放大矩阵QeC^x~的处理。因此，可以得到在时刻n，中继接收端的信息 f[H]e C、"1 为：
[0027] γ [n] =HRRVs[n]+HRRXR[n]+nR[n] (1)
[0028] 其中，s[n]~CN(0，Id)是由d个传输符号组成的向量，n,?~表示加性 高斯白噪声(AWGN)，Hrrxr [η ]表示中继输出端到中继输入端的回环自干扰。在中继发送端的 发送信号xR[n]可以表示为：
[0029] XR[n] =Q γ [η-τ] (2)
[0030] 结合式(1)和式(2)，中继输出信号可以表不为：
[0031]
(3)
[0032]在式(3)中的QHRRQY[n-2T]是关于Q的复杂函数，同时令系统设计变得非常复杂。 为了简化设计，本发明利用迫零法使从中继输出端到中继输入端的回环自干扰(SI)为零， 也即：
[0033] QHrrQ = 0 (4)
[0034] 将(4)代入(3)可以得到：
[0035] XR[n] =QHRRVs[n-T]+QnR[n-T] (5)
[0036] 最后，终端D所收到的信号即可表示为：
[0037]
[0038] 其中，11"^]~€7\\〇力_〗>丨)为复性加性高斯白噪声。根据式(6)，系统速率可以表示 为：
[0039]
(7；
[0040] 另外，在中继的功率消耗可以表示为：
[0041]
(8)
[0042]同时,Tr (VVH)表示信源的功率消耗。
[0043] 本发明的目标是通过信源-中继联合预编码设计来优化在信源/中继功率约束以 及回环自干扰(SI)约束下的系统速率。该速率最大化问题用数学语言来表示即：
[0044]
(9) \ τ ' I _ ?
[0045] 其中，问题(9)中目标函数表示整个系统的传输速率，PR、Ps分别表示为支持整个系 统正常运营所需的最小中继消耗功率和最小信源消耗功率。
[0046] 问题(9)表现为非凸问题，在数学上很难求解。针对该复杂问题，本发明设计了一 种在中继放大矩阵秩为1情况下的基于梯度上升算法的信源中继联合预编码方法，实现了 在达到各节点消耗功率要求的同时提高系统信息传输速率的目的。该方法核心思想为:利 用中继放大矩阵= 将原非凸的系统传输速率最大化问题转化一个易解无约束问题， 再利用梯度上升法迭代求解近似问题，最终得到在信源和中继的预编码矩阵。根据该方法 的思想，具体实施如下：
[0047] 本发明是在中继端放大矩阵秩为1的情况下提出的，当放大矩阵Q秩为1时，具有 的结构，其中，xt，Xr都是一维向量，本发明将利用这个特殊关系来解决原速率最大 化问题(9)。
[0048] 当Q = X,xf时，可以得到：
[0049]
(13)
[0050] 利用以上两个关系，R(V，Q)又可以转化为：
[00511
[0052] 其中，第二个式子运用了（Ι+ΑΒΓ?ζΑα+ΒΑΓ1这个性质。同样的，根据Tr(AB) = Tr(BA)，pR(V，Q)又可以表示为
[0053]
(15)
[0054] 另外，当QHRRQ = 0时，可以得到xf = 0。因此，根据以上转化，就可得到一个新 的目标函数和相应的约束条件，原问题也可以等价为：
[0055]
[0056] 虽然问题（15)比问题(9)简单了很多，但因为回环自干扰和功率的限制，这个问题 依然很难解决。但对于这个放大矩阵Q秩为1的特殊结构，本发明继续将问题简化为一个只 关于xr的无约束问题。
[0057] 本发明首先定夕
'同时定义Amax(Xr)是矩阵H fi/,丨| 的最大 特征值。
[0058] 1)问题（16)可表示成以下无约束性问题
[0059]
[0060] 2)本发明通过式（17)得到xr的最优解之后，就可以得到相应的如以下 所示：
[0061]
[0062] 具体过程如下：
[0063]第一步:可以发现，当给定xr时，V的最优解应该在使最大化的同时受到 功率的约束。这两者是相互矛盾的存在。但可以反过来假设问题（16)取到最优解时，关于预 编码矩阵v的||xfHffiv| 2并没有取到最大值。因此，可以通过选择合适的放大矩阵V，来使 |xfHffiV『增加。同时，可以令^减小一些，使（16)中目标函数的因式『及 中继功率约束中因式|xfHwV|2|x,f均保持一个常数不变。这也就意味着可以取一个合适 的(v，Xt)来增加目标IxfHuVl2的值，虽然此时V和以并不是相应的最优取值，但目标函数 因此可以得到关于v的最优解。同时，可以得到，IxfHiVl2的最优解等于e |χ=η5Λ|2，所以， 问题（16)可以重新写成以下优化问题：
[0064] (19)
[0065] 第二步:在这个步骤中，本发明将利用Xt与Xr之间存在的关系消除回环自干扰的约 束。首先，假设(xt，xr)是最优取值，注意到对于任意实数α，都有(《χ,,?χ,.)是问题（19)的最 a 优取值。
[0066] 因此，若是不考虑对Xt、Xr的最优取值，本发明可以假设令| |xt| |=1(该假设并不 影响目标函数的最优解），所以，问题（19)可以重新写为：
[0067]
[0068] 另一方面，从第二个约束条件看到xt存在于Η^χ,.的零空间，根据前文对
&的定义，在xfH#x, = 0约束下的每一个以都可以表示为Xt=IIu，其中 u为任意向量。将Xt=IIu代入(20)，又可以得到一个新的表达式：
[0069] (21)
[0070] 除此之外，可以注意到现在的目标函数随着| |HRDnu| |2的增加而增加，所以，当关 于U的因式| IHrdIIuI |2取得最大值时，就可以得到目标函数的最大值。根据以上分析，首先 假设為^代表1 |HRDnu| I2的最优解，有
[0071]
(22)
[0072] 其中，本发明用到了 Π 这样一个性质：Π = Π 2。根据拉格朗日乘子法可以得到，问 题(22)中，U为矩阵ΠΗ^Ημ Π 的最大特征值，也即?，利用这个性质，本 发明可以将原优化问题进一步简化
[0073]
(23)
[0074]第三步：在这一步骤中，本发明将（23)够转化为一个无约束问题。从上式可以发 现，目标函数是随着11 Xr 11的增加而增加的，因此，当原不等式约束的取到边界值时，相应 的目标函数也能取到其最优解。可以得到：
[0075]
(24)
[0076] 通过对原系统模型的简化后，得到了一个简单易解的数学问题，在本发明中，将利 用梯度上升法来求解相应的目标函数。
[0077] 首先，Amax(Xr)是矩阵Η^ΠΗ:的最大特征值。而对于随机产生的信道矩阵Hrd、 HRR/H^niC^非零特征值λ(Χ?)总是不同的，因此，矩关于Xr的最大特征值 也不同，所以，本发明令m作为该矩阵最大特征值所对应的特征向量，可以得到A max(Xr)的梯 度表达式：
[0078]
发射波束成形向量和一维接收波束成形向量，将ΜΙΜΟ全双工中继系统的速率最大化问题等 价为仅关于的无约束优化问题：
[0083]
%
[0084] 其中λ腹(Xr)为矩阵H Η;的最大特征值，
[0085] (3)利用梯度上升法得到步骤2所述无约束问题的最优解Xr，并按如下公式计算信 源
[0086] 波束成形V、中继发送波束成形Xt，中继接收波束成形xr:
[0087]
[0088] (4)在信源端根据V对传输符号进行预编码，同时根据Xt，Xr计算得到中继放大矩阵 Q = ，利用Q对中继接收端的接收信号进行放大处理，完成MBTO全双工中继系统的预编 码设计，最终实现信源与信宿之间以全双工的方式进行信息传输。
[0089 ]进一步地，所述步骤3中梯度上升方法，具体包括以下步骤：
[0090] (3.1)初始化迭代次数t = 0,迭代次数上限tmax，初始化Xr为矩阵Η_Η孟最大特征 值所对应的特征向量；
[0091] (3.2)更新迭代次数4 =丨+1;令111为1：^)111^0最大特征值\_,所对应特征向量， 由下式计算得到矩阵Π Η&最大特征值的梯度^^^^:^^即：
[0092]
[0093] 继而根据V4nax )计算得到步骤2所述无约束问题目标函数关于 Xr的梯度grad_ f;
[0094] (3.3)判断是否满足grad_f彡ε或t彡tmax，其中ε表示判定阈值，其值在〇 . 〇〇 1~ 0.00000 1之间，若满足，则输出最终解Xr，否则根据线搜索方法确定步长β，更新Xr = Xr+β X grad_f，并重复步骤3 · 2~3 · 3 〇
[0095] 符号说明：对于方阵4乂，，人^'1>(4)，以1^(4)分别表示对矩阵厶做转置 (transpose)运算、共辄车专置（conjugate transpose)运算、伪逆（seudoinverse)运算、逆 (inverse)运算、求迹（trace)运算、求秩（rank)运算，另外det(A)为求矩阵A的行列式 (determinant)运算等。对于复向量X，复矩阵X，| |x| |和| |X| |分别表示相应的Euclidean范 数和frobenius范数;CMXN表示Μ X N的复矩阵集合;~CN( X，y)表示为循环对称复高斯随机变 量且满足均值为X，协方差为y。
[0 09 6 ]图3 - 4是本发明通过M a 11 a b对所设计方案的仿真验证。由于传统的半双工中继的 性會κ已经在文南犬"Suraweera H A,Krikidis I,Zheng G,et al ·Low-Complexity End-to- End Performance Optimization in ΜΙΜΟ Full-Duplex Relay Systems[J]. IEEE Transactions on Wireless Communications,2014，13(2) :913-927." 中被讨论过了，本发 明在这里只分析在不同的预编码方法下全双工中继系统的性能。参数具体设置为:信源和 信宿的天线数为Ns = Nd = Nsd，中继接收端和发送端的天线数为Nt = Nr=Ntr，设定信息传输中 的噪声功率为σ!=砣同时，令Ps = PR=P，&V/^ lObgd.除非另有规定，一般情况 CT 下，设定P=10dB，〇2 = 0dB。进一步假设信源-中继和中继信宿信道经历独立瑞利平坦衰落。 另外，在回环信道HRR中每个元素定义为高斯分布的随机复变量，其平均值为0,方差为-20dB。请注意，所有的仿真结果为进行100个独立信道实现后的平均值。
[0097]图3在NSD = NTR = 2的条件下，比较了在多种方法下的随SNR增加而变化的系统传输 速率。当利用中继放大矩阵Q = x,Xf将原多约束非凸问题简化为可以用梯度上升法解决的 无约束问题后，可以从图中观察到，此时梯度上升法比全局搜索法和发射/接收迫零法能实 现更优的传输性能。
[0098] 图4展示了平均传输速率与在信源中继的传输天线的数量N之间的关系，此时Nsd = Ntr = N。同样的，也可以从这幅图中看出梯度上升法性能优于发射/接收迫零法。然而，值得 注意的是梯度法和迫零法之间的性能增益差随着在信源/中继端的天线数量N增加而慢慢 消失。可以得到这样的结论:在小规模的MMO全双工中继系统中，梯度上升法实现的平均速 率优于迫零法所实现的平均速率，而考虑系统传输性能和成本后可以得到，迫零法（即低复 杂度次优方案)适用于大规模的ΜΙΜ0全双工中继系统。
[0099] 本发明不仅局限于上述【具体实施方式】，本领域一般技术人员根据本发明公开的内 容，可以采用其它多种具体实施方案实施本发明。因此，凡是采用本发明的设计结构和思 路，做一些简单的变化或更改的设计，都落入本发明保护范围。
【主权项】
1. 一种ΜΙΜΟ全双工中继系统信源中继联合预编码方法，其特征在于，包括以下步骤： (1) 初始化:系统确定源节点的天线数Ns、中继节点接收端的天线数Nr、中继节点发射端 的天线数Ντ、信宿节点的天线数Nd、信源的功率上限Ps、中继的功率上限Pr;设定中继处噪声 功率cr〗，信宿处噪声功率; (2) 利用导频方法进行信道估计得到信源-中继、中继发射端-中继接收端、中继-信宿 的信道矩阵，分别记为HSR、H RR、HRD;设定中继放大矩阵(3 = *^，其中^、&分别为一维发射 波束成形向量和一维接收波束成形向量，将MMO全双工中继系统的速率最大化问题等价为 仅关于Xr的无约束优化问题：(3) 利用梯度上升法得到步骤2所述无约束问题的最优解Xr，并按如下公式计算信源预 编码矩阵V、中继发送波束成形 xt，中继接收波束成形Xr:(4) 在信源端根据V对传输符号进行预编码，同时根据Xt，Xr计算得到中继放大矩阵 Q = XiXf，利用Q对中继接收端的接收信号进行放大处理，完成MMO全双工中继系统的预编 码设计，最终实现信源与信宿之间以全双工的方式进行信息传输。2. 根据权利要求1所述的一种MIMO全双工中继系统信源中继联合预编码方法，其特征 在于，所述步骤3中梯度上升方法，具体包括以下步骤： (3.1) 初始化迭代次数t = 0、迭代次数上限tmax，初始化Xr为矩阵最大特征值所 对应的特征向量； (3.2) 更新迭代次数:t = t+Ι;令U1为Π Η:最大特征值Amax所对应特征向量，由下式 计算得到矩阵Had riHfD最大特征值的梯度Vitaax ，BP :继而根据计算得到步骤2所还尤约束问题目标凼数夫fXr的梯度grad_f; (3.3 )判断是否满足grad_f彡ε或t彡tmax，其中ε表示判定阈值，其值在〇 . OOl~ 0.00000 1之间，若满足，则输出最终解Xr，否则根据线搜索方法确定步长β，更新Xr = Xr+β X grad_f，并重复步骤3 · 2~3 · 3 〇
【文档编号】H04B7/04GK105897319SQ201610371146
【公开日】2016年8月24日
【申请日】2016年5月27日
【发明人】陈志勇, 史清江, 徐伟强
【申请人】浙江理工大学