一种基于判决统计量设计的信号检测方法及系统

技术领域

本发明涉及声纳信号处理领域，特别涉及一种基于判决统计量设计的信号检测方法及系统。

背景技术

信号检测在雷达、声纳和语音等领域中占有重要的地位。因此，如何在复杂背景中准确检测感兴趣的信号一直是声纳、雷达和语音等领域中信号处理所关心的研究内容之一。目前，比较经典和常用的信号检测方法有互相关检测法、平方律检波器法、能量累积检测法等。其中，能量累积检测法由于其复杂度低，被广泛应用。但是该类检测法需要一定的信噪比，在信噪比较低的情况下，对未知信号检测效果较差，需要采取一定的滤波手段对数据进行增强处理，以提高信噪比。

为了对处理数据中期望信号实现增强，提升信号检测效果，研究者提出了一些信号的增强方法，并取得了一定的应用效果，如参数方法、非参数方法以及统计方法。在参数方法中主要代表性方法为：梳状滤波器、维纳滤波和卡尔曼滤波方法等；非参数方法主要代表性方法为：自适应抵消和谱减法等；在统计方法中主要代表性方法为：映射方法，参数估计和隐马尔科夫模型等。但这些方法在实际应用中还存在一些问题，如参数方法对所使用的模型依赖性强，需要提取模型参数；非参数方法不需要从传感器拾取数据中估计模型参数，由于没有利用潜在的统计信息，不适合处理含非平稳背景噪声数据；统计方法充分利用了统计特点，以概率计算方法实现对信号的增强，但需要依赖一定准则实现对最优结构的选择。

综上，如何应对低信噪比情况下、对信号进行检测仍然是富有挑战性的难题。

发明内容

本发明的目的在于解决上述问题，提出一种基于判决统计量设计的信号检测方法及系统。

为了实现上述目的，本发明提出了一种基于判决统计量设计的信号检测方法，所述方法包括：

采用伽马通滤波器组对传感器拾取数据进行分频带处理得到各频带数据；

根据信号与背景噪声在强度和稳定性上的差异性，对各频带数据进行方差归一化处理，建立判决统计量；

根据判决统计量对各频带数据进行判决处理，实现对信号的增强和检测。

作为上述方法的一种改进，所述方法具体包括：

步骤1)采用伽马通滤波器组对传感器拾取数据x(n)进行处理，得到I组频带数据，第i组频带数据y

y

其中，n为采样点，n＝1,2,…,N，N为一次处理采样长度，g

g

其中，A为伽马通滤波器组增益，L为伽马通滤波器组的阶数，L＝4，f

其中，m表示时域信号；则μ(n)＝1；

b

b

其中，ERB(f

步骤2)根据信号与背景噪声在强度和稳定性上的差异性，对每组频带数据计算强度归一化方差，消除不同滤波器等效矩形带宽对判决统计量的影响；第i组强度归一化方差Dδ

其中，

步骤3)建立各频带的判决统计量，第i组频带的判决统计量T

步骤4)根据各频带判决统计量，对各频带数据进行判决加权处理，并进行合成，得到增强信号z

步骤5)采用能量累积检测方法对z

一种基于判决统计量设计的信号检测系统，所述系统包括伽马通滤波器组、归一化及判决统计模块和信号增强检测模块；其中，

所述伽马通滤波器组，用于对传感器拾取数据进行分频带处理得到各频带数据；

所述归一化及判决统计模块，用于根据信号与背景噪声在强度和稳定性上的差异性，对各频带数据进行方差归一化处理，建立判决统计量；

所述信号增强检测模块，用于根据判决统计量对各频带数据进行判决处理，实现对信号增强和检测。

作为上述系统的一种改进，所述伽马通滤波器组的具体处理过程包括：

对传感器拾取数据x(n)进行处理，得到I组频带数据，第i组频带数据y

y

其中，n为采样点，n＝1,2,…,N，N为一次处理采样长度，g

g

其中，A为伽马通滤波器组增益，L为伽马通滤波器组的阶数，L＝4，f

其中，m表示时域信号；则μ(n)＝1；

b

b

其中，ERB(f

作为上述系统的一种改进，所述归一化及判决统计模块的具体处理过程包括：

根据信号与背景噪声在强度和稳定性上的差异性，对每组频带数据计算强度归一化方差，消除不同滤波器等效矩形带宽对判决统计量的影响；第i组强度归一化方差Dδ

其中，

建立各频带的判决统计量，第i组频带的判决统计量T

作为上述系统的一种改进，所述信号增强检测模块的具体处理过程包括：

根据各频带判决统计量，对各频带数据进行判决加权处理，并进行合成，得到增强信号z

采用能量累积检测方法对z

与现有技术相比，本发明的优势在于：

利用信号与背景噪声在强度和稳定性上的差异特性，通过对各频带数据进行方差归一化处理，建立判决统计量，对各频带数据进行判决处理，在无需搜索单个滤波频带的情况下，提高了合成数据信噪比，实现对信号增强和检测。

附图说明

图1为本发明实施例1的基于判决统计量设计的信号检测方法流程图；

图2(a)为本发明实施例1的伽马通波器组第12个滤波器时域波形；

图2(b)为本发明实施例1的伽马通波器组各滤波器频响；

图2(c)为仿真实例所造数据时域波形；

图2(d)为仿真实例所造数据频域波形；

图3(a)为仿真实例采用FIR滤波后数据时域波形；

图3(b)为仿真实例采用FIR滤波后数据频域波形；

图4(a)为仿真实例采用本发明的方法处理后的时域波形；

图4(b)为仿真实例采用本发明的方法处理后的频域波形；

图5(a)为仿真实例本发明所得归一化方差；

图5(b)为仿真实例本发明所得判决统计量；

图6(a)为仿真实例未采用本发明方法处理前信号检测结果的时域波形；

图6(b)为仿真实例未采用本发明方法处理前信号检测结果的时域包络；

图7(a)为仿真实例采用本发明方法处理后信号检测结果的时域波形；

图7(b)为仿真实例采用本发明方法处理后信号检测结果的时域包络。

具体实施方式

本发明提出一种基于判决统计量设计的信号检测方法，该方法以人耳耳蜗时频分解原理为基础，设计伽马通滤波器组，对传感器拾取数据进行分析处理，并根据信号与背景噪声在强度和稳定性上的差异，建立判决统计量，对各频带数据实现判决加权处理，在无需搜索单个滤波频带的情况下，提高了合成数据信噪比，实现了对信号的增强和检测，进而提升了能量累积检测方法对信号的检测效果。最后通过数值仿真进一步验证了本发明所述方法的有效性和可行性。

下面结合附图和实施例对本发明的技术方案进行详细的说明。

实施例1

如图1所示，本发明的实施例1提供了一种基于判决统计量设计的信号检测方法。

在对本发明的方法做详细说明前，首先对本发明的方法所述信号检测流程图加以描述。传感器拾取数据首先通过已设计的滤波器组进行分解，将其分解成与耳蜗非线性频率相对应的子带数据；然后采用强度和稳定性信息形成判决统计量提取组成因子，构建判决统计量；最后采用非线性判决加权方法对信号增强。

下面对本发明方法做进一步的说明。

伽马通(Gammatone)滤波器滤波器组基本原理

在人耳听觉系统中，不同频率的信号可在人耳耳蜗基底膜的不同位置产生最大响应。对此，可依据人耳耳蜗基底膜频响特性，设计相应滤波器组模拟耳蜗基底膜频响特性，对传感器拾取数据进行分解处理，增强未知多频带组合信号所在频带在整个处理频带中的比重。

在模拟人耳耳蜗基底膜时频分解中，本发明采用一组相互交叠的伽马通滤波器组模拟实现，滤波器组的时域响应可表示为：

式中，n为采样点，A为滤波器组增益，L为滤波器组的阶数，b

衰减因子b

b

式中，ERB(f

对于采样率为f

判决处理

首先，采用Gammatone滤波器组对传感器拾取数据x(n)(n＝1,2,…,N，N为一次处理采样长度，)进行处理，可得I组分解后数据y

y

式中，g

其次对每个频带数据y

式中，

再次建立判决统计量。

最后根据提到的各频带判决统计量T

步骤5)可采用能量累积检测方法对z

仿真实例分析

下面结合实例，对本发明方法的效果进行比较。

为了进一步验证本发明方法可以有效实现对信号检测，进行如下数值仿真分析，见下表：

表1数值仿真参数

仿真中，假定数据处理前对脉冲信号频带、脉宽长度均未知，所以在处理数据前只采用了带宽为100Hz～8kHz的128阶FIR滤波器对处理数据进行初始滤波处理。

图3至图4给出了采用本发明方法处理前后的数据对比，图5为本次处理中所得的判决统计量，所用滤波器组即为图2(a)-(d)所示Gammatone滤波器组。

数值仿真中采用图2所示Gammatone滤波器组。

图3(a)、(b)为FIR滤波后数据时频波形；

图4(a)、(b)为本发明方法处理后时频波形；

图5(a)、(b)为本发明所得判决统计量；

图6(a)、(b)为未采用本发明方法处理前信号检测结果；

图7(a)、(b)为采用本发明方法处理后信号检测结果。

由图6和图7结果可知，采用能量累积检测法无法对原始数据进行信号检测，无法对第2000个采样点附近未知脉冲信号实现检测；而经本发明方法处理后，再采用能量累积检测法对其进行检测，可有效实现对第2000个采样点和第4000个采样点附近未知的脉冲信号检测。该结果进一步证实了本发明方法可提升信号检测方法对未知多频带组合信号的检测效果。

实施例2

本发明的实施例2提出了一种基于判决统计量设计的信号检测系统。该系统包括：一种基于判决统计量设计的信号检测系统，所述系统包括伽马通滤波器组、归一化及判决统计模块和信号增强检测模块；具体处理方法同实施例1。

所述伽马通滤波器组，用于对传感器拾取数据进行分频带处理得到各频带数据；

所述归一化及判决统计模块，用于根据信号与背景噪声在强度和稳定性上的差异性，对各频带数据进行方差归一化处理，建立判决统计量；

所述信号增强检测模块，用于根据判决统计量对各频带数据进行判决处理，实现对信号增强和检测。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。