一种混合噪声干扰下微弱脉冲信号检测方法技术

本发明专利技术公开了一种用于检测微弱脉冲信号的方法，其先对实测信号进行一维多尺度小波分解，即通过一维小波变换将实测信号分解为N层，并利用小波分解系数对信号的高频部分进行重构，得到N个高频重构信号dj和1个低频重构信号a；然后对重构信号dj分别进行滑动峰态处理，并得到相应的滑动峰态时间序列cj；最后对滑动峰态时间序列cj进行叠加合成和归一化处理，得到最终的滑动峰态时间序列c。通过一维多尺度小波分解预处理，实测信号在不同的分解尺度上表现出了不同的时频域特征，在对其进行滑动峰态处理时，可以实现多滑动窗同步作用的检测效果，有利于对实测信号中混合的噪声、干扰和待检测脉冲等实现差异化分析和处理。

【技术实现步骤摘要】

本专利技术涉及信号处理方法研宄领域，特别涉及一种在混合噪声干扰下检测微弱脉 冲信号的方法。
技术介绍
微弱信号通常是指深埋在噪声干扰中极其微弱的有用信号，这里的微弱主要表示 有用信号的能量或幅值因相对于噪声干扰能量或幅值较小而淹没在噪声干扰中。目前，微 弱信号处理技术在通信工程、电子对抗、生物医学等领域的微弱周期信号检测中已得到广 泛应用，且取得了显著成效。而在实际工程应用领域，还有大量的脉冲信号需要检测，比如 地震勘探、故障诊断、无损检测等工程应用中产生的脉冲回波信号、放电脉冲信号、振动冲 击信号、超声脉冲信号等，在上述脉冲信号的检测过程中，受实测环境下背景噪声干扰、信 号传播衰减以及测试系统灵敏度等因素的影响，均会导致待检测脉冲信号被噪声干扰淹没 的情况出现，即实测信号的信噪比显著降低，无法直接判断待检测脉冲信号是否存在，从而 严重影响其特征信息的准确提取。在信号处理领域，高阶统计量在抑制高斯噪声、反映非高 斯信号特性方面具有显著优势，并得到广泛关注和应用。比如，零均值、零时滞的四阶累积 量即峰态对非高斯信号就非常敏感，能够有效描述信号偏离高斯分布的程度。由于脉冲信 号具有典型的非高斯特性，因此有科研人员将基于四阶累积量的滑动峰态算法应用在机械 故障弱冲击信号的处理和特征提取中。例如，《振动与冲击》第28卷第4期中《基于滑动峰 态算法的信号弱冲击特征提取及应用》给出了滑动峰态算法原理，同时，杨富春等人认为机 械故障振动信号中往往含有故障引起的弱冲击成分，冲击信号具有显著的非高斯特性，而 零时滞四阶累积量即峰态能够描述信号偏离高斯分布的程度；基于峰态这一特性，提取一 种基于滑动峰态算法的弱冲击特征提取方法，首先对原信号进行滑动峰态计算，获得一个 新的峰态时间序列，然后对该峰态时间序列进行傅里叶变换，提取出信号中冲击成分的频 率特征。通过强背景信号及噪声环境下弱冲击特征提取的仿真研宄，证明了该方法具有很 好的冲击特征提取能力。 在这些现有滑动峰态算法中，滑动窗宽度的取值对微弱脉冲信号的检测结果有直 接影响，而滑动窗宽度值的最优取值可以根据峰态值最大化原则计算得到，且其在计算信 号的滑动峰态时间序列时是固定不变的。但在实际应用中发现，背景噪声干扰对最优滑动 窗宽度的取值确定影响较大，不同类型、甚至同类型不同强度的噪声干扰信号，均会影响最 优滑动窗宽度值的计算结果，从而直接影响微弱脉冲信号检测结果。因此，现有滑动峰态算 法在混合噪声干扰下对微弱脉冲信号的检测性能会显著劣化，无法达到预期的检测效果。
技术实现思路
本专利技术要解决的技术问题是提供一种能够在混合噪声干扰下提取微弱脉冲信号 的检测方法，以提高测试系统在低信噪比条件下对微弱脉冲信号的检测能力。 为解决上述技术问题，本专利技术提供了一种混合噪声干扰下检测微弱脉冲信号的方 法，其关键技术在于：其包括下述步骤： 第一步：对实测信号进行一维多尺度小波分解，即通过一维小波变换将实测信号 分解为N层，并利用小波分解系数对信号的高频部分进行重构，得到N个高频重构信号dj(j =1，2, . . .，N)和1个低频重构信号a; 其中，在进行多尺度小波分解时，小波基函数的选择根据待检测脉冲信号的先验 知识来确定；小波分解层数N的选择根据待检测脉冲的提取结果来确定，具体来说，就是当 待检测脉冲的提取结果随分解层数N的增加不再发生变化时，此时的N就是满足要求的最 小且最优分解层数；N个高频重构信号dj(j= 1，2，...，N)中，dl的分布频段最高，dN的分 布频段最低，其余各重构信号按照分布频段从高到低依次排序； 第二步，对重构信号dj(j= 2, 3, ...，N)分别进行滑动峰态处理，并得到相应的滑 动峰态时间序列cj(j= 2, 3,. . .，N); 其中，滑动峰态处理方法为：对于一个具有n个采样点的实测信号x(i) (i= 1，2,. . .，n)，峰态值c4x定义为c4x=E -3{E}2,滑动峰态定义为【主权项】1. ，其特征在于：其包括下述步骤： 1) 对实测信号进行一维多尺度小波分解：通过一维小波变换将实测信号分解为N 层，并利用小波分解系数对信号的高频部分进行重构，得到N个高频重构信号dj(j= 1，2,...，脚和1个低频重构信号a; 2) 对重构信号dj(j= 2, 3,...，脚分别进行滑动峰态处理，并得到相应的滑动峰态时 间序列cj(j= 2, 3,...，脚； 3) 对滑动峰态时间序列cj(j= 2, 3,...，脚进行叠加合成和归一化处理，得到最终的 滑动峰态时间序列C，该时间序列中的脉冲个数及其峰值所对应时刻与实测信号中待检测 微弱脉冲信号的个数及其发生时刻相对应，从而实现了对微弱脉冲信号的检测。2. 根据权利要求1所述的混合噪声干扰下微弱脉冲信号检测方法，其特征在于：所述 步骤1)中，在进行多尺度小波分解时，小波基函数的选择根据待检测脉冲信号的先验知识 来确定；小波分解层数N的选择根据待检测脉冲的提取结果来确定；当待检测脉冲的提取 结果随分解层数N的增加不再发生变化时，此时的N就是满足要求的最小且最优分解层数； N个高频重构信号dj(j=1，2，...，脚中，dl的分布频段最高，dN的分布频段最低，其余各 重构信号按照分布频段从高到低依次排序。3. 根据权利要求1所述的混合噪声干扰下微弱脉冲信号检测方法，其特征在 于：所述步骤2)中，滑动峰态处理方法为；对于一个具有n个采样点的实测信号x(i) (i=l，2，...，n),峰态值C4鹿义为C4,=E-3{E}2,滑动峰态定义为 r^.v(/)，A-(/ +I-l)]，该里C。.；!表示对xa)到x(i+kl)的L个点计算峰态值并取绝对值，利 用滑动峰态遍历实测信号X(i)后，就得到与X(i)相对应的滑动峰态时间序列C(i);L是计算滑动峰态时选取的滑动窗宽度，它的取值根据峰态值最大化原则来确定：就是对应 不同的L值，当滑动峰态时间序列中计算得到的峰态值最大时所对应的L值就是最优的滑 动窗宽度。【专利摘要】本专利技术公开了一种用于检测微弱脉冲信号的方法，其先对实测信号进行一维多尺度小波分解，即通过一维小波变换将实测信号分解为N层，并利用小波分解系数对信号的高频部分进行重构，得到N个高频重构信号dj和1个低频重构信号a；然后对重构信号dj分别进行滑动峰态处理，并得到相应的滑动峰态时间序列cj；最后对滑动峰态时间序列cj进行叠加合成和归一化处理，得到最终的滑动峰态时间序列c。通过一维多尺度小波分解预处理，实测信号在不同的分解尺度上表现出了不同的时频域特征，在对其进行滑动峰态处理时，可以实现多滑动窗同步作用的检测效果，有利于对实测信号中混合的噪声、干扰和待检测脉冲等实现差异化分析和处理。【IPC分类】G01R31-00【公开号】CN104849590【申请号】CN201510226197【专利技术人】刘卫东, 刘尚合, 胡小锋, 魏明, 王雷 【申请人】中国人民解放军军械工程学院【公开日】2015年8月19日【申请日】2015年5月6日本文档来自技高网...

【技术保护点】
一种混合噪声干扰下微弱脉冲信号检测方法，其特征在于：其包括下述步骤：1)对实测信号进行一维多尺度小波分解：通过一维小波变换将实测信号分解为N层，并利用小波分解系数对信号的高频部分进行重构，得到N个高频重构信号dj(j＝1,2,...,N)和1个低频重构信号a；2)对重构信号dj(j＝2,3,...,N)分别进行滑动峰态处理，并得到相应的滑动峰态时间序列cj(j＝2,3,...,N)；3)对滑动峰态时间序列cj(j＝2,3,...,N)进行叠加合成和归一化处理，得到最终的滑动峰态时间序列c，该时间序列中的脉冲个数及其峰值所对应时刻与实测信号中待检测微弱脉冲信号的个数及其发生时刻相对应，从而实现了对微弱脉冲信号的检测。

【技术特征摘要】

【专利技术属性】
技术研发人员：刘卫东，刘尚合，胡小锋，魏明，王雷，
申请(专利权)人：中国人民解放军军械工程学院，
类型：发明
国别省市：河北;13