非平稳随机过程功率谱密度估计的小波方法 Power spectrum estimation of non-stationary processes via wavelet.pdfVIP

第26卷第3期 振 动 工 程 学 报 V01．26No．3 2013年6月 ofVibration Jun．2013 Journal Engineering 非平稳随机过程功率谱密度估计的小波方法 孔 凡1，李 杰1’2 (1．同济大学土木工程学院，上海200092；2．同济大学土木工程防灾国家重点实验室，上海200092) 摘要：讨论了已有文献中基于一般非正交小波以及广义谐和小波的非平稳随机过程演变功率谱密度(Evolutionary Power Sta— SpectralDensity，EPSD)估计的问题。在一种新的非平稳随机过程模型(局部平稳小波过程，Locally wavelet tionary Process，Lsw)的基础上，提出了一种新的估计非平稳随机过程时变功率谱密度的方法。所建议 的新方法能与估计非平稳随机过程EPSD的经典方法统一起来，当以上两种方法均使用广义谐和小波时，二者退 化为同一形式。为了验证所建议方法的有效性，给出了基于广义谐和小波的多变量均匀调制下非平稳随机地震动 互／自功率谱估计的算例。并以汶川8．o级地震中某近场地及远场地上的地震加速度为例，计算得到了其能量在 时一频域上的不同分布。 关键词：随机过程；非平稳；功率谱密度；地震动；小波分析 中图分类号：U461．5+6；TU352．1十1 文献标识码：A 文章编号：1004—4523(2013)03—0418—11 引 言 用了这一方法估计由非平稳过程谱表达方法合成的 非平稳随机过程的功率谱密度。值得指出的是，Ba— 工程结构在役期间，会受到各种非平稳随机动 su关系式只适用于频域内互不重叠的小波族n引。 力的作用，如地震、强风以及海浪等。其非平稳性不 Spanos等人发展了利用正交且频域内不重叠的谐 仅表现在依赖于时间的幅值上，而且其频谱特性也 和小波(Harmonic Harmonic 是随时间变化的[1]。由平稳随机过程经过确定性函 波(GeneralizedWave【et，GHW)估计非 数调制所得的非平稳过程并不能合理反映频域非平 平稳随机过程功率谱密度的方法(简称spanos— 稳性质o]。获取非平稳随机过程的时间一频率特性， 成为描述随机动力激励、计算结构随机动力响应的 非平稳随机过程[27’2 重要一环。 一般非正交小波(简称为“一般小波”)估计非平稳随 最早的频率分析工具为Fourier变换：它在频 域内分辨率可视为6函数，没有时间分辨率口“]，或 方法)[2明；Huang等人把这种方法推广于多维随机 者精确地讲，其时间信息湮灭在相位信息之中口]。 为了改进Fourier变换的不足，人们提出了Gabor 等人采用这种方法计算了大跨桥梁的随机非平稳响 TimeFourier 展开／短时Fourier变换(Short Transform，STFT)以及Wigne卜ville分布(Wign—非平稳随机过程小波谱的思想[31|，并且证明此方法 e卜Ville Distribution，WVD)方法口~9]。这些方法 能推广应用于一般非Crame卜wold模型的非平稳 虽然能体现非平稳过程的某些时间一频率特征，但它