HHT变换在地震波谱分析中应用.pdfVIP

现代地震工程进展 373 HHT变换在地震波谱分析中的应用+ 罗奇峰 石春香 (同济大学结构工程与防灾研究所，上海200092) 摘要：本文对HHT变换的理论及其在地震波谱分析中的应用进行了阐述。说明 它是对非平稳非线性时程进行分析的有效工具。 关键词：HHT变换；固有模态函数分离；Hilbert变换 1概述 近年来发展起来的Hilbert—Huang变换理论(以下简称HHT变换)，是～种适合分析非平 稳、非线性过程的谱分析方法…。本文首先介绍HHT变换，而后以地震波的HHT变换为 例，对比分析HHT变换和傅里叶谱分析方法。 2 HHT变换 HHT变换由Huang变换和Hilbert谱分析组成。 2．1 Huang变换 Huang变换的关键是经验模态分离法，该方法认为任何复杂的时间序列都是由～些相 互不同的、简单的、并非正弦函数的固有模态函数组成，基于此可从复杂的时间序列直接分 离成从高频到低频的若干阶固有模态函数，即基本时间序列。固有模态函数需满足以下2 个条件‘“： 1)对整个时间序列来说，极值的个数与穿过零点的个数相同或其差值为1； 2)在任何一点最大值的包线与最小值的包线的均值为零。 图1为一条原始的地震记录，图2为一个典型的经Huang变换后得到的固有模态函数， m o o 一～＼吡)毯嘲曩 o：o以““矗Ⅲ 田】一条原始的地震记录 ·国家自害!{科学基金资助项目(批准。 374 罗奇峰等：HI-IT变换在地震波谱分析中的应用 从曲线可看出其极值的个数与曲线穿过零点的个数相同，且该时程的上下包线相对零线对 称。固有模态函数看起来像Fourier变换得到的正弦波，但它不是单一频率的波，且每一时≮ 刻的振幅也不相同。对任意一条地震波都可按下述方法进行Huang变换““。 0 4 一 毫 0‘2 意 。 型．02 日一O4 u ) 10 15 删 23 30 35 40 45 50 时间(s) 圈2经Huang变换后得到的一个典型的固有模态函数 一陋—一i 首先，确定地震波x(I)的所有峰值点，然后用三次样条函数曲线循序连接所有的最大 值，得到Ⅳ(t)的上包络线，采用同样的方法循序连接所有的最小值，得到x(￡)的下包络线， 循序连接上、下包络线的均值可得一条均值线叭(￡)，于是，可从下式计算得到h。(￡)： hl(￡)=X(f)一B1,，(t) (1) 如果h，(#)满足固有模态函数所需的条件，则h，(z)即为第一阶固有模态函数。～般地说，它 并不满足条件，此时，将h。(f)看成新的时间曲线，重复上述方法，可得 h1．(￡)=hI(￡)一mll(f) (2) 式中，m1](f)是h．(1)的上、下包络线。按上述方法重复k次，直到h，。(￡)满足固有模态函数 的条件为止，^，。(￡)由下式计算： hIt(￡)=hm一1)(f)一m1(I_(￡) (3) 式中的^。+(￡)即为第一阶固有模态函数，记作c，(t)．即 cI(￡)=五t}(￡) (4) 第～阶固有模态函数c。(￡)包含着时程x(￡)的频率最高的成分。从x(￡)中减去高频成分 c。