基于数学形态学的旋转机械振动信号降噪方法.docx
- 文档编号:26424436
- 上传时间:2023-06-19
- 格式:DOCX
- 页数:17
- 大小:24.58KB
基于数学形态学的旋转机械振动信号降噪方法.docx
《基于数学形态学的旋转机械振动信号降噪方法.docx》由会员分享,可在线阅读,更多相关《基于数学形态学的旋转机械振动信号降噪方法.docx(17页珍藏版)》请在冰豆网上搜索。
基于数学形态学的旋转机械振动信号降噪方法
第42卷第4期2006年4月
机械工程学报
v01.42No.4CHINESEJOURNALOF
MECHANICAL
ENGINEERING
Apr.
2006
基于数学形态学的旋转机械振动信号降噪方法木
胡爱军1唐贵基1
安连锁2
(1.华北电力大学机械工程学院保定071003;2.华北电力大学动力工程学院保定071003)
摘要:
基于数学形态学实现振动信号降噪。
研究了数学形态滤波器对振动信号在不同类型、不同强度噪声干扰下
的降噪能力,提出了采用开一闭和闭一开组合数学形态滤波器实现旋转机械振动信号降噪处理的方法。
通过仿真
计算及实例,检验了形态滤波器的滤波效果,表明数学形态滤波器可以有效剔除脉冲、降低随机噪声干扰,提高振动信号的信噪比。
对强烈噪声干扰采用傅里叶变换与形态滤波器结合的处理方法可以取得明显的滤波效果。
并具有算法简单、运算速度快的特点。
关键词:
数学形态学旋转机械滤波器振动噪声中图分类号:
TN911
0前言
振动信号分析是旋转机械状态监测与故障诊断中应用最广泛的方法。
在实际工程测量中,现场采集的振动数据往往被各种噪声污染,在某些情况下噪声干扰甚至大于实际的真实信号,信号降噪成为动态信号测试和设备故障诊断研究的重要内容。
近年来基于小波分析等对信号奇异性特征提取uJ的故障早期诊断取得了一定进展,然而这些研究多是在无噪声情况下进行的,由于小波对噪声和微弱信号同样敏感,降噪也成为其工程应用的重要内容。
数字滤波器是振动信号预处理的常用手段,大多数场合已代替了传统的模拟滤波器。
常用的数字滤波器有时域平均法、IIRfFIR滤波器及小波滤波
剁2】等。
时域平均方法在具体实施过程中需要对大
量的数据进行处理,且要求有时标信息的支持;小波降噪技术的降噪效果则在很大程度上取决于滤波器性能的优劣,即选择不同的滤波器所得的降噪效果也有所区别。
另外,数字滤波器由于基于时域、频域或时频域(如小波)构建,存在着诸如时滞、相移等缺点;对于信号频率和噪声干扰的频率重叠在一起的情况,常用滤波器都无法将两者区分开来。
数学形态学p1(Mathematicalmorphology)是基于积分几何和随机集论建立起来的有别于基于时域、频域的数学方法。
该方法进行信号处理时只取决于待处理信号的局部形状特征,通过数学形态变换将一个复杂的信号分解为具有物理意义的各个部分,将其与背景剥离,同时保持信号主要的形状特
・华北电力大学博士教师基金资助项目(92104392)。
20050525收到初稿20051110收到修改稿
征,要比传统的线性滤波更为有效。
利用构成数学
形态滤波器[31,即使原始信号伴随较强的噪声,甚
至发生了严重的畸变,其基本形状仍可以被识别和重构及增强。
数学形态学在图像处理【4J、形状分析、模式识别和计算机视觉等方面己得到广泛应用,在电力系统的信号分析【5J中也得到了一定应用,而应用于振动信号处理有待深入研究。
根据旋转机械振动信号的特点,引入形态学方法,采用数学形态滤波器处理含噪声的振动信号,消除振动信号噪声,提取出信号特征,重构原始信号。
研究了数学形态滤波器对含有不同类型、不同强度噪声信号的降噪能力,使数学形态滤波器对旋转机械振动信号处理取得了较好的滤波效果。
l
数学形态滤波的基本原理
1.1数学形态基本变换
数学形态学是一种非线性图像(信号)处理和分析工具。
数学形态学以集合来描述目标信号,其思想是设计一个称作结构元素的“探针”(相当于滤波窗),收集信号的信息,通过该探针在信号中不断移动,对信号进行匹配,以达到提取信号、保持细节和抑制噪声的目的。
形态变换一般分为二值形态变换和多值形态变换(灰度变换)。
全部形态变换包括7种基本运算【3],即腐蚀、膨胀、开运算、闭运算、击中、细化和粗化等。
由于在振动信号分析中一般只涉及一维图像信号,只讨论一维离散情况下的灰度值形态变换,包括上述前四种变换及开、闭的级联组合。
基本形态运算在一维信号处理中的定义为设原始信号厂(以)为定义在F=(0,l,…,N—1)上的离散函数,定义序列结构元素g(n)为G=(0,1,…,
万方数据
】28
机械工程学报第42卷第4期
M-1)J2的离散函数,且N≥M,则厂(,z)关于g(聆)波效果,其中以三角形结构元素处理效果最好。
的腐蚀和膨胀分别定义为
(厂融)(刀):
lIlin[厂(n+m)一g(m)]
(1)
2仿真计算
,..、,?
∈u,1'…,朋-1、
,、,
为考察形态滤波器对振动信号中噪声的滤波能
(厂。
g)@)=max[f(/./-m)+g沏)]
(2’
力,葚存孚茹苄若i孬釜:
≤赢菇啬善为‘
“
m∈0,1,…,M一1
z(f)=sin(2XpiX50Xf)+sin(2×piXl00Xf)(8)
f(n)关于g(,2)的开运算和闭运算分别定义为
(/og)(n)=(f0909)(玎)(3)
采样频率取2kHz。
在原始信号中加入不同的(厂・g)0)=(厂ogpg)(以)
(4)
噪声s(f),形态滤波得到滤波后信号厂(f),信噪比为
这里符号。
和・分别表示开运算和闭运算。
开运算使目标轮廓光滑,并去掉毛刺和孤立点,它可以抑制信号中的峰值(正脉冲)噪声;闭运算则填平了
小沟、弥合了孔洞和裂缝,可以滤除信号中的低谷(负脉冲)噪声。
开、闭运算均具有低通特性,它们的组合可以用做形态学滤波。
1.2形态开、闭组合滤波器
为了同时去除信号中正、负两种噪声,通常采用形态开、形态闭的级联形式。
MARAGOS采用相同尺寸的结构元素,通过不同顺序级联开、闭运算,定义了形态开一闭和闭一开滤波器【6~1。
Foc(厂(n))=(厂og・g)(,z)(5)最D(厂(n))=(/・gog)0)
(6)
由以上定义知,形态开一闭和闭一开滤波器具有开闭运算的所有性质,可以同时滤除信号的正负脉冲噪声,但是,由于开运算的收缩性导致开一闭滤波器的输出偏小,闭运算的扩张性导致闭一开滤波器的输出偏大,因而存在统计偏倚现象,直接影响到滤波器的噪声抑制性能,所以单独使用它们并不能取得较好的滤波效果。
为了有效地抑制信号中的各种噪声,采用级联开、闭运算,构造开一闭和闭一开组合形态滤波器【8】,用于振动信号的降噪
1
y(n)=亡[‰(厂(甩))+Foc(厂(以))】
(7)
二
1.3结构元素
除了运算方式的组合外,结构元素的选择对信号处理结果也有很大影响。
常用的结构元素有直线、曲线(如二次,三次等)、三角形、圆形和其他多边形(如钻石形,六角形等)及其组合等。
结构元素的设计决定于要处理的信号的形状,其结构要尽可能接近待分析信号的图形特点。
郭景峰等[4】采用组合形状结构元素处理图像边缘噪声。
李刚等[9】采用直线形、三角形、正弦形状的结构元素处理心电信号取得相近效果。
对于振动信号的处理,三角形、圆形、直线形三种形状的结构元素均可取得较好的滤
以NR
2
lO
J.g(p,/5)
(9)
式中,竹为信号功率,只为噪声功率。
定义太NR,为滤波前信号的信噪比,冬NR,为滤波后的信噪比。
滤波后的噪声由厂(幻减去x(0得到。
2.1脉冲噪声形态滤波分析
在旋转机械振动信号现场采集中经常存在随机脉冲干扰,工程上一般采用测量相邻采样值,取其差值作为判断依据的限幅滤波法来剔除脉冲,但阈值选择及剔除脉冲后相应数据点的替代值的选择均会给数据带来很大误差。
为考察形态滤波器对脉冲干扰的滤波效果,在x(f)中间隔一定时间,加入幅值为lO的正负脉冲干扰信号,未滤波的信号波形及频谱如图1a、lb所示。
从图1中可以看出,由于脉冲有一定的周期性,频谱图中出现了x(f)中没有的谐波分量,基频幅值增大,已不能反映原始信号的真实状态。
对于此情况,常用的滤波器己无法得到好的滤波效果。
g
÷
飞
迥馨
1
呈-飞
器。
时间r/点数
(a)
0
100
200
300
400
500
600
700
800
频率f/Hz
(b)
图1脉冲干扰信号的时域波形及频谱
选择三角形结构元素,采用式(7)构造开一闭和
闭~开组合形态滤波器,对信号进行形态滤波处理后的时域波形及频谱如图2a、2b所示,时域图中实线为滤波后的波形,虚线为原始信号x(f)。
对比图1、
必
撇
万方数据
2006年4月胡爱军等:
基于数学形态学的旋转机械振动信号降噪方法
129
2可以看出,滤波后的波形完整再现了原始信号的时域及频域特征。
说明形态滤波取得了满意效果。
g
三
飞
j鹫肇
时间f/点数
(a)
频率f/Hz
(b)
图2脉冲干扰经形态滤波后波形及频谱
2.2对随机噪声形态滤波分析
在z(幻中加入均值为0,标准差为艿的均匀分布的白噪声,并通过改变的万大小来改变噪声强度。
6=0.3时对信号进行形态滤波处理结果如图3所示,图3a为未滤波的信号时域波形,图3b为形态滤波处理后的时域波形。
时间t/点数
(a)
ⅢM八忡
时间t/点数
(b)
图36=0.3白噪声形态滤波结果
从图3中可以看出,滤波后的波形与原始信号
波形非常接近。
滤波前后的信噪比太Ⅲ和忍N配分
别为20.86dB及29.88dB。
在仿真计算中还发现,
随着蹦加大,滤波后信噪比的提高值有增加的趋
势,但滤波后信噪比的绝对值也在降低,如万=0.6时,太NR,和太NIl2分别为8.86dB及20.76dB。
这是由形态滤波器的波形特征保持性所决定的,因为在强噪声干扰下,原始波形部分特征已被噪声淹没。
2.3
强噪声干扰下的频域形态滤波分析
虽然强烈噪声淹没了时域信号的原始特征,但
在频域,旋转机械的振动信号一般仍会保持对应信
号频率的幅值及相位特征,考虑形态滤波器对脉冲的良好滤波性能,提出了采用频域形态滤波处理强噪声干扰的方法。
首先对含噪声的信号进行傅里叶变换,保留相位信息,同时对幅频特性曲线进行形
态滤波处理,提取峰值信号,再用傅里叶逆变换得到滤波后时域信号,为与第2.2节中方法区分,称之为频域形态滤波法。
在x(0中加入白噪声,选择三角形结构元素,采用如上方法进行滤波处理。
万=1时的滤波效果如图4所示,2sNRl=--0.02
dB,2SNR2=33.30dB。
可
以看出,即使存在强噪声干扰,原始波形无法分辨,形态滤波仍能够较好的恢复原始波形。
g
{
飞
坦馨
g
毒
飞
邋馨
0
50
lOO
150
200
250
300
350
400
时间r/点数
(a)
时间t/点数
(b)
图4万=l白噪声频域形态滤波结果
3
应用实例
在水泵试验台上测得的振动信号时域波形及频
谱如图5所示,测量时转速为2
944
dmin,采样频
率为2k。
为清楚显示图形特征,时域信号均给出部分波形,频谱图上限频率取400Hz。
信号显示,时域信号含有噪声成分,脉冲干扰明显且无规律,振动峰值较高,频谱图中以基频分量为主,含有无规律的高次谐波。
对停机后的噪声测试显示,噪声含有高峰值的脉冲干扰,噪声频谱中各尖峰的峰值及频率分布与原始信号中的谐波成分非常接近。
对信号进行形态滤波处理后的波形及频谱如图6所示。
由图6看出,噪声得到了很好的抑制,时域信号较光滑,频谱图中主要为基频分量,其他谐波分量峰值很小,说明形态滤波取得了良好效果。
万方数据
130
机械工程学报
第42卷第4期
吕
{
弋
遥馨
g
÷
飞
遥馨
100
0
100
2010
0
E
÷
飞
j四馨
0
O.05
O,lO
O.15
时间t/s
频率,/Hz
图5滤波前振动信号时域波形及频谱
20
莹
运10
骧
0
频率f/Hz
图6滤波后振动信号时域波形及频谱
4结论
针对旋转机械振动信号数据采集的特点,提出了采用开一闭和闭一开组合形态滤波器处理振动信号的新方法,研究了数学形态滤波器对不同类型、不同强度噪声下振动信号的降噪能力,仿真计算及试验结果表明数学形态滤波器对旋转机械振动信号取得了良好的滤波效果,可以有效提高振动信号的信噪比,特别是对于不明原因的脉冲干扰处理,比其他数字滤波器有明显的优势。
对于强烈噪声干扰采用傅里叶变换与形态滤波器结合的处理方法可以取得明显的滤波效果。
形态变换的实现一般只包含布尔运算、加减法运算而不需要做乘法,计算简单,运行速度快,这是数学形态学的又一突出优点。
参考文献
[1]彭志科,何永勇,卢青,等.小波多重分形及其在振动
信号分析中应用的研究[J].机械工程学报,2002,38(8):
59-63.
[2】徐敏强,王日新,张嘉钟.基于梳状小波的旋转机械振
动信号降噪方法的研究[J].振动工程学报,2002,15(1):
89.91.
[3]崔屹.图像处理与分析——数学形态学方法及应用
【M】.北京:
科学出版社,2000.
[4】郭景峰,申光宪,郑绳楦,等.数学形态学在数字滤波
中的应用研究[J].机械工程学报,2002,38(10):
144—147.【5】岳蔚,刘沛.基于数学形态学消噪的电能质量扰动检测
方法[J】.电力系统自动化,2002,26(7):
13.17.[6]MARAGOS
P,SCHAFERR
W.Morphologicalfilters-一
Parth
Theirsettheoreticanalysisandrelationtolinear
shiftinvariant
filters[J].IEEETram.on
ASSP,1987,
35(8):
1
153—1
169.
[7]MARAGOSP,SCHAFERRW.Morphologicalfilters—
PartII:
Theirrelationto
median,order.statistic,andstack
filters[J】.IEEETrans.onASSP,1987,35(8):
1
170・l184.
[8】赵春晖.数学形态滤波器理论及其算法研究【D].哈尔
滨:
哈尔滨工业大学,1998.
[9】李刚,何宏献,刘巍,等.数学形态滤波器及其在心电
图机中的应用【J】.仪器仪表学报,1999,20(4):
335.339.
DE—NoISINGTECHNIQUEFoRVIBRATIoNSIGNALSoFRoTATING
MACHINERYBASEDoN
MATHEmTICAL
MoRPHoLoGY
FILTER
HUAijuni
TANGGuijii
ANLiansu02
(J『.Mechanical
EngineeringDepartment,
NorthChinaElectricPowerUniversity,
Baoding071003;
2.PowerEngineeringDepartment,NorthChina
ElectricPowerUniversity,Baoding071003)
Abstact:
De—noisingthevibrationsignalofrotatingmachinery
realizedbasedon
mathematicalmorphology.Newnoisereduc—
tiontechniquebased
on
open・-closing
andclose・-openingcorn--
bined
morphologicalfiltersproposed.Theabilityofmorpho-
logicalfilter
dealswithdifferenttypesofnoisediscussed.The
methodofFTcombinedmorphologicalfilteralsoproposedsto
reduce
intensity
noise.Simulation
andexperimentexample
demonstratetheeffectivenessofthealgorithm.Itiseffectiveto
filter
thenoisesin
vibrationsignaland
improve
S/Nratio
prominently.111ealgorithmisfastandpropitioustoengineering.Keywords:
Mathematicalmorphology
Rotating
machinery
Filter
Vibration
Noise
作者简介:
胡爱军,男,1971年出生,博士研究生。
主要研究方向为旋转机械振动测试与故障诊断。
E—mail:
hua日llll316@sohu.corn
万方数据
基于数学形态学的旋转机械振动信号降噪方法
作者:
胡爱军,唐贵基,安连锁,HUAijun,TANGGuiji,ANLiansuo
作者单位:
胡爱军,唐贵基,HUAijun,TANGGuiji(华北电力大学机械工程学院,保定,071003,安连锁,ANLiansuo(华北电力大学动力工程学院,保定,071003刊名:
机械工程学报
英文刊名:
CHINESEJOURNALOFMECHANICALENGINEERING年,卷(期:
2006,42(4被引用次数:
13次
参考文献(9条
1.徐敏强;王日新;张嘉钟基于梳状小波的旋转机械振动信号降噪方法的研究[期刊论文]-振动工程学报2002(012.彭志科;何永勇;卢青小波多重分形及其在振动信号分析中应用的研究[期刊论文]-机械工程学报2002(083.李刚;何宏献;刘巍数学形态滤波器及其在心电图机中的应用[期刊论文]-仪器仪表学报1999(044.赵春晖数学形态滤波器理论及其算法研究1998
5.MARAGOSP;SCHAFERRWMorphologicalfiltersPartⅡ:
Theirrelationtomedian,order-statistic,andstackfilters[外文期刊]1987(08
6.MARAGOSP;SCHAFERRWMorphologicalfiltersPartⅠ:
Theirsettheoreticanalysisandrelationtolinearshiftinvariantfilters1987(08
7.岳蔚;刘沛基于数学形态学消噪的电能质量扰动检测方法[期刊论文]-电力系统自动化2002(078.郭景峰;申光宪;郑绳楦数学形态学在数字滤波中的应用研究[期刊论文]-机械工程学报2002(109.崔屹图像处理与分析--数学形态学方法及应用2000
引证文献(13条
1.杨杰.郑海起.刘晓平.王彦刚多尺度柔性形态滤波在轴承故障诊断中的应用[期刊论文]-轴承2011(3
2.李兵.张培林.刘东升.任国全.米双山基于形态梯度解调算子的齿轮故障特征提取[期刊论文]-振动、测试与诊断2010(1
3.张君.王金平.朱波旋转机械振动信号处理中的滤波及特征提取技术[期刊论文]-南京工程学院学报(自然科学版)2009(3
4.杜必强.唐贵基.石俊杰旋转机械振动信号形态滤波器的设计与分析[期刊论文]-振动与冲击2009(95.沈路.周晓军.张文斌.张志刚广义数学形态滤波器的旋转机械振动信号降噪[期刊论文]-振动与冲击2009(96.曹冲锋.杨世锡.杨将新大型旋转机械非平稳振动信号的EEMD降噪方法[期刊论文]-振动与冲击2009(97.杜必强.王松岭.唐贵基形态学降噪在转子分形故障状态识别中的应用[期刊论文]-中国电机工程学报2008(268.刘占生.窦唯旋转机械振动参数图形边缘纹理提取的数学形态学方法[期刊论文]-振动工程学报2008(39.张文斌.周晓军.林勇广义形态滤波器在振动信号处理中的应用研究[期刊论文]-农业工程学报2008(610.窦唯.刘占生.王政先.何鹏旋转机械状态参数图形识别的免疫-模糊形态学方法[期刊论文]-航空动力学报2008(6
11.章立军.徐金梧.阳建宏.杨德斌自适应多尺度形态学分析及其在轴承故障诊断中的应用[期刊论文]-北京科技大学学报2008(4
12.杜必强.王松岭.唐贵基转子系统振动信号形态滤波器的设计[期刊论文]-中国工程机械学报2008(213.章立军信号的数学形态学分析方法及其应用研究[学位论文]博士2007
本文链接:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 数学 形态学 旋转 机械振动 信号 方法