一种基于循环相关熵的载频估计方法
【技术领域】
[0001] 本发明涉及循环平稳信号处理技术领域，涉及到调制信号载频估计方法，特别涉 及到一种基于循环相关熵的载频估计方法。
【背景技术】
[0002] 传统的载波频率估计方法主要是基于平稳信号的频率谱或者功率谱来实现的，但 是当存在同频带干扰时，这类方法就无法准确地实现载频估计的功能。为了解决该问题，研 究工作者提出了基于循环相关谱的载频估计方法，首次部分解决了同频带干扰条件下的载 频估计问题，但是该类方法由于其严重依赖高斯噪声的假设而存在缺陷，当噪声模型不服 从高斯分布时，此类方法的性能严重退化。因此，本发明首先提出了循环相关熵的概念，然 后提出了一种基于循环相关熵谱的载频估计方法，为解决非高斯噪声与同频带干扰并存条 件下的非平稳信号的载频估计提出了切实可行的方法。
【发明内容】
[0003] 本发明的主要目的是为了解决现有技术问题，提供一种基于循环相关熵谱的具有 抑制非高斯噪声和同频带干扰的非平稳信号载频估计方法。
[0004] 本发明是基于循环相关熵的载频估计方法，首先利用公式计算信号的相关熵和对 应的循环相关熵，然后利用公式计算循环相关熵谱，并获得循环相关熵投影，最后根据循环 相关熵谱投影峰值所对应的循环频率与载频的关系，估计信号的载波频率，技术方案包括 以下具体步骤：
[0005] 第一步，计算待估计信号的循环相关熵。
[0006] ⑴由公式 vx(t，τ) =E[K G(x(t)-x(t+T))]计算信号的相关熵 Vx(t，τ)，其 中X(t)是信号，τ表示信号的时延，E表示数学期望操作符，K。表示高斯核函数，它满足
[0007] (2)按照公另
计算信号的 循环相关熵巧id，循环相关熵是相关熵的傅里叶级数的系数，其中<Τ>表示被积分的区间 长度为Τ，ξ表示循环频率，〈· >t表示求时间平均。
[0008] 第二步，根据该信号的循环相关熵计算循环相关熵谱在循环频率域的投影。
[0009] (1)根据公式
计算信号的循环相关熵谱允?/)，循环相 关熵谱是循环相关熵的傅里叶变换。
[0010] (2)将循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影。
[0011] 第三步.根据循环相关熵谱投影估计载波频率。
[0012] (1)根据循环相关熵谱投影中的峰值，获得对应的循环频率ξ。
[0013] (2)根据循环频率和载频的关系ξ =af。，估计信号的载波频率f。，其中a为整数， 由不同的调制方式决定。
[0014] 本发明首先提出了循环相关熵以及循环相关熵谱等一系列新的概念，解决了脉冲 噪声和同频带干扰并存条件下的非平稳信号的载波频率估计。
【附图说明】
[0015] 附图1是本发明中的载频估计方法的总体流程图。
[0016] 附图2是本发明中以BPSK信号为例的循环相关熵谱。
[0017] 附图3是本发明中以BPSK信号为例的循环相关熵谱投影。
【具体实施方式】
[0018] 为使本发明实施例的目的、技术方案及其优点更加清楚，下面结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚完整的描述，整体算法流程图如图1所 示：
[0019] 第一步，计算待估计信号的循环相关熵。
[0020] ⑴按照公式Vx(t，τ) =E[k G(x(t)-x(t+T))]计算信号的相关熵，其中：
[0021] x(t)是信号，τ表示信号的时延，E表示数学期望操作符；
[0022] κ。表示高斯核函数，
[0023] ⑵按照公式
计算信号 的循环相关熵，其中<Τ>表示被积分的区间长度为Τ，ξ表示循环频率，〈· >t表示求时间 平均。
[0024] 第二步，根据该信号的循环相关熵计算在信号的循环相关熵谱投影。
[0025] (1)根据公式
计算信号的循环相关熵谱，它是循环相关 熵的傅里叶变换。如图2所示，是以BPSK信号为例的循环相关熵谱的示意图。
[0026] (2)将循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影。如图3所 示，是将图2中BPSK信号的循环相关熵谱投影至循环频率域而获得的讯相关熵谱投影
[0027] 第三步.根据循环相关熵谱投影估计载波频率。
[0028] (1)根据循环相关熵谱投影中的峰值，获得对应的循环频率ξ。由图3可见除了 代表直流的中心分两外，在ξ = ±2000Ηζ处有脉冲成分，而在ξ = ±4000Ηζ处则是其谐 波成分。
[0029] (2)根据循环频率和载频之间的公式ξ = af。，估计信号的载波频率f。，其中a为 整数，不同调制方式a取值不同。针对BPSK信号，a = ±2,推导过程如下：由于BPSK信号 的循环相关熵谱中对应二阶矩的成分可以写成
[0030]
[0031] 其中，表示2-范数操作符，δ ( ·)是单位冲激函数，Sa(f)表示BPSK信号包络 a(t)的频率谱，*表示共辄操作符。由公式可见，在ξ = ±2f。处有脉冲成分，即说明载频 为fc的BPSK信号在循环频率ξ = ±2fc处存在可以用来估计载频的成分。
【主权项】
1. 一种基于循环相关熵的载频估计方法，其特征在于，包括以下具体步骤： 第一步，计算待估计信号的循环相关熵 (1) 由公式vx(t，τ) =E[^(x(t)-x(t+T))]计算信号的相关熵Vx(t，τ)，其 中，x(t)是信号，τ表示信号的时延，Ε表示数学期望操作符；κ。表示高斯核函数，计算信号的循 环相关熵其中<Τ>表示被积分的区间长度为Τ，ξ表示循环频率，〈*>t表示求时间 平均； 第二步，根据该信号的循环相关熵计算循环相关熵谱在循环频率域的投影计算信号的循环相关熵谱; (2) 将循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影； 第三步，根据循环相关熵谱投影估计载波频率； (1) 根据循环相关熵谱投影中的峰值，获得对应的循环频率ξ; (2) 根据循环频率和载频之间的公式ξ=af。，估计信号的载波频率f。，其中a为整数， 不同的调制方式a取值不同。
【专利摘要】本发明属于循环平稳信号处理技术领域，提供了基于循环相关熵的载频估计方法。本发明首先通过公式计算信号的相关熵和循环相关熵；然后，通过公式计算信号的循环相关熵谱，并通过投影获得信号的循环相关熵谱投影；最后，根据调制方式和载波频率与循环频率之间的关系，通过循环相关熵谱投影峰值所对应的估计信号的载波频率。实验证明本发明算法性能良好，能够在非高斯噪声和同频带干扰并存的条件下对信号进行载频估计。
【IPC分类】G06F17/50
【公开号】CN105302940
【申请号】CN201510591312
【发明人】邱天爽, 栾声扬, 王鹏, 李景春, 谭海峰, 鲍尧, 杨文翰, 朱永杰, 戚寅哲, 马济通, 史益新
【申请人】大连理工大学
【公开日】2016年2月3日
【申请日】2015年9月16日