第40卷第4期 昆明理工大学学报(自然科学版) Vo1.40 No.4 2015年8月 Journal of Kunming University of Science and Technology(Natural Science Edition) Aug.2015 doi:10.16112/j.cnki.53—1223/n.2015.04.009 基于压缩感知的滚动轴承振动信号压缩方法 刘 畅 ,伍 星 ,毛剑琳 ,柳小勤 (1.昆明理工大学机电工程学院,云南昆明650500;2.昆明理工大学信息工程与自动化学院,云南昆明650500) 摘要:压缩感知是一种基于信号稀疏性的信号采集与处理框架,能够在信号采集的同时对信号 进行压缩.本文提出一种基于压缩感知的滚动轴承振动信号压缩方法,在信号的变换域内,使用 幅值百分比作为阈值对变换系数进行稀疏处理,采用高斯随机矩阵对信号进行降维观测,实现数 据的压缩.通过研究稀疏处理对信号压缩比和逼近误差的影响,分析阈值选择与信号稀疏比和逼 近误差的关系,分析不同变换基对稀疏处理的影响.实验数据分析表明,相对于未经过稀疏处理 的信号来说,该方法能有效地提高信号的压缩效果,且保持较好的逼近误差. 关键词:压缩感知;滚动轴承;稀疏表示;信号压缩 中图分类号:TH133.33 文献标志码:A 文章编号:1007—855X(2015)04—0046—05 Rolling Bearing Signal Compression Using Compressive Sensing LI U Chang ,WU Xing ,MAO Jian—lin ,LI U Xiao—qin (1.Faculty of Mechanical and Electircal Engineering,Kunming University of Science and Technology,Kunming 650093,China; 2.Faculty of Information Engineering and Automation,Kunming University of Science and Technoloyg,Kunming 650093,China) Abstract:Compressed sensing is a signal acquisition and processing method based on signal sparsity.This paper presents a compression method of vibration signal of the rolling bearing based on compressed sensing.In trans— form domain,the amplitude as a percentage threshold for sparse treatment is used to process the transform eoefif— cients.The Gauss random matrix is used to measure signal and compress data.The relationship between sparse signal compression ratio and approximation error is studied.Different thresholds are then selected to examine the relationship between sparse signal ratio and the approximation error.The influence of different bases on sparsity is also analyzed.The resuhs show that this method can improve the compression effect,and reduce the approxi— mation error. Key words:compressive sensing;rolling bearing;sparse representation;signal compression O引言 在滚动轴承的状态监测与故障诊断中,振动信号采集都是采用基于Nyquist—shannon采样定理方法, 即通过预先估计信号中最高频率成分,或使用抗混滤波器来限定信号的最高频率成分,根据Nyquist—shan— non采样定理确定合适的采样频率后进行数据采集.这种采集方式会带来两个方面的问题:一是由于轴承 的特征频率出现在中高频段,为满足采样定理所需要的采样频率就会很高,这就对数据采集设备提出较高 的要求;另一个问题在保证同样的分辨率的情况下,采样频率越高,采样得到的数据量就会越大,这样就对 采样数据的传输、存储带来巨大的压力.能否找到一种方法,既能够保证采集到足够的信息用于信号分析, 收稿日期:2014—03—10.基金项目:国家自然科学基金项目(51265018). 作者简介:刘畅(1979一),男,博士研究生,工程师.主要研究方向:机械设备状态监测与故障诊断技术、智能诊断、压缩 感知.E-mail:lxD385@163.com 通信作者:伍星(1973一),男,博士,教授,博士生导师.主要研究方向:机械信号处理与故障诊断、振动分析等.E-mail: xingwu@aliyun.com 第4期 刘畅,伍星,毛剑琳,等:基于压缩感知的滚动轴承振动信号压缩方法 47 又能够减少采集的数据量,利于采样数据的存储和传输.压缩感知理论为本问题提供了一个新的思路. 压缩感知理论由Cand ̄s,Donoho等人于2006年提出 ,该理论是建立在矩阵分析、泛函分析、稀疏 表示、最优重构等理论基础上的一种全新的理论框架.文献[1]给出了压缩感知的描述:当信号在某个变 换域是稀疏或可压缩的,可以利用与变换矩阵非相干的测量矩阵将信号投影到一个低维空间.由于这种投 影保持了重建信号所需要的信息,那么通过求解稀疏最优化问题就能够从低维的观测数据中高概率地重 构出原始信号.压缩感知理论在提出后不久,就在雷达探测 J、医疗成像 j、无线通信 j、数据采集 、语 音信号处理 等领域开展研究. 本文提出一种基于压缩感知的滚动轴承振动信号压缩方法,在信号的变换域内,使用幅值百分比作为 阈值对变换系数进行处理,采用高斯随机矩阵对信号进行降维观测,实现数据的压缩.相对于未经过处理 的信号来说,该方法能有效地提高信号的压缩效果,且同时获得较低的稀疏比和较好的逼近误差. 1压缩感知基本理论 压缩感知基本理论是:设长度为n的信号 在某个正交基或紧框架 上的变换系数 是稀疏的,如果 用一个与变换基 不相关的m×n维观测基 (m<n)对信号 进行线性变换,得到m×1维观测集合Y, 利用最优化求解方法从观测集合中精确或近似精确的重构出原始信号_9 J. 设{ 是R 空间的正交基向量,由其作为列向量组成正交基矩阵 =【 , ,…, 】.在R 空间中 任意向量 都可以表示为: =∑ i 或 =伽 (1) 式中: = ’ 是由投影系数 =x, )构成的n×1维列向量.如果 ∈R 中仅仅有k个非零项,且 k<<n,则表明信号是稀疏的或可压缩的,那么就可以用k个大系数较好的逼近信号 ,称 是k一稀疏, 或 在基 上是k项稀疏 . 对信号进行观测得到降维的观测集合Y∈R ,观测过程可以表述为: Y= = 蝴 (2) 由于m<n,该问题是一个欠定问题.文献[10]指出在满足信号 是稀疏或可压缩的前提下,且观测 矩阵 满足RIP条件,可以通过求解一个最优化的z。范数问题来重构原始信号 . 信号的重构是指由m×1维观测向量Y中重构出n×1维信号的过程,重构问题可以由下式描述: min l I1 s.t.Y= 蜘 (3) 通过求解z,范数最小问题来重构信号¨ .由于式(3)是凸优化问题,可以很方便地转化成线性规划 问题进行求解. 2基于压缩感知的滚动轴承振动信号压缩 2.1信号压缩与解压缩的基本流程 基于压缩感知的信号压缩本质是利用信号在变换域上的稀疏性,基于信号中的信息进行采样,以 远低于Nyquist—shannon采样频率的速率实现采样的同时进行压缩,得到降维的观测数据(压缩数据). 然后通过求解最优化问题,从少量的观测数据中精确地重构出原始信号(数据解压).该方法的基本流 程如图1所示. 变换域 稀疏系数 观 阵 重构算法 图1基于压缩感知的信号压缩与解压缩流程 Fig.1 Signal compression and decompression process based on compressed sensing 48 昆明理工大学学报(自然科学版) 第40卷 对滚动轴承振动信号应用压缩感知方法进行压缩的前提是振动信号必须是稀疏或可压缩的,且信号 越稀疏,压缩效果越好.文献[12]指出光滑信号的Fourier系数、小波系数等都具有足够的稀疏性.此外离 散余弦变换DCT具有实信号的变换系数仍是实信号,且具有非自适应性质,也是一种较好的稀疏表示方 法.但由于振动信号是非严格稀疏信号,因此在变换域内存在大量的非零、小值出现,导致信号重构得到的 恢复信号与原始信号存在一定的重构误差. 信号压缩在压缩感知的观测过程中进行.文献[13]证明当观测矩阵 是随机矩阵时, 能以较大 概率满足RIP条件,即观测矩阵能够与大多数的变换基不相关,以保证能够精确地重构出原始信号.文献 [14]给出了在压缩感知框架下精确重构信号观测值需要满足的条件. m≥C・ ( .qt)・k・log凡 (4) 公式(4)同时给出了观测数m与信号稀疏度k之间的关系,即信号稀疏程度越高,观测时所需要的观 测值越少,信号的压缩效果越好. 在基于压缩感知的信号压缩中,信号的重构将从观测信号(压缩信号)中重构出原始信号(信号解 压缩).常用的重构算法包括OMP算法、STOMP算法、BP算法、GPSR算法、ISTA算法等.其中以OMP为 代表的贪婪类算法具有计算速度快的特点,以BP为代表的凸松弛类算法则重构精度高,但相对来说计 算复杂度也较高. 2.2变换域内阈值处理 针对轴承振动信号在变换域内非严格稀疏的特点,通过在变换域内对变换系数进行阈值处理,将 小于阈值的非零小值置零,提高变换系数的稀疏化程度,能够减少信号的观测数,进而提高信号的压缩 率,这种处理方法称为稀疏处理.在稀疏处理过程中,如果阈值选择不当,将导致信号中关键信息损失, 从而导致重构信号的误差增大,使得信号的压缩失去意义.因此阈值的选择对于信号的稀疏度和重构 信号的误差有很大的影响,阈值取值越大,丢失的信息越多,导致重构信号的误差越大;而阈值取值越 小,信号的稀疏度也就越大,但是信号的压缩效果也就越差.此外使用不同的变换基,得到的信号稀疏 度和压缩效果也不同. 本文分析不同阈值对信号稀疏度和重构误差的影响,在变换域内根据变换系数的幅值百分比来自适 应地确定阈值,根据该阈值进行稀疏化处理,以期能够在稀疏度和重构误差之间保持较好的平衡.设有n ×1维信号 ,对 进行变换得到0= ’ , 是n x 维变换基,0是 在 基上的变换系数,确定阈值8: s=A I max( )一min(0)l (5) 式中A为幅值百分比.为了验证本方法的可行性,在不同的变换基下,对稀疏处理和未经稀疏处理的信号 压缩率和逼近误差进行对比分析,并得到阈值与信号稀疏比和逼近误差的关系曲线. 3滚动轴承振动数据的分析与验证 本文使用的实验数据来自凯斯西储大学轴承数据中心的滚动轴承振动数据 ,轴承型号为SKF 6205-2RS型深沟球轴承,转速为1 797 r/min,采样频率为12 kHz.观测矩阵选择高斯随机测量矩阵,采用 OMP进行信号重构. 为了更好地对压缩过程和压缩效果进行表述,定义稀疏比 =k/n,凡为信号元素的个数(信号长度), k为信号中非零项的个数,显然信号中非零元素的数目越少,信号稀疏程度越高,稀疏比越小;定义压缩比 C=m/n,m为信号观测所需要的观测值数目,显然观测数越少,压缩比越小,信号的压缩效果就越好;定义 逼近误差 =_ lXk— ll/ ll用来衡量重构信号与原始信号的差异,逼近误差越小,说明重构信号与原 始信号的差异越小,重构信号越能逼近或代替原始信号. 3.1稀疏处理与压缩比与逼近误差的关系 为了验证稀疏处理方法能够有效地提高信号的压缩效果,首先对比分析稀疏处理和未经稀疏处 理的信号逼近误差.选择变换基为FFT,数据长度1 024,在保证观测数满足m一4k的条件下,分析压 缩比与逼近误差err的关系.图2是压缩比与逼近误差err的关系曲线.显然,在相同的压缩比下,稀疏 第4期 刘畅,伍星,毛剑琳,等:基于压缩感知的滚动轴承振动信号压缩方法 49 处理后的逼近误差比未经过稀疏处理的信号小,表明对信号进行稀疏处理能够有效地提高信号的压 缩效果. 3.2幅值百分比A与稀疏比和逼近误差的关系 阈值的选择对稀疏比和逼近误差有直接的 影响.根据公式(5),A与占成线性对应关系,因 此下文使用A代替对 的分析.取数据长度为 4096,采用FFT作为变换基,分析A与稀疏比 和逼近误差err的关系.图3是A与稀疏比OL和逼 近误差err的关系曲线.A与稀疏比 成近似指数 对应关系,提高A值能够大幅提高信号的稀疏 度,但当超过某个值后(如A≥0.1),提高A值对 于信号稀疏度的提高作用就不再明显;§ 而A与逼 近误差err成近似线性关系,提高A值直接导致 压缩比 逼近误差成线性增加.因此选择合适的A值对权 图2压缩比与逼近误差的关系曲线 衡稀疏比和逼近误差相当重要.太大的A值将导 Fig・2 The curve of compress ratio and approximation error 致逼近误差的增加,而A值太小时,信号的稀疏 比较大,信号的压缩效果不明显.通过实验方法,当A≥0.1能够保证既能够得到较大的稀疏比,同时也能 保证较好的压缩效果. 3.3不同变换基下 与稀疏比和逼近误差的关系 使用不同的变换基得到的变换系数具有不同的稀疏程度.使用本文提出的方法对比分析使用FFT、 DCT和DWT(采用8阶Daubechies离散正交小波) 基时A与稀疏比 和逼近误差err的关系.图4为不 同变换基下A与稀疏比Ol和逼近误差err的关系曲线.在相同的阈值下,使用DCT得到的稀疏比最低,但 是逼近误差也最高,而FFT和DWT在逼近误差上相差不大,但是FFT具有更好的稀疏性,因此针对滚动 轴承的振动信号,在相同的阈值下,使用FFT作为变换基能够同时保持较好的稀疏度和较低的逼近 误差. 0.8 0.7 0.2 0.1 O 幅值百分比/2 幅值百分比/f (a)幅值百分比 与稀疏比 的关系 (a)幅值百分比 与稀疏比 的关系 0 0 j}{11 0 0 0 0 0 0 幅值百分比/2 幅值百分比 (b)幅值百分比/2与逼近误差P 的关系 (b)幅值百分比 与逼近误差P 的关系 图3 与稀疏比和逼近误差的关系 图4不同变换基] 与稀疏比和逼近误差的关系 Fig.3 The curve of2,sparsity and approximation error Fig.4 The curve of2,sparsity and approximation error under different bases 昆明理工大学学报(自然科学版) 第4O卷 4结论 ’本文提出一种基于压缩感知的滚动轴承振动信号压缩方法,在信号的变换域内,使用幅值百分比作为 阈值对变换系数进行稀疏处理,采用高斯随机矩阵对信号进行降维观测,实现数据的压缩.研究稀疏处理 对信号压缩比和逼近误差的影响,分析阈值选择与信号稀疏比和逼近误差的关系,分析不同变换基对稀疏 处理的影响.实验数据分析表明,相对于未经过稀疏处理的信号来说,该方法能有效地提高信号的压缩效 果,且保持较好的逼近误差.但是稀疏处理增加信号处理的步骤,且对于非严格稀疏的振动信号来说,要想 同时获取较低的稀疏程度和较好的逼近误差,还需要进一步开展研究. 参考文献: [1]David L Donoho.Compressed sensing[J].Information Theory,IEEE Transactions,2006,52(4):1289—1306. [2]Emamnuel J Candes.Compressive sampling[R].Proceedings on the International Congress of Mathematicians:Madrid,Au— gust 22—30,2006:invited lectures. [3]Emmanuel Candes,Justin Romberg,Terence Tao.Robust uncertainty principles:Exact signal reconstruction from highly in— complete frequency ifnormation[J].Information Theory,IEEE Transactions,2006,52(2):489—509. [4]Baraniuk R,Steeghs P.Compressive radar imaging[R].IEEE Radar Cofnerence,Waltham,Massachusetts,April 2007. [5]Lustig M,Donoho D L,Pauly J M.Sparse MRI:The application of compressed sensing for rapid MR imaging[J].Magnetic Resonance in Medicine,2007,58(6):1182—1195. [6]Quer G,Masiero R,Munaretto D,etc.On the Interplay Between Routing and Signal Representation for Compressive Sensing in Wireless Sensor Networks[R].Ifnormation Theory and Applications Workshop(ITA 2009),San Diego,CA. [7]Mishali M,Eldar Y C.Xampling:Compressed sensing of naalog signals[M].Cambridge,UK:Cambridge University Press, 2012. [8]Griifn A,Hirvonen T,Tzagkarakis C,et 1a.Single—channel and muhi—channel sinusoidal audio coding using compressed sensing[J].Audio,Speech,and Language Processing,IEEE Transactions on,2011,19(5):1382—1395. [9]石光明,刘丹华,高大化,刘哲,林杰,王良君.压缩感知理论及其研究进展[J].电子学报,2009,37(5):1070—1081. [10]Emamnuel J Cand ̄s.The restircted isometry property and its implications for compressed sensing[J].Comptes Rendus Mathematique,2008,346(9):589—592. [1 1]Scott Shaobing Chen,David L Donoho,Michael A Saunders.Atomic decomposition by basis pursuit[J].SIAM journal on scientific computing,1998,20(1):33—61. [12]Emamnuel J Cand ̄s,Terence Tao.Near—optimal signal recovery rfom random projections:Universal encoding strategies? [J].Ifnormation hTeory,IEEE Transactions,2006,52(12):5406—5425. [13]Emamnuel J Cand ̄s,Justin Romberg,Terence Tao.Stable signal recoveyr from incomplete and inaccurate measurements [J].Communications on pure and applied mathematics,2006,59(8):1207—1223. [14]Emamnuel J Cand ̄s,Justin Romberg.Spamity and incoherence in compressive sampling[J].Inverse problems,2007.23 (3):969. [15]Loparo K.Bearing vibration data set[EB/OL].(2003—10—23).http://csegroups.case.edu/bearingdatacenter/home.