HHT法边界问题的处理及其在时变结构动力特性识别中的应用研究.pdfVIP

HHT法边界问题的处理及其在时变结构 动力特性识别中的应用 张永利李鸿晶 (南京工业大学土木工程学院 江苏南京210009) 摘要 在研究基于经验模态分解(EMD)算法的基础上，提出了～种处理EMD边界问题的新方法一 —拟正弦边界延拓法．为检验该方法的适用性和优越性-给出了一个数值算例。最后，将其应用于线 性时变SDOF体系的动力特性识别中，并与采用特征波法的识别结果进行比较．对比表明．应用拟正 弦边界延拓法处理ENID的端点效应具有一定的优越性，值得推广应用．同时，本文也为FIHT在时变 体系动力特性识别中的应用提供一些参考． ， 关键词 经验模态分解；固有模态函数；希尔伯特变换；瞬时频率；线性时变体系 0引言 1998年，美国国家宇航局的N．E．Huang等人提出了一种新的非线性、非平稳信号处理方法， Transform，简称HHT)。该方法由两部分组成：首先是 被称为希尔伯特一黄变换(Hilbert-l-luang Mode Decomposition)将一个时问序列信号分解成有 用基于经验的模态分解方法EMD(Empirical Mode 限个不同时间尺度的固有模态函数IMF(Intrinsic 换构造解析信号，得出信号的瞬时频率和振幅，进而得到能量分布谱图(Hilbert Spectrum)【l’21, Hilbert谱在幅值一频率一时间分布图上对各分量的瞬时频率、瞬时振幅都进行了比较精确的刻画， I-IHT对于处理非线性、非稳态数据有着强大的生命力。HHT一经提出便在各领域得到广泛的应 用。目前，该方法在地震工程中主要应用于地震信号处理14j和结构分析、结构健康诊断、损伤检 测及系统识别p1等领域，并且取得了较好的结果。但由于该方法刚提出不久，还存在一些问题亟 待解决，如包络线拟合、边界问题的处理、筛法及模态混叠等问题，其中边界问题的处理是应用 EMD方法的瓶颈问题。EMD是用极值点拟合信号包络线进行分解的，在信号内部，极值点总是 能获得的，局部均值也就很容易确定；在边界处，拟合包络线需用到边界外的极值点或信号，这 些极值点或信号需要经过分析人为给定，以确定边界处信号的局部均值，这便是EMI)的边界问 题。EMD的每次分解都要经过多次筛选，每次筛选又是用信号的极值点拟合信号的均值曲线，其 边界处理的影响不只局限在边界的有限范围内，会向信号内部传播，这就使得边界问题的处理显 得尤其重要o。目前，研究人员针对EMD的边界问题提出了多种处理方法，如特征波法、边界波 形匹配预测法【6J、“边筛分、边延拓”法口】、极值延拓法及镜像闭合延拓法口】等。特征波法和极值 延拓法这两种处理方法比较简单，有时效果也很好，但没有考虑边界的非完全波；边界波形匹配 ①钟佑明．希尔伯特一黄变换局瞬信号分析理论的研究．重庆：重庆大学博士学位论文，2002 预测法考虑了边界的非完全波，但所采用的具体算法应用范围比较有限，适用于无低频叠加波和 调幅信号幅度相差不大的信号。：“边筛分、边延拓”法和镜像闭合延拓法适用范围较广，效果也 比较好，但比较复杂，编程实现较困难。本文在研究EMD算法的基础上，提出了一种新的边界 处理方法一拟正弦边界延拓法，并通过仿真试验和比较，验证了此方法的可行性及优越性。 1基于经验的模态分解算法 在EMD方法中，Huang将满足如下两个条件的信号定义为固有模态函数(IMF)； (1)极值点数量与过零点数量相等或至多差一个； (2)在任意点，由局部极大值点和局部极小值点分别构成的包络线的均值为0。 这种方法的本质是通过特征时间尺度获得固有振动模态，然后用有限个固有振动模态来分解 时间序列数据。其实现过程如下； (1)找出原始数列z(f)的所有局部极大值点和极小值点，将其用三次样条函数拟合为x(0的 上下包络线； (2)求上下包络线的均值嘲(f)； (3)将原始数列减去rth(0可得到一个去掉低频