热门关键词:

基于独立分量分析的欠定盲源分离方法

  • 该文件为pdf格式
  • 文件大小:265.12KB
  • 浏览次数
  • 发布时间:2014-10-12
文件介绍:

本资料包含pdf文件1个,下载需要1积分

Underdetermined blind source separation method based on independent component analysisYANG Jie ,YU Wen-wen ,TIAN Hao ,GUAN Zhen-zhen(1.Ordnance Engineering Colege,Shijiazhuang 050003,China;2.Huayin Ordnance Center,Huayin 714200,China)Abstract: The underdetermined blind source separation(UBSS)algorithm at present can only separate sparsesignals,but cant separate non-sparse signals successfully.Classical ICA algorithms,such as,extended Infomax,canseparate both super-Gaussian and sub-Gaussian signals,but it is used only in an over-determined BSS(OBSS).Combinedwith an extended Infomax,an underdetermined ICA algorithm was proposed here.By generating hidden data,an UBSSproblem was transform ed into an OBSS one,and then an extended Inform ax algorithm was used to analyze the signals。

This method could separate both super-Gaussian and sub-Gaussian signals in an UBSS problem. Through analyzingtransient signals of a gearbox by use of the underdetermined ICA combined with the order envelope spectral analysis,itsfault features were fuly detected and the effectiveness of the proposed method was verified。

Key words:independent component analysis(ICA);extended infomax;underdetermined BSS;fault diagnosis齿轮箱发生故障时常出现多种故障并发情形,用诊断单-故障的频谱诊断法、轴心轨迹法、小波诊断法已不再有效,而盲源分离技术为较有效的机械复合故障诊断方法。但已有研究基本针对超定盲分离而言,即要求观测器个数多于或等于源信号个数。少于源信号个数的欠定盲分离问题,已成为盲信号处理领域中挑战性问题。

源数多于传感器数时,机械故障源分离主要考虑稀疏源,利用源信号稀疏性估计混合矩阵;源稀疏性较好时,通过聚类方法可达到较好的盲分离效果 。

源信号稀疏性较差时,虽可利用稀疏性得到混合矩阵估计,但通过聚类方法得不到较好的分离效果。

独立分量分析(ICA)是解决盲源分离问题 (BSS)基金项 目:国家 自然科学基金资助项目(50775219);军械工程学院院基金项目(yjxml0019)收稿日期:2011-04-08 修改稿收到日期:2012-04-18第-作者 杨 杰 女 ,博士生,1983年生的主要算法,其可使输出信号各分量统计独立以实现信号分离。已有多种不同算法 J,如扩展 Infomax算法。由于其既能分离超高斯信号,也能分离亚高斯信号,因而得到广泛应用。但该算法只能用于超定盲分离。本文提出的欠定 ICA算法,通过生成隐藏数据将欠定盲分离转化为超定盲分离问题,再用扩展 Infomax算法对信号进行分析。该方法可分离欠定情形下超高斯与亚高斯混合信号♂合包络阶次方法对齿轮箱实测瞬态信号进行分析,找到故障特征,验证该方法的有效性。

1 基于独立分量分析的盲源分离算法1.1 独立分量分析问题描述图1为 (t)由m个传感器获得m路观测信号,每个观测信号由 个独立源信号 (t)(i1,2,, )的线性混合,即:( )As(t) (1)式中: (t)[ 。(t), :(t),, (t)r为m个混合第7期 杨 杰等:基于独立分量分析的欠定盲源分离方法 31信号,s(t)[s (t),s (t),,s (t)] 为n个未知源信号,A[a ,a ,,a ]为未知的 m×/-t维混叠矩阵,ICA算法目的为从多路观测信号 中分离出隐含其中的源信号 s,即寻找分离矩阵 w,有:uWx (2)式中: (t):[ ,(t), :(t),, (t)]T为对隐含源信号分量 s的逼近。

- - - 1、、 束知混合系统 /~ 图 1 盲源分离原理图Fig.1 The principle sketch of blind source separation1.2 扩展的 Infomax算法Infomax算法L4 处理过程为:经矩阵 w解混后对所得 的每个分量 分别用-个非线性函数g ( )进行处理 ,得:rg( )[g1(1),,gⅣ(Ⅳ)] , , r ]T (3)自适应处理 目标函数可调节 w 使 r的总熵量(r, )极大,由信息理论知,日(r, )极大表示u的各分量间互信息,(u -,u )极校传统Infomax算法中分离矩阵 w迭代式为:W(n1)W(n) [,- ( )M ]W(n) (4)式中: 为时变调节步长。 ( )- 等,,- ,理论上g (M )应取信源s 的概率密度函数,( i)即p(Si)的累积概率分布,而 (i)- 。

源信号中既有超高斯信号又有亚高斯信号时,采用 (/t)的选择将无法得到满意分离,Lee等 在传统Infomax基础上提出扩展的Infomax算法。在信号分离过程中,根据信号峰度变化,转换动态模型系数。分离矩阵调整公式为:W( 1)W(n) [,-Ktanh(u)M - ]W(n)(5)式中:Kdiag[sign( (u。)),,sign(k ( ))],k (·)为信号峭度,超高斯情况 k >0,亚高斯情况k <0,输出信号 i峭度可计算如下:k (U )E[sech ( )]E(u )-E[tanh( )-M ](6)式中:E(·)表示期望,处理过程中可用样本均值进行估计。

2 欠定ICA算法对欠定盲分离问题,大多数算法均在假设源信号稀疏前提下采用聚类或优化算法。该算法对不稀疏信号无能为力;对亚高斯信号缺乏有效欠定盲分离方法。

本文通过生成隐藏数据 ,将欠定盲分离问题转化为传统超定盲分离问题,从而应用传统的ICA算法进行分离。隐藏数据 z的生成通过最大化已知混叠信号和解混矩阵 w条件概率实现,有:arg maxlogp(z I ,W).drg l0g arg maxlogp (Wy)I detWlarg m.ax∑log P 。w y) (7)生成隐藏数据z后,源信号分离问题转化为传统盲源分离问题,如图2所示,新的混叠信号为:y: (8)式中:Y∈ , ∈ 。

图2 欠定 ICA算法流程图Fig.2 Block diagram of underdetermined ICA为较好估计隐藏数据 ,需较准确估计源信号的概率密度。概率密度估计方法分参数估计与非参数估计。参数估计方法,需假定作为判别依据的、随机取值的数据样本在各个可能的类别中均服从特定分布经验及理论说明,参数模型的基本假定与实际物理模型间常存在较大差距,因此该方法并非总能取得满意结果。

鉴于此,Rosenblat等 提出非参数估计方法,即核密度估计方法。由于核密度估计方法不用有关数据分布的先验知识,对数据分布不附加任何假定,从数据样本本身出发研究数据分布特征,因此非常适用于信号盲处理。

给定独立同分布随机变量 的-组实现 , ,, , 的概率密度可估计为:1 n( )÷∑Kh(x-X ) (9) i 1式中:K(x)称为概率密度函数估计核函数,要求K(x)满足条件:(1)K( )对称,且p K( )K(- );- ∞(2)J K( )dx1;32 振 动 与 冲 击 2013年第 32卷(3) 1,2,,k-1时有 J K( )dx0,而J-∞r十∞I K( )dx≠0。

J-∞通倡满足此三个条件的连续函数 K( )称为 k阶核函数。在概率密度函数估计的核函数法中,最常用的核函数为高斯核函数 (b( ) e 。

√基于扩展 Infomax的欠定 ICA算法为:(1)初始化分离矩阵 W;(2)给定观测信号 与分离矩阵w后,据式(10)及当前估计源信号 u的核密度估计式(12)由极大似然估计法生成隐藏数据 z;(3) 由扩展 Infomax算法,式 (5)更新分离矩阵2O广-------------------------10 .H- 叫 -20 - - -- -- - - 垂。

Oj翌馨 .0乏 0j型·00点数 ×102/n(a)三个源信号w,式(2)估计当前源信号 ;(4)返回步骤(1)直到算法收敛。

3 仿真信号分析选-标准的 Laplace信号作为仿真超高斯源信号 51(f);选 S2(t):sin(407rt)sin(600wt),S (t)sin(200a'rt)作为仿真两亚高斯源信号,采样点数为1 200,用 Matlab产生的随机混合矩阵A对源信号进行混合产生混合信号XA [s。;s ;S ],三个源信号与两路混叠信号波形如图 3所示。图4为分离后信号 ,可以看出,除 ICA算法本身固有的分离信号排列顺序及幅值不确定性,信号分离效果较好。

4> 2迪 0坚 -2. 40 2 4 6 8 1O 12点数 ×102/n(b)两路混叠信号图3 源信号与混叠信号Fig.3 Source signals and mixed signals图5 实测振动信号Fig.5 Vibration signals阶次图6 实测振动信号的包络阶次谱图Fig.6 Order envelope spectrumof vibration signals4 齿轮箱复合故障信号分析用某型单级齿轮箱进行试验验证。该系统由电机带动输入轴,输出轴带动负载。主动齿轮齿数Z 30,被动齿轮齿数 Z 50,输入轴轴承为 6206型单列向心滚动轴承,输出轴轴承为 7207型单列向心滚动轴承,通过在齿轮啮合面切去-定厚度材料模拟齿面磨损故障。轴承故障是在输入轴轴承内圈和滚动体通过线切割实现,对振动信号在时域里进行等时间间隔采样,输MM 州- - - - -- 2r-- 。。。。。。。。。。。。 ------- 。 。。。。。。。。 1馨 -o5 - v1v v vv-v wl2.: !!图4 分离信号Fig.4 Separated signals图7 分离信号时域图Fig.7 Separated signals入轴转速由0均匀加速 1 200 r/min,对加速过程的振动信号进行分离实验。内圈故障特征阶次为 X 5.42阶,外圈故障特征阶次为X 3.58阶,滚动体故障特征阶次为 2.35阶。

由于机械结构振动信号通过结构传递产生延迟,机械系统瞬态特征不会保持不变。但当考虑~刚度很大、位移很小的结构时,其振动传递延迟可忽略不计,可保持瞬态模型。为此,认为齿轮箱振动信号满足线性瞬时混合模型。

5 O 5 2 O 2 2 O 2 O O O O O O . A,覃梁 -第 7期 杨 杰等:基于独立分量分析的欠定盲源分离方法 33在输入轴两端轴承座上各安装-个型号为 B&K4508加速度传感器,两个传感器采集的振动信号时域如图5所示,随转速的提高,信号幅值逐渐增大,此为非稳态过程。由于阶次包络谱 能有效提取旋转机械升降速过程中的故障特征阶次,因此本文采用阶次包络谱方法进行分析。实测振动信号阶次包络谱如图6所示,从图中只能看到轴承内圈故障的特征阶次,而外圈及滚动体故障特征却因噪声及混合故障存在变得不明显。

采用本文所提欠定ICA算法对实测振动信号进行处理,分离信号时域图如图7所示,从时域上并不能直观看出分离效果,图8为分离信号包络阶次谱图,从图8(a)、(b)、(c)三幅图中分别可看出内圈、外圈和滚动体故障的特征阶次,说明欠定 ICA算法可将三种故障成功分离。

萼 工 s: -- ]i5 ]: 二0 2 4 6 8 l0 l2 l4 l6 l8 2O阶次图8 分离信号的包络阶次谱图Fig.8 Order envelope spectrum of separated signals5 结 论将独立分量分析引入欠定盲分离问题中,解决了欠定情形下超高斯与亚高斯混合信号的分离问题,为欠定盲分离提供了新方法,同时扩展了独立分量分析的应用范围。通过对实测齿轮箱故障信号分析,表明该方法可有效分离齿轮箱故障特征,验证了方法的有效性。

正在加载...请等待或刷新页面...
发表评论
验证码 验证码加载失败