使用扩展点阵列的map解码方法

文档序号:10541112阅读:374来源:国知局
使用扩展点阵列的map解码方法
【专利摘要】本发明涉及用于通过噪声信道接收到的信号的MAP解码方法。所述信号是由具有非均匀概率分布的码表符号组成。所述符号通过来自点阵列(Λ)的点来表示。所述符号的概率分布通过高斯分布来近似。根据所述点阵列(Λ)以及噪声方差与符号的高斯分布中的方差之间的比值(β)来构造扩展的点阵列(Λaug)。然后所提供的MAP解码方法回到根据ML标准通过从扩展阵列搜索最接近于表示接收到的信号(yaug)的点来执行解码。
【专利说明】
使用扩展点阵列的MAP解码方法
技术领域
[0001] 本发明涉及数字电信领域,尤其涉及信息源的解码,该信息源的符号不遵循均匀 分布。本发明尤其适用于在使用物理层网络编码(Physical Layer Network Coding)技术 的网络中进行解码。
【背景技术】
[0002] 数字通信通常基于能够发送由符号组成的消息的发射机和能够复制来自接收到 的信号的消息的接收机之间的离散信道模型。
[0003] 离散信道可看作是接受属于码表(alphabet)Ax的符号Xl作为输入并且提供属于码 表Αγ的符号 yi作为输出的随机系统,信道输入通过由以下概率所定义的概率模型与信道输 出相关:
[0004] P(Yl = yi, . . . ,Ym = ym|Xl = Xl, . . . ,Xn = Xn) (1)
[0005] 在实践中信道是有噪声的,换言之,信道输出端处的符号受到通常被认为是加性 高斯白噪声(,BBAG)的噪声的影响。接收到的信号采样的采样记作Ri = Yi+Zi,i = l,. . .,m, 其中,Zi为噪声采样。
[0006] 接收机搜索给定的接收信号的符号Xl,...,&的最可能的序列,换言之,接收机搜 索使下列条件概率最大化的输入符号的序列:
[0007] P(Xi = xi, . . . ,Xn = xn| Ri = ri, . . ,Rm=rm) (2)
[0008] 我们将上式简单地记作P(X I R),其中,X为输入符号的待估计向量,R为接收到的信 号采样的向量。
[0009] 换言之,我们搜索使下列概率最大化的向量Λ'
[0010] P(X|R)=P(R|X)P(X)/P(R) (3)
[0011] 因为项P(R)不依赖于X,需要做的就是搜索使下列概率最大化的向量1:
[0012] P(R|X)P(X) (4)
[0013] 在该阶段,两类解码方法之间存在区别。第一类解码在假设码表Αχ中的符号是等 概率的基础上使用被称为最大似然(Maximum Likelihood,ML)标准的标准。然后所需要的 就是将P(R IX)最大化,换言之,当信道没有记忆时:
[0014]
(5) ^ i-Ι .'χ?
[0015] 当输入符号实际上等概率时,也就是当码表Αχ具有均匀概率分布时,该解码方法 是最佳的。
[0016] 被称为最大后验(Maximum A Posteriori,ΜΑΡ)解码的第二类解码不进行这种简 化假设。因此,条件概率必须根据表达式(5)通过输入符号的概率进行加权。当信道没有记 忆时,上述表达式变为:
[0017]
( 6 )
[0018] 虽然大量的信息源具有其码表的符号的非均匀概率分布(以下用非均匀信息源来 表示),尽管ML解码具有次佳(sous-optimal)属性,但是由于其较简单,所以ML解码通常是 优选的。
[0019] 由于包含在非均匀信息源中的连续符号(比特)之间的强相关性,非均匀信息源包 括大部分的多媒体信息源(文本、语音、音频、视频)。虽然信息压缩能够显著地减少冗余性, 但是并不能完全消除冗余性。
[0020] 待解码符号的概率分布的非均匀性并不是信息源的必然结果,但这可能是由所使 用的通信策略造成的。因此,使用物理层网络编码(Physical Layer Network Coding)技术 的协作网络中的中继终端对源于不同源终端的信号的组合进行接收。S.Zhang等人发表在 2006年ACM移动通信会刊(Proc.Of ACM Mobicom)第358-365页的论文《物理层网络编码》 (《Physical Layer Network Coding》)给出了物理层网络编码的介绍。虽然由源终端发送 的符号的概率分布是均匀的,但是这些符号在中继终端处的组合不是均匀的。因此,在图1 示出的简单示例中,两个源终端Si,S 2&送符号χ#ΡΧ2到两个中继终端h,R2,两个中继终端 R!,办在对符号χ#ΡΧ2进行重编码(如果有必要的话)以及将其重发到目标终端D之前解码符 号xdP X2;在该示例中,如果源和中继之间具有单元信道的增益,则两个中继将接收到和符 号^+幻。如果源终端使用具有均匀概率的BPSK码表,则和符号_2、0、2的概率将分别为1/4、 1/2和1/4。更一般地,如果中继终端接收到来自K个源终端的BPSK符号,则从中继端接收到 的和符号将遵循二项式概率分布B(K,2)。相同的结论对QPSK符号是有效的,此时概率分布 在每个同相轴和正交轴上是二项的。
[0021] 在上述信息源最初的非均匀性、由网络引起的非均匀性)的两种情况下,接收机需 要对属于具有非均匀概率分布的码表的信息符号进行解码。
[0022] 如上所述,最佳解码方法为MAP解码,但是由于其高度的复杂性,这种解码方法在 实践中很难实现。特别地,在网络编码情况下,要在来自不同源终端的符号的所有可能的组 合上进行穷举搜索。此外,ML解码就误码率来说是次佳的。此外,使用网络编码的电信系统 中的现有ML解码方法也非常复杂。
[0023] Neva等人发表在2010年4月IEEE通信学报(IEEE Trans.on Com.)第58卷第4期第 1151-1160页上的论文《在不完全CSI下MUM0系统中高斯星座的探测》(《Detection of Gaussian constellations in MUM0 systems under imperfect CSI》)公开了一种ΜΙΜΟ系 统中的MAP检测器,由于被设置为用于传输的格式,所以所接收到的符号具有拟高斯振幅分 布。
[0024] 因此,本发明的目的是提供一种对属于具有非均匀概率分布的码表的信息符号进 行解码的方法,该方法比通过穷举搜索的MAP解码简单得多,并且与通过最大似然(ML)的传 统解码相比具有更低的误码率。

【发明内容】

[0025] 本发明被限定为通过有噪声信道接收到的信号的MAP解码方法。所述信号由属于 预设码表的符号组成并且受到加性高斯白噪声的影响。所述码表中的所述符号的概率分布 是非均匀的,所述符号通过由NXN矩阵Μ生成的格A = {X|x = M:a, aeZj中的点来表示,其 中:
[0026] 确定表示接收的信号的N维向量y;
[0027] 构造增广向量Υ ,其中On是N维0向量; 9 L- Af_
[0028] 使用生成元矩网
杉成增广格
, 其中# = ^是噪声的方差与对所述码表中所述符号的非均匀概率分布进行建模的高斯概 率分布的方差之间的比值;
[0029] 在所述增广格中的点中进行搜索以找出最接近于所述向量的表示接收到的信号 的点L
[0030] 通过生成矩阵Μ和最接近的近邻对接收到的所述信号进行估计,即免= ΛΜ。
[0031] 有利地,对最接近的近邻的搜索被限制在所述格的表示信号能量上的约束的有限 子集Ω a上。
[0032] 根据本发明的第一实施例,有噪声信道是由源终端和中继终端之间的多条K基本 信道组成的信道,并且符号中的每一个符号是由通过对应的源终端发送到所述中继终端的 基本符号的组合来组成,每个所述基本符号属于源终端的基本码表。
[0033] 在一个示范应用中,源终端的基本码表都是相同的,并且通过& =厂〇-|得到对码 表中所述符号的非均匀分布进行建模的高斯概率分布的方差2,其中是基本码表中基 本符号的概率分布的方差。
[0034] 在这种情况下,通过σ〗=ΛΤ_得到所述方差4,其中Pmax为源终端的最大传输功 率。
[0035]有利地,使用球形解码对最接近的近邻进行搜索,并且球形解码根据Pohst算法使 用点的列举法。
[0036] 可选地,可通过使用具有球形界的堆栈解码器进行解码来搜索最接近的近邻。
【附图说明】
[0037] 以下将结合附图通过优选实施例对本发明的其他特征和有益之处给出更详细的 说明,其中:
[0038] 图1为使用物理层网络编码的协作网络的简化示例;
[0039] 图2示出了针对本发明一个实施例的MAP增广格解码方法的示例性应用的随信噪 比变化的误码率曲线;
[0040] 图3为根据本发明一个实施例的MAP增广格解码方法的流程图。
【具体实施方式】
[0041] 以下我们将考虑具有非均匀概率分布的信息符号的码表。换言之,该码表中的符 号不是等概率的。
[0042] 例如该码表可以是非均匀信息源所使用的码表或者可以是采用物理层网络编码 的协作网络中的中继端处的组合结果。
[0043] 为了说明的目的而且又不降低普遍性,我们将公开本发明的作为协作网络的一部 分的解码方法。为此,我们将假设K个源终端S1; . . .,SK分别传输K个符号X1,. . .,χκ,并且假 设中继终端接收与这些符号的组合相对应的信号,也就是:
[0044]
(7)
[0045]其中gk,k=l,...,Κ为源终端和中继终端之间的信道中的每条信道的增益。在不 降低普遍性的前提下,我们假设这些增益都等于1。
[0046] 符号Xk为被称为基本码表的有限码表Αχ= Λ Π Ω的元素,其中Λ为Ν维格,Ω为 的凸面(convexe),其包括原点并且限定了源终端的传输功率限制(假设所有终端的传输功 率限制相同)。
[0047]例如,如果符号X#QAM码表中的元素,则格的维数将为2。注意,通常N维格是由N XN矩阵Μ定义的,其中矩阵的列向量称为格生成元向量(vecteurs g6n6rateur du r6 seau),格中的由x = Ma ( aeZ〃)以向量形式定义的点χ是Ν维向量,其中元素为整数。然 而,可以理解的是对于QAM码表(a矣On)来说并没有考虑格的原点。然而,可以理解的是如果 矩阵Μ能够修改,我们可以一直认为增益g k等于1。
[0048] 凸面Ω典型地为的中心在原点的球面,球面的半径由源终端的最大传输功率 给出。因此,如果pmax是源终端的最大传输功率,则凸面Ω可由点X进行定义,使得:
[0049]
( B )
[0050] 其中E{.}是均值。
[0051] 表达式(7)中的向量z为N维噪声向量,其分量为具有零均值和方差为〇2的随机高 斯变量。
[0052] 中继终端设置有解码器以对和ν = 进行解码。解码后如果适用的话,中继终端 k 二 1 将其传输到目标终端。
[0053] 和符号v = 属于由基本码表Ax的叠加所产生的码表Αν。假定格Λ通过加法是 k'=l 稳定的(具有加法的格有群结构),和符号v仍然属于格Λ。结果是定义了码表Αν的格Λν为Λ 的子集()。和符号的码表Αν使得Αν= ΛνΠ Ων,其中Ων为R#的反映 Κ个源终端的传 输功率限制的凸面。
[0054] 在中继处的MAP解码将导致对和向量七^的搜索,使得:
[0055] .(. 9' ).
[0056]
[0057] ( 10 )
[0058] 本发明的基本构思是使用具有零均值和方差为4 =1^4 = 的随机高斯变 量来对和符号ν进行建模。通过将中心极限定理应用到随机变量处上来证明这种近似法是 正确的。
[0059] 基于该近似法,和符号的概率分布定律如下:
[0060]
.(11)
[0061] 将对P(v)的估计带入到表达式(10)中,我们得到新的解码标准:
[0062]
C 12 )
[0063] 考虑到第一个项不取决于V,解码由以下确定组成:
[0064]
( 13 )
[0065]增加通过向量y和N维零向量的垂直串接得到的2N维向量yaug(称为增广向量),即
,以及满秩矩¥
,其中In为NXN单位矩阵,β为系数,其反映了噪声方差 和将被解码的和符号的方差之间的比值,即# =
[0066] 则解码标准(13)可表示为:
[0067]
(⑷
[0068] 其中通过将向量a限制为属于的子集来示出传输功率极限上的条件,使得传输 功率限制得到满足,换言之
L i V 、 ^ J
[0069] 如果我们通过
g义出2NXN增广矩阵Maug,则解码标准最终表示 为以下形式:
[0070]
(15)
[0071] 换言之,这等同于在由矩阵Maug生成的增广格Aaug中搜索最接近于点yaug的点,点 yaug表示所接收到的信号。因此,可以理解的是这等同于2N维增广空间中的传统ML解码,其 中格Λ aug和表示接收到的信号的点yaug两者都是2N的。
[0072] 典型地可使用球形解码来搜索增广格Aaug中yaug的最接近的近邻。应该注意到,球 形解码可将对最接近的近邻的搜索限定到属于噪声球的格,该噪声球以表示接收到的信号 的点为中心Viterbo等人发表在IEEE信息理论学报(IEEE Transactions on Information Theory) 1999年7月第45卷第1639-1642页上的文献《用于衰落信道的通用格 码解码器》(《A universal lattice code decoder for fading channels》)给出了球形解 码的描述。根据G.Rekaya Ben-Othman发表在IEEE无线&移动计算、网络&通信国际会议录 (Proc. of IEEE International Conference on ffireless&Mobile Computing, Networking&Communi cat ion)第322-327页上的文献《球形界堆栈解码器》(《The spherical bound stack decoder》)中的描述,还可以使用球形界堆栈解码器(SB)。
[0073]使用球形解码器搜索最接近的近邻就其本身而言尤其可以使用Pohst列举法。 [0074]球形解码器从以表示接收到的信号的点yaug为中心且具有半径C的噪声球面开始。 格Λ aug中属于该球面的点a满足下列关系式:
[0075] | |yaug-Mauga| |2<C2。 (16)
[0076] 搜索最接近的近邻还必须考虑以下事实:点a必须满足功率限制,即ae Ω a。
[0077] 在准备步骤中,增广矩阵Maug遵从QR分解,即Maug = QR,其中Q为2NXN矩阵,Q中的列 向量是正交的并且R为NXN上三角矩阵。若记凡,$以及| = p-a,则关系式(16)可写 为:
[0078] |RC | | < C2 (17)
[0079] 或者:
(18 )
[0082]用a的第N个分量开始进行搜索:
[0083]
(19)
[0084] 考虑到馬£泛这一事实,也就是:
[0085]( 20 ) v ^ yy
[_ V r 期 _
[0086] 类似的,对其他分量的搜索被限制在以下区间上:
[0087]
( 21 )
[0088] 其中Λ = ρ,. + Σ馬A以及Ti = Ti-rpii (Sii-ai),Ti中TN = C2。以从最后一个分量到 /=1 第一个分量的顺序进行搜索,选择具有减小了对较低下标分量进行搜索的区间的给定下标 的分量。
[0089] 可使用以下简化的关系式来表示搜索区间(20)和(21):
[0090] b~ / = 1,..., A/ ( 22 )
[0091]
[0092] 搜索区间还必须限制为考虑最大传输功率限制ae 的情况。该限制规定a中的 每个分量ai应该保持在两个界之间:
[0093] (〇~ <α? <ω: , / = 1,..., ¥ ( 23:)
[0094] 因此,最后对于每个分量ai的搜索区间减小为:
[0095] supift;; ,h, ) < a, < inf Jr ) ( 24 )
[0096] 然后继续进行以下解码过程:第一步是选择区问/、=卜叩内 的分量aN,第二步是搜索区间心广卜叩丨丨,内的候选分量aN-!,区间 Ih上的边界从aN开始计算。如果值aN-A在这个区间内,贝峨们回到N阶以选择另一候选分 量aN。逐步继续该过程直到达到1阶。一旦找到向量S .,其中用于S的分量满足所有的条件 (24),则更新球面的半径并且通过扫描所有的区间来重复进行搜索过程,直到找到最接 近于y aug的向量?。符号^.P*由兔=??给出。
[0097] 图2示出了中继端处的依据信噪比的符号解码误码率。
[0098]在所示的情况下,已经假设N=4,K = 2并且最大传输功率Pmax=l。
[0099] 210涉及对应于经典最大似然(ML)解码方法的曲线,而220涉及本发明的对应于 MAP增广格解码方法的曲线。从该图中可以确认增广格解码方法能够对于给定的信噪比达 到较低的符号误码率,或者对于给定的目标符号误码率达到较好的信噪比。还应该注意到 当格具有高的维数N时,增益也更高。
[0100] 本领域技术人员应理解的是,在不超出本发明的保护范围的条件下,可以使用对 于格内最接近的近邻的其他搜索算法,以及特别是Schnorr-Euchner解码来替代球形解码 器方法。
[0101 ]在总体框架下,我们将公开MAP增广格解码方法。
[0102] 解码方法的目的在于对通过有噪声信道接收到的信号进行解码,所述信号由属于 预设码表的信息符号组成。该码表中信息符号的功率分布是非均匀的,换言之,一些信息符 号与其他信息符号相比更为频密。影响信道的噪声假设为白的和高斯的。
[0103] 假设码表中的符号可通过由关系{x|x = Ma,定义的N维格Λ中的点来表 示,其中Μ为格生成元矩阵(matrice g6n6ratrice du r6seau)。当符号属于QAM型调制星座 时,码表是Λ的有限子部分,其被定义为所述格与JR"的凸面Ω的交集。然而,在QAM码表的 情况下,应该注意到格的原点是被排除在外的。通过由码表中符号的概率分布的不均匀性 导致的不同的功率分布对格中不同的点进行加权。
[0104] 根据本发明,码表中符号的概率分布通过具有方差为的高斯分布来建模。
[0105] 要被解码的信号以y = x+z的形式来表达,其中X是Λ Π Ω中的点,ζ是Ν维噪声向 量,ζ的分量是具有方差为σ2的随机变量。
[0106] 图3概略地示出了本发明的一般实施例的增广格解码方法的流程图。
[0107] 在步骤310中,形成向量y,该向量y表示Ν维空间中接收到的信号,其中Ν为格Λ的 维数。例如该向量可通过对接收到的信号进行解调制来获得。
[0108] 在步骤320中,构造增广向量少胃=^ 。 L iY'_
[0109] 在步骤330中,使用生成元矩罔
形成增广格Λaug,其中Μ为格Λ的生成 元矩阵,而彦为噪声的方差与对符号在码表中的非均勾概率分布进行建模的高斯概率 σ·丨- 分布的方差之间的比值。
[0110]在步骤340中,在增广格Λ aug中搜索增广向量yaug的最接近的近邻,g卩,的使距 离I Iyaug-Mauga| I2最小化的向量i,其中〇3是£?的一部分,其反映了接收到的信号的最大能 量限制。
[0111] 在步骤350中,从MAP标准的意义上讲,通过1_* = Μ?来估计接收到的符号。
[0112] 可以理解的是,如上所述,可以使用球形解码器或者具有球形界的堆栈解码器来 搜索格中最接近的近邻。在这种情况下,已知的不同的球形解码变体是适用的,尤其是使用 Pohst算法对点实施的列举法。
[0113] 本领域技术人员可以理解的是,因为这些和符号形成了通过不同源终端的基本码 表的叠加而得到的码表的一部分,并且该码表中的和符号的分布是非均勾的,所以对于协 作网络的信道的和符号的解码是图3中特殊情况的解码。
[0114] 最后,如上所述,不超出本发明的保护范围的条件下,可以使用对于格内最接近的 近邻的其他搜索算法,以及特别是Schnorr-Euchner解码。
【主权项】
1. 通过有噪声信道所接收到的信号的MAP解码方法,所述信号由属于预定码表的符号 组成并且受附加的高斯白噪声影响,符号在所述码表中的概率分布是非均匀的,所述符号 通过由N X N矩阵M生成的#P的点来表示,其特征在于: 确定(310)表示所接收到的信号的N维向量y;构造(320)增广^ 其中,On是N维零向量; 使用生成元矩P 杉成(330)增广格,其中,I噪声的方差与对所述符号在所述码表中的非均匀概率分布进行建模的高斯概率 分布的方差之间的比值; 在所述增广格中的点中进行搜索(340)以找出与表示所接收到的信号的向量最接近的 点S; 根据生成矩阵M和最接近的近邻通过兔=MS来对所接收到的符号进行估计(350)。2. 根据权利要求1所述的MAP解码方法,其特征在于,对所述最接近的近邻的搜索被限 制到所述格的表示信号能量上的约束的有限子集Ω a。3. 根据权利要求1或2所述的MAP解码方法,其特征在于,所述有噪声信道是由网络的多 个源终端和中继终端之间的K条基本信道组成的信道,并且在于,所述符号中的每个符号由 通过对应的源终端发送到所述中继终端的基本符号的组合来组成,每个所述基本符号属于 源终端的基本码表。4. 根据权利要求3所述的MAP解码方法,其特征在于,所述多个源终端的基本码表是相 同的,并且通过σ,! = 得到对符号在所述码表中的非均匀分布进行建模的所述高斯概率 分布的方差,其中,是基本符号在基本码表中的概率分布的方差。5. 根据权利要求4所述的MAP解码方法,其特征在于,通过4 得到所述方差?τ?, 其中,Pmax为源终端的最大传输功率。6. 根据前述权利要求中任一项所述的MAP解码方法,其特征在于,使用球形解码器来进 行对所述最接近的近邻的搜索。7. 根据权利要求6所述的MAP解码方法,其特征在于,球形解码根据Pohst算法使用点的 列举法。8. 根据权利要求1至5中任一项所述的MAP解码方法,其特征在于,通过使用具有球形界 的堆栈解码器进行解码来进行对所述最接近的近邻的搜索。
【文档编号】H04L1/00GK105900364SQ201480054792
【公开日】2016年8月24日
【申请日】2014年9月29日
【发明人】阿斯马·梅杰里, 伽雅·哈卡亚-本·奥特曼
【申请人】矿业电信学校联盟
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1