振动信号的分形滤波法.pdfVIP

第11卷增刊 振动 工程学报 V01．1lNoS 1998年{月 Jou】．nalofVibration Sep．1998 Engineering 振动信号的分形滤波法 侯祥林虞和济 (东北太学理学院 沈阳．110006) 摘要以一定采样频率采得振动信号也括反映机械设捂状态的真宴信号和噪声信号，为了用真 实信号实项监测和故障诊断，就要去除噪声信号。多数噪声为白噪声或有色噪声，具有很宽颡带 或全频带。这样浞有噪声的信号就可以看成柑对较高频率的多种振动信号的叠加，具有维数大 于1的丹形结构。奉丈提出采并{分形力法消除噪声信号。使用维数函数作为组合参数逐步减少 噪声信号的影响，通过实例验证了遘种丹形滤波方法的可行性。 关键词：信号处理；分形滤渡；振动；故障诊断 f年田分类号：THll3．C“ L l’ 1分形滤波原理 l_1滤波一般过程 设采样时nd间隔为工f，采得样率为∞，jz，…，z。，五一M十￡．，M为真实信号，=．为噪声信 号，滤波是去掉噪声信号≈．获得真实信号且的过程，实际只能获得M的近似估计值多，， 器，滤波输出值与观测值和前步滤波输出相关。可建立递推关系，滤波输出值是上步滤波输 出值和本步信号观测值构成本步信号估计值的线性纽合： Y，=kz．一如Y“，∥l=j，，i=2．…，月 (1) 矗t+矗z=1．^，，七。为滤波输出值和信号观测值的权。五。，五?是观测值的函数。对于低通滤波来 说，如果无噪声信号，邯么七t=l，凫，=0．而有噪声信号0尼。，志：1，噪声信号越强，^，越小， 女2越大，噪声信号越弱，^。越大．k：遗小， 1．2分形滤波原理 用网格法来计算信号分形盒维数，当出在无标度区内，可获得不变的维数d，实际对带 有噪声的采样信号所计算出的维数d为近似维数，这里称为计算维数。对于低通滤波器来 说，噪声信号处i二高频率，结果反映采样信号的计算维数较高。击除噪声信号的结果是滤波 后信号向一维低频转化，使其反映真实佑号。对于高通滤波器来说，噪声信号处于低频率，结 果反映采样信号的计算维数较低。这要使处理后信号向高维高频转化。 收稿日期：Is98—06，{o 增刊 候样林等：振动信号∞分形滤波法 Infy{(tF d— ㈣ (2) ㈣熙‘—2￡一i．i茏1， 这里主要采用低通滤波器。在(1，式中，取序。一厂l(J)．^。一^(0)，自，。^：成为计算维数的函 数，取屉z=1,／d，k。=1一是，一l一1／d作为它们的一种形式。^．，正：是相关的。这种作法，可见高 频噪声信号越强，计苒维数d越大，亡，越小，正!越九：高频噪声信号越弱，计算维敬d越小．五． 越大，七。越小。存在噪声信号时，满足Ok，，^：1．通过递推，达到采用分形方法逐步去噪达 到滤波的目的。 2分形滤波程序框图 K一定粕阜 图1分彤浊披框吲 3分形滤波实例 本文通过上面分析对W眦-rsl 25jo-21t,sin(Akt)。这里f以秒汁√ff)以IDIII|{这是个分形函数，它具有处处连续不可 … 微，盒维数是j的特性。实际处理t：蔺数只能将无阳项取成有限项．可将函数表示为．^(，) Ⅳ ±一j o““．图象如图2所示。 fl。。一1．5L00／2”一4．065*10“．7’juu=2．461—1 1s数量级，得到表1。 果称计