地震工程作业：求解地震波P、P及卓越频率.pdfVIP

工程第一次作业

结32吕2013010203

1分别求两条波的PGA与PGV

1.1求解PGA

利用Matlab编写程序分别求解ElCentro波与San波的PGA。绘制二者加速度时程曲

线，并利用最大值函数得到PGA。

图1Elcentro波加速度时程曲线图2San波加速度时程曲线

由运行结果可得：ElCentro波PGA=341cm/s22

；San波PGA=185.87cm/s

1.2求解PGV

利用matlab编写代码，利用trapz函数积分得到速度时程曲线。

程序代码

San波速度时程曲线，PGV=-49.433cm/s

ElCentro波速度时程曲线PGV=38.18cm/s

注意到ElCentro波速度时程曲线存在基线漂移现象，即结束时速度不为0。这主要因

为数据时强震仪记录的信息不仅是时纯粹的地面信息，还包含复杂的噪音，

其中的低频噪音会导致加速度时程出现基线漂移，积分求速度时，基线漂移进一步扩大。

因此，需要对Elcentro波进行修正。采用最小二乘法进行时域优化，修正加速度基线。

查阅有关文献[1]，四次多项式拟合程度较好。故选用4次多项式作为拟合均值线。

程序代码：

2

修正后ElCentro波速度时程曲线，PGA=339.74cm/sPGV=33.963cm/s

注意到尽管最小二乘法并不能得到精确速度时程，但基线漂移现象已得到很好的修正。

2分别求两条波的卓越频率

ElCentro波与San波反应谱

图3ElCentro波反应谱图4San波反应谱

Elcentro波卓越频率：1.468Hz

San波卓越频率：0.3479Hz

3分别求两条波的持时

按抗震规范规定，有效持时：一般从首次达到该时程曲线最大峰值的10%那一点算起，到

最后达到最大峰值的10%为止。

ElCentro波：PGA=38.18cm/s，开始时间0.86s，结束时间30.22s，持时29.36s

San波：PGA=185.87cm/s2，开始时间1.9s，结束时间326.73s，持时324.83s

参考文献

[1]陈灯红,彭刚,姚艳华,张微微.波时域数值优化研究及应用[J].世界