当前位置：文档之家 › 基于MATLAB的地震反应谱与傅里叶谱计算分析

基于MATLAB的地震反应谱与傅里叶谱计算分析

  • 格式：pdf
  • 大小：472.41 KB
  • 文档页数：5

下载文档原格式

  / 5

合集下载

基于MATLAB地震反应谱数值算法的稳定性和精度分析

基于MATLAB地震反应谱数值算法的稳定性和精度分析

基于MATLAB地震反应谱数值算法的稳定性和精度分析摘要：地震反应谱是进行结构抗震分析与设计的重要工具，反应谱的计算在反应谱法和时域逐步积分方法中有重要地位，引起了学者的重视和广泛研究。

而对计算方法优劣的评定常取决于其计算的耗时、稳定性和精度等因素。

本文基于数值算法的相关研究及应用现状，以MATLAB为平台，建立数值算法在不同影响因素下的三维图形，并结合理论进行对比分析。

通过算例进一步分析验证，得出不同数值算法在实际计算中的表现，为工程实际计算中选取哪种积分算法更为合适提供参考。

关键词：地震反应谱；时域逐步积分算法；稳定性和精度；MATLAB1、地震反应谱的基本假定地震反应谱基于的三个基本假设[1]：(1)结构物所处的地面假定为刚性面，认为体系各质点的运动是完全一致的。

(2)强震观测仪的记录为地面运动的过程。

(3)结构体系不能是双或多质点体系，必须是单质点体系；同时应是弹性体系状态。

这里所谓的单自由度体系结构，就是用无量刚的弹性杆件支承于地面上，将结构体系中参与振动的质量用一点表示。

同时，假定结构振动和地面运动不发生扭转，只是水平平移运动并且是单方向的。

2、基于MATLAB地震反应谱数值算法的稳定性和精度分析2.1 概述目前MATLAB地震反应谱数值理论算法主要有中心差分法[2]、Wilson-法、Houbolt法、线性加速度法及Newmark-法等，理论算法主要是以求解线性结构体系动力方程时所表现出的特性作为数值算法优劣的评价依据[3]，但是在实际工程运用中，人们常常凭借经验来判定选取较为合适的积分方法。

随着工程问题越来越复杂，在对大型复杂结构的结构动力反应分析更为复杂，要求高效率计算情况下获得较精确地计算结果。

然而各计算方法的精度和稳定性对结构动力反应分析的发————————————E-mail:skyuanyan@展引起了很大的影响和制约[4]。

2.2数值算法的稳定性分析基于上述情况，本文对上述几种常用数值算法的稳定性方面通过图形进行比较分析，并结合算例进一步验证分析。

matlab对地震波进行傅里叶变换

matlab对地震波进行傅里叶变换

matlab对地震波进行傅里叶变换地震波是指地震时由地震源产生的机械波，它在地球内部传播并在地球表面或近表面造成振动。

对于研究地震波的特性和分析其成因机制，傅里叶变换是一种非常重要的数学工具。

在matlab中，我们可以使用fft函数来对地震波进行傅里叶变换。

1. 准备数据首先需要准备一组地震波数据。

这里我们可以使用matlab自带的load函数加载一个示例数据文件，该文件包含了一个从南极到北极的走时曲线：load seismictest.mat;2. 绘制时域图像利用plot函数可以绘制出该走时曲线的时域图像：plot(seismictest);可以看到该图像呈现出明显的周期性振动。

3. 进行傅里叶变换接下来，我们可以使用fft函数对这组数据进行傅里叶变换：Y = fft(seismictest);其中Y为变换后得到的频域信号。

4. 绘制频域图像利用abs函数和fftshift函数可以将频域信号转化为幅度谱，并通过plot函数绘制出频域图像：f = (-length(Y)/2:length(Y)/2-1)/length(Y);Y_shift = fftshift(Y);plot(f, abs(Y_shift));可以看到该图像呈现出多个峰值，这些峰值对应着不同的频率成分。

5. 分析结果通过傅里叶变换，我们可以将地震波信号从时域转化为频域，进而分析地震波的频率成分和振幅。

在上面的例子中，我们可以看到该地震波信号包含了多个频率成分，这些成分对应着不同的振幅。

通过进一步的分析和处理，我们可以更深入地研究地震波的特性和成因机制。

总之，matlab提供了强大的工具来进行地震波信号处理和分析。

通过使用fft函数对地震波进行傅里叶变换，我们可以将时域信号转化为频域信号，并对其进行进一步的分析和处理。

这对于研究地震学和相关领域具有非常重要的意义。

基于MATLAB的地震数据的分析

基于MATLAB的地震数据的分析

基于MATLAB的地震数据的分析地震数据的分析是地震科学研究中的重要环节之一，可以帮助地震学家了解地震的特征、预测地震的趋势以及评估地震的影响程度。

MATLAB作为一种功能强大的数据处理和分析工具，在地震数据分析中也扮演着重要的角色。

本文将介绍基于MATLAB的地震数据分析方法和应用。

首先，地震数据通常是通过地震仪器采集到的地震波形数据，以地震波形数据为基础进行地震分析是地震学研究中的常见方法。

MATLAB提供了丰富的信号处理函数和工具箱，可以用于地震波形数据的预处理和分析。

通过MATLAB可以对地震波形数据进行滤波、降噪、去趋势、去仪器响应等操作，减少噪声对地震数据分析的影响。

其次，地震数据的谱分析也是地震学研究中的一项重要内容。

谱分析可以帮助地震学家了解地震数据在不同频率上的能量分布情况，揭示地震波的频谱特征。

MATLAB提供了多种谱分析函数和工具箱，如快速傅里叶变换（FFT）、功率谱密度估计、波谱比等，可以用于地震数据的频谱分析。

地震学家可以通过MATLAB计算地震波的功率谱密度，绘制地震波的频谱图，进一步了解地震数据的频率特征。

此外，地震数据的时频分析也是地震学研究中的重要内容之一、时频分析可以揭示地震波的时变特征，对地震波形的瞬态信号进行分析。

MATLAB提供了时频分析函数和工具箱，如小波变换、短时傅里叶变换等，可以用于地震数据的时频分析。

地震学家可以通过MATLAB计算地震波形的时频谱，提取地震波形的瞬态特征，进一步分析地震的发展过程。

最后，MATLAB还可以用于地震数据的可视化分析。

通过MATLAB的绘图函数，可以将地震数据以图形的形式展示出来，直观地反映地震数据的变化趋势和特征。

地震学家可以通过MATLAB绘制地震波形图、频谱图、时频图等，辅助地震数据的分析和研究。

在应用方面，基于MATLAB的地震数据分析方法已经广泛应用于地震学研究和地震监测预警等领域。

例如，在地震预测方面，研究人员可以通过分析历史地震数据，利用MATLAB对地震数据进行模式识别和预测建模，从而提高地震预测的准确性和可靠性。

利用MATLAB,plot等软件来实现地震资料的处理

利用MATLAB,plot等软件来实现地震资料的处理

５ ｌ
上 述过 程得 到 的数据 实 际上是 一些 能够 反 映地 震 信号 振 幅值 的 ｓｇ ｅ ｙ格式 的 二维 数组 ， 是 由于 不 但 能 对 ｓｇ ｅ ｙ文 件进 行观 察 ， 以可 以利用 ｓｇ ．ｘ 所 ｅ ｙ ｅ ｅ软 件将 ｓｇ ｅ ｙ文 件转 化成可 以用 Ｍａｌｂ进 行处理 的ｄ ｔ ｔ ａ ａ 文件 ， 然后 考 虑 用 Ｓ ｉ ｏｍａ． ｘ ｅｓ ｒ ｔｅ ｅ软 件 来对 转 化 后 Ｆ 的ｄｔ ａ 文件进行转化 ， 转化成可以被Ｐｏ．ｘ 软件显 ｌ ｅｅ ｔ 示的ｄｔ 件 ， ａ文 再用 Ｐｏ．ｘ 软 件对 转化 好 的 ｄｔ ｌｔｅ ｅ ａ 文 件 进行 显示 ， 图 ３ 如
开发 运 用 ， 是研 究人 员如 果一 些新 的思维 和想法 ， 但 要在大 的软件上应用起来不是很方便 ， 针对这个问 题 我们 可 以通 过 利用 Ｍａｌｂ进 行 算法 的编 写 ， ｔ ａ 然后 利 用 各 种小 的软 件 来 呈现 出地 震记 录 ， 且对 于 算 并
法 的修 改和 优化 也很 方便 。下 面我 们就 以在一个 三 套 水平 地层 上运 用 自己的 算法进 行 处理 的实例来 说
一
图 ２
பைடு நூலகம்
或者 不 用 抓 图 ， 得 到 的 记录 由于 波 在传 播 的 但 过程 中的波 的 衰减 、 吸收 和 发 散导 致 浅 层 和深 层 的 能量 差 别很 大 ， 且在 同一 层 的近 道 和 远道 的差 别 而 也很 大 ， 以得 到 的离散 记录 数量 级会 差别 很大 ， 所 从 而使 得到 的视 图效果 不好 。于 是进 行 自动 增益控 制 （ ＡＧＣ 处 理 ， ） 因为 在 Ｔ ｓｅａ 上 直接 进 行 的增 益 处 ｅｓｒｌ 理 只会改变 视觉 上 的效果 而 不会 使记录 的 离散记 录 值改 变 ， 对后 续 的 Ｍａｌｂ编 程处 理 没 有 用 处 ， 是 ｔ ａ 于 考 虑 将记 录 的 ｓｇ ｅ ｙ文 件转 到 Ｐｏ ｘ中进 行 Ａ ｒＭａ ＧＣ 处理这样才会使记录的离散值彻底 的进行增益 改

matlab地震反应谱

matlab地震反应谱

matlab地震反应谱摘要：I.引言- 介绍Matlab 在地震工程领域的应用- 介绍本文的主要内容II.地震反应谱的计算方法- 加载地震波数据- 预处理地震波数据- 计算地震波数据的频谱- 计算地震反应谱III.地震反应谱的分析方法- 评估结构的地震响应特性- 为地震工程设计提供依据IV.结论- 总结Matlab 在地震反应谱计算中的应用- 展望未来可能的研究方向正文：Matlab 是一种功能强大的数学软件，可以用于解决各种工程和科学问题。

在地震工程领域，Matlab 可以用于计算地震反应谱，从而评估结构的地震响应。

本文将介绍如何使用Matlab 计算地震反应谱，以及如何根据反应谱分析结构的地震响应。

首先，我们需要加载地震波数据。

可以使用Matlab 内置的函数fiducial() 从文件中读取地震波数据，也可以使用自定义函数从其他来源获取数据。

接下来，我们需要对地震波数据进行预处理，以去除噪声和异常值。

可以使用Matlab 内置的函数removeoutliers() 对数据进行去噪处理。

然后，我们可以使用Matlab 内置的函数freqz() 计算地震波数据的频谱。

通过频谱，我们可以了解地震波的频率分量以及振幅。

最后，我们可以使用Matlab 内置的函数responsespectrum() 计算地震反应谱。

该函数可以根据结构的阻尼比和自振周期计算反应谱。

通过反应谱，我们可以了解结构在不同地震波作用下的加速度响应。

综上所述，Matlab 可以用于计算地震反应谱，从而评估结构的地震响应。

基于matlab的地震反应谱计算方法的比较【matlab源码】

基于matlab的地震反应谱计算方法的比较【matlab源码】

毕业论文（设计）题目学院学院专业学生姓名学号年级级指导教师教务处制表基于MATLAB的地震反应谱计算方法的比较一、程序说明本团队长期从事matlab编程与仿真工作，擅长各类毕业设计、数据处理、图表绘制、理论分析等，程序代做、数据分析具体信息联系二、写作思路与程序示例地震反应谱是进行结构抗震分析与设计的重要工具，反应谱的计算在反应谱法和时域逐步积分方法中有重要地位，引起了学者的重视和广泛研究。

而对计算方法优劣的评定常取决于其计算的耗时、稳定性和精度等因素。

目前，计算反应谱的方法有很多，以往常规方法主要有中心差分法、Newmark-法、线性加速度法及Wilson-法等，这些方法虽然存在一些弊端，但是其计算精度能够满足一般实际工程计算的需要，仍被广泛应用于工程实践。

因此有必要对这些方法进行深入的比较和探讨。

目前，我国现行建筑抗震设计规范中对阻尼比规定：混凝土结构通常取0.05，钢结构通常取0.02，阻尼比均比较小，因此利用拟反应谱是可靠的。

然而，随着高层和超高层的出现，建筑新形式、新材料的使用、隔震、减震耗能机构的研究以及抗震减灾新要求的不断提出，传统理论也越来越不能适应新要求，使用拟反应谱可能带来较大误差，这应当引起我们的注意。

同时，通过反应谱理论分析得到：当周期超过3s以后，结构地震反应已不是由地面加速度控制，而可能是由速度甚至是位移控制。

然而，现代的超大跨、超高层和巨型结构的自振周期大都达到了10s以上，如果仍然采用加速度谱来对其进行分析将不合适。

我们应当采用能更好反应三个周期段的反应谱来对长周期结构进行分析，而三联反应谱就具有这一特点，它是将位移、速度以及加速度联合反应谱同时表示在一张图上的四坐标对数图，能够使我们更直接地掌握上述三种物理量与结构自振周期之间的控制关系。

因此，采用MATLAB的GUI编程三联反应谱图形界面，并将其应用于工程实践，将是一项很有意义的工作。

论文基于以上内容，主要进行了以下具体工作：1.基于数值算法的相关研究及应用现状，本文以MATLAB为平台，建立数值算法在不同影响因素下的三维图形，并结合理论对比分析。

基于matlab的地震活动性分析

基于matlab的地震活动性分析Matlab在地震活动性图像分析中的应用1），李红光2）1）河北省地震局2）中国地震应急搜救中心摘要：地震活动性分析是地震预测、地震工程的一个重要依据，地震活动性分析又多是通过图像来表现。

Matlab是一种简单易学、强大的计算功能和编程可视化的计算机语言。

本文用Matlab语言编程，实现了地震统计区内地震的快速选取，并根据这些选中的地震进行地震活动性分析。

关键词：Matlab语言；地震活动性引言地震活动性研究就是通过分析一定震级区间内的地震时间、空间的分布特征，探讨其物理含义，进而对地震发生的规律进行科学总结。

通过地震活动性研究，可对地壳介质非均匀性和运动形态有宏观的了解和总体把握，因此可服务于地震预测和地球动力过程等研究[1]。

在地震安全性评价中，通过地震活动性分析，为工程场地一定时间内的地震活动性趋势和地震环境做出评价，为合理划分潜源区和确定其地震活动性参数提供依据[2]。

地震活动图像的分析方法很多，有简单的图像描述法，如地震震中分布、蠕变曲线、M-T图等；也有采用统计参数表征地震活动时空图像特征的方法，如b值、地震活动度S、地震能流密度、地震强度因子MF分布等。

Matlab具有强大的计算能力、计算结果可视化和编程效率高的优势，它是地震活动性分析的一个有力工具。

Matlab是1984年由美国MathWorks公司推出的荣誉产品。

早在20世纪80年代中期，Matlab就在我国出现，大规模流行时再90年代中期以后。

现在Matlab已被广泛应用在科学研究、工程计算等方面。

M atlab采用全新的数据类型和面向对象编程技术，采用了新控制流和函数结构，特别是包含很多常用的子函数，非计算机专业人员非常容易用Matlab来实现很复杂的计算程序。

并且Maltab提供了图像处理功能，可以很方便的生成图形。

在地震活动性数字图像分析中，用Matlab可以很简单、方便的实现研究人员的思想。

基于MATLAB的地震反应谱与傅里叶谱计算分析

图5 Fig.5 傅里叶相位谱
参考文献(References)：
[1] 郭晓云. 汶川地震反应谱研究分析[D]. 哈尔滨： 中国地震局工程力 学 研 究 所 ， 2011.(DUO xiaoyun. Study on response spectrum of WenChuan Earthquake.[D]. Harbin ： Institute of Engineering Mechanics，China Earthquake Administration ，2011.(in Chinese)) [2] 骆剑锋. 框架结构静力与动力弹塑性抗震分析对比研究[D]. 上海： 同济大学，2007. [3] 张晓志，谢礼立，于海英. 地震动反应谱的数值计算精度和相关问 题[J]. 地震工程与工程振动，2004，24(6)：15-20. [4] 赵凤新，胡聿贤. 地震动非平稳性与幅值谱和相位差谱的关系[J]. 地震工程与工程振动，2003，23(1)：1-5. [5] 金 星， 廖振鹏. 地震动强度包线函数与相位差谱频数分布函数的
[2]
的方法，因此对Δt 时段内的系统反应或加速度的变 化规律未做出假定。 事实上， 采用 Fourier 变换方法 求解动力平衡方程时，对输入荷载采用了一个全局 性基本假定，即将非周期的加速度荷载假定为一个 具有有限带宽的周期荷载。在此假定下，加速度荷 载被表达为有限个连续简谐波的线性叠加。因此， 频域方法也是近似方法，不能作为度量其他方法精 度的准则。 傅立叶谱不仅包含幅值信息，而且包含相位信 息，这是地震反应谱所不具备的。所以，当傅里叶 谱与地震反应谱进行比较时，只能采用幅值谱。
3 强震记录资料
北岭大地震， 1994 年 1 月 17 日凌晨 4 时 31 分， 洛杉矶地区发生里氏 6.6 级地震，震中位于圣费兰 多峡谷， 北纬 34°12′54″， 西经 118°32′16.8″， 在洛杉矶西北方向 20 英里处， 属浅源地震， 震源深 度 为 18 英 里 。 当 地 震 级 ML=6.6 ， 面 波 震 级 为 MS=6.7（NEIC），地震矩震级为 MW=6.7（CIT）， 持续时间约 30 秒。在持续 30 秒的震撼中，大约有 11000 多间房屋倒塌， 震中 30 公里范围内的高速公 路、高层建筑或毁坏或倒塌，煤气、自来水管爆裂， 电讯中断，火灾四起，直接和间接死亡 58 人，受伤 600 多人，财产损失 300 多亿美元。 本文以北岭大地震中，某县医院台站强震动仪 所记录到的强震记录为例，通过 MATLAB 编程实 现傅里叶谱和地震反应谱的计算， 并进行比较分析。 该台站位于北纬 34°19′33.6″，西经 118°26′ 38.4″，所记录的振动持续时间为 39.98s，峰值加 速度为-1007.942cm/s²，记录时间点在 4.328s；峰值 速度为-124.896cm/s，记录时间点在 3.700s；最大位 移为-31.078cm，记录时间点在 4.200s；初始速度为 -3.694cm/s；初始位移为 298cm。

基于MATLAB的地震数据的分析资料

基于MATLAB的地震数据的分析孙玉柱冯光房桂梅摘要：地震波原始数据中存在的干扰信号，会影响震相分析的准确性。

为了滤除干扰信号，对地震波原始信号进行了频谱分析，给出了一种基于MATLAB的FIR数字滤波器的优化设计方案，将其用于地震波数据的分析中，并进行了仿真分析。

仿真结果表明，FIR数字滤波器对地震波原始信号进行滤波处理后，提高了震相分析的准确性，得到了理想的效果，达到了预期的目的。

关键词：MATLAB；FIR数字滤波器；优化；滤波the Analysis of Earthquake Data Based on MATLAB SUN Yuzhu,FENG Guang,FANG Guimei Abstract: The interference that existed in the earthquake data will affect the accuracy of the seismic phase analysis. In order to filter the disturbance signal, this paper carries out spectrum analysis of the earthquake data, proposes an optimum design method for FIR digital filter based on MATLAB and applies it to the analysis of earthquake data. After the filter of the noise jamming, the true information of the earthquake wave is clearly reflected. The simulation results manifest that it canimprove the accuracy of seismic phase analysis and arrive at the purpose desired.Key words: MATLAB；FIR digital filter；optimization；filter1 引言地震带给人类的损失是巨大的，汶川大地震依旧在我们的记忆深处清晰存在。

matlab地震反应谱

在MATLAB中，地震反应谱的计算可以帮助我们更好地理解地震波动的特征和规律。

具体来说，地震反应谱可以表示地震动强度特性和频谱特性的关系，而傅里叶谱则可以用来检出时间过程中所含的频率分量并进行时域到频域的变换。

在MATLAB编程计算分析中，地震反应谱的特征参数包括平台值和特征周期。

平台值表示地震动的强度特性，特征周期则反映了地震动的频谱特性。

通过MATLAB计算分析，我们可以得到标准加速度反应谱峰值、相对加速度谱峰值的数值，这些数据可以用来确定地震动的相关模型及其关键参数。

此外，傅里叶谱表示地震波的重要意义还体现在两个方面：一是检出时间过程中所含的频率分量，二是进行时域到频域的变换。

通过傅里叶振幅谱的最大振幅值和其所对应的频率，我们可以进一步研究地震波动的特征和规律。

总的来说，MATLAB中的地震反应谱和傅里叶谱分析工具对于研究地震波的特征和规律具有重要意义，可以广泛应用于地震工程、结构抗震分析和设计等领域。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于 MATLAB 的地震反应谱与傅里叶谱计算分析
XXX
摘要：地震反应谱特征参数包括平台值和特征周期，其分别表征了地震动的强度特性和频谱特性；傅里叶谱表示 地震波的重要意义有两个，一是检出时间过程中所含的频率分量，另一是进行时域到频域的变换。通过 MATLAB 编程计算分析，标准加速度反应谱峰值 3517.7gal，相对加速度谱峰值为 244.61m/s，析主要影响原因为震中距小， 震源深度浅，震级大；傅里叶振幅谱的最大振幅值为 665.6gal*s，其所对应的频率为 2.5Hz。 关键词：地震反应谱；傅里叶谱；MATLAB；北岭大地震
3 强震记录资料
北岭大地震， 1994 年 1 月 17 日凌晨 4 时 31 分， 洛杉矶地区发生里氏 6.6 级地震，震中位于圣费兰 多峡谷， 北纬 34°12′54″， 西经 118°32′16.8″， 在洛杉矶西北方向 20 英里处， 属浅源地震， 震源深 度 为 18 英 里 。 当 地 震 级 ML=6.6 ， 面 波 震 级 为 MS=6.7（NEIC），地震矩震级为 MW=6.7（CIT）， 持续时间约 30 秒。在持续 30 秒的震撼中，大约有 11000 多间房屋倒塌， 震中 30 公里范围内的高速公 路、高层建筑或毁坏或倒塌，煤气、自来水管爆裂， 电讯中断，火灾四起，直接和间接死亡 58 人，受伤 600 多人，财产损失 300 多亿美元。 本文以北岭大地震中，某县医院台站强震动仪 所记录到的强震记录为例，通过 MATLAB 编程实 现傅里叶谱和地震反应谱的计算， 并进行比较分析。 该台站位于北纬 34°19′33.6″，西经 118°26′ 38.4″，所记录的振动持续时间为 39.98s，峰值加 速度为-1007.942cm/s²，记录时间点在 4.328s；峰值 速度为-124.896cm/s，记录时间点在 3.700s；最大位 移为-31.078cm，记录时间点在 4.200s；初始速度为 -3.694cm/s；初始位移为 298cm。
[2]
的方法，因此对Δt 时段内的系统反应或加速度的变 化规律未做出假定。 事实上， 采用 Fourier 变换方法 求解动力平衡方程时，对输入荷载采用了一个全局 性基本假定，即将非周期的加速度荷载假定为一个 具有有限带宽的周期荷载。在此假定下，加速度荷 载被表达为有限个连续简谐波的线性叠加。因此， 频域方法也是近似方法，不能作为度量其他方法精 度的准则。 傅立叶谱不仅包含幅值信息，而且包含相位信 息，这是地震反应谱所不具备的。所以，当傅里叶 谱与地震反应谱进行比较时，只能采用幅值谱。
Hz，因此可推知 该条地震波对自振周期在 0.4s 左右的建筑物，也就 是说对 5~6 层的钢筋混凝土结构，会产生严重的影
通过对该台站地震记录的计算分析，发现该县
响。
本文通过 MATLAB 编程计算得到了反应谱与 傅里叶谱，由分析总结得出： （1）标准加速度反应谱峰值加速度峰值为 3517.7gal，特征周期短，分析主要影响原因为震中 距小，震源深度浅，震级大。 （2）相对加速度谱的速度峰值为 244.61m/s， 分析其原因也是由于震级大，震源深度浅、震中距 小对其造成的影响。 （3）傅 里 叶 振 幅 谱 的 最 大 振 幅 值 为 665.6gal*s，其所对应的频率为 2.5Hz，可推知该条 地震波对自振周期在 0.4s 左右的建筑物会产生严重 的影响。
图 1 加速度时程曲线 Fig.1 Acceleration time-histories curve
震源深度为 18 公里。 图 3 为速度反应谱，通过最小二乘大拟合利用 MATLAB 编程得到，由计算结果可知，该条强震记 录所对应的相对速度峰值为 244.61m/s， 分析其原因 是由于震级大，震源深度浅、震中距小对其造成的 影响。 4.1 傅里叶谱计算分析 地震波的分量重，哪些分量的振幅大，哪些振 幅小，对地震波的性质有重要的影响。 通过 MATLAB 编程计算分别得到傅里叶振幅 谱、傅里叶相位谱、功率谱和自相关系数谱。
4
基于 MATLAB 的傅里叶谱与反应
谱的实现
本文通过 MATLAB2010b 计算软件编程，实现 反应谱与傅里叶谱的计算，进而得出相关参数并进
行对比分析。 4.1 反应谱计算分析
医院台站加速度峰值为 1007.942gal（1.029g），记 录时间在 4.328s，此峰值也是县医院台站本次地震 记录到的最大加速度值。 图 1 为该县医院台站所记录到了加速度时程曲 线，通过 MATLAB 编程计算得出加速度峰值在 1G 左右，且在 4.328 秒位置，该时程曲线记比较完整， 振动衰减清晰可见。 图 2 为台站强震记录通过最小二乘法拟合计算 得到的标准加速度反应谱，由图中可得，该条地震 记录的标准加速度峰值将近 3500gal （3.57g），且 特征周期很短，曲线为高瘦型，分析其影响原因有 二：一是距离震中很近，该县医院台站距离震中仅 有 15.1 公里； 二是震源浅， 属浅源地震， 且震级大，
the peak flat values denote the strength characteristic of ground motions and the characteristic period reflects the stong motion’s characteristic. Fourier spectrum indicates that there are two important meanings of seismic wave, one is the frequency component in the process of detection, the other is to carry out the time domain to the frequency domain. Through MATLAB programming, the peak value of the standard acceleration response spectrum is 3517.7gal, the relative acceleration spectrum is 244.61m/s, and the main influence factors are the small epicenter distance, the shallow source depth, and the large magnitude. The maximum amplitude value of the Fourier amplitude spectrum is 665.6gal*s, and the corresponding frequency is 2.5Hz.
点是其与时程的对应性， 有时程可以获得傅里叶谱， 由傅里叶谱也可以反演出时程 。对于建筑结构而 言，通过地震危险性分析和已有资料的统计回归， 可以确定地震动的相关模型及其关键参数，从而进 行结构的抗震设计，因此实现地震反应谱、傅里叶 谱等计算并对其有较清楚的认识尤为重要。
2 反应谱与傅里叶谱的关系
对于同一条地震动，地震反应谱与傅里叶谱描
Keywords：Seismic response spectrum；Fourier spectrum；MATLAB；Northridge Earthquake
的、平稳的还是非平稳的。傅里叶谱的另外一个优
1 引 言
地震反应谱是进行结构抗震分析与设计的重 要工具，反应谱的计算在反应谱法和时域逐步积分 方法中有重要地位，引起了国内外学者的重视和广 泛研究。而对计算方法优劣的评定常取决于其计算 的耗时、稳定性和精度等因素。反应谱利用一个单 自由度体系的反应间接反映地震的频谱特性；傅里 叶谱则是一种信号分析工具，对于任意时程信号都 可以获得其傅里叶谱，而不论其是随机的还是确定
Analysis of seismic response spectrum and Fourier spectrum based on MATLAB
XXXXXXX
Abstract：The characteristic parameters of seismic response spectrum include peak flat value and characteristic period,
图5 Fig.5 傅里叶相位谱
参考文献(References)：
[1] 郭晓云. 汶川地震反应谱研究分析[D]. 哈尔滨： 中国地震局工程力 学 研 究 所 ， 2011.(DUO xiaoyun. Study on response spectrum of WenChuan Earthquake.[D]. Harbin ： Institute of Engineering Mechanics，China Earthquake Administration ，2011.(in Chinese)) [2] 骆剑锋. 框架结构静力与动力弹塑性抗震分析对比研究[D]. 上海： 同济大学，2007. [3] 张晓志，谢礼立，于海英. 地震动反应谱的数值计算精度和相关问 题[J]. 地震工程与工程振动，2004，24(6)：15-20. [4] 赵凤新，胡聿贤. 地震动非平稳性与幅值谱和相位差谱的关系[J]. 地震工程与工程振动，2003，23(1)：1-5. [5] 金 星， 廖振鹏. 地震动强度包线函数与相位差谱频数分布函数的
图2 Fig.2
标准加速度反应谱
Standard acceleration response spectrum
图4
傅里叶振幅谱
Fig.4 Fourier amplitude spectrum
图 4 为傅里叶振幅谱，通过计算结果可知，该 条强震记录所对应的傅里叶最大振幅值为
图 3 相对速度反应谱 Fig.3 Relative velocity response spectrum
Wilson - 法等，从精度、稳定性和计算效率三个
方面进行了一个简单的定性对比。 2.2 傅里叶谱 傅立叶谱是利用傅立叶分析对地震动时程进行 分解的结果，其基本原理是任意非周期性的复杂函 数均可以表示成简单周期函数的组合，由此可以获 得傅氏幅值谱和傅氏相位谱， 两者合称傅立叶谱[4]。 日本学者大崎顺彦最早强调了加速度时程的相 位谱对时程的影响。同以往研究者不同的是，大崎 顺彦不是直接研究相位谱，而是研究相位差谱对地 震动包线的作用。他发现地震加速度记录的相位差 的频数分布曲线与加速度时程的强度包线形状在一 定程度上具有相似性。 Nigam 研究了均匀调制的高斯白噪声随机过程 的确定性包线同相位谱导数之间的关系； Sawada 推导了相位差谱在窄带内随机分布时窄带波包的到 时同相位差谱的关系；金星等[5]根据群速度与相速 度的关系， 阐明了相位差谱这一综合物理量的含义， 讨论了它与等效群速度的联系。 2.3 两者之间的关系 反应谱计算方法是从时域的角度出发的，而 Fourier 变换方法则是一种不同于时域逐步积分的 方法，而是从频域计算反应谱的方法。由于 Fourier 变换方法不是一个在时域逐步递推求解系统反应量