[发明专利]一种针对无频率跳变宽带信号的分段降噪方法

说明书

技术领域

本发明涉及一种信号降噪方法，尤其是一种针对无频率跳变宽带信号的降噪方法。

背景技术

目前，无频率跳变的宽带信号(如线性调频信号，S型和多项式相位型非线性调频信号等)广泛应用于雷达、通讯及声呐等领域，而工程应用中，噪声总不可避免，在某些情况下，瞬时噪声甚至可以完全淹没信号，因此，高效、易于工程实现的降噪处理方法是此类信号处理领域的重要研究内容。

国内外针对无频率跳变宽带信号的降噪方法开展了大量研究，基于统计模型和时频分析等方法不断被提出(如小波变换和分数阶傅里叶分析等)，并且降噪的效果也不断提升。然而，这些方法基本都需要进行时频变换和二维搜索，计算量非常大，对处理系统的运算能力和功耗都带来了巨大挑战。此外，随着现代雷达、通讯及声呐系统向高度集成化方向发展，系统要处理的任务越来越繁重，不可能分配过多的时间和资源用于降噪处理。因此，如何在保证降噪效果的前提下，研究计算量小、易于实时实现的针对无频率跳变信号的降噪方法是当务之急。

发明内容

本发明要解决的技术问题是现有的无频率跳变宽带信号的降噪方法计算量较大，影响信号降噪处理效率。

为了解决上述技术问题，本发明提供了一种针对无频率跳变宽带信号的分段降噪方法，包括如下步骤：

步骤1，对含有噪声的信号序列x(n),n＝0,…,N-1进行分段，形成分段后的信号集为：

x(n)＝{x₁(m),x₂(m),...,x_i(m),…,x_q(m)},i＝1,2,…,q,m＝0,1,…,M-1

式中，q为分段数，M为每段信号的样本数，N为含有噪声的信号序列总的样本数；

步骤2，对每个分段信号x_i(m)分别进行快速傅里叶变换，并分别搜索获得每个分段信号x_i(m)中最大模的谱线，同时确定各个最大模的谱线所对应的索引号k_i0；

步骤3，根据步骤2确定的各个索引号k_i0分别计算频偏δ_i，然后根据各个频偏δ_i将各个分段信号x_i(m)进行频移得到相应的频移信号x_i(m)；

步骤4，对各个频移信号x_i(m)进行快速傅里叶变换，在频域搜索获得各个频移信号x_i(m)种最大模的谱线，保留各个最大模的谱线以及与该最大模的谱线前后相邻的各两条谱线，其他谱线置0，得到各个处理后的频域样本，然后对各个处理后的频域样本做逆快速傅里叶变换，得到相应的重构信号

步骤5，对各个重构信号做逆向频移得到逆频移信号并对各个逆频移信号进行首尾拼接，得到降噪后的原始信号x(n)。

采用将无频率跳变信号进行分段方法，构建分段样本，然后通过计算分段样本的频偏进行频移，使每段信号的中心频率靠近量化频率点，其目的是使信号能量在频域中尽量集中在少数几条谱线，从而有效提高频域降噪的效果；本发明的计算量主要集中在快速傅里叶变换，相比较现有技术中采用时频分析法的二维搜索，本发明的计算量小、易于实现，有效地满足了现代雷达、通信等领域信号处理实时性的要求，因此有效地解决了相应领域信号处理负荷量大和功耗高的问题，具备较高的工程应用价值。

作为本发明的进一步改进方案，在步骤1中，若信号序列总的样本数N不能被每段信号的样本数M整除，则在信号序列的尾部补0，以满足q个分段中均有M个样本；在步骤5中剔除步骤1中末端补零形成的样本。采用在尾部补0能够在N不能被M整除时仍确保q个分段中均有M个样本，方便后期整体计算。

作为本发明的进一步限定方案，步骤3中频偏δ_i的计算方法为：