专利名称：：自适应随机共振微弱信号检测方法自Mi^随机共振微弱信号检测方法駄领域本发明属于低信噪比^i牛下微弱特征频率信号的检测方法，特别涉及一种利用自适应随机共振检测微弱周期信号的方法。
背景技术：
：随机共振技术于上世纪八十年代由意大利学者巴兹等人在研究地球古气候变化时提出。它指一个非线性双稳系统，当仅在噪声或仅在小周期信号作用下都不足以使系统输出在两个稳态之间瑕歐，而在噪声和小周期信号的共同作用下，系统输出的功率谱中，在信号的频率处出现一峰值，当噪声强度达到某一适当值时输出功率谱的峰值达到最大，该峰值的出现是由于产生了噪声能量向信号能量的转移。随机共振利用噪声增强微弱信号传输的机制，使其与其它微弱信号检观仿法相比具有独特的优势，因而在生物信号处理、视觉图像与听觉识别、电磁系统及光信号处理等领域得到广泛重视。在禾拥随机共振检测微弱信号的过程中，涉及到噪声、信号及非线性系统参数三者的最佳匹配问题，只有三者达到最佳匹配关系，随机共振效果最好，微弱信号检测效果也相应最佳。在工程实践中，我们遇到的微弱信号，往往信号频率、噪声强度的大小事先是不知道的，信号频率可以根据相关先验知识判断出大致的频率范围，噪声的强度有大有小，对于不同强度的噪声干扰及不同频率的信号通常采用调节双稳系统参数的办法来使三者之间达到最佳匹配关系，如何调节系统参数，自适应地使双稳系统达到随机共振状态，来增强微弱信号的传输，是工程实践中面临的急需要解决的问题。
发明内容本发明目的在于解决随机共振技术在工程实际中检测微弱特征频率信号面临的问题，提供一种自适应随机共振微弱信号检测方法。该方法特别适合于机械、电子、通信、机械故障检测等领域低信噪比条件下微弱特征频率信号的检测。为了实现上述发明目的，本发明的自适应随机共振微弱信号检测方法，包括以下步骤1)根据随机共振理论，对周期信号与噪声共同作用的双稳系统i:=ax-6+w(O+r(f)进行数值仿真，得出极低频率范围内，小信号条件下非线性双稳系统最佳匹配随机共振的噪声强度与频率的对应关系，式中义为双稳系统输出，仏6为系统结构参数，rw是均值为0，噪声强度为D的高斯分布白噪声，w(O为微弱周期信号，当"(O=^cos(2兀,+p)时；2)利用归一化尺度变换z=jc、,一,r=W、or将双稳模型变换为归一化形式+4)C。s(2;r/0r+p)+cr0"r)用于处理工程中遇到的高频微弱信号;利用最佳匹配随机共振的噪声强度与频率的对应关系，根据公式求出最佳匹配条件下的双稳系统参数"》；式中z是变换后双稳系统的输出，r是变换后时间变量，為是变换后周期信号的幅度，/是变换前周期信号的频率，/o是变换后周期信号的频率，P是周期信号的相位，C7是变换前噪声的均方根，C7。是变换后噪声的均方根，""是均值为O，噪声强度为1的高斯分布白噪声；3)将待检观微弱信号输入最佳匹配条件下的双稳系统，系统输出即为增强的微弱周期信号；针对信噪比极低的微弱周期信号检测问题，将外加信号i^(;c)施加于双稳系统模型i(0=ox—fcc3+w(o+r")+f尸o)用以改变势垒的高度，帮助微弱特征频率信号实现势垒间的跃迁，实现外加信号诱发随机共振，实现微弱周期信号的自适应检测。本发明方法是一种非线性信号检测方法，其特点是它利用非线性双稳系统作为检测模型，将匹配随机共振与外加信号诱发随机共振相结合来检测微弱周期信号，根据不同频率、不同信噪比的信号自适应地调节系统参数，达到最佳匹配随机共振状态，以增强微弱周期信号的传输。该方法为解决随机共振技术在工程实际中检测微弱特征频率信号面临的问题，提供了一种新的有效的技术手段。该方法适合于物理、化学、生物、机械故障早期检测等领域低信噪比条件下微弱周期特征频率信号的检测。下面结合附图进一步说明本发明的技术方案。图1是周期信号频率".001Hz时，双稳系统产生随机共振时信号幅值与最优匹配噪声强度的关系图。图2是工频60Hz时转子系统早期碰摩故障的振动信号波形图。图3是图2所示信号的功率谱图图4是图2所示信号加入强度/>=2.5的高斯白噪声后的信号波形图。图5是图4所示信号的功率谱图。图6是图4所示信号经自适应随机共振处理后的输出波形图。图7是图6所示自适应随机共振输出的功率谱图。本发明自适应随机共振微弱信号检测方法是利用最佳匹配随机共振与外加信号诱发随机共振相结合的方法来实现微弱特征频率信号自适应检测的。其步骤如下一、根据双稳系统随机共振理论，进行数值仿真，得出极低频率范围内，小信号条件下双稳系统最佳匹配随机共振的噪声强度与频率的对应关系。受周期信号与白噪声驱动的非线性双稳系统实质上描述的是一个质点同时受到外力和噪声驱动时，在对称双势阱中的运动。双稳系统可以由朗之万方程描述<formula>formulaseeoriginaldocumentpage7</formula>其中K》表示对称双势阱<formula>formulaseeoriginaldocumentpage7</formula>，得到<formula>formulaseeoriginaldocumentpage7</formula>(2)上式中x为系统输出，a、6为非线性系统结构参数，r(/)是均值为0，噪声强度为D的高斯分布白噪声，当"(0=Jcos(27i,+p)时，输入外力为高斯噪声驱动的余弦信号，当输入(噪声+信号)为零时，±#为系统的两个势阱。当v4X)时，势低点相对于垒高交替地升降。当j24=^/^7^(临界值)时，系统丧失双稳性。因此，在没有输入激励04=0，D=0)时，系统状态局限在两势阱之一中，且与初始条件有关。当输入周期信号频率/1他，幅值爿处于临界值以下04<0.3)时，对模型(2)进行数值仿真，图l是当频率戶O.001Hz，采样频率义=1他时，^以0.01为步长在区间(0.08，0.3)内取值，对输出采用10次平均，通过数值仿真得出的输出信噪比最大所对应的噪声强度D随信号幅值^的变化曲线，其中的直线为最小二乘拟合得到的结果，从结果可知输出信噪比最大时对应的噪声强度D与爿无关。根据随机共振绝热近似理论以及线性响应理论可知，双稳随机共振系统仅有噪声输入时，经M稳系统后输出的频谱能量主要集中在低频段，当信号频率落在该频段时，噪声能量就会向信号转移，从而使该周期成分在频域凸显出来。频率处于0.0014.009Hz之间的信号完全符合系统的低频要求，一旦频率确定，即便幅值爿在一定范围内变化，总存在某个强度相对稳定的噪声能将该周期信号、WJ显现出来，达到最佳匹配的要求。最佳匹配时，噪声强度大小主要由信号的频率决定。有了以上结论，利用数值仿真方法得出输入周期信号频率在0.001~0.009Hz之间变化时，信号幅值小于0.3，双稳系统发生随机共振时，最佳噪声强度与不同信号频率对应的关系。表1给出了不同信号频率在输出信噪比最大时对应的噪声方差(取十次平均)。表1劍圭匹EB寸不同信号频率/^应的麟方差<table>tableseeoriginaldocumentpage8</column></row><table>二、归一化尺度变换。由于以上规律成立的频率范围非常有限，这一步就要解决对于设计任意频率的输入信号来利用以上规律实现最佳匹配随机共振。首先通过变量替换，使随机共振理论能够应用于任意频率的信号。模型(2)中，引入变量替换<formula>formulaseeoriginaldocumentpage9</formula>(3)可以将模型变换为归一化形式<formula>formulaseeoriginaldocumentpage9</formula>(4)<formula>formulaseeoriginaldocumentpage9</formula>其中<formula>formulaseeoriginaldocumentpage9</formula>(5)(4)式即为模型(3)的归一化形式，两者是等价的，且归一化变换后信号频率变为原来信号频率的1/"。因此对于高频信号，可以通皿择较大的参数"，使之归一化为低频信号，来满足随机共振理论的要求。总之，归一化变换以后相当于对信号和噪声同时乘上比例因子^J。如果归一化变换前已经是一个最佳匹配随机共振系统，那么变换以后仍然满足最佳匹配条件，根据表l的结果，如果已知待检信号的频率/，总可以通过变量替换将其对应到表1的频段，由(5)式可以计算出对于任意频率信号最佳匹配随机共振系统的参数"、&由此即可设计出针对任意频率信号的最佳匹配随机共振系统。公式中Z是变换后双稳系统的输出，r是变换后时间变量，Jo是变换后周期信号的幅度，/是变换前周期信号的频率，力是变换后周期信号的频率，伊是周期信号的相位，是变换前噪声的均方根，是变换后噪声的均方根，《0")是均值为0，噪声强度为1的高斯分布白噪声(归一化尺度变换方法见申请人的实用新型专利ZL200520052219.2)。三、对于极低信噪比的周期信号禾l」用外加信号来诱发随机共振。通过以上两步已能够实现最佳匹配随机共振，但是在极低信噪比条件下，即使双稳系统达到了最佳匹配随机共振，仍不能有效地检测出周期信号的频率。为此，需要增强由式(1)所描述的振子对单频信号的响应，给式(1)右边加入可控信号&("，得到丈(0=—rr(/)+i^(jc)=—j^(;c)+爿cos2兀方+r(/)(7)外加信号&(x)是一个关于位置的函数，所以当振子处于双阱中的左边势阱时，可以调节信号使&(^<0)=+4，使振子向右边势阱移动；如果当振子处于右边势阱时，调节信号为&(^0)=-4，使振子向左边势阱跳动。这就使得势垒高度来回地(非周期地)摆动，最终可能使振子越过中间势垒。实际上，输入外加信号后振子的势能表达式为^=K-;ci^，此时的势垒高度低于输入^W信号之前的。其中4>4，4为原双稳系统跃迁临界值。采用该方法可以降低势垒高度，微弱周期信号就可以实现势阱间跃迁，从而检测出该周期信号的存在。下面用一个具体实例来说明本发明用于微弱信号检测的具体过程与有益效果动静4牛碰摩是转子系统中较常见且难以检测和捕捉的一类故障J见象。主要发生在机组动静叶片密封、转子轴封以及滑动轴承等处。碰摩发生的时间短且有时位置不确定，如果碰摩Jim经常发生，可能引发其它严重故障的发生，造成重大经济损失。因此，设法在转子系统碰摩故障发生早期将其检测出来，对于转子系统避免碰摩及继发性故障以及机组的健康运行具有现实意义。通过试验研究转子系统中动静件间尖锐碰摩时的振动特征规律可知在早期碰摩阶段，发现有工频的1/3X、2/3X等分量稳定存在，这一特性可以为这类故障的早期诊断提供依据。下面据此特性，采用夕卜加信号随机共振方法，对碰摩故障进行早期检测。在转子系统试验装置上，我们采集了转子、定子间发生早期尖锐碰摩的径向振动位移信号(以电压信号表征)，图2是工频为60Hz时转子系统早期碰摩故障的振动信号波形图。上述i^数据是在实验理想环境中获得的，表征早期碰摩故障的各种微弱特征还能在其功率谱图3中隐约可见，如果在现场的旋转机械运行环境中，这样微弱的特征信号将被背景噪声所完全淹没。为分析方便，假设背景噪声为高斯白噪声，将图2中信号经一定处理(滤掉1/3工频以上频率成分，采样频率为乂=1000Hz，数据长度为^=3000)再加入同样长度且强度为》2.5的高斯白噪声数据，结果如图4所示。图3、图5分别为上述两种信号的功率谱图。由图3可知，图2所包含的数据存在频率为20Hz、40Hz两个明显的分量。由图5根本分辨不出图4信号中包含有弱周期信号。但采用基于外加信号增强随机共振原理微弱信号检测方法却育&对上述的进行辨识。首先对20Hz频率成分进行检测，按照上述设计方法，首先确定基本模型中的各个参数值根据参数"、6变4t^稳系统输出信噪比增益的影响规律分析，确定『6=1;根据基本模型中频率应该满足随机共振对低频的要求，选定_^0.1Hz;根据前述方法原理中^p的取值原则，确定^=0.9*^0=0.3464。实际检测信号频率为20Hz，根据归一化尺度变换原贝ij，此时^6=200，放大倍数为200(将基本模型中的信号、噪声以及外加信号放大200倍)，采样为50倍周期采样，采样频率为1000Hz，数据长度^=3000，检测结果如图6、图7所示，20Hz周期信号成分清晰可见。由此可以实现低信噪比条件下转子系统早期碰摩故障检测。权利要求1、一种自适应随机共振微弱信号检测方法，包括以下步骤1)根据随机共振理论，对周期信号与噪声共同作用的双稳系统全文摘要本发明公开了一种自适应随机共振微弱信号检测方法。它是通过将匹配随机共振与外加信号诱发随机共振相结合来实现微弱周期信号的自适应检测。根据随机共振理论，首先进行数值仿真得出了极低频率范围内，小信号条件下非线性双稳系统最佳匹配随机共振的噪声强度与频率的对应关系；利用归一化尺度变换，处理工程中遇到的高频微弱信号。对于信噪比极低、噪声强度很大微弱信号的检测，利用诱发随机共振来改变双稳系统势垒的高度，来帮助微弱特征频率信号实现势垒间的跃迁，使双稳系统达到随机共振状态，从而实现微弱周期信号的自适应检测。文档编号G01R23/16GK101191804SQ20071019248公开日2008年6月4日申请日期2007年12月3日优先权日2007年12月3日发明者刘冠军,杨定新,杨拥民,温熙森,秦国军,政胡,胡海峰,胡茑庆,静邱,敏陈申请人:中国人民解放军国防科学技术大学