专利名称：：随机共振-混沌微弱信号检测方法
技术领域：
：本发明涉及一种强噪声背景下的微弱信号检测方法，特别适用于低信噪比条件下，机械系统早期故障诊断中对相关故障是否发生以及发生程度的全面检测。
背景技术：
：随机共振技术于上世纪八十年代由意大利学者巴兹等人在研究地球古气候变化时提出。随机共振利用噪声增强微弱信号传输的机制，使其与其它微弱信号检测方法相比具有独特的优势，因而在生物信号处理、视觉图像与听觉识别、电磁系统及光信号处理等领域得到广泛重视。研究表明，利用随机共振方法能把微弱特征频率信号是否存在有效检测出来，但由于随机共振发生机制中包含噪声的能量向信号能量的转移，且系统产生随机共振时噪声不是某个确定的值，而是某个区间，故无法直接用随机共振方法进行微弱信号幅值估计。机器实际运行过程中，需要全面掌握故障是否发生以及发生的程度，以便进行合理维修。而传统的微弱信号幅值估计方法(如最大似然法等)，在信噪比较低时，估计效果不理想。由此，关于随机共振在机械系统早期故障微弱特征频率信号检测中的应用存在以下亟待解决的问题低信噪比条件下，如何在利用随机共振方法实现微弱特征频率信号是否存在的基础上，对其幅值进行估计？
发明内容本发明目的在于解决随机共振技术在实际微弱特征频率信号检测中无法估计信号幅值大小的问题，提供一种低信噪比条件下，既能快速、准确判断出特征频率信号有无，又能对其幅值大小进行估计的微弱信号检测方法。为了实现上述发明目的，本发明的随机共振-混沌微弱信号检测方法包括如下步骤首先通过一种自适应参数调节随机共振方法实现微弱特征频率信号是否存在的快速、准确检测，得到实际待检信号最可能的信噪比大小；在此基础上，设计出相应的Duffmg混沌振子，并计算得出不同信噪比条件下混沌振子的输入周期激励与输出大尺度周期运动幅值间的变化关系，再根据随机共振方法估计得出的待检信号信噪比大小即可估计得出微弱特征频率信号的幅值大小。从而实现微弱特征频率信号的"有无"和"大小"的全面检测。根据上述方法，本发ji进一步的技术方案是采用如下步骤1)根据随机共振理论，对周期信号与噪声共同作用的双稳系统模型i=ax-fex3+w(o+r(y),，、<r(o>=o,<r(o,r(,')>=2賤h')进行数值仿真，得出模型(1)中"=6=1、周期信号w(/)频率/^0.1Hz时，非线性双稳系统最佳匹配随机共振的噪声强度A)(即噪声均方根值为&，=凤)与输入信噪比57Win(SNRin=201og(i))间的对应关系；-(1)式中X为双稳系统输出，fl、6为系统结构参数，f是时间变量，rw是均值为0、噪声强度为ZK即噪声均方根值为c7，a=V^)的高斯分布白噪声，"w是幅值为j、频率为y;的微弱周期信号；2)滤掉待检信号中的工频成分，估计其噪声均方根值cr,;3)改变待检信号信噪比值，同时对双稳系统模型(1)进行数值仿真，得到实际待检信号最可能的信噪比大小；根据步骤l)所得结果，即可得到该信噪比条件下，频率/。=0.1^所对应的最佳匹配随机共振噪声均方根值CT。，再按照归一化尺度反变换<formula>formulaseeoriginaldocumentpage7</formula>(2)计算出模型(1)中的"》值；4)将"，6值代入模型(1)，得到对应于该输入信号的随机共振检测模型，根据该模型输出频谱图即可判断输入信号中是否存在特征频率为/(角频率为q)的信号；5)如果待检信号中存在特征频率为y;的信号，即可针对乂(即^)，设计出相应的Duffing混沌振子i<formula>formulaseeoriginaldocumentpage7</formula>(3)(3)式中，5为阻尼系数，^为根据Melnikov方法计算得出的混沌振子混沌解与周期解的分叉阈值，"(O和rW在实际信号检测中往往混和在一起，构成实际待检信号；6)采用四阶Rung-Kutta算法对混沌振子(3)进行数值求解，在此基础上估计其大尺度周期解的幅值4;7)通过对混沌振子(3)进行数值仿真，得到待检信号信噪比在某区间内变化时，混沌振子(3)的激励幅值大小y与相应的振子大尺度周期解幅值^,间的变化关系为4(力U+P2(4)(4)式中，y为混沌振子(3)中分叉阈值^与周期信号w(/)幅值J之和，4由步骤6)得出，Pt和P2分别为线性拟合系数；8)根据步骤3)所估计的待检信号信噪比值，结合歩骤7)所得到的该信噪比条件下的振子大尺度周期解幅值^和振子激励幅值H司的拟合经验公式(4)，即可求出y值；9)用y减去混沌振子(3)中的分叉阈值^，即可估计出所检测的微弱特征频率信号幅值。本发明是一种非线性信号检测方法，其特点在于将随机共振与混沌两种方法巧妙地结合起来所设计的参数调节随机共振方法(前4个步骤)在快速、准确检测出微弱特征频率信号是否存在的同时，还为基于混沌振子的微弱信号幅值估计方法(后5个步骤)提供所需的待检信号信噪比大小，再设计相应的Duffing混沌振子估计出微弱信号幅值大小，估计效果优于传统的最大似然法。本发明将以上微弱信号检测方法命名为随机共振-混沌方法。该方法回避了混沌振子微弱信号检测中系统是否发生相变的判断这一难点，弥补了随机共振无法检测信号大小的不足，为微弱信号检测提供了一条新的途径。该方法适合于物理、化学、生物、机械故障早期检测等领域低信噪比条件下微弱特征频率信号的检测。下面结合附图进一步说明本发明的技术方案。图1是原始碰摩信号时域波形。图2是原始碰摩信号频谱图。图3是处理后的碰摩信号时域波形。图4是处理后的碰摩信号频谱图。图5是W^H)时，信噪比增益与信号幅值关系曲线。图6是自适应参数调节随机共振模型输出时域波形。图7是自适应参数调节随机共振模型输出局部频谱图。具体实施例方式下面用一个具体实例来说明本发明用于微弱信号检测的具体过程与有益效果动静件碰摩是转子系统中较常见且难以检测和捕捉的一类故障现象。主要发生在机组动静叶片密封、转子轴封以及滑动轴承等处。碰摩发生的时间短且有时位置不，定，如果碰摩现象经常发生，可能引发其它严重故障的发生，造成重大经济损失。因此，设法在转子系统碰摩故障发生早期将其检测出来，对于转子系统避免碰摩及继发性故障以及机组的健康运行具有现实意义。通过试验研究转子系统中动静件间尖锐碰摩时的振动特征规律可知在早期碰摩阶段，发现有工频的1/3X、2/3X等分量稳定存在，这一特性可以为这类故障的早期诊断提供依据。下面据此特性，采用随机共振-混沌方法，对碰摩故障进行早期检测。图1是工频为60Hz时的一组典型数据，采样频率为/^1000Hz，数据长度为^=2000。上述获得的试验数据是在实验理想环境中获得的，表征早期碰摩故障的各种微弱特征还能在谱图中隐约可见，如果在现场的转子系统运行环境中，这样微弱的特征信号将被背景噪声所完全淹没。为分析方便，假设背景噪声为高斯白噪声，将图1中数据进行滤波处理(此操作同时完成了步骤2)中的滤除工频信号影响)，仅考虑包含被检测频率的窄带宽的情况，所采用的带通数字滤波器(ChebyshevI型6阶滤波器)的通带截止频率分别为16Hz、30Hz，阻带截止频率分别为5Hz、40Hz，通带波纹0.5dB，阻带衰减-45dB。再加入同样长度且强度为》2.5的高斯白噪声数据，结果如图3所示。图2、图4分别为上述两种信号的功率谱图。由图1可知，图l所包含的数据存在频率为20Hz、40Hz两个明显的分量。由图3、图4根本分辨不出该信号中包含有弱周期信号。首先采用参数调节随机共振方法对其特征频率信号进行检测。假设针对20Hz信号设计出的双稳系统随机共振检测模型为^i=(3X-ZX3+力)(5)其中，^)为图3描述的待检信号，相当于模型(1)中的周期信号"(0和噪声信号r(O之和。步骤l):通过对模型(1)进行数值仿真，得到"=6=1、周期信号t^)频率/。-O.lHz时，非线性双稳系统最佳匹配随机共振的噪声强度A)(即噪声均方根值为C7。)与输入信噪比间的对应关系。模型(1)中^庐1、/(H).lHz时，对模型(1)进行数值仿真，采用四阶Rung-Kutta(参见"黄明游编.数值计算方法.北京科学出版社，2003")算法进行数值求解，取计算步长/^0.05，数据长度TV=2000。实际操作中采用噪声均方根值与信号幅值的比值cr"来描述输入信噪比<S7Win。下面以=10时双稳系统最佳匹配随机共振噪声(以下简称最佳噪声)的计算为例，来说明不同C7"时最佳噪声的计算过程。令模型(1)中周期信号"W幅值^以0.05的步长在[O.l2.5]区间内变化，保持CT/JH0不变，对模型(1)进行数值仿真，容易得到信噪比增益(输出信噪比与输入信噪比的比值)与信号幅值之间的变化曲线如图5。由图可知^"=10时，双稳系统最佳匹配随机共振(图5中信噪比增益最大值处)时周期信号幅值^。=0.4，即最佳噪声均方根值cr。=4。为进一步获得各个信噪比条件下的最佳噪声，让输入信噪比^"值在5到20之间(大量实验结果表明，随机共振方法在该区间内检测效果较好)，以0.5为歩长递增，按照上述cr"-10时的计算方法，依次计算得出各个信噪比下的最佳噪声，计算结果见表l。表1最佳匹配随机共振噪声大小与输入信噪比间的对应关系<table>tableseeoriginaldocumentpage11</column></row><table>步骤2):滤掉待检信号中的工频成分，估计其噪声均方根值^。首先对待检信号进行滤波去除工频信号的影响(此操作在对信号进行预处理时已经完成)，再将待检信号近似为纯噪声估计出其噪声均方根A=0.2505。步骤3):计算实际待检信号最可能的输入信噪比，根据归一化尺度反变换关系式(2)计算实际检测模型(5)中的",M直。'步骤2)已经估计得出了待检信号的噪声均方根值o",，假设待检信号的输入信噪比按照表1所列各值变化，并对模型(5)进行数值仿真，依次计算各个输入信噪比条件下的信噪比增益，找出信噪比增益最大值所对应的输入信噪比大小，即为该实际待检信号最可能的信噪比。针对实际转子碰摩信号，计算得出该信号最可能的信噪比为cr"=10。由表1，对应模型(1)中"=6=1时，双稳系统的最佳噪声0"。=4.0。至此，式(2)中的各个参数值均以得到:/尸20Hz、/cKUHz、CTl=0.2505、cr。=4.0，代入式(2)，即可计算得出模型(5)中的参数!a=200，6=2.048xl07。步骤4):将a》值代入实际检测模型(5)，得到对应于该输入信号的随机共振检测模型，根据该模型输出频谱图即可判断输入信号中是否存在特征频率为20Hz的信号。将"=200,&=2.048><107代入模型(5)，采用四阶Rung-Kutta算法对其进行数值计算，计算步长》=0.001、数据长度为A^2000。图6、图7是其输出的时域和频域谱图，从图7中可以清晰分辨出20Hz信号，即可确定待检信号中1/3分频信号的存在。在此基础上，进一步采用混沌方法对其故障特征频率信号幅值大小进行估计，以便对其故障发生程度进行分析。步骤5):针对20Hz信号，设计出相应Duffing混沌振子为1dx2<5CX,/、，、/八-+---x+x=yccos(叫)+s(0其中，5=0.5，w,407r，0.8260，/c通过Melnikov方法(参见"刘曾荣编.混沌的微扰判据.上海上海科技教育出版社，1994")计算得到，^)为待检信号，相当于模型(3)中的周期信号W0和噪声信号r(o之和。步骤6):应用四阶Rung-Kutta算法对混沌振子(6)进行数值求解，在此基础上估计其大尺度周期解的幅值4;采用四阶Rung-Kutta算法对混沌振子(6)进行数值计算，步长/=0.005，计算长度A^2000，得到其大周期解序列，在此基础上，根据最大似然法估计出其大周期解幅值4=1.4202。步骤7):数值计算得出待检信号信噪比在某区间内变化rf，混沌振子(3)的激励幅值大小/(K与K之和)与相应的振子大尺度周期解幅值^间的变化关系。讨论《=1时的Duffing振子为i+5i-x+x3=yccos"w(0+厂(0其中，5=0.5;"(^为待检微弱特征频率信号，i"(0=^cc^;^为根据Mdnikov方法计算得出的混沌振子混沌解与周期解的分叉阈值，且;；=0.8260。由基于混沌的微弱信号检测相关理论可知，当待检信号中存在《=1微弱周期信号时，振子处于大尺度周期运动状态。设已经处于大尺度周期状态的Duffing振子周期运动振幅大小为A，振子激励的振幅大小为y(K与K之和)。大量仿真计算结果表明，当^以0.002为步长在区间内取值，即y在区间内变化，待检信号信噪比j化值在区间[5，20]内变化时，^与/基本保持正比关系。用4("=p,"p2(8)表示各个c7化时拟合的4随;k变化的经验公式，表2列出了^t"以0.5为步长在区间[5,20]内变化时的各个p,，p2值。仿真实验结果表明，以下表格对混沌振子(3)中的频率《*1时都能成立。表2不同信噪比对应的a与z间经验公式拟合系数<table>tableseeoriginaldocumentpage14</column></row><table>步骤8):计算混沌振子(7)中频率为20Hz周期信号的激励幅值大小P根据步骤3)所估计的信号输入信噪比c"40，步骤6)计算得出的4-l.機，以及表2所列出的c7"-10时，p产1.388，p2=0.198，将以上数值代入公式(8)，即可计算得出混沌振子(7)中频率为20Hz的激励幅值大小y=0.88055。步骤9):估计所检测的微弱特征频率信号幅值。由于混沌振子(7)中的激励幅值y相当于分叉阈值^和实际待检信号微弱周期信号幅值&。之和，故碰摩信号中20Hz成分幅值大小&He=0.05455，这与其实际大小^。=0.05538非常接近。按照同样方法可以实现故障信号中40Hz和80Hz频率成分的检测，在应用随机共振方法检测出存在基础上，估计出幅值大小分别为八。=0.0636，&。=0.093，实际幅值为y4。=0.0645，y8。=0.0958。可见，应用随机共振-混沌方法在成功检测出碰摩故障特征信号的同时，还实现了其特征频率信号幅值的估计，能够为实际故障检测中故障发生程度分析提供一定依据。1权利要求1、一种随机共振-混沌微弱信号检测方法，包括如下步骤首先通过一种自适应参数调节随机共振方法实现微弱特征频率信号是否存在的检测，得到实际待检信号最可能的信噪比大小；在此基础上，设计出相应的Duffing混沌振子，并计算得出不同信噪比条件下混沌振子的输入周期激励与输出大尺度周期运动幅值间的变化关系，再根据随机共振方法估计得出的待检信号信噪比大小即可估计得出微弱特征频率信号的幅值大小。2、根据权利要求1所述的随机共振-混沌微弱信号检测方法，其特征是包括以下步骤；1)根据随机共振理论，对周期信号与噪声共同作用的双稳系统模型<formula>formulaseeoriginaldocumentpage2</formula>进行数值仿真，得出模型(1)中"=6=1、周期信号w(Z)频率/广0.1Hz时，非线性双稳系统最佳匹配随机共振的噪声强度A)(即噪声均方根值为&，cr。=V^)与输入信噪比SA^n(SNRin=201og(^))间的对应关系；(1)式中X为双稳系统输出，a、6为系统结构参数，f是时间变量，r(O是均值为0、噪声强度为"(即噪声均方根值为C7，CT=V^)的高斯分布白噪声，t^)是幅值为j、频率为y;的微弱周期信号；2)滤掉待检信号中的工频成分，估计其噪声均方根值(T,;3)改变待检信号信噪比值，同时对双稳系统模型(1)进行数值仿真，得到实际待检信号最可能的信噪比大小；根据步骤l)所得结果，即可得到该信噪比条件下，频率/。-0.1Hz所对应的最佳匹配随机共振噪声均方根值cr。，再按照归一化尺度反变换a=A,(2)计算出模型(1)中的",6值；4)将a,M直代入模型(1)，得到对应于该输入信号的随机共振检测模型，根据该模型输出频谱图即可判断输入信号中是否存在特征频率为/(角频率为q)的信号；5)如果待检信号中存在特征频率为乂的信号，即可针对y;(即^)，设计出相应的Duffing混沌振子JT^_+_^_^—x+JC3二乙cos(6^)+"(0+r(0(3)ft),Aqi/,G)式中，3为阻尼系数，^为根据Mdnikov方法计算得出的混沌振子混沌解与周期解的分叉阈值，t^)和rW在实际信号检测中往往混和在一起，构成实际待检信号；6)采用四阶Rung-Kutta算法对混沌振子(3)进行数值求解，在此基础上估计其大尺度周期解的幅值4;7)通过对混沌振子(3)进行数值仿真，得到待检信号信噪比在某区间内变化时，混沌振子(3)的激励幅值大小/与相应的振子大尺度周期解幅值^间的变化关系为4("U+P2(4)(4)式中，y为混沌振子(3)中分叉阈值^与周期信号w(0幅值J之和，4由步骤6)得出，p,和P2分别为线性拟合系数；8)根据步骤3)所估计的待检信号信噪比值，结合步骤7)所得到的该信噪比条件下的振子大尺度周期解幅值^和振子激励幅值y间的拟合经验公式(4)，即可求出/值；9)用y减去混沌振子(3)中的分叉阈值^，即可估计出所检测的微弱特征频率信号幅值。全文摘要本发明公开了一种随机共振-混沌微弱信号检测方法，首先通过一种自适应参数调节随机共振方法实现微弱特征频率信号是否存在的检测，得到实际待检信号最可能的信噪比大小；在此基础上，设计出相应的Duffing混沌振子，并计算得出不同信噪比条件下混沌振子的输入周期激励与输出大尺度周期运动幅值间的变化关系，再根据随机共振方法估计得出的待检信号信噪比大小即可估计得出微弱特征频率信号的幅值大小。该方法适合于物理、化学、生物、机械故障早期检测等领域低信噪比条件下微弱特征频率信号的检测。文档编号G01H17/00GK101561314SQ20091004336公开日2009年10月21日申请日期2009年5月12日优先权日2009年5月12日发明者潘中银,秦国军,胡茑庆,敏陈申请人:中国人民解放军国防科学技术大学