论文部分内容阅读
工程中设备运行状态千变万化,存在着大量的非平稳动态信号,机械设备在运行过程中当故障发生或发展时也将导致动态信号非平稳性的出现,因此,非平稳性可表征某些故障的存在。对于机械设备运行中的非平稳信号,其频谱特性是随时间而变化的,单纯的时域或频域分析不能够充分描述非平稳信号,而时频分析提供了信号的频谱内容随时间变化的信息,是非平稳信号分析的有力工具。本文结合国家自然科学基金(No.50675209)、河南省杰出人才创新基金(No.0621000500)和河南省重大科技攻关资助项目(No.0122022000),针对目前时频分析方法的局限,基于不同的构造方法对时频分辨率的改善进行了深入研究。其主要内容包括以下几方面:1.本文在对二次型时频分布及自适应核进行深入研究的基础上,分析了固定核函数的缺陷,基于信号的自项、交叉项在模糊域的特点,以及模糊域和时频域的关系,提出了自适应抛物线核时频分析方法,给出了其定义和算法。通过对几种不同类型的仿真信号进行分析,并与其他时频分布结果进行对比,表明本文提出的自适应抛物线核时频分析方法克服了固定核函数的缺陷,使得核函数能够基于信号特征而自适应变化,避免了0/1截断自适应矩形核的振荡缺陷,减小了核通带和阻带之间的距离,从而实现了更高的时频分辨率。自适应抛物线核时频分析方法由于核函数对信号自项的选通作用,不仅较好地抑制了交叉项,同时对噪声起到了很好的抑制效果。自适应抛物线核时频分析方法由于对信号不作任何要求和假设,这种方法所能适应的信号范围更广泛。2.利用谐波小波具有更普遍意义上的正交性、优异的时频分解能力和消除频带重迭与能量泄漏的优点,本文提出了谐波小波包分解方法,并给出了其具体算法。然后,根据匹配追踪原理,在谐波小波包分解的基础上,应用最大投影分解理论,进一步提出了基于谐波小波包的自适应分解算法,根据谐波小波包分解得到的原子字典,将信号分解为一系列谐波小波基原子的线性组合。谐波小波包自适应分解方法可以完全消除二次型时频分布中无法避免的交叉项的干扰,充分利用谐波小波理想的频带滤波特性,消除了能量的泄漏,从而消除了小波包分解及小波包匹配追踪中经常出现的虚假频率。谐波小波包自适应分解方法根据信号结构特点,在谐波小波包原子字典中自适应选择与其局部结构最佳匹配的时频原子,避免了小波包分解中信号信息被冲淡的缺陷,从而使得对信号的描述更为准确、细致,因而提高了时频分辨率。3.为了更好的描述工程中大量存在的非平稳信号的时变特征,研究了基于Chirplet的自适应分解方法。基于原子字典,使用匹配追踪方法进行分解的Chirplet方法,由于随着参数类型的增加,字典中所含的原子个数将大为增加,字典将变得庞大,若信号不能完全与字典中的原子相吻合时,则需要字典中的多个原子进行拟合,这样就会将信号的特征冲淡,因而不能更好地表征信号,此外,在低信噪比的情况下,因为噪声的影响也使得信号特征提取不好。由于Chirplet分解方法的关键是基函数参数的选择,本文针对自适应Chirplet时频分析,提出了一种Chirplet参数估计的算法来自适应分解信号,自适应设计最佳的基函数。首先研究了分数阶傅立叶变换的算法,利用分数阶傅立叶变换进行信号分量的调频率搜索,这种调频率估计与Randon-Wigner变换相比,既减少了搜索量,又可以降低噪声对估计值的影响。虽然基于分数阶Fourier变换的信号调频率检测方法在检测单分量信号和强度相近的多分量信号时很有效,但在检测强度相差较大的多分量信号时,则可能存在强分量掩盖弱分量的问题。而本文的自适应分解方法每分解出一个信号分量,则从残余信号中减去该分量,利用这种逐次消去思想,正好解决了检测信号时强分量掩盖弱分量的问题,因而得到了更好的估计精度。在此基础上,提出构造函数对信号进行变换,利用短时傅立叶变换进行时间中心、频率中心和尺度因子的估计,利用拟牛顿和期望最大化(EM)方法进一步提高参数估计的精度,本文将EM方法的原理引入匹配追踪的参数估计,考虑了完备数据中噪声的存在,使得每次分解结果收敛于全局最大值,因而得到了更好的参数估计。仿真试验也表明,当信号具有较低的信噪比时,本文方法仍然给出了分辨率良好的时频分布,体现出了本文方法的鲁棒性。本文提出的自适应Chirplet参数估计及信号分解方法保证了每一个基函数都是按照与信号最匹配的最大投影分解原则来选择,因而能够以更少的基函数、更稀疏而精确的描述信号,准确地刻画信号的时频特征,不仅可以得到良好的时频分布,更为故障诊断中的智能诊断提供了极为有利的条件。当信号类型与基函数类型相同时,信号分解次数与信号分量个数相同就能完全将信号有效的表示。当信号类型与基函数类型不相同时,信号分量可以用多个基函数来逐段近似,因而仍能得到反映信号特征的时间-频率特性。4.从参数化建模的角度出发,研究了用系数参量表征非平稳信号的时变参数模型方法,通过将模型的时变参数展开在一系列随时间变化的时间基上的方法,将时变参数的估计转化为时不变参数的估计,用带遗忘因子的递推最小二乘法进行参数估计,从而实现了时变信号的自适应时频分析。提出了小波基时间函数的TVAR自适应时频分析方法,研究了该分析方法的理论,推导了其算法。给出了其计算步骤。用仿真算例进行了TVAR自适应时频分析方法的可行性验证及其性能研究,结果显示,基于TVAR建模的自适应时频分析方法对时变信号的跟踪及时频分辨率均优于传统的时频分析方法,与假设在一段时间间隔上信号是平稳的参数估计方法相比,时变参数模型法可以进一步提高参数估计的精度和对时变信号的跟踪能力。基于TVAR建模的自适应时频分析方法可以克服时频分布大数据量的问题,参数化建模时频分析又避免了分辨率受采样频率和采样点数的限制,具有分辨率高的优点。基于TVAR的时频分析方法对时间基函数有一定的依赖作用,不同的基函数对时频分布的结果有一定影响,文中提出的小波基函数由于具有良好的时频聚集能力,因而比傅立叶基、勒让德基等对信号的适应能力更广,对信号的拟合性能更优。采用递推最小二乘算法避免了矩阵的求逆运算,从而避开了病态矩阵对算法结果的影响,同时由于观测数据的递推进行,从而节约了存储量和运算时间,更利于适时、在线的进行信号分析与处理。但是,基于TVAR的时频分析方法由于在参数估计过程中以误差平方最小为目标,当信号染噪时,这种方法不可避免地受到较大的影响,因而这种方法的抗噪性能虽然比常用的Wigner分布好,但总体来说对染噪信号分析的鲁棒性有待于进一步提高。5.由于旋转机械的绝大多数故障都与转子及其组件(齿轮、轴承)直接相关,其它位置的故障相对较少,而滚动轴承属于故障易发性零件,并且轴承出现故障时,其信号属于典型的非平稳信号,因而本文进行了试验的构建,获取了转子和滚动轴承的故障实验信号。应用本文所提出的几种自适应时频分析方法对故障信号进行了分析,获得了其时频分布特征。与常用时频分布相比,本文所提出的自适应时频分析可以更好地刻画信号的时频特征,尤其是自适应抛物线核时频分析和自适应Chirplet分解方法不仅具有良好的时频分辨率,而且对信噪比较低的信号仍能较好地提取故障特征,显示了其优良的抗噪性能。根据滚动轴承信号的实验结果对几种自适应时频分析方法进行了分析对比,指出了各自的特点和应用范围,验证了本文算法的有效性。