一种双向协作中继信道计算转发编码系数向量搜索方法及通信方法
【专利摘要】本发明公开了一种双向协作中继信道计算转发编码系数向量搜索方法及通信方法，该方法针对双向协作中继信道中计算转发编码工作方式的特点，将中继节点的编码系数向量优化问题建模为带有二次约束的整数二次规划模型；针对该优化问题的特点，通过提升、凸松弛、生成割平面等步骤将原优化问题转化为较易求解的新松弛规划问题，通过对松弛规划问题的求解以有效获取原问题的最优解。编码系数向量的选取对网络可达性能指标有着重要影响，本发明提出的方法能有效获取当前信道状态下的最优系数向量组合，为计算转发编码的应用提供了良好的优化基础。
【专利说明】
一种双向协作中继信道计算转发编码系数向量搜索方法及通 信方法
技术领域
[0001] 本发明属于无线通信技术领域，涉及物理层网络编码下的协作通信技术，提出了 一种双向协作中继信道计算转发编码系数向量搜索方法。
【背景技术】
[0002] 无线通信系统中，发射信号通过电磁波在物理层进行传输。在无线多源通信网络 中，发射节点的广播特性使得接收节点可能会在同一时隙内接收到来自于多个不同源节点 的发射信息，这将引起不同传输信号间的相互干扰，从而影响整个网络性能。因此，接收端 如何有效处理多个接收信号之间的相互干扰问题是无线通信技术研究的一个重大挑战。
[0003] 近年来，线性网络编码技术在有线网络应用中已取得了令人瞩目的研究成果。网 络编码具有较强的兼容能力及信息提取能力，这使得解决上述多用户信号间的干扰问题成 为可能。传统的网络编码方案大都运行在MAC层，为了减少对现有无线通信系统软硬件设备 和协议的相应修改，一般采用MAC层资源及用户调度算法尽量减少干扰。但在发送多源数据 时，传统的网络编码方法仍然效率不高。在无线网络中，如何有效利用发射节点的广播特性 来提升无线信道容量显得更加重要。
[0004] 基于嵌套Lattice的计算转发网络编码方案，不但可以解决高阶调制下中继节点 处的解码问题，而且能接近AWGN双向中继信道容量。Lattice编码的结构特点使得叠加后的 信号矢量仍然是一个码字，中继节点只需解码各个码字的线性组合。目的节点通过获取各 中继节点转发的线性组合信息，即可有效解码源节点的发送信息。
【发明内容】
[0005] 发明目的：为了克服现有技术中存在的不足，本发明提供一种双向协作中继信道 计算转发编码系数向量搜索方法及通信方法，本发明能够有效获取当前信道状态下的最优 系数向量组合，为计算转发编码的应用提供了良好的优化基础。
[0006] 技术方案:为实现上述目的，本发明采用的技术方案为：
[0007] -种双向协作中继信道计算转发编码系数向量搜索方法，在双向协作中继信道计 算转发编码方案中，将中继节点的计算转发编码系数向量优化问题建模为带有二次约束的 整数二次规划模型。根据中继节点的转发编码系数向量引入辅助变量将该整数二次规划模 型中的的二次项提升到新的高维空间。通过凸化和松弛过程，将提升到新的高维空间的整 数二次规划模型转化为新的松弛问题，求解该松弛问题，如果得到的最优解满足整数要求， 则为原问题的最优解。否则，根据引入辅助变量及上一步的凸化和松弛约束条件定义一个 割平面约束条件，增加到松弛问题的约束集中，以切掉一部分不满足要求的可行解，缩小可 行域，然后，求解新的松弛规划问题。重复以上过程，直至求出整数最优解。
[0008] 具体包括以下步骤：
[0009]步骤1:在双向协作中继信道的计算转发编码策略中，获取中继节点的系数向量a 的优化目标函。
[0010] 步骤2:根据中继节点的系数向量a引入辅助变量^将中继节点的系数向量a的优 化目标函数的二次项提升到新的高维空间，得到新的中继节点的系数向量a的优化目标函 数。
[0011] 步骤3:采用PSD松弛方法将步骤2中原问题中的非凸约束进行凸化松弛，将步骤2 得到的新的中继节点的系数向量a的优化目标函数改写为松弛问题。
[0012] 步骤4:求解步骤3中的松弛问题，若所得最优解满足整数要求，则为原问题的最优 解，并跳转步骤6。否则跳转到步骤4。
[0013] 步骤5:在步骤3中得到的松弛问题的约束集中根据引入辅助变量及上一步的凸化 和松弛约束条件增加割平面约束条件，以切掉一部分不满足要求的可行解，缩小可行域，然 后，求解新的松弛规划问题。
[0014] 步骤6:返回，输出最优解。
[0015] 所述步骤1中获取的中继节点的系数向量a的优化目标函数如下：
[0016] a = arg min(aTGa)
[0017] s.t. I |a| |2^Ξ1+Ρ| |h| I2
[0019] ai^O ,a2^0
[0020] 其中，a表示中继节点的系数向量，也就是网络编码线性组合的整数系数， a = [at，fl2f eZ2>表示节点51网络编码的整数系数，助表示节点&网络编码的整数系数，E 表示矩阵转置，
.斗表示约束功率，11=[111，112]^為￡肢为节点5 1、32和中 继间的实值信道增益，I表示单位矩阵。
[0021] 所述步骤2中的新的中继节点的系数向量a的优化目标函数为：
[0022] a = arg min<G,A>
[0026]其中，I为2X2的单位矩阵，为2X2的对称矩阵集合，b=l+P| |h| I2。对称矩阵A 可提升表示为A = aaT<XA，B>表示对称矩阵A与B的Frobenius内积，即tr(ATB)。则原问题中的 二次项a TGa可表示为<G，aaT>。
[0027]所述步骤3中得到的松弛问题为：
[0031] 其中，A彡0表示A为对称半正定矩阵。
[0033] Vv.e. I3 ?ν表示任意实数向量。
[0034]所述步骤5中得到的割平面约束条件表示为：
[0036] 中，Vk分别表示的特征值对应的特征向量，A'表示▲的变量空间内的任意一点， k= 1,2〇
[0037] -种双向协作中继信道通信方法，包括以下步骤：
[0038] 步骤(1)，在双向协作中继信道计算转发编码方案中，源节点将各自的信息从有限 域映射到一个嵌套Lattice码字。
[0039] 步骤(2)，源节点同时将映射后的码字信息发送至中继节点。
[0040] 步骤(3)，中继节点接收到来自各源节点的复合信号，并将接收到的信息解码成 Lattice码字的线性组合方程。
[0041 ] 步骤(4)，中继节点广播Lattice方程至第一源节点Si和第二源节点S2,各源节点将 Lattice码映射回有限域，利用自身存储的信息完成解码。
[0042] 有益效果:本发明提供的一种双向协作中继信道计算转发编码系数向量搜索方法 及通信方法，相比现有技术，具有以下有益效果：
[0043] 本发明针对双向协作中继信道中计算转发编码工作方式的特点，将中继节点的编 码系数向量优化问题建模为带有二次约束的整数二次规划模型;针对该优化问题的特点， 通过提升、凸松弛、生成割平面等步骤将原优化问题转化为较易求解的新松弛规划问题，通 过对松弛规划问题的求解以有效获取原问题的最优解。编码系数向量的选取对网络可达性 能指标有着重要影响;经过仿真验证，本发明提出的方法能有效获取当前信道状态下的最 优系数向量组合，为计算转发编码的应用提供了良好的优化基础。
【附图说明】
[0044] 图1为本发明中双向中继信道计算转发策略框图；
[0045] 图2为本发明中编码系数向量对可达计算速率的影响；
[0046] 图3为本发明中不同系数向量搜索方案的计算速率比较；
[0047]图4为本发明中P = 10dB时，信道系数对可达计算速率的影响；
[0048]图5为本发明中P = 20dB时，信道系数对可达计算速率的影响；
[0049]图6为本发明中P = 30dB时，信道系数对可达计算速率的影响；
[0050] 图7为双向协作中继信道通信方法的示意图。
【具体实施方式】
[0051] 下面结合附图和具体实施例，进一步阐明本发明，应理解这些实例仅用于说明本 发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种 等价形式的修改均落于本申请所附权利要求所限定的范围。
[0052] -种双向协作中继信道计算转发编码系数向量搜索方法，如图1所示，在双向协作 中继信道计算转发编码方案中，将中继节点的计算转发编码系数向量优化问题建模为带有 二次约束的整数二次规划模型。根据中继节点的转发编码系数向量引入辅助变量将该整数 二次规划模型中的的二次项提升到新的高维空间。通过凸化和松弛过程，将提升到新的高 维空间的整数二次规划模型转化为新的松弛问题，求解该松弛问题，如果得到的最优解满 足整数要求，则为原问题的最优解。否则，根据引入辅助变量及上一步的凸化和松弛约束条 件定义一个割平面约束条件，增加到松弛问题的约束集中，以切掉一部分不满足要求的可 行解，缩小可行域，然后，求解新的松弛规划问题。重复以上过程，直至求出整数最优解。
[0053] 如图1所示，中继节点将多址接入阶段获取的接收信号yeir的线性组合方程的估 计值广播至各源节点。各源节点通过解码器q -Ff将接收码字映射回有限域，减去自 身存储的信息，即可恢复出所需的估计信息&。
[0054]具体包括以下步骤：
[0055] 步骤1:在双向协作中继信道的计算转发编码策略中，获取中继节点的系数向量a 的优化目标函，中继节点恢复的线性方程系数向量不同，则各节点可达计算速率也不同，为 了最大化系统容量，将中继节点的计算转发编码系数向量优化问题建模为带有二次约束的 整数二次规划模型。其中，中继节点的系数向量a的优化目标函数如下：
[0056] a = argmin(aTGa)
[0058] 其中，a表示中继节点的系数向量，也就是网络编码线性组合的整数系数， 3=扣1，《2/ €^2，&1表示节点&网络编码的整数系数，&2表示节点&网络编码的整数系数， 1 表示矩阵转置：
，？表示约束功率，h = [hi，h2 ]τ ΛΛ e R为节点Si、S2和中 继间的实值信道增益，i表示单位矩阵。
[0059] 步骤2:根据中继节点的系数向量a引入辅助变量^将中继节点的系数向量a的优 化目标函数式（1)的二次项提升到新的高维空间，。令Αυ = &^(1彡i，j<2)，对称矩阵A可 提升表示为A = aaT，〈A，B>表示对称矩阵A与B的Frobenius内积，即tr(ATB)。则原问题（公式 (1))中的二次项aTGa可表示为<G，aaT>。因此得到新的中继节点的系数向量a的优化目标函 数：
[0063] 其中，I为2X2的单位矩阵，沪为2X2的对称矩阵集合，b=l+P| |h| I2。对称矩阵A 可提升表示为A = aaT<XA，B>表示对称矩阵A与B的Frobenius内积，即tr(ATB)。则原问题中的 二次项a TGa可表示为<G，aaT>。
[0064] 由上式(2)可见，提升后的最小化问题由原来的二次函数和二次约束转变为关于 (a，A)的线性整数优化问题，降低了求解难度。
[0065] 步骤3:采用PSD松弛方法将步骤2中原问题中的非凸约束进行凸化松弛，将步骤2 得到的新的中继节点的系数向量a的优化目标函数改写为松弛问题。
[0066] 在原问题式（1)中，二次约束为凸约束。但提升过程中引入的约束条件A = aaT# 凸，采用适当的松弛方法使优化问题凸化。原问题式（1)可以写为如下等效形式 ζ·ν.
[0067] a = arg min {(G, Λ): (a, Λ) e conv/'j
[0068] 其中，(；〇(1、..歹表不/·的凸包络。
[0069] 不考虑A = aaWPa e f两个约束条件，则问题约束可以表示如下：
[0070] 2 := {(a, A) e IT X 义：〈I，A〉<
[0071] 令厶多0表示A为对称半正定矩阵，由A = aaT可得A-aaT>0,因此可以在此基础上采 用PSD松弛方法将原问题中的非凸约束进行凸化松弛。
[0072] 对于aeR"，满足：
[0074] 因此，线性半正定不等式1:匕〇适用于￡；()1^，#0约束定义如下：
[0075] PSD= ((a,A): A >1-0}
[0076] 在增加 PSD松弛约束后，conv^ = Ζ Π PSD .:
[0077] 令
[0079]则松弛后的目标函数式(2)可表示为：
[0082] 由于X为半正定矩阵，即：
[0083] yr&>0,VveM3
[0084] 可将式(3)中PSD约束替换为如下形式(松弛问题）：
[0085] a = argmin(G, A)
[0086] s,t, (6^,^)2 0
[0087] v'AvSO
[0088] 其中，v表示任意实数向量。
[0089] 步骤4:求解步骤3中的松弛问题，若所得最优解满足整数要求，则为原问题的最优 解，并跳转步骤6。否则跳转到步骤4。
[0090] 步骤5:在步骤3中得到的松弛问题的约束集中根据引入辅助变量及上一步的凸化 和松弛约束条件增加割平面约束条件，以切掉一部分不满足要求的可行解，缩小可行域，然 后，求解新的松弛规划问题。
[0091] 令Γ表示蓋的变量空间内的任意一点。可以通过I：的特征分解来判别Γ是否位 于PSD圆锥(PSD cone)的约束区域内。令AdPvk*别表示f的特征值及其对应的特征向量，k =1，2。为不失一般性，假设λ: <λ2，At〈(X λ?+1，te〇，l，2。若t = 0,说明所有特征值均为非 负，A为半正定矩阵;若技〇,则vA__ = <0 ,k=l，. . .，t，无法满足A为正定 矩阵的要求。
[0092] 因此，可将下式作为割平面约束条件，增加到原来的约束集中，求解新的松弛规划 问题。
[0093] yj Awk > 0 (4)
[0094] 其中，Vk分别表示f的特征值对应的特征向量，瓦''表示$的变量空间内的任意一 点，k=l，2。
[0095] 步骤6:返回，输出最优解。
[0096] -种双向协作中继信道通信方法，如图7所示，包括以下步骤：
[0097] 步骤(1)，在双向协作中继信道计算转发编码方案中，源节点将各自的信息从有限 域映射到一个嵌套Lattice码字。
[0098] 步骤(2)，源节点同时将映射后的码字信息发送至中继节点。
[0099] 步骤(3)，中继节点接收到来自各源节点的复合信号，并将接收到的信息解码成 Lattice码字的线性组合方程。
[0100] 步骤(4)，中继节点广播Lattice方程至第一源节点Si和第二源节点S2,各源节点将 Lattice码映射回有限域，利用自身存储的信息完成解码。
[0101] 图2为本发明中编码系数向量对可达计算速率的影响。码字线性组合系数向量的 搜索问题是计算转发编码策略的关键问题。在信道系数固定的条件下，不同的编码系数向 量将会对网络可达速率指标产生重要影响。假设双向协作中继网络各源节点发送速率对 称，gPRFfe。令信道系数向量11=[1，112]7，112￡[0，1]，功率约束?为1〇(^。不同码字线性组合 的系数向量对信息速率的影响如图2所示。本文分别比较了 a=[0，l]T、a=[-l，-l]T、a=[l， 2]τ三种系数向量情况下的可达计算速率情况。由图中可以看出，可达计算速率与编码系数 向量的选取密切相关。当编码系数向量与当前信道系数相匹配时，网络可达计算速率才能 实现最大化。当h 2 = l时，a=[-l，_l]T与信道系数的匹配状态最佳，故其可达计算速率最大。 因此，为了保证系统计算速率的最大化，需要对编码系数向量做出最优选择。
[0102] 图3为本发明中不同系数向量搜索方案的计算速率比较情况。采用Monte-Carlo方 法，进行1000次随机实验，并将所得结果进行统计平均。节点ShSs和中继间的实值信道增益 111，112相互独立且服从1(0，1)分布，功率约束？的变化范围为0(^至30(^。图\3所示为本文所 提方法、可达计算速率上界、穷搜算法及无系数向量搜索PNC几种不同方案下的可达计算速 率比较情况。计算转发策略的可达计算速率上界可表示为凡= minM :丨/ 21〇以丨_- |/?/1: P)。 穷搜算法通过穷搜约束范围内编码系数向量的所有整数组合，满足当前计算速率最大的解 即为当前信道状态下的最优系数向量。LRCP算法表示本文所提出的系数向量搜索算法。无 编码系数向量搜索PNC方案定义为中继节点不进行系数向量搜索而直接译码并向各源节点 转发信息幻+^。从图中可以看出，不同编码策略所得可达计算速率的差异明显。由于穷搜算 法获取的解为最优系数向量，因此最接近理论上界。LRCP算法所得结果与穷搜算法一致，表 明了本文所提出算法能有效获取网络的最优系数向量组合。无系数向量搜索PNC方案由于 没有针对信道系数的变化进行系数向量的对应调整，因此性能最差，在高信噪比条件下存 在较大的性能差距。
[0103] 图4、5、6分别给出了功率约束P分别为10dB、20dB、30dB时，信道系数改变对可达计 算速率的影响。令信道系数向量11=[1，11 2]7，112￡[0，1]，首先通过本文所提方法获取最优系 数向量，然后计算最优系数向量下的可达计算速率。由图中可以看出，在低信噪比条件下， 本文所提方法、穷搜算法及无系数向量搜索PNC方案的性能较为接近，且均接近可达计算速 率的上界。但随着信噪比的增加，本文所提方法及穷搜算法的性能将显著优于无系数向量 搜索PNC方案。需要注意的是，无编码系数向量搜索PNC方案默认的编码系数向量为a =[ 1， 1 ]τ，当h2>0.8时，该方案与信道系数的匹配程度提高，因此其性能与本文所提方法及穷搜 算法相当。同时，当h2 = 0.5时，本文所提方法可以达到计算速率的理论上界。
[0104] 上述仿真结果表明，编码系数向量的选取对网络可达性能指标有着重要影响，本 文所提方法能有效获取当前信道状态下的最优系数向量组合，为计算转发编码的应用提供 了良好的优化基础。
[0105] 以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种双向协作中继信道计算转发编码系数向量搜索方法，其特征在于：在双向协作 中继信道计算转发编码方案中，将中继节点的计算转发编码系数向量优化问题建模为带有 二次约束的整数二次规划模型;根据中继节点的转发编码系数向量引入辅助变量将该整数 二次规划模型中的的二次项提升到新的高维空间；通过凸化和松弛过程，将提升到新的高 维空间的整数二次规划模型转化为新的松弛问题，求解该松弛问题，如果得到的最优解满 足整数要求，则为原问题的最优解;否则，根据引入辅助变量及上一步的凸化和松弛约束条 件定义一个割平面约束条件，增加到松弛问题的约束集中，以切掉一部分不满足要求的可 行解，缩小可行域，然后，求解新的松弛规划问题;重复以上过程，直至求出整数最优解。2. 根据权利要求1所述的双向协作中继信道计算转发编码系数向量搜索方法，其特征 在于，包括以下步骤： 步骤1:在双向协作中继信道的计算转发编码策略中，获取中继节点的系数向量a的优 化目标函； 步骤2:根据中继节点的系数向量a引入辅助变量^将中继节点的系数向量a的优化目 标函数的二次项提升到新的高维空间，得到新的中继节点的系数向量a的优化目标函数； 步骤3:采用PSD松弛方法将步骤2中原问题中的非凸约束进行凸化松弛，将步骤2得到 的新的中继节点的系数向量a的优化目标函数改写为松弛问题； 步骤4:求解步骤3中的松弛问题，若所得最优解满足整数要求，则为原问题的最优解， 并跳转步骤6;否则跳转到步骤4; 步骤5:在步骤3中得到的松弛问题的约束集中根据引入辅助变量及上一步的凸化和松 弛约束条件增加割平面约束条件，以切掉一部分不满足要求的可行解，缩小可行域，然后， 求解新的松弛规划问题； 步骤6:返回，输出最优解。3. 根据权利要求1所述的双向协作中继信道计算转发编码系数向量搜索方法，其特征 在于:所述步骤1中获取的中继节点的系数向量a的优化目标函数如下：其中，a表示中继节点的网络编码整数系数向量>ai表示节AS1网络编 码的整数系数，助表示节点&网络编码的整数系数，Z表示矩阵转置3表 示约束功率，11=[111，112]^^]?为节点51、&和中继间的实值信道增益，1表示单位矩阵。4. 根据权利要求1所述的双向协作中继信道计算转发编码系数向量搜索方法，其特征 在于:所述步骤2中的新的中继节点的系数向量a的优化目标函数为：其中，I为2 X 2的单位矩阵，S2为2 X 2的对称矩阵集合，b = 1+P||h||2;对称矩阵A可提升 表示为A = aaT;〈A，B〉表示对称矩阵A与B的Frobenius内积，即tr(ATB); 则原问题中的二次项aTGa可表示为〈G，aaT>。5. 根据权利要求1所述的双向协作中继信道计算转发编码系数向量搜索方法，其特征 在于:所述步骤3中得到的松弛问题为：WeR3 ;:v表示任意实数向量。6. 根据权利要求1所述的双向协作中继信道计算转发编码系数向量搜索方法，其特征 在于:所述步骤5中得到的割平面约束条件表示为： 其中，Vk分别表示玉^的特征值对应的特征向量，表示A的变量空间内的任意一点，k =1，2〇7. -种基于权利要求1至6任一所述的双向协作中继信道计算转发编码系数向量搜索 方法的双向协作中继信道通信方法，其特征在于，包括以下步骤： 步骤（1)，在双向协作中继信道计算转发编码方案中，源节点将各自的信息从有限域映 射到一个嵌套Lattice码字；； 步骤(2)，源节点同时将映射后的码字信息发送至中继节点； 步骤（3)，中继节点接收到来自各源节点的复合信号，并将接收到的信息解码成 Lattice码字的线性组合方程； 步骤（4)，中继节点广播Lattice方程至第一源节点Si和第二源节点S2,各源节点将 Lattice码映射回有限域，利用自身存储的信息完成解码。
【文档编号】H04L1/00GK106027206SQ201610616663
【公开日】2016年10月12日
【申请日】2016年7月29日
【发明人】衡伟, 梁天
【申请人】东南大学