一种跳频通信中的迭代检测译码方法及装置的制造方法
【专利摘要】本发明公开了一种跳频通信中的迭代检测译码方法及装置，所述方法包含：步骤101)根据发送端的信息比特序列b、发送序列x、接收端的接收序列r、随机相位偏转θ以及表示多普勒频偏的变量之间的约束关系建立因子图，并依据因子图得到向上，向下，向前和向后的消息更新公式；步骤102)采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散化的频偏值，在步骤101)中建立的因子图上迭代更新相位信息的消息，每次迭代前根据现有信息重新调整频偏设置；且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译码器，进而与译码器联合迭代，迭代结束时，输出译码结果。
【专利说明】
-种跳频通信中的迭代检测谭码方法及装置
技术领域
[0001] 本发明设及跳频通信检测领域，具体设及一种跳频通信中的迭代检测译码方法及 装置。
【背景技术】
[0002] 跳频通信系统因其良好的抗干扰性和低截获概率，广泛应用于现代军事通信领 域。然而在实际的跳频通信系统中，收发晶振的相位噪声使得接收信号存在随机相位偏转， 而系统的高动态特性产生多普勒频移，都对接收端恢复原始信号造成严重干扰。因此必须 对接收信号的相位偏转进行补偿，即载波同步。传统的载波同步方法，如锁相环，收敛速度 太慢，难W应用于快速跳频系统中。其他开环的算法，比如对QPSK信号的四次方的方法，只 在信噪比高时性能良好。
[0003] 基于因子图及和积算法框架的联合检测译码算法近十年来得到广泛关注。因子图 用来表示各变量的联合分布，包含相位噪声和频偏检测部分W及译码部分，和积算法用来 求解该联合分布的边缘分布，W便接收方根据最大后验概率准则恢复信号。每一次迭代， LDPC译码模块与相位恢复模块之间互相传递消息，不断挖掘有用信息，因此系统能工作在 低信噪比下，且需要极少导频。然而，和积算法求解时会出现积分困难，为此现有技术多采 用相位噪声离散化的方法来解决运个问题。然而离散化方法只有在量化精度足够高时才能 达到较好的性能，复杂度过高。而有些现有技术还提出Ti化onov-Single算法，采用单个 Ti化onov分布近似因子图中有关相位的消息，只需传递Ti化onov分布的参数，计算复杂度 大大降低。在此基础上，现有技术还对近似进行优化，得到更为精确的估计。然而由于跳频 通信系统每跳序列较短，上述运些算法性能表现不佳。
[0004] W上运些技术方案都是针对多普勒频偏可忽略的情况研究的。而当飞行器高速运 动、载波频率较高等情况下，多普勒频偏并不能忽略，此时载波相位估计中还应包括对频偏 的估计。现有的算法采用的是将频偏离散化和单个Ti化onov近似，然而运样的处理方式同 样在跳频通信中性能不佳。
【发明内容】
[0005] 本发明实施例提供了一种跳频通信中的迭代检测译码方法及系统，能够利用译码 器和相位检测模块的消息传递和反复迭代来补偿相位噪声和多普勒频偏，进而提高接收机 性能。为了解决上述技术问题，本发明实施例公开了如下技术方案：
[0006] 第一方面，提供的一种跳频通信中的迭代检测译码方法，所述方法包含：
[0007] 步骤101)根据发送端的信息比特序列b、发送序列X、接收端的接收序列r、随机相 位偏转0 W及表示多普勒频偏的变量口之间的约束关系建立因子图，并依据因子图得到向 上，向下，向前和向后的消息更新公式；
[0008] 步骤102)采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散 化的频偏值，在步骤101)中建立的因子图上迭代更新相位信息的消息，每次迭代前根据现 有信息重新调整频偏设置;且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译 码器，进而与译码器联合迭代，迭代结束后输出译码结果；
[0009] 其中，将因子图中的相位消息近似为若干个吉洪诺夫Ti化onov分布的和，且构成 每个消息的吉洪诺夫Ti化onov分布的个数根据实际情况中对复杂度和估计性能折中要求 进行设定;根据相位估计均值划分区域，在每个区域内将多个吉洪诺夫Ti化onov分布近似 为单个吉洪诺夫Ti化onov。
[0010] 结合上述第一方面，在第一种可能的实现方式中，上述步骤102)进一步包含：
[00川步骤102-U将多普勒频偏连续变量巧近似为[-狗，狗化L个均匀分布的离散值 口W，作为初始迭代的L条路径；
[0012] 针对每条路径进行如下处理：
[0013] 步骤102-1-1)将所有表示Ti化onov分布权值大小的W设为1，表示相位目分布的参 数。受为0,表示频偏路径权重的T设为1;
[0014] 初始化因子图上和发送序列X相对应的边X向下的消息和随机相位偏转变 量对应的边0向下消息的参数^"，作为初始迭代的值；
[0015] 步骤102-1-2)计算向下消息的权值;
[0016] 步骤102-1-3)在每条路径上，从左到右计算前向消息；
[0017] 步骤102-1-4)在每条路径上，从右向左计算反向消息；
[0018] 步骤102-1-5)计算向上传递的消息；
[0019] 步骤102-2)当达到对应信噪比下设定的迭代次数或者译码器相邻两次输出的外 信息变化值小于一个设定值时，迭代计算结束，迭代译码输出的结果即为恢复后的信息序 列；
[0020] 若不满足上述迭代结束的条件，则计算向下的消息选择另一条路径并返 回步骤102-1-2)。
[0021] 结合上述第一方面，和第一种可能的实现方式，在第二种可能的实现方式中，当采 用PSK调制时，
[0022] 所述步骤102-1-3)进一步包含:首先，向下传递的消息和向左传递的消息，经过等 号结点，相乘得到MQ个吉洪诺夫Ti化onov分布的和；
[0023] 其次，用MinKU ?)将相乘得到的MQ个吉洪诺夫Ti化onov分布变为Q个；
[0024] 然后，再用如下公式计算如下参量:前向消息中每个吉洪诺夫Ti化onov分布的权 值苗;;、参数立3化及频偏路径权值品> .
[0025]
[0026]
[0027]
[0028] 其中，角标k表示一跳中第k个符号，角标q表示一条消息中第q个吉洪诺夫 Ti化onov分布;上角标《'表示第1条频偏路径;1?々和表示第"k-1"个符号的消息经过 MinKU ?)后合并得到的消息中各吉洪诺夫Ti化onov分布的权值和参数；
[0029] 参数表示消息更新后各吉洪诺夫Ti化onov分布的权值和参数；表示 利用前向消息计算得到的该频偏路径的权值，且越大表示该频偏越越接近实际频偏;〇? 完 表示随机游走相位模型中参数A的方差;Q为前向和反向相位消息中Ti化onov分布的个数， M为PSK的调制阶数。
[0030] 结合上述第一方面，和/或第一种可能的实现方式，和/或第二种可能的实现方式， 在第=种可能的实现方式中，步骤102-1-4)进一步包含：
[0031] 首先，向下传递的消息和向左传递的消息经过等号结点相乘得到MQ个吉洪诺夫 Ti Idiono V分布的和；
[0032] 其次，用最小化距离近似手段MinKU ?)将相乘得到的MQ个吉洪诺夫Ti化onov分 布变为Q个；
[0033] 然后，再用如下公式求得如下参数:反向消息中每个吉洪诺夫Ti化onov分布的权 值每;;、参数衣；W及频偏路径权值之().
[0034]
[0035]
[0036]
[0037] 其中，参数:设4和話,)。表示第"k+1"个符号的消息相乘并经过MinKU ?)后合并得 到的消息中各了；[化0]10￥分布的权值和参数;参数1_1，|/"^表示消息更新后各了1化〇11〇￥分布 的权值和参数；^')表示利用反向信息后计算得到的该频偏路径的权值。
[0038] 结合上述第一方面，和/或第一种可能的实现方式，和/或第二种可能的实现方式， 在第四种可能的实现方式中，步骤102-1-5)进一步包含：
[0039] 首先，向左传递的消息和向右传递的消息经过等号结点，相乘得到Q2个Ti^onov 分布的和；
[0040] 其次，在每条路径上用MinKU ?)将相乘得到的妒个Ti化onov分布变为M个；
[0041 ]然后，再用如下公式计算yt(xk):
[0042]
[0043] 其中，yt(xk)为边Xk上向上传递的消息，该消息将送给译码器；；；的4，/''.>和 J做，巧巧)分别表示输入同一个等号结点的向右和向左的消息，和幻：T为经过MinKL (?)合并后向上相位消息中各Ti化onov分布的权值和参数。
[0044] 第二方面，本发明实施例还提供了一种跳频通信中的迭代检测译码装置，所述装 置包含：
[0045] 存储器，用于存储有关数据；
[0046] 处理器，用于执行所述程序，进而进行如下处理：
[0047] 根据发送端的信息比特序列b、发送序列X、接收端的接收序列r、随机相位偏转0 W 及表示多普勒频偏的变量與之间的约束关系建立因子图，并依据因子图得到向上，向下，向 前和向后的消息更新公式；
[0048] 采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散化的频偏 值，在建立的因子图上迭代更新相位信息的消息，每次迭代前根据现有信息重新调整频偏 设置;且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译码器，进而与译码器联 合迭代，迭代结束后输出译码结果；
[0049] 其中，将因子图中的相位消息近似为若干个吉洪诺夫Ti化onov分布的和，且构成 每个消息的吉洪诺夫Ti化onov分布的个数根据实际情况中对复杂度和估计性能折中要求 进行设定;根据相位估计均值划分区域，在每个区域内将多个吉洪诺夫Ti化onov分布近似 为单个吉洪诺夫Ti化onov。
[0050] 结合上述第二方面，在第一种可能的实现方式中，当采用PSK调制时，所述处理器 进一步包含：
[0051] 相位检测模块，用于估计频偏的取值和相位噪声的分布，进而进行补偿，计算出关 于发送序列的软信息，送给译码模块进行迭代译码;所述软信息为相位检测计算得到的发 送序列的概率分布；
[0052] 逆映射模块，用于PSK解调，连接译码模块和相位检测模块；
[0053] 译码模块，用于进行迭代译码，每次迭代完成将软信息回送给相位检测模块;满足 迭代条件时，输出译码结果;其中，所述软信息为译码计算得到的序列取值的概率分布；
[0054] 迭代判决模块，用于判断是否满足迭代结束的条件，当满足结束条件时则停止迭 代，由译码模块输出译码;否则返回相位检测模块继续上述相位检测模块、逆映射模块和译 码模块的迭代处理，直至满足迭代结束条件。
[0055] 结合上述第二方面，和/或第一种可能的实现方式中，在第二种可能的实现方式 中，相位检测模块进一步包含：
[0化6]频率离散化处理子模块，用于将接收的多普勒频偏连续变量0近似为[-挪，骑]上L 个均匀分布的离散值口W，作为初始迭代的L条路径；W及
[0057] 相位估计子模块，用于估计相位噪声的分布；
[0058] 补偿子模块，用于根据得到的相位噪声的分布计算发送序列X的软信息，即所述软 信息为通过相位检测获得的发送序列X的概率分布；
[0059] 输出子模块，用于向逆映射模块输出软信息，进而进行PSK逆映射。
[0060] 总之在发明中，提出一个针对跳频通信系统改进的联合检测译码算法。本发明提 供的算法采用频偏离散化，并且将因子图中的相位消息近似为Ti化onov分布的和，其中构 成每个消息的Ti化onov分布的个数可W人为设定，来平衡计算复杂度和性能。为解决和积 算法带来的计算爆炸问题，我们根据相位估计均值划分区域，在每个区域内将多个 Ti化onov分布近似为单个Ti化onov。本文采用随机游走模型作为每跳内部的相位模型。
[0061] 与现有技术相比，本发明的技术优势在于：
[0062] 本发明提出了一种适合快速跳频系统的联合迭代解调译码算法，相比现有的技 术，对信息挖掘更充分，在低信噪比下拥有优越的估计性能。虽然现有某些算法能达到相同 性能（如对相位和频偏都进行离散化），但复杂度远远高于本发明。因此本发明提供了一种 性能优越同时低复杂度的联合迭代解调译码实现方案。第二方面，本发明技术方案应用于 通信系统接收方基带解调模块，将相位检测和迭代译码联合，利用译码器和相位检测模块 的消息传递和反复迭代来补偿相位噪声和多普勒频偏，进而提高接收机性能。
【附图说明】
[0063] 图1是本发明实施例提供的因子图；
[0064] 图2是本发明实施例提供的算法流程图；
[0065] 图3是本发明实施例提供的系统的组成框图。
【具体实施方式】
[0066] 下面结合附图对技术方案进行详细阐述。
[0067] 实施例
[0068] 1、本发明技术方案的系统模型：
[0069] 在实际的跳频系统中，假设传输的信息比特序列为b，在发送端先经过LDPC编码变 为序列C，再经过MPSK映射，并添加导频后，形成基带发送序列X。设N为LDP邱马字长度，H为序 列X每帖包含的跳数，L为序列X每跳包含的信息符号个数，P为每跳添加的导频符号，则每帖 实际发送的符号数为"N+pH"。假定信道为AWGN信道，接收信号经解跳、匹配滤波、采样等后 得到的数字基带信号为：
[0070]
[0071] 其中，nk为接收机引入的独立同分布的复高斯噪声，均值为0,实部和虚部的方差 均为〇2D0k为未知的相位偏转。由于射频跳频单元的影响，各跳之间的随机相位噪声是相互 独立的。
[0072] 假定每跳第一个符号的随机相位偏转服从[0,231)上的均匀分布，且在各跳内服从 如下随机游走横巧：
[0073]
[0074] 其中，Ak为独立同分布的实高斯随机变量，均值为0,方差夫
其中
Ts为符号间隔，Af为多普勒频移，假定它服从[-Af" Af"1 巧匀分布，并且在一帖之内 保持不变。因此，巧服从[-錢，掷]上的均匀分布，其中
[0075] 为了恢复原始序列b，根据最大后验概率准则，接收机的目标是：
[0076]
口）
[0077] 显然，P(bi Ir)可W通过求解联合后验概率密度函数P化代0内的边缘概率得到。 由概率论的巧关知巧易知；
[0079] 其中I[?访指示函数。此外，[0080]
[007引 (4)
[0081]
[0082] 根据W上内容得到图1所示的因子图，该因子图表示了信息比特序列b、发送序列 X、接收序列rW及随机相位偏转0的约束关系。
[0083] 2、下面将SPA应用于该因子图。得到消息更新公式如下：
[0084]
[0085]
[0086]
[0087]
[0088] 如前所述，上式在实际计算中会出现积分困难，现有技术中采用离散化或近似来 解决。为了进一步解决上述问题，本发明还提供了一种基于Ti化onov-Sum算法的改进方法 计算上述公式的过程。具体为：
[0089] 首先，介绍Tildionov-S皿近似：
[0090] Ti化onov分布T( 0 ; O也被称作循环正态分布，是一种圆上概率分布模型。参数为 复数，ZC表示随机变量0的均值，I引反比与0的方差。
[0091] Tildionov分布有如下性质：
[0092]
[0093]
[0094]
[00M]基于上述算法计算向下传递的消息，易得：
[0096]
[0097]
[009引
[0099] 再计算前向和反向传递的消息，由因子图和概率论的相关知识可得：
[0100]
[0101]
[0102]
[0103]
[0104]
[0105]
[0106] W前向为例J 首先将连续变量梦近似为[-辦，稱] O 上L个离散值口= 1，..Z >可视作L条并行计算的路径。在每条路径上，每跳内部经过n个 符号的前向传递后，得到的消息为r个Ti化onov分布的和，运显然会导致计算爆炸。因此， 本专利申请将每个J辦)和的J都近似为Q个Ti化onov分布的和。运里Q是可变的，用来平 衡复杂度和性能的需求。则：
[0107]
[0108] 首先来看前向的消息传递情况，
[0109]
(21)
[0110] 即//(马1),7/0? I,/I)为MQ个Tildionov分布的和。我们将运MQ个Tildionov分布根 据其ZC按[-n+23i(q-l)/Q，-3i+23iq/Q)划分为Q个区间（q=l，...，Q)，在每个区间内将属于 该区间的若干个Ti化onoV分布近似为一个Tikhonov分布，近似方法采用最小化距离 化uUback-Leibler Divergence)法。KL距离是相同事件空间内两个概率分布差异的度量。 定义如下：
[0111]
[0112] 我们的目标是用单个Ti化onov分布近似多个Ti化onov分布的和，即求解：
[0113] )
[0114] 略去中间推导的详细步骤，可得：
[0115]
[0116]
[0117] 其中
Io(X)和Ii(X)分别为零阶和一阶的第一类贝塞尔函数。
[011引上述近似方法记做MinKU .)。[0119] 现在，我们有
[0120]
[0121]
[0122]
[0123]其中，
[0124]
[0125]
[0126]
[0127]
[0128] 前向和反向的消息更新后，就可W计算向上的消息了。运里需要将Q2个分布的和 近似为M个，得到：
[0129] )
[0130] 需要注意的是，第一次迭代时所选择的L个离散值《於是均匀分布[-狗，狗]上的，因 为巧服从[-卿，狗]上的均匀分布。定义
[0131]
[0132] 当第一次迭代结束时，丫 最大的路径被认为是频偏估计最准的路径，后续迭代 时只需保留3条路径进行计算。下次迭代时只需在该条路径对应的频偏附近进行计算，频偏 捜索半径每次迭代时减少为上次的一半，因此估计精度随着迭代次数的增加而增加。
[0133] 结合上述公式说明本发明详细的相位检测和迭代的具体处理过程：
[0134] (1)初始化:分为=步。首先将将连续变量^^近似为[-狗，％]上1个均匀分布的离散 值巧W >作为初始迭代的L条路径;其次，在每跳路径上，将所有表示Ti化onov分布权值大小 的W设为1，表示0分布的参数C设为0,表示频偏路径权重的丫都设为1;最后，初始化irkXk,m) 和每跳中未知符号的概率取值yl(Xk,m)在各可能的取值上等概分布，而导频符号的概 O 率取值在某个取值处为1，在其他取值处为0。而(t,,可根据式(15)计算。
[0135] (2)根据式（15)计算向下消息的权值初始迭代时(Xk,m)如步骤（1)所述设 O 置，后续迭代时yl(xk,m)为LDPC迭代译码、MPS即央射模块计算后向下传递的消息。
[0136] (3)在每条路径上，从左到右计算前向消息。首先用MinKU ?)将相乘得到的MQ个 TiWionov分布变为Q个，再用式(28)求得1? ^及^') a
[0137] (4)在每条路径上，从右向左计算反向消息。首先用MinKU ?)将相乘得到的MQ个 Ti化onov分布变为Q个，再用式(30)求得:记、
[0138] (5)计算向上传递的消息。首先，在每条路径上，用MinKU ?)将相乘得到的Q2个 Ti化onov分布变为M个，再用式(31)计算Jit(Xk)。
[0139] (6)MPS即央射、LDPC迭代译码。LDPC迭代译码已经非常成熟，运里不再寶述。
[0140] (7)迭代计算是否结束，可W根据迭代次数是否达到设定次数，或LDCP译码器相邻 两次迭代输出的外信息的变化值小于某个设定值。若计算结束，迭代译码输出的结果即为 恢复后的信息序列;若计算没有结束，译码器和PSK映射模块计算向下的消息将消 息传递给信道估计模块。执行步骤(2)~(7)。
[0141] 图3是本发明所研究点对点跳频通信系统的等效基带模型图。在发射端，传输的信 息比特序列为b，在发送端先经过LDPC编码变为序列C，再经过MPSK映射，并添加导频后，形 成基带发送序列X。接收端收到信号r，在相位检测模块，和迭代译码模块之间不断迭代，最 终实现解调和译码。
[0142] 将接收的序列r输入相位检测模块，相位检测模块用来补偿相偏0和频偏口，即该 相位检测模块的输入为已知的接受序列r，输出为序列X的软信息yt(xk)。所述软信息yt(xk) 为相位检测模块计算的X的取值概率。
[0143] 实施例2
[0144] 本发明实施例还提供了一种跳频通信中的迭代检测译码装置，所述装置包含：
[0145] 存储器，用于存储程序；
[0146] 处理器，用于执行所述程序，进而进行如下处理：
[0147] 根据发送端的信息比特序列b、发送序列X、接收端的接收序列r、随机相位偏转0 W 及表示多普勒频偏的变量0之间的约束关系建立因子图，并依据因子图得到向上，向下，向 前和向后的消息更新公式；
[0148] 采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散化的频偏 值，在建立的因子图上迭代更新相位信息的消息，每次迭代前根据现有信息重新调整频偏 设置;且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译码器，进而与译码器联 合迭代，迭代结束后输出译码结果；
[0149] 其中，将因子图中的相位消息近似为若干个吉洪诺夫Ti化onov分布的和，且构成 每个消息的吉洪诺夫Ti化onov分布的个数根据实际情况中对复杂度和估计性能折中要求 进行设定;根据相位估计均值划分区域，在每个区域内将多个吉洪诺夫Ti化onov分布近似 为单个吉洪诺夫Ti化onov。
[0150] 结合上述装置，在第一种可能的实现方式中，当采用PSK调制时，所述处理器进一 步包含：
[0151] 相位检测模块，用于估计频偏的取值和相位噪声的分布，进而进行补偿，计算出关 于发送序列的软信息，送给译码模块进行迭代译码;所述软信息为相位检测计算得到的发 送序列的概率分布；
[0152] 逆映射模块，用于PSK解调，连接译码模块和相位检测模块；
[0153] 译码模块，用于进行迭代译码，每次迭代完成将软信息回送给相位检测模块;满足 迭代条件时，输出译码结果;其中，所述软信息为译码计算得到的序列取值的概率分布；
[0154] 迭代判决模块，用于判断是否满足迭代结束的条件，当满足结束条件时则停止迭 代，由译码模块输出译码;否则返回相位检测模块继续上述相位检测模块、逆映射模块和译 码模块的迭代处理，直至满足迭代结束条件。
[0155] 结合上述装置和/或第一种可能的实现方式，在第二种可能的实现方式中，相位检 测模块进一步包含：
[0156] 频率离散化处理子模块，用于将接收的多普勒频偏连续变量0近似为[-妍，妍]上L 个均匀分布的离散值口W，作为初始迭代的L条路径；W及
[0157] 相位估计子模块，用于估计相位噪声的分布；
[0158] 补偿子模块，用于根据得到的相位噪声的分布计算发送序列X的软信息，即所述软 信息为通过相位检测获得的发送序列X的概率分布；
[0159] 输出子模块，用于向逆映射模块输出软信息，进而进行PSK逆映射。
[0160] W上所述仅为本发明的一个【具体实施方式】，并非用于限定本发明的保护范围，本 领域的技术人员应当理解，在不脱离发明原理的前提下，对本发明的技术方案进行修改或 者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的保护范围中。
【主权项】
1. 一种跳频通信中的迭代检测译码方法，所述方法包含： 步骤101)根据发送端的信息比特序列b、发送序列X、接收端的接收序列r、随机相位偏 转9W及表示多普勒频偏的变量资之间的约束关系建立因子图，并依据因子图得到向上，向 下，向前和向后的消息更新公式； 步骤102)采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散化的 频偏值，在步骤101)中建立的因子图上迭代更新相位信息的消息，每次迭代前根据现有信 息重新调整频偏设置;且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译码器， 进而与译码器联合迭代，迭代结束时，比如迭代次数达到设定次数，或译码器相邻两次迭代 输出的外信息变化值小于设定值时，输出译码结果； 其中，将因子图中的相位消息近似为若干个吉洪诺夫Ti化onov分布的和，且构成每个 消息的吉洪诺夫Ti化onov分布的个数根据实际情况中对复杂度和估计性能折中要求进行 设定;根据相位估计均值划分区域，在每个区域内将多个吉洪诺夫Ti化onov分布近似为单 个吉洪诺夫Ti化onov。2. 根据权利要求1所述的跳频通信中的迭代检测译码方法，其特征在于，所述步骤102) 进一步包含： 步骤102-1)将多普勒频偏连续变量0近似为[-挪，觀i上L个均匀分布的离散值口('>，作为 初始迭代的L条路径； 针对每条路径进行如下处理： 步骤102-1-1)将所有表示Ti化onov分布权值大小的W设为1，表示相位Θ分布的参数ζ设 为0,表示频偏路径权重的γ设为1; 初始化因子图上和发送序列X相对应的边X向下的消息和随机相位偏转变量对 应的边Θ向下消息的参数，作为初始迭代的值； 步骤102-1 -2)计算向下消息的权值; 步骤102-1-3)在每条路径上，从左到右计算前向消息； 步骤102-1-4)在每条路径上，从右向左计算反向消息； 步骤10 2 -1 - 5)计算向上传递的消息； 步骤102-2)当达到对应信噪比下设定的迭代次数或者译码器相邻两次输出的外信息 变化值小于一个设定值时，迭代计算结束，迭代译码输出的结果即为恢复后的信息序列； 若不满足上述迭代结束的条件，则计算向下的消息选择另一条路径并返回步 骤102-1-2)。3. 根据权利要求2所述的跳频通信中的迭代检测译码方法，其特征在于，所述步骤102- 1-3)进一步包含:首先，向下传递的消息和向右传递的消息，经过等号结点，相乘得到MQ个 吉洪诺夫Ti化onov分布的和； 其次，用MinKL( ·)将相乘得到的MQ个吉洪诺夫Ti化onov分布变为Q个； 然后，再用如下公式计算如下参量：前向消息中每个吉洪诺夫Ti化onov分布的权值 i、参数瓦料及频偏路径权值完句.其中，角标k表示一跳中第k个符号，角标q表示一条消息中第q个吉洪诺夫Ti化onov分 布;上角标< 表示第1条频偏路径;v^i)4和^|^表示第"k-r个符号的消息经过Min化（·）后 合并得到的消息中各吉洪诺夫Ti化onov分布的权值和参数;参数"表示消息更新后 k.q 3皮句 各吉洪诺夫Ti化onov分布的权值和参数；^，)表示利用前向消息计算得到的该频偏路径的 权值，且越大表示该频偏越越接近实际频偏；<7!表示随机游走相位模型中参数A的方 差;Q为前向和反向相位消息中Ti化onov分布的个数，Μ为PSK的调制阶数。4.根据权利要求2或3所述的跳频通信中的迭代检测译码方法，其特征在于，所述步骤 102-1-4)进一步包含： 首先，向下传递的消息和向左传递的消息经过等号结点相乘得到MQ个吉洪诺夫 Tikhonov分布的和； 其次，用最小化距离近似手段MinKU ·)将相乘得到的MQ个吉洪诺夫Ti化onov分布变 为Q个； 然后，再用如下公式求得如下参数：反向消息中每个吉洪诺夫Ti化onov分布的权值 1?、参数云2 W及频偏路径权值茲.其中，参数:乎'；和^',>表示第"k+1"个符号的消息相乘并经过Min化（·）后合并得到的 .克+l，g ^克 Α?/ 消息中各Ti化onov分布的权值和参数;参数wf嘴…表示消息更新后各Ti化onov分布的权 值和参数;^ >表示利用反向信息后计算得到的该频偏路径的权值。5. 根据权利要求2或3所述的跳频通信中的迭代检测译码方法，其特征在于，所述步骤 102-1-5)进一步包含： 首先，向左传递的消息和向右传递的消息经过等号结点，相乘得到Q2个Ti化onov分布的 和； 其次，在每条路径上用MinKL( ·)将相乘得到的妒个Ti化onov分布变为Μ个； 然后，再用如下公式计算yt(xk):其中，yt(xk)为边化上向上传递的消息，该消息将送给译码器;;和；分 别表示输入同一个等号结点的向右和向左的消息，w'['r和结Γ为经过MinKU .)合并后向上 相位消息中各Ti化onov分布的权值和参数。6. -种跳频通信中的迭代检测译码装置，其特征在于，所述装置包含： 存储器，用于存储相关数据； 处理器，用于执行所述程序，进而进行如下处理： 根据发送端的信息比特序列b、发送序列X、接收端的接收序列r、随机相位偏转ΘΚ及表 示多普勒频偏的变量口之间的约束关系建立因子图，并依据因子图得到向上，向下，向前和 向后的消息更新公式； 采用频偏离散化，即假定频偏为若干个离散化的值，对每个确定的离散化的频偏值，在 建立的因子图上迭代更新相位信息的消息，每次迭代前根据现有信息重新调整频偏设置； 且每次迭代完成后，将相偏和频偏估计结果的软信息传输给译码器，进而与译码器联合迭 代，迭代结束时输出译码结果； 其中，将因子图中的相位消息近似为若干个吉洪诺夫Ti化onov分布的和，且构成每个 消息的吉洪诺夫Ti化onov分布的个数根据实际情况中对复杂度和估计性能折中要求进行 设定;根据相位估计均值划分区域，在每个区域内将多个吉洪诺夫Ti化onov分布近似为单 个吉洪诺夫Ti化onov。7. 根据权利要求6所述的跳频通信中的迭代检测译码装置，其特征在于，所述处理器进 一步包含： 相位检测模块，用于估计频偏的取值和相位噪声的分布，进而进行补偿，计算出关于发 送序列的软信息，送给译码模块进行迭代译码;所述软信息为相位检测计算得到的发送序 列的概率分布； 逆映射模块，用于PSK解调，连接译码模块和相位检测模块； 译码模块，用于进行迭代译码，每次迭代完成将软信息回送给相位检测模块;满足迭代 条件时，输出译码结果;其中，所述软信息为译码计算得到的序列取值的概率分布； 迭代判决模块，用于判断是否满足迭代结束的条件，当满足结束条件时则停止迭代，由 译码模块输出译码；否则返回相位检测模块继续上述相位检测模块、逆映射模块和译码模 块的迭代处理，直至满足迭代结束条件。8.根据权利要求6或7所述的跳频通信中的迭代检测译码装置，其特征在于，所述相位 检测模块进一步包含： 频率离散化处理子模块，用于将接收的多普勒频偏连续变量與近似为[-妍，妍]上L个均 匀分布的离散值，作为初始迭代的L条路径；W及 相位估计子模块，用于估计相位噪声的分布； 补偿子模块，用于根据得到的相位噪声的分布计算发送序列X的软信息，即所述软信息 为通过相位检测获得的发送序列X的概率分布； 输出子模块，用于向逆映射模块输出软信息，进而进行PSK逆映射。
【文档编号】H04L1/00GK105978661SQ201610301037
【公开日】2016年9月28日
【申请日】2016年5月9日
【发明人】倪祖耀, 尹悦, 匡麟玲, 吴胜
【申请人】清华大学