基于多传感器时频分布的机械故障信号欠定盲源分离方法

  • 格式:pdf
  • 大小:10.25 MB
  • 文档页数:6

第41卷第6期武汉科技大学学报Vol.41,No.62018年12月JournalofWuhanUniversityofScienceandTechnology췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍췍Dec.2018

收稿日期:2018-08-20 基金项目:国家自然科学基金面上项目(51875416,51805382). 作者简介:李小彪(1994-),男,武汉科技大学硕士生.E-mail:1157620502@qq.com 通讯作者:易灿灿(1989-),男,武汉科技大学讲师,博士.E-mail:meyicancan@wust.edu.cn

DOI:10.3969/j.issn.1674-3644.2018.06.009基于多传感器时频分布的机械故障信号欠定盲源分离方法

李小彪,吕 勇,易灿灿(1.武汉科技大学冶金装备及其控制教育部重点实验室,湖北武汉,430081;2.武汉科技大学机械传动与制造工程湖北省重点实验室,湖北武汉,430081)

摘要:结合多传感器时频分布(multisensortime-frequencydistributions,MTFD)和盲源分离(blindsourceseparation,BSS)的特点,提出一种针对机械复合故障信号的欠定盲源分离方法。首先利用Wigner-Ville分布

将观测信号转化为MTFD矩阵;然后对该矩阵进行白化处理和噪声阈值处理,并对其自动项进行选择,对其特征向量进行集群处理,从而得到源信号TFD的估计;最后对源信号进行重建,得到源信号的估计。仿真及试验结果表明,本文所提出的方法在处理非平稳复合信号的欠定盲源分离方面具有很好的效果。关键词:多传感器;时频分布;欠定盲源分离;故障诊断;Wigner-Ville分布;故障信号

中图分类号:TH165+.3 文献标志码:A 文章编号:1674-3644(2018)06-0453-06

在实际工业生产中,多种设备故障有可能同时出现,其振动信号表现为故障特征的相互耦合,给故障诊断带来很大困难。盲源分离(blindsourceseparation,BSS)是从可观测的混合信号中恢复不可观测的源信号的方法,为机械复合故障诊断提供了一种有效途径,相关研究成果也有不少。焦卫东等[1]结合独立分量分析和主分量分析的各自优势,将隐藏于多通道观测信号中的独立源信号分离出来,从而实现不同源信号的识别。Li等[2]研究了基于两阶段稀疏表示的欠定盲源分离方法。董绍江等[3]在形态学滤波的基础上结合匹配跟踪算法和盲源分离算法各自的特点,提出了一种基于最优匹配追踪信号分解的单通道盲源分离方法。然而这些机械故障源分离方法要求源信号满足非高斯、平稳且相互独立的假设,但在复杂工况下,实际采集到的机械振动信号并不一定能满足这些前提条件。针对以上不足,研究者提出了一种基于时频(time-frequency,TF)分析的设备故障源盲分离方法,由于时频分布(time-frequencydistribu-tion,TFD)提供了信号的频谱内容随时间变化的信息,因此该方法充分利用了信号的非平稳特征。目前,基于Cohen类时频分布[4]

、分数傅里叶变

换[5]、局域均值分解[6]和Wigner-Ville分布(WVD)[7]

等时频分析工具的机械故障信号盲源

分离研究成果有不少,但仍有进一步拓展的空间。多传感器时频分布[8](multisensortime-fre-

quencydistribution,MTFD)是从单传感器时频

分布[9]发展而来的,其结合了多通道高分辨率的时频分析和阵列信号处理方法,可用于源定位的信号估计和非固定源的多组分信号分离,本文将其应用于机械故障欠定盲源分离,提出一种基于MTFD的盲源分离方法(简称为MTFD-BSS),并

通过仿真分析和试验研究对其有效性进行验证。1 MTFD-BSS方法设计

MTFD-BSS方法的基本思路是利用Wigner-

Ville分布将观测信号转化为MTFD矩阵,对其

进行白化处理和噪声阈值处理,然后通过自动项选择和对特征向量的集群处理得到源信号TFD

的估计,最后对源信号进行重建,从而得到源信号的估计。具体步骤如下。(1)计算观测信号。

在假定源位于远场的条件下,给出观测信号武汉科技大学学报2018年第6期的线性数据模型:z(t)=As(t)+η(t)(1)式中:z(t)=[z1(t),…,zm(t)]T为观测信号向量;s(t)=[s1(t),…,sn(t)]T为源信号向量;η(t)为噪声向量;A为m×n的混合矩阵。(2)计算观测信号的MTFD矩阵。首先,引入非固定的单传感器TFD模型。对于给定的观测信号z(t),其Wigner-Ville分布Wz(t,f)定义为信号z(t)的瞬时自相关函数Kz(t,f)的傅里叶变换:Wz(t,f)=∫RKz(t,τ)e-j2πfτdτ(2) 瞬时自相关函数Kz(t,f)定义为:Kz(t,τ)=z(t+τ2)z*(t-τ2)(3) 为了适用性更加广泛,将WVD与一个相关的二维TF内核进行卷积运算,得到二次TFD矩阵ρz(t,f):ρz(t,f)=Wz(t,f)**tfγ(t,f)(4)式中:**tf表示t和f的二重卷积;γ(t,f)为TF平滑的内核。然后,将式(4)扩展到多传感器的情形下。对于一个观测信号向量z(t)=[z1(t),…,zm(t)]T,其MTFD矩阵ρzz(t,f)为:ρzz(t,f)=Wzz(t,f)**tfγ(t,f)(5)在空间WVD矩阵Wzz(t,f)中单个元素的每一项上都应用双卷积**if,即Wzizj(t,f)=∫Rz(it+τ)2z*(jt-τ)2e-j2πfτdτ(6)式中:i,j=1,2,…,m。因此有Wzz(t,f)=Wz1z1(t,f)Wz1z2(t,f)…Wz1zm(t,f)Wz2z1(t,f)Wz2z2(t,f)…Wz2zm(t,f)︙︙︙Wzmz1(t,f)Wzmz2(t,f)…Wzmzm(t,féëêêêêêêùûúúúúúú)(7)观测信号与源信号之间的关系则转变为:ρzz(t,f)=Aρss(t,f)AH+σ2Im(8)式中:ρss(t,f)为源信号的时频分布矩阵;σ2Im代表噪声信号,其中σ2为所添加噪声的方差,Im为单位矩阵。(3)数据白化预处理。首先,考虑一个n×m矩阵H,可以证明(HA)(HA)H=UUH=In,其中U=HA是一个n×n白化酉矩阵(它将观测信号部分白化)。然后通过以下公式可计算得到白化的MTFD矩阵ρ-zz

(t,f):

ρ-zz(t,f)=Hρzz

(t,f)HH(9)

使用白化矩阵H和式(8),能在无噪声的情

况下表示ρ-zz

(t,f):

ρ-zz(t,f)=Uρss

(t,f)UH(10)

可以采用不同的方法估计白化矩阵H

。一种

方法是通过自相关矩阵的反平方根来估计H,本

文则使用MTFD矩阵来计算H。(4)对观测信号的

MTFD

矩阵进行噪声阈值

处理。对信号TFD应用噪声阈值处理可以降低计算成本。阈值ε1用于只保留有足够能量的点(ts

,

fs

),通常取ε1=5%。这个步骤可以表示为:

Keep(ts,fs)if‖ρzz(ts,fs)‖>ε1(11) (5)对观测信号的MTFD矩阵进行自动项选

择。将自动项与交叉项分开需要一个适当的检测标准。若给定源信号时频图不相交的条件,则每个自动项MTFD矩阵的秩都是1,或者是自动项

MTFD矩阵的特征值比其他矩阵的特征值大。

于是可以使用秩作为选择标准来区分(t,f)点,即

选择所有MTFD矩阵中秩为1的点作为自动项点。这个步骤可以表示为:

ifλmaxρzz(t,f{})

normρ

zz

(t,

f{})

-1>ε

2

,

then点(t,f)是一个交叉项点(12)式中:参数ε2是一个可以忽略的正标量,通常取ε2=10-4;λmaxρzz(t,f{})为ρzz(t,f)的最大相对

特征值。(6)向量集群和源信号

TFD

估计。

具体步骤包括:①STFD特征向量分解;②向量集群;③对TF特征进行排序;④对相交点使用子空间投影分离技术;⑤对源信号TFD矩阵进行估计。(7)对源信号进行重建。

本文提出的MTFD-BSS方法既适用于欠定型盲源分离问题,也适用于超定型盲源分离问题。而传统的基于时频分析(如Cohen类时频分布、分数傅里叶变换等)的盲源分离方法往往具有一定的局限性,只能处理单一情形。2 仿真分析

为了验证MTFD-BSS方法的有效性,构造如下3个非平稳的线性调频信号:

454