地震动瞬时谱估计的Unscented Kalman滤波方法.pdfVIP

维普资讯 应用数学和力学，第28卷 第 l1期 2007年 11月 15日出版 Ap~iedMathematicsandMechanics Vo1．28，No．11，Nov．15，2007 文章编号：100042887(2007)11-1370-09 @应用数学和力学编委会，ISSN1000-0887 地震动瞬时谱估计的UnscentedKalman 滤 波 方 法 李英民 ，董银峰 ，赖 明。 (1．重庆大学 土木工程学院，重庆400045； 2．重庆大学 山地城镇建设与新技术教育部重点实验室，重庆 400045， 3．建设部 科技司，北京 100835) (陈山林推荐) 摘要 ： 用时变ARMA模型描述地震动时程，提出了采用UnscentedKalman滤波技术实现地震动瞬 时谱估计的思路．算例分析表明，UnscentedKalman滤波方法较Kalman滤波方法适用范围广，具有 较高的时间和频率分辨率，能够更好地跟踪地震动的局部特性，适合处理非线性模型或有突变特 性的模型的辨识问题．不同阶数ARMA模型的估计结果还表明，以往被忽略的ARMA模型的理论 频率分辨力对地震动瞬时谱估计精度有重要影响，应作为一个参考指标在ARMA模型的判阶中加 以考虑． 关 键 词：UnscentedKalman滤波； 地震动； 瞬时谱； 时变ARMA模型 中图分类号： 0211，64；P315．3 文献标识码： A 引 言 近年来的研究表明，用来描述地震动工程特性三要素的传统参数并不总是足以构成工程 结构的控制地面运动特征，而其它一些特性(如强度和频率非平稳特性)对于结构地震反应特 别是非线性反应在一定情况下同样具有不可忽略的影响．对于这些特性，研究者从不同的角 度相继提出了一些指标来描述，如用穿零率描述频率非平稳特性，用瞬时谱同时从时域和频域 描述强度和频率非平稳特性等 J．其中，瞬时谱是描述强度和频率非平稳特性非常有效的概 念，可以直观地反映地震动能量在时域和频域的分布．因而，地震动的瞬时谱估计成为 目前 地震动工程特性研究的一个重要方向，瞬时谱的精确估计及其衰减规律研究是这一研究的难 点． 概括研究者所提出的各种地震动瞬时谱估计方法，基本可分为非参数化方法和参数化方 法两类．非参数化方法主要有时频分析(分布)方法 (如短时傅氏变换、谱图、Wigner-Ville分 布、Choi—Williams分布以及基于高阶累积量的双谱或多谱等高阶谱分析方法)、演变谱方法、小 收稿 日期： 2006-09．20；修订 日期： 200709．04 基金项 目： 国家 自然科学基金资助项 目；教育部新世纪优秀人才支持计划资助项 目 (NC口 O6．Or765) 作者简介： 李英民(1968一)，男，山东人 ，教授，博士生导师 (联系人 ．rel：+86-23 E-mail：liyingmin@cqu．edu．cn)． 】370 维普资讯 李 英 民 董 银 峰 赖 明 波变换方法以及HHT方法等ll4；参数化方法主要指基于 ARMA模型的现代谱估计方法．参 数化方法因分辨率高、谱曲线光滑、便于谱特性模型化及衰减规律研究等优点，受到更多的重 视 [1，51． 参数化方法一般包括以下4个步骤：1)地震动的时变ARMA模型建模；2)时变 ARMA模 型的状态．空间模型表示；3)通过状态．空间模型的估计得到时变ARMA模型的参数；4)根据 时变ARMA模型的参数计算地震动瞬时谱．其中，状态．空间模型