CN105067025A

CN105067025A - 一种利用单稳系统随机共振效应检测微弱信号的方法

Info

Publication number: CN105067025A
Application number: CN201510458769.2A
Authority: CN
Inventors: 郭锋; 朱宝龙
Original assignee: Southwest University of Science and Technology
Current assignee: Southwest University of Science and Technology
Priority date: 2015-07-31
Filing date: 2015-07-31
Publication date: 2015-11-18

Abstract

本发明提供了一种利用单稳系统随机共振效应检测微弱信号的方法。该方法基于随机共振理论，利用噪声，而非抑制噪声，通过多次调节一个系统参数和注入系统的噪声强度，使系统输出信噪比最佳，可以最大化地提高系统输出信噪比，从而实现对微弱信号的最佳检测。本发明特别适合机械故障、电子故障检测等低信噪比、低频微弱特征信号检测。

Description

一种利用单稳系统随机共振效应检测微弱信号的方法

技术领域

本发明属于低信噪比条件下微弱信号的检测方法，特别涉及一种利用单稳系统随机共振效应检测微弱信号的方法。

背景技术

微弱信号检测的首要任务是提高信噪比，以便从强噪声中检测出有用的微弱信号。由于微弱信号的检测能提高测量灵敏度和可检测下限，因此在物理、化学、生物以及许多工程技术领域都得到了广泛应用。常采用的微弱信号检测方法都是抑制噪声来提取微弱信号，在强噪声背景和信号的频率非常低的情况下，显得无能为力。

当系统的非线性、输入信号与噪声三者之间存在某种匹配时，系统输出信噪比最优，这种非线性现象即随机共振。随机共振的概念是由意大利物理学家RobertoBenzi、美国物理学家AlfonsoSutera和意大利物理学家AngeloVulpoiani等人1981年在研究古气象冰川问题时提出的。当有噪声的系统发生随机共振时，部分噪声能量会转化为有用信号的能量，从而使系统输出信噪比大大提高,即给特定系统加入一定强度的噪声，不但不会阻碍反而会提高信号检测的性能。

目前基于随机共振检测信号模型通常为双稳系统或分段线性系统模型。其中，

(1)双稳模型满足方程

其中为输入信号，为噪声，,,势能函数为

(2)分段线性模型满足方程

应用双稳模型检测信号时模型参数与系统特性有关，改变系统参数则系统特性会随着改变，需要调节系统a、b两个参数；而分段线性系统需要调节a、b、c三个参数，构建模型复杂，不利于工程应用。

发明内容

为解决上述问题，本发明提供了一种利用单稳系统随机共振效应检测微弱信号的方法，所述方法包括以下步骤：

步骤一：初始化单稳系统中的b、；在该系统中，为系统输出信号，为系统结构参数，是加性高斯白噪声，s(t)为周期信号，,是注入系统的乘性高斯分布白噪声。

步骤二：选择固定系统参数b或者乘性高斯分布白噪声中的一个，调节另一个，直到系统输出信号的信噪比达到最大；

步骤三：固定上一步中被调节的对象最后的取值，调节另一个，直到系统输出信号的信噪比达到最大。

步骤四：返回执行步骤三，直到系统输出信噪比不再增大为止。

进一步的，系统输出信号进行傅里叶变换(FFT)后计算信噪比。

进一步的，系统结构参数调节范围为（0,5）。

进一步的，乘性高斯分布白噪声强度D的调节范围为（0,1）。

本发明的有益效果为：

本发明针对现有双稳随机共振和分段线性随机共振检测微弱信号需要调节参数较多，不利于工程适用的不足，提供一种基于单稳态随机共振的微弱周期信号检测方法，该方法基于随机共振理论，利用噪声，而非抑制噪声，通过多次调节一个系统参数和注入系统的噪声强度，使系统输出信噪比最佳，可以最大化地提高系统输出信噪比，从而实现对微弱信号的最佳检测。本发明特别适合机械故障、电子故障检测等低信噪比、低频微弱特征信号检测。

附图说明

图1为本发明单稳系统结构框图。

图2为本发明方法的微弱低频周期信号检测流程图。

图3为仿真信号时域波形图。

图4为仿真信号幅度谱图。

图5为调节系统结构参数后单稳系统输出信号幅度谱图。

图6为调节乘性高斯白噪声的强度后单稳系统输出信号幅度谱图。

图7为多次调节参数和，系统输出信噪比不再增大时的输出信号幅度谱图。

具体实施方式

本发明针对现有双稳随机共振和分段线性随机共振检测微弱信号需要调节参数较多，不利用工程适用的不足，提供一种基于单稳态随机共振的微弱周期信号检测方法。下面进行详细说明。

为实现本发明所述方法，先建立一个单稳模型系统为：，系统的框图表示如图1所示。

其势能函数为单稳态形式，各个参数的含义为：为系统输出信号，为系统结构参数，是加性高斯白噪声，为被噪声淹没的周期信号，是注入系统的乘性高斯白噪声，其均值为0，强度为D。

根据随机共振理论及申请人数次实验研究，参数调节范围为（0,5），乘性高斯白噪声强度D的调节范围为（0,1）。

利用上述模型检测微弱信号的方法包括如下步骤：

步骤一：初始化单稳系统中的b、。

步骤二：选择固定系统参数b或者乘性高斯白噪声中的一个，调节另一个，直到系统输出信号的信噪比达到最大。

步骤四：重复步骤三，直到系统输出信噪比不再增大为止。

采用上述方法，只需要调节一个系统参数，采用向单稳系统注入乘性高斯白噪声的方法，通过调节这个系统参数和注入系统的噪声强度，实现对微弱低频周期信号的检测。

本系统的技术原理为：单稳系统只有一个稳态，当向该单稳系统注入乘性高斯白噪声后系统变为，设乘性噪声强度为D，则此时系统有两个工作状态，一个是稳态，一个是非稳态。当系统输入含噪声的微弱信号时，通过调节系统参数和注入的乘性噪声强度，使该非线性系统发生随机共振现象，将一部分噪声能量转换为信号的能量，从而大大提高系统输出信噪比，检测出微弱周期信号。

需要强调的是，选择先固定b调节或者选择先固定调节b都是可行的，具体根据实际使用选择。这种可选择性的设计与本单稳模型单一调节方式使得本发明具有更强的操作性，提高应用范围，降低设计难度。

优选的，系统输出信号进行傅里叶变换(FFT)后计算信噪比。

下面结合图2到图6对本发明方法具体应用进行举例阐述。图2为信号检测流程示意图。取输入信号，其中为信号幅度，代表时间，代表信号频率。与信号一起输入单稳系统的加性高斯白噪声是零均值的高斯分布白噪声，其强度为。本实例中取，，=0.1Hz。图3为加性高斯白噪声与输入信号相加后的仿真信号时域波形图，可以看出，输入信号已经被噪声完全淹没。图4为加性高斯白噪声与输入信号相加后的仿真信号幅度谱图，从中不能分辨出信号的频率，所以不能检测出。

实例的具体实现包括以下步骤：

步骤一：初始化单稳态模型的参数，包括系统参数，以及注入的乘性高斯白噪声的强度。本实例中初始化值取=0.2，=0.4。

步骤二：仿真信号加入单稳系统。

基于随机共振理论，本系统的输出信噪比为

其中,.

其中为无信号输入时系统的两个状态间的特征转移速率，为系统两个状态间的特征转移速率对信号的导数。

步骤三：固定注入的乘性高斯白噪声强度，调节系统参数。当=0.6时，单稳系统发生随机共振，输出信噪比最大为7，输出信号幅度谱如图5所示，可以明显分辨出在0.1Hz处存在信号。

步骤四：固定系统参数为=0.6，调节乘性高斯白噪声强度。当=0.35时，单稳系统发生随机共振，输出信噪比最大为12，输出信号幅度谱如图6所示，从图中可以更明显分辨在0.1Hz处的信号。

步骤五：固定步骤四中乘性高斯白噪声的强度最后的取值，调节系统参数，使系统输出信噪比最大；然后固定系统参数为，调节乘性高斯白噪声的强度。重复该步骤，直到输出信噪比不再增大为止。

本实例中当=0.05，=0.75时，系统输出信噪比不再增大。此时系统输出信噪比为25，输出信号幅度谱见图7。

仿真结果表明：通过本发明的单稳系统，通过调节系统参数和注入的乘性高斯白噪声强度，利用随机共振的特性，可以很好地把淹没在噪声中的微弱低频周期信号提取出来。本系统调节的参数少，可以减少系统成本，易于工程实现。本实例中输出信号中还含有噪声，可以再级联一个本单稳系统以提高输出信噪比，也可以采用滤波平滑等措施进一步处理。

以上所述，仅为本发明的一个实例实施方式，但本发明的保护范围不局限于此，任何在本发明涉及的技术范围内进行的变化或替换，都属于本发明的保护范围

本发明的有益效果为：

Claims

1.一种利用单稳系统随机共振效应检测微弱信号的方法，其特征在于，所述方法包括以下步骤：

步骤一：初始化单稳系统中的b、；在该系统中，为系统输出信号，为系统结构参数，是加性高斯白噪声，为周期信号，是注入系统的乘性高斯分布白噪声。

步骤二：选择固定系统结构参数b或者乘性高斯分布白噪声中的一个，调节另一个，直到系统输出信号的信噪比达到最大；

步骤三：固定上一步中被调节的对象最后的取值，调节另一个，直到系统输出信号的信噪比达到最大;

2.如权利要求1所述的利用单稳系统随机共振效应检测微弱信号的方法，其特征在于，系统输出信号进行傅里叶变换(FFT)后计算信噪比。

3.如权利要求1所述的利用单稳系统随机共振效应检测微弱信号的方法，其特征在于，系统结构参数调节范围为（0,5）。

4.如权利要求1至3任一项所述的利用单稳系统随机共振效应检测微弱信号的方法，乘性高斯分布白噪声强度D的调节范围为（0,1）。