一种基于循环相关熵的载频估计方法

摘要

本发明属于循环平稳信号处理技术领域，提供了基于循环相关熵的载频估计方法。本发明首先通过公式计算信号的相关熵和循环相关熵；然后，通过公式计算信号的循环相关熵谱，并通过投影获得信号的循环相关熵谱投影；最后，根据调制方式和载波频率与循环频率之间的关系，通过循环相关熵谱投影峰值所对应的估计信号的载波频率。实验证明本发明算法性能良好，能够在非高斯噪声和同频带干扰并存的条件下对信号进行载频估计。

主权项

一种基于循环相关熵的载频估计方法，其特征在于，包括以下具体步骤：第一步，计算待估计信号的循环相关熵(1)由公式V_x(t,τ)＝E[κ_σ(x(t)‑x(t+τ))]计算信号的相关熵V_x(t,τ)，其中，x(t)是信号，τ表示信号的时延，E表示数学期望操作符；к_σ表示高斯核函数，$K_{\sigma }(·)=\frac{1}{\sqrt{2\pi }\sigma }\exp (-{(·)}^2/(2{\sigma }^2\left)\right);$(2)按照公式$V_x^{\xi }\left(\tau \right)=\frac{1}{T}{\int }_{<T>}{V}_x(t,\tau )e^{-j2\pi \xi t}dt=<{\kappa }_{\sigma }\left(x\right(x)-x(t+\tau \left)\right)e^{-j2\pi \xi t}{>}_t$计算信号的循环相关熵其中<T>表示被积分的区间长度为T，ξ表示循环频率，<·>_t表示求时间平均；第二步，根据该信号的循环相关熵计算循环相关熵谱在循环频率域的投影(1)根据公式计算信号的循环相关熵谱(2)将循环相关熵谱投影至循环频率域，获得信号的循环相关熵谱投影；第三步，根据循环相关熵谱投影估计载波频率；(1)根据循环相关熵谱投影中的峰值，获得对应的循环频率ξ；(2)根据循环频率和载频之间的公式ξ＝af_c，估计信号的载波频率f_c，其中a为整数，不同的调制方式a取值不同。