一种色噪声中的载波频率获取方法技术

本发明专利技术公开了一种色噪声中的载波频率获取方法，涉及通信技术领域，解决N次方谱算法的效果差的技术问题，所述方法包括：对信号进行时域采样，截取总长为N的信号；估计所述信号的SNR；对所述信号做EMD分解，得到若干个低阶的IMF分量；根据所述SNR，选取对应的IMF分量，对该IMF分量做FFT变换，得到所述信号的时频谱图F

【技术实现步骤摘要】
一种色噪声中的载波频率获取方法


[0001]本专利技术涉及通信
，更具体地说，它涉及一种色噪声中的载波频率获取方法。

技术介绍

[0002]信号的载波频率(Symbol Rate frequency,SRF)是通信系统中的一个非常重要参数，在非协作通信中，SRF估计就是根据接收信号的过采样，通过非数据辅助的方法，得到SRF的估计值。
[0003]考虑幅相类调制信号MQAM/MPASK，信息由信号的幅度与相位承载，则接收信号可以用式(1)表示。
[0004][0005]其中，n(t)为0均值且方差已知的加性高斯白噪声，A
m
与θ
m
分别表示调制信号的第M个符号所对应的幅度与相位，g(t)为发送的低通波形脉冲，f
c
为载波频率。SRF估计的任务就是依据接收信号r(t)，在事先对信号没有先验信息的条件下，对载波频率f
c
做出精确的估计。假设发送信号为单极性脉冲，则其信号的功率谱如式(2)所示。
[0006][0007]其中P表示发送传号的概率，G(ω)为g(t)的幅度谱。式(2)的第一项为连续谱，包含通信的信息，第二项为离散谱，通常包含载波频率的信息，即在f
b
的整数倍上可能存在冲激分量。若存在离散分量，则可以直接对SRF做出估计。
[0008]然而，遗憾的是幅相类调制的接收信号r(t)没有离散谱，只有连续谱，即不能直接通过对接收信号做频谱分析，获得对SRF的估计。需要对接收信号做非线性变换，使之产生离散谱。通常的做法是，根据信号的二次平稳特性，对接收信号做N次方运算，即r
N
(t)，再做频谱分析。此时信号有可能存在与SRF相关的离散谱，然后再做频谱分析，则在SRF的N倍处会出现冲激信号，即幅度有一个明显的峰值，如图1所示。根据这个峰值，可以得到N倍SRF的估计，进而得到SRF的估计。
[0009]上述算法在白噪声中可以获得SRF的较为精确的估计，然而，由于通信信号在发送时，受低通波形脉冲g(t)的影响，在做类似于N次方的非线性变换时，不仅可以产生SRF，同时还会产生大量的低频成分，即背景噪声并不是白噪声，而是色噪声，这在通信系统中是普遍存在的。这使得SRF所对应的频点，往往并不是全局的最大值，而仅仅是个局部最大值，如图2所示。这使得对SRF的估计，不能简单地通过搜索频域极值点，来获得对SRF的估计，而要通过复杂的算法研究SRF对应频点满足的条件。不仅算法复杂，且经常发生误判，这对SRF的估计提出了严峻的挑战。而且在SNR(Signal Noise Ratio信噪比)较低时，算法经常失效，这是SRF估计面临的又一个严峻挑战。

技术实现思路

[0010]本专利技术要解决的技术问题是针对现有技术的上述不足，本专利技术的目的是提供一种可以有效获取载波频率的色噪声中的载波频率获取方法。
[0011]本专利技术的技术方案是：一种色噪声中的载波频率获取方法，包括：
[0012]步骤1.对信号进行时域采样，截取总长为N的信号；
[0013]步骤2.估计所述信号的SNR；
[0014]步骤3.对所述信号做EMD分解，得到若干个低阶的IMF分量；
[0015]步骤4.根据所述SNR，选取对应的IMF分量，对该IMF分量做FFT变换，得到所述信号的时频谱图F
m
(ω)；
[0016]步骤5.搜索所述时频谱图F
m
(ω)的最大值，所述最大值点对应的频率即为SRF的估计。
[0017]作为进一步地改进，在步骤3中，所述信号表示为X(t)，对所述信号X(t)做EMD进行K次分解后，信号X(t)分解为K个IMF
k
分量与余量r
(K)
的和，
[0018][0019]其中，每个IMF
k
分量占据着由高到低的不同的频域空间。
[0020]进一步地，具体分解过程如下：
[0021]步骤31.找到原信号X(t)的所有局部极大值点，然后用三次样条插值函数拟合得到原信号的上包络线，找到原信号X(t)的所有局部极小值点，然后用三次样条插值函数拟合得到原信号的下包络线，上包络线和下包络线的均值用m
l
(t)表示；
[0022]步骤32.原信号序列X(t)与该平均m
l
(t)作差，得到新信号序列，
[0023]h1(t)＝X(t)
‑
m1(t)；
[0024]步骤33.若新信号序列h1(t)不存在负的局部极大值且不存在正的局部极小值，则令IMF1＝h1(t)；否则，将h1(t)作为新的X(t)，重复步骤31～步骤33；
[0025]步骤34.令r
(1)
＝X(t)
‑
IMF1；将r
(1)
作为新的X(t)，重复步骤31～步骤34，依次得到IMF2……
IMF
k
，直到r
(K)
不可以被分解。
[0026]进一步地，所述SRF的估计为：
[0027][0028]进一步地，所述估计的性能可以通过归一化均方根误差来表征：
[0029][0030]其中f
e
为f
c
的估计。
[0031]有益效果
[0032]本专利技术与现有技术相比，具有的优点为：
[0033]本专利技术通过低阶的EMD分解，获得包含SRF的IMF分量，然后通过频谱分析可得对SRF的估计，把复杂的时域信号分解为若干个在频域上体现出极高的规律性的分量，从频域地面来看观察结果，所得分量排列规律，几乎没有能量遗漏，可自动适应信号的变化，本发
明不仅可以在色噪声中取得良好的估计效果，且在低信噪比条件，算法仍然有效。
附图说明
[0034]图1为白噪声中的通信信号频谱图(22M)；
[0035]图2为色噪声中的SRF不能保证能量最大；
[0036]图3为本专利技术的流程图；
[0037]图4为SNR＝
‑
10dB时EMD的时域与频域性能图；
[0038]图5为SNR＝0dB时EMD的时域与频域性能图；
[0039]图6为SNR＝6dB时EMD的时域与频域性能图；
[0040]图7为SRF估计性能图。
具体实施方式
[0041]下面结合附图中的具体实施例对本专利技术做进一步的说明。
[0042]参阅图3～7，一种色噪声中的载波频率获取方法，包括：
[0043]步骤1.对信号进行时域采样，截取总长为N的信号；
[0044]步骤2.估计信号的SNR；
[0045]步骤3.对信号做EMD(Empirical Mode Decomposition基于经验模态分解)分解，得到若干个低阶的IMF(Intrinsic Mode Function本征模函数)分量；...

【技术保护点】

【技术特征摘要】
1.一种色噪声中的载波频率获取方法，其特征在于，包括：步骤1.对信号进行时域采样，截取总长为N的信号；步骤2.估计所述信号的SNR；步骤3.对所述信号做EMD分解，得到若干个低阶的IMF分量；步骤4.根据所述SNR，选取对应的IMF分量，对该IMF分量做FFT变换，得到所述信号的时频谱图F
m
(ω)；步骤5.搜索所述时频谱图F
m
(ω)的最大值，所述最大值点对应的频率即为SRF的估计。2.根据权利要求1所述的一种色噪声中的载波频率获取方法，其特征在于，在步骤3中，所述信号表示为X(t)，对所述信号X(t)做EMD进行K次分解后，信号X(t)分解为K个IMF
k
分量与余量r
(K)
的和，其中，每个IMF
k
分量占据着由高到低的不同的频域空间。3.根据权利要求2所述的一种色噪声中的载波频率获取方法，其特征在于，具体分解过程如下：步骤31.找到原信号X(t)的所有局部极大值点，然后用三次样条插值函数拟合得到原信号的上包络线，找到原信号X(t)的所有局部极小值点，然后用三次样条插值...

【专利技术属性】
技术研发人员：王勤民，汪洋，刘世刚，崔世友，魏淑稳，赵珂，张灿，黄宁，
申请(专利权)人：河南炬讯信息技术有限公司，
类型：发明
国别省市：