跳频通信干扰条件下的信噪比估计方法
【技术领域】
[0001 ] 本发明属于跳频自适应传输领域，具体涉及一种跳频通信干扰条件下的信噪比估 计方法。
【背景技术】
[0002] 信噪比是通信系统中的一项重要参数指标，它反映了信号与噪声之间的相对关 系，与系统误码率性能有直接的对应关系，可为数据解调、译码等提供必要的依据，是通信 系统中经常用到的一种信道状态信息。链路自适应技术就是根据信道的状态对传输信号参 数进行自适应的调整，使通信质量达到要求。在慢跳频系统中一种典型的干扰就是部分频 带干扰，由于每一跳频时隙中包含多个符号，在这一情况下的干扰检测可视为对每跳频时 隙内的信噪比估计。在信噪比估计方面前人已经进行了大量研究，这些研究可以分为两大 类：基于数据辅助的方法（DA)和基于非数据辅助的方法（NDA)。数据辅助的方法要求在数 据流中插入接收端已知的导频信号或者训练序列。显然，采用这种方法会降低系统传输效 率。在很多应用场合，从提高频谱效率的角度来考虑，非数据辅助方法是一种更好的选择。
[0003] 由于非数据辅助的信噪比估计技术在节省系统资源上显现出的巨大优势，让其在 众多学者中得到了更多的青睐。目前常用的NDA信噪比估计算法有一阶二阶矩和二阶四 阶矩（M2M4)估计算法。其中一阶和二阶来估计信噪比的方法比M2M4算法精度高，在低信 噪比时估计偏差更小。但是实际应用中受查找表大小限制，高信噪比时估计性能会随着 信噪比增加而下降。文献"Non-data-aid SNR estimation method for APSK exploting rank discrimination test"采用一致最大功效检验的方法针对APSK信号的NDA方法进 行了研究，在高信噪比时其性能要优于M2M4估计算法。文献"Signal-to-noise ratio estimation algorithm for adaptive coding and modulation in advanced digital video broadcasting - radar cross sectionsatellite systems" 针对 MPSK 调制系统提 出了一种短突发SNR估计算法，高信噪比时性能优于M2M4估计算法，但其需要准确的相位 信息，对相位估计提出了较高的要求。以上是针对高信噪比时的信噪比估计进行的研究，且 未与编码系统相结合。针对低信噪比估计文献"Non-data-aided signal-to-noise-ratio estimation"提出了一种迭代似然函数期望最大化NDA算法。文献"How to use a priori information of data symbols for SNR estimation" 提出了如何将数据符号 的先验信息用于 SNR 估计。文献"Robust frequency hopping for interference and fading channels"在解调器和译码器之间增加一个信道估计器，它采用期望最大（EM - Expectation maximization)算法，利用每次译码反馈的后验信息得到信道估计值。这种 方法具有良好的性能，在每跳包含一定数量符号的条件下与已知信道边信息的情况十分接 近，但是由于EM算法的计算量很大，这种方法的复杂度较高，不利于实际应用。
[0004] "Robust frequency hopping for interference and fading channels，'（IEEE Transactions on Communications, vol. 56, no. 8, pp. 1343-1351，2008.)在解调器和译码 器之间增加一个信道估计器，它采用期望最大（EM - Expectation maximization)算法，利 用每次译码反馈的后验信息得到信道估计值。这种方法具有良好的性能，在每跳包含一定 数量符号的条件下与已知信道边信息的情况十分接近，但是由于EM算法的计算量很大，这 种方法的复杂度较高，不利于实际应用。
[0005] 综上所述，传统的信噪比估计方法估计精度低、信噪比估计范围窄、在低信噪比情 况下估计偏差很大，在信号长度很短且信噪比很低的情况下，很多信噪比估计方法都会出 现较大估计偏差。
【发明内容】
[0006] 本发明的目的在于提供一种跳频通信干扰条件下的信噪比估计方法，用于解决在 估计精度高、信噪比估计范围广且在低信噪比情况下估计偏差很小条件下的信噪比估计。
[0007] 实现本发明目的的技术解决方案为：一种跳频通信干扰条件下的信噪比估计方 法，方法步骤如下：
[0008] 步骤1、在跳频通信干扰条件下，将接收到的信号通过匹配滤波器进行匹配滤波， 再将其传入能量检测器进行能量检测，并预判出干扰跳，再将其解调；
[0009] 步骤2、将解调后的信号分别送入低密度奇偶校验码（LDPC)迭代译码器和信噪比 估计器，经LDPC迭代译码器迭代译码后，将信号的迭代译码信息送入信噪比估计器，结合 解调后的信号信息，在信噪比估计器中通过估计算法进行性噪比估计，将估计值送入LDPC 迭代译码器，进行下一次迭代译码，然后将迭代译码信息送入信噪比估计器，再结合解调后 的信号信息，通过估计算法进行性噪比估计，再将估计值送入LDPC迭代译码器，重复上述 过程，直至达到迭代次数，并将最终估计值输出。
[0010] 上述步骤2中，估计算法为通过最大似然算法估计出信号幅度A和噪声方差2 〇2， 算法步骤如下：
[0011] 匹配滤波后信号yn的函数为
[0012] yn = Aang〇+wn (1)
[0013] 信号yn的对数似然函数111?(7/^，〇2)为
[0015] 其中，Μ为调制阶数，K为估计符号长度，在慢跳频中每跳的符号数应能被K整除， Pr(s(1))为第i个发送符号序列s(1)的先验概率，i = 0，1，2...11<-1，&是抽样脉冲的峰值， ak是发送符号，wk表示经匹配滤波采样后的噪声，p表示概率密度函数；
[0016] 分别对A和〇 2求偏导
[0019]条件似然函数 In p(y/s(1)，Α，σ 2)为
[0021] 将式（6)分别代入式（4)和式（5)得
[0024] 其中nk为译码信息期望值，％为nk的共轭，p k为后验均方值；
[0025] 令式（7)和式⑶分别等于零，即可得到A和0 2的估计值，并将所得到的估计值 代入基于最大似然的非数据辅助（ML-NDA)估计式，得到估计值f
[0027] 再乘以因子（Κ_3)/Κ以降低估计偏差，得到修正后的估计值，并将其输出
[0029] 上述译码信息期望值η k取值如下：
[0030] 对于BPSK调制信号，译码信息期望值
是迭代译码器 迭代计算1次后输出的对数域译码信息；
[0031] 对于QPSK调制信号，1为
[0034] 上述[ak(1)，bk(1)]是映射到第k个QPSK调制符号的两个比特，L(a k(1))和L(bk(1)) 是
[0035] 映射到第k个QPSK调制符号的两个比特对应的迭代译码器迭代计算1次后输出 的对数域译码信息。
[0036] 步骤2中的迭代次数范围为[8,40]。
[0037] 本发明与现有技术相比，其显著优点：（1)本发明解决了在估计精度高、信噪比估 计范围广且在低信噪比情况下估计偏差很小条件下的信噪比估计；（2)通过分析估计算法 的的偏差特性，当采用SNV-SRT算法时，虽然用于估计的数据个数仅有64个，但在信噪比低 至-2dB依然可以保持很小的估计偏差；(3)在慢调频通信中，每一跳的数据个数往往较少， 因此这种修正方法特别适用于慢调频通信，对于提高慢跳频通信自适应传输时信道状态估 计准确性有积极作用。