傅立叶讲义谱与反应谱
傅里叶光学讲义
傅⾥叶光学讲义傅⾥叶光学实验傅⾥叶光学原理的发明最早可以追溯到1893年阿贝（Abbe ）为了提⾼显微镜的分辨本领所做的努⼒。
他提出⼀种新的相⼲成象的原理，以波动光学衍射和⼲涉的原理来解释显微镜的成像的过程，解决了提⾼成像质量的理论问题。
1906年波特（Porter ）⽤实验验证了阿贝的理论。
1948年全息术提出，1955年光学传递函数作为像质评价兴起，1960年由于激光器的出现使相⼲光学的实验得到重新装备，因此从上世纪四⼗年代起古⽼的光学进⼊了“现代光学”的阶段，⽽现代光学的蓬勃发展阶段是从上世纪六⼗年代起开始。
由于阿贝理论的启发，⼈们开始考虑到光学成像系统与电⼦通讯系统都是⽤来收集、传递或者处理信息的，因此上世纪三⼗年代后期起电⼦信息论的结果被⼤量应⽤于光学系统分析中。
两者⼀个为时间信号，⼀个是空间信号，但都具有线性性和不变性，所以数学上都可以⽤傅⽴叶变换的⽅法。
将光学衍射现象和傅⽴叶变换频谱分析对应起来，进⽽应⽤于光学成像系统的分析中，不仅是以新的概念来理解熟知的物理光学现象，⽽且使近代光学技术得到了许多重⼤的发展，例如泽尼克相衬显微镜，光学匹配滤波器等等，因此形成了现代光学中⼀门技术性很强的分⽀学科—傅⾥叶光学。
实验原理：我们知道⼀个复变函数f(x,y)的傅⽴叶变换为：+-=?=dxdy vy ux 2i y x f y x f v u F )](exp[),()},({),(π ( 1 )F (u,v)叫作f(x,y)的傅⽴叶变换函数或频谱函数。
它⼀般也为复变函数，f(x,y)叫做原函数，也可以通过求 F(u,v)逆傅⽴叶变换得到原函数f(x,y)：+=?=-dudv vy ux 2i v u F v u F y x f 1)](exp[),()},({),(π（2）在光学系统中处理的是平⾯图形，当光波照明图形时从图形反射或透射出来的光波可⽤空间两维复变函数（简称空间函数）来表⽰。
在这些情况下⼀般都可以进⾏傅⾥叶变换或⼴义的傅⾥叶变换。
反应谱和傅里叶谱+地震波选取
【拓展知识1-2】功率谱，反应谱和傅里叶谱，地震波选取，地震持续时间确定功率谱功率谱是功率谱密度函数的简称。
对于一般情况的随机振动，其时间历程具有明显的非周期性，具有连续的多种频率成分，每种频率有对应的功率或能量，用图像来表示这种关系，称为功率在频率域内的函数，简称功率谱密度。
加速度功率谱是对地震动加速度时程进行快速傅里叶变换（FFT）得到的[1]。
对于非平稳随机过程，功率谱密度的单位是G的平方/频率。
G指的是随机过程。
对于加速度功率谱，加速度的单位是m/s2，则功率谱密度的单位是（m/s2）2/Hz，Hz的单位是1/s，故加速度功率谱密度的单位为m2/s3。
加速度功率谱密度函数曲线下方的面积代表随机加速度的总方差，即加速度功率谱可以理解为“随机加速度方差的密度分度”。
参考文献[1] 庄表中. 随机振动入门.科学出版社，1981.反应谱和傅里叶谱反应谱（earthquake response spectrum），是单自由度弹性系统对于某个实际地震加速度的最大反应（可以是加速度、速度和位移）和体系的自振特征（自振周期或频率和阻尼比）之间的函数关系。
反应谱是地震工程中分析结构和设备在地震中的性能的非常有用的工具，因为许多主要表现为简单的振荡器(也称为单自由度系统)。
因此，如果能找出结构的固有频率，那么建筑的峰值响应可以通过从地面响应谱中读取相应频率的值来估计。
在地震区域的大多数建筑规范中，这个值构成了计算结构必须抵抗的力的基础(地震分析)。
如前所述，地面响应谱是在地球自由表面所做的响应图。
如果建筑物的响应与地面运动(共振)的组成部分“协调”，可能会发生重大的地震破坏，这些成分可以从响应谱中识别出来。
傅里叶谱，全称为傅里叶振幅谱。
地震波是在时间上连续的随机过程，地震动记录仪是按照一定的采样频率得到该连续曲线上离散的点，想要还原这个曲线，可以通过解N 元1次方程组，更简洁有效的方式是采用有限傅里叶级数来近似原始的时间历程。
地震工程学地震动特性
8
3.3 地震动特性 3.3.2 地震动频谱特性
• 地基状况相同，震级与震中距不同
• 震级越大，地震动记 录中长周期（低频） 分量越显著；
• 震中距越远，地震动 记录中的长周期（低 频）分量越显著；
9
3.3 地震动特性 3.3.2 地震动频谱特性
F
m
xg
x max
✓ 意义：数学上 当 x xm时ax
Fourier变换： F( )
xg(t)e i tdt
F( )
T
0 xg(
)e i
d
T
0 xg( )cos d
T
i 0 xg( )sin d
F( )
T
2
0 xg( )cos d
T
2
0 xg( )sin d
无阻尼的相对速度反应谱： Sv( )
t
0 xg ( ) cos[ (t
T
2
0 xg( )cos d
傅里 叶变 换对 复数 表达
x(t)
C ei(2 kt/T ) k
k
k
1 T
df ,Tk
f
x(t)
TCk ei2 (k/T )t
1 T
F f ei2 ftdf
1
2
F w ei tdw
连续傅里 叶变换对
F f TCk
x t e i(2 ft)dt
x t e i tdt
傅里叶谱
Ff
TCk
T 2
Ak2
Bk2
✓ 初始条件： x(0) 0 x(0) 0
✓ 运动方程的解： 令 D 1 2
地震工程学-傅里叶变换-反应谱计算
一、得到地震波数据 (1)二、地震波分析 (2)时程曲线绘制 (2)傅里叶谱的绘制 (3)反应谱结果分析 (8)[参考文献] (18)附件说明 (18)一、得到地震波数据访问网站，下载相应的附录说明文件，在这些文件中可以找到相关地震的信息，便于有效的规划检索词，可以看到题目需要的地震波RSN 编号为6,；图查询相应的编号根据所查到的RSN编号直接检索，如图所示输入RSN号即可检索图检索条件得到检索结果下载即可得到所需的地震波，如图所示图检索结果具体地震波数据见附件1.二、地震波分析时程曲线绘制根据得到的地震波数据，进行MATLAB程序编制，绘制的到相应的竖向地震、180度地震程曲线，270度地震的加速度时程曲线，分别为图，图竖向地震加速度时程曲线图 180度地震加速度时程曲线傅里叶谱的绘制根据离散傅里叶的变换准则可以得到A k （k=0,1，…,N/2）、B k (k=1,….N/2-1)，进而计算得到相应的参数； 式 ; 式式式式中N=1024，T=;在具体计算时采用了MATLAB 中的FFT 函数，并且对参考文献中[1]中的例题波进行了试算得到了理想的结果；故可证明程序中所使用的算法是没有问题的；得到的结果见下。
图竖向地震波的傅里叶幅值曲线图 180度地震波的傅里叶幅值曲线图 270度地震波的傅里叶幅值曲线图竖向地震波的傅里叶相位曲线图 180度地震波的傅里叶相位曲线图 270度地震波的傅里叶相位曲线图竖向地震波的功率谱曲线图 180度地震波的功率谱曲线图 270度地震波的功率谱曲线反应谱结果分析根据地震波是可以求得相应的不同周期的单自由度体系的各种反应谱的，这包括：位移反应谱，速度反应谱，准速度反应谱，加速度反应谱，绝对加速度反应谱，具体的原理是根据杜哈梅积分对单自由度进行积分获得，相应的公式如下[2]式式其中，,式中ε为阻尼比，在本题中分别取为0和。
此处对MATLAB计算程序的相应算法进行说明，在实际计算反应谱时使用的是复合的辛普森积分公式[3]式因此在计算中是根据△t=进行数值积分运算的，这样选择的理由有二，1.辛普森公式具有高阶的精度，2由于时间间隔较小是可以捕捉到最大值的，不至于引起太大误差又可以减小计算量，大大较少程序运行时间。
地震动高频衰减参数(kappa)模型及汶川M_S8.0地震地震动模拟
第1期中国地震局地球物理研究所2016届博士论文摘要（I )47块体和小江断裂带这两个低速异常区之间，即四川盆地西南部至攀枝花地区的地壳速度表现为相对高速特征；滇西地区的地壳速度整体偏低，但不均匀连续。
地壳厚度自北向南逐渐变 薄，川滇块体的地壳厚度由北部的60 km减薄至南部的40 km，且川滇块体的地壳厚度比两侧块体厚，北部的川西北块体在东边界厚度变化剧烈，南部的滇中块体在东西边界厚度变化比较平缓。
推测四川盆地西南部至攀枝花地区的高速特征与峨眉山玄武岩的分布有关，该高速物质对川西北次级块体的低速体的运移具有阻挡作用，导致地壳增厚和地形隆升；小江断裂带北部的低速带范围减小可能是由于峨眉山火成岩区的高速物质侵入导致。
青藏高原东南缘地壳整体表现为大范围的低速弱强度区，但低速异常区并非预想的那么均匀，呈现出侧向不均匀性。
如果青藏高原物质是以地壳流模式运移，那么地壳流在青藏高原东南缘会受到高速异常体的影响，从而不能连贯地发生。
关键词青藏高原东南缘；接收函数；各向异性；S波速度结构；联合反演(作者电子信箱，蔡妍：caiyan@)地震动高频衰减参数（kappa)模型及汶川Ms8. 0地震地震动模拟傅幕(中国地震局地球物理研究所，北京100081)中图分类号：P315.8; 文献标识码：A; doi：10. 3969/j. issn. 0235-4975. 2017. 01. Oil在地震动记录稀疏或缺乏的地区，模拟地震动的随机性方法在建立地震动预测方程和估计设定地震的地震动场中得到了广泛且成功的应用，证明了该方法具有很高的可信度和发展前景。
在地震动模拟中，高频衰减参数O c)是一个重要的参数，它控制了地震动的傅里叶振幅谱和反应谱的高频部分（>1Hz)，对地震动时程的峰值有显著的影响。
但是，/C的物理起源、计算方法和经验模型的建立，至今存在较多的争议。
在对2008年5月12日汶川Ms8.0地震和2013年4月20日芦山Ms7.0地震的强震动记录进行预处理之后，我们选取龙门山地区69个强震动观测台站记录到的震中距小于300 km、面波震级在3.3〜8.0之间的500多个主余震的3分量强震动记录共1700多组，对该地区的高频衰减特性及高频衰减参数(/c)模型进行了研究，取得以下进展：(1)对强地面运动的模拟方法以及国内外的研究现状进行了系统的总结，并详细介绍了 傅里叶振幅谱高频衰减特性的研究现状和面临的问题。
傅立叶变换及幅值谱
对于零均值化的信号，第一项为零。而在满足采样定理的情 况下最后一项为零：
mn mn jm jm x(n) X (m) e exp( j 2 ) X (m) e exp( j 2 ) N N m 1
N 1 2
x ( n)
N 1 2
N 1 2 m 1
N 1
（ n N-1） N 1 n X (N-1) x(n) exp( j 2 )= x( n) exp( j 2 n j 2 ) N N n0 n0 N 1 n x(n) exp( j 2 ) N n0
N 1
K=1与N-1的X(k)是共轭对，与k=2与N-2也是,….. 。再看
x(nt ) X (k f ) exp( j 2 k fnt )f
k 0
N 1
1 1 f T N t
N 1
X (k )
X (k f ) ts
kn 1 x(n) X (k ) exp( j 2 ) N N k 0
在Matlab中，假设有时域信号 x(n) （n=0,1,…N-1），对应的 时间是0,dt,2dt,….(N-1)dt。对x(n)进行DFT，即 fft(x) 得到
mn mn X (m) exp( j 2 +m + exp( j -2 -m N N
mn = 2 X (m) cos 2 + m N m 1
可见与第m个频率分量对应的幅值是
2 X (m)
X (m) x(n) exp( j 2
kn 1 x(n) X (k ) exp( j 2 ) N N k 0
K=1与N-1的
kn exp( j 2 ) N n exp( j 2 ) N n( N 1) n) kn exp( j 2 ) exp( j 2 n j 2 ) exp( j 2 ) N N N
傅里叶谱的简单概述
傅⾥叶谱的简单概述
【转载问题】地震波经过傅⾥叶变换，得到的频谱，主要有哪些应⽤？
先来介绍⼀下Fourier谱：基本思想是把⼀个复杂的地震动过程按傅⽒级数展开⼀系列具有不同频率的周期函数，包括傅⽒幅值谱和傅⽒相位谱（详地震⼯程相关书籍）。
幅值谱与相位谱从两个不同⾓度描述了地震动的频谱特性，与反应谱、功率谱相⽐，傅⽒谱对地震动的描述描述更加全⾯，包括了各频率分量的相位分布信息（更接近于实际，规范利⽤反应谱来进⾏设计存在⼀定的弊端，但便于操作），即频率间的相互影响也考虑了。
幅值谱与相对速度反应谱存在⼀定联系，幅值谱是相对速度反应在地震动终⽌时的值，⽆阻尼相对速度反应谱是整个地震动过程中相对速度反应的最⼤值（反应谱是绝对值概念），即⼀般来讲，傅⽒幅值谱⼩于等于相对速度反应谱值，两者⼤体相当。
楼主图中给出的是地震动傅⽒幅值谱，该谱透露的信息类似于反应谱，但较反应谱更加真实，因为其考虑了频率间的相互影响，由幅值谱可以⼤概判断处地震动的卓越周期（频率），反应幅值谱可以判断结构反应的峰值周期（频率）。
补充：频谱分析⼀般设计到功率谱、反应谱及傅⽒谱，后者与前两者的最⼤区别在于傅⽒谱考虑了相位因素（让我想到前不久的⼀个问题，关于楼板应⼒设计时采⽤振型反应谱法的弊端，即反应谱法未考虑反应的相位因素，⽽楼板主拉应⼒是有相位因素的），三者之间存在着⼀定的联系，功率谱是⼀簇样本样本函数（可以是地震动、外荷载或结构反应）的傅⽒幅值谱的平⽅的平均值，反应谱与功率谱的联系稍复杂，⼆者通过反应谱的超越概率建⽴关系。
傅里叶谱和反应谱之间的关系及其对经验地震动预测方程(GMPE)调整的影响研究(Ⅰ)
傅里叶谱和反应谱之间的关系及其对经验地震动预测方程(GMPE)调整的影响研究(Ⅰ)Sanjay Singh Bora;Frank Scherbaum;Nicolas Kuehn;Peter Stafford【摘要】经验反应谱地震动预测方程(GMPEs)的函数形式通常由地震动傅里叶谱模拟的概念导出, 随后这些GMPEs由经验观测数据校准, 所以, 对特定地震场景的地震动预测不会构成重大问题. 然而, 当调整反应谱GMPEs来计算未被原始经验数据集覆盖的条件时, 傅里叶谱模拟导出反应谱的假设可能会产生意想不到的结果. 因此, 几个问题出现了. 例如, 地震动傅里叶谱和反应谱的区别和相似之处是什么? 如果它们是不同的, 什么机制可以解释这样的差异? 对傅里叶振幅谱(FAS)所做的调整在反应谱中怎么表示? 本文利用随机振动理论(RVT)探讨了地震动傅里叶谱和反应谱之间的关系. 借助简单的Brune震源模型[1-2], 在固定震级和距离的情况下, 用RVT生成了加速度谱. RVT分析表明, 反应谱低频谱值标定可以视为与相应傅里叶谱标度值等同. 然而, 反应谱高频谱值受到了带宽很大的傅里叶谱的控制. 实际上, 地震动峰值加速度不能反映高频地震动特征, 它受到了整个地震动傅里叶谱的控制. 此外, 本文说明了对FAS做的调整如何相似或区别于对反应谱坐标所做的相同的调整. 为此, 我们研究了调整应力参数(Δσ) (震源项)和调整反应场地响应的属性(VS-κ0)这两种情况.【期刊名称】《国际地震动态》【年(卷),期】2017(000)003【总页数】8页(P25-32)【作者】Sanjay Singh Bora;Frank Scherbaum;Nicolas Kuehn;Peter Stafford【作者单位】【正文语种】中文【中图分类】P315.9在当前概率地震危险性分析(PSHA)中，最常用的地震动强度指标/参数(GMIM)是用特定阻尼水平的(如，临界阻尼的5%)单自由度(SDOF)结构反应谱坐标进行的。
傅里叶红外光谱样品调制及图谱解析技巧
制备高聚物薄膜常用溶剂适合的溶剂高聚物苯 聚乙丁烯、聚丁二烯、聚苯乙烯等 甲》聚醋酸乙烯酯、乙基纤维董 二甲基甲酰胺 聚丙烯m氯仿或丙《 聚甲基丙烯酸甲酩 甲酸尼龙6 ;二氯乙烷聚碳酸ffii 丙n纤维素 四氯乙烷 涤纶 四氢咲响\聚氯乙烯二甲亚» 聚酰亚胺、聚甲醛（热》 甲苯、四氢蔡聚乙烯（热〉、聚丙烯（热》 1水_«乙烯静（热）、甲基纤维素常用的反射配件•固体样 品 -粉末样 晶 ・定性及 定*分析漫反射附件、7、-•a可用于，单层膜分析LB般镀膜■薄层分子取向研究黑色样品单次反射ATR3、液体样品的制备对于沸点较高且粘度较大的液体样品, 或一滴样品直接涂在KBr 窗片上进行测试对于沸点较低的样品及粘度小.流动性较大的 高沸点液体样品放在液体池中测试液体池是由两片KBr 窗片和能产生一定厚度的 垫片所组成切记不得有水-—•固体 •液体•不规则的 样品•非破坏性取 2ing气体池以及气体池架将气体池放在气体池架上即可，气体池的 两边由KBr 窗片或其它类型的盐片密封，要特 别注意防止盐片受潮。
.o®二、红外光谱解析技巧1、分子结构对基团吸收谱带位置的影响衽双原子分子中，其特征吸收谱带的位2由键力常数和原 子质fi 决定-在复杂的有机化合物分子中，某一基团的特征 吸收频率同时还要受到分子结构和外界条件的彭响。
同一种 基团由于其周围的化学环境不同，使其吸收频率会有所位移， 而不在同一位置出峰。
即基团的吸收不是固定在某一个频率 上，而在一定范围内波动. 如：C-H 的伸缩振动频率受到与这个碳原子邻接方式的影响C-C-H ： C=C-H ： C=C-H ：外部条件对吸收的影响有！物态效应、晶体状态和溶剂效应• 主要讨论分子结构的影响因素有以下7个方面：3000—2850cm-1 3100-3000cm-1 3300 cm-1 附近-二元酸两个按基之间只有1一2个碳原子时，会出现两个UC=0相隔3个碳原子以上则没有这种偶合。
傅里叶红外光谱分析课件ppt
2/17/2024
经营者提供商品或者服务有欺诈行为 的，应 当按照 消费者 的要求 增加赔 偿其受 到的损 失，增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
1.2、红外吸收光谱产生的条件
condition of Infrared absorption 满足两个条件： spectroscopy
(1)辐射应具有能满足物质产生振动跃迁所需的能量； (2)辐射与物质间有相互偶合作用。
对称分子：没有偶极矩，辐 射不能引起共振，无红外活性。 如：N2、O2、Cl2 等。
非对称分子：有偶极矩，红 外活性。
偶极子在交变电场中的作用示 意图
(动画)
2/17/2024
经营者提供商品或者服务有欺诈行为 的，应 当按照 消费者 的要求 增加赔 偿其受 到的损 失，增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
2c
12 / 2
正己烯中C=C键伸缩振动频率实测值为1652 cm-1
2/17/2024
经营者提供商品或者服务有欺诈行为 的，应 当按照 消费者 的要求 增加赔 偿其受 到的损 失，增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
1.3、分子中基团的基本振动形式
basic vibration of the group in molecular
2.2.2. 傅里叶变换红外光谱仪结构（框动图画）
干涉仪
样品室
检测器
光源
计算机
显示器 绘图仪
干涉图 FTS
光谱图
2/17/2024
经营者提供商品或者服务有欺诈行为 的，应 当按照 消费者 的要求 增加赔 偿其受 到的损 失，增 加赔偿 的金额 为消费 者购买 商品的 价款或 接受服 务的费 用
《傅立叶谱与反应谱》课件
图像处理
傅立叶谱在图像处理中有着重 要应用，可用于图像增强、边 缘检测等领域。
通信系统
傅立叶谱在通信系统中广泛应 用，用于频谱分析、信号调制 解调等关键技术。
反应谱的计算方法
反应谱的计算是衡量结构抗震性能的重要步骤。它涉及多个参数和计算方法，如地震动输入、结 构响应计算与频响函数的估算。不同的方法可以得到不同的反应谱曲线。
转换关系
傅立叶变换可以将时域信号转换为频域信号，而反应谱则是对结构在不同频率上响应的表征。
频率成分
傅立叶谱展示了信号在不同频率上的成分，而反应谱则描述了结构对地震波的频率敏感性。
抗震效果
通过分析傅立叶谱和反应谱的关系，可以评估结构对不同频率地震波的抗震效果。
傅立叶谱的应用
傅立叶谱在信号处理与频域分析中有着广泛的应用。它不仅用于分析信号的频谱特征，还在许多领域中 起到重要作用，如音频处理、图像处理、通信系统等。
1 地震动输入
根据设计地震动参数，选择合适的地震波输入作为反应谱计算的基准。
2 结构响应计算
通过动力分析等方法，计算结构在不同地震波作用下的响应。
3 频响函数估算
根据结构的动力性能参数和振型，估算结构的频响函数，用于计算反应谱。
结论与展望
傅立叶谱与反应谱是频域分析中重要的工具和方法。它们在信号处理和结构抗震设计中发挥着重要作用。 未来，随着技术的不断进步，傅立叶谱和反应谱的应用将更加广泛，为解决实际问题提供更多可能性。
1
频域分析
通过傅立叶谱，可以将信号从时域转换为频域，揭示信号的频率成分。
2
频谱特征
傅立叶谱能够展示信号在不同频率上的能量分布情况，为信号特征的分析提供了依 据。
反应谱
) x y 此时结构加速度为 (
x m
kx
k
c
) m( x y cx
(t ) y
绝对位移：x+y，相对位移：x
) cx kx 0 m( x y
地震反应振动方程
cx kx m m x y
先看没有地面运动，称为阻尼自由振动
cx kx 0 m x
改写方程，引入自振频率和阻尼常数
0 k / m, c / 20m c / 2 km
得到
(c / m) x (k / m) x x
2 x 2 0 x 0 x 0
代入解
x Ce
得到
位移反应谱； 速度反应谱； 加速度反应谱（绝对）
反应谱计算公式：
u t 1
0 t u e sin t d g 0 t t
0 t u u e g 0 0 t u g t u
EPV Rv / 2.5
Ra为阻尼5%的加速度反应谱在2～10Hz频率范围内 平滑化的平均值 ，Rv为1Hz附近平滑化的速度反应
谱的幅值，
反应谱的数值计算方法（不讲）：
单步逐步递归法； 两步逐步递归法； 中心差分法； 纽马克广义加速度法； 威尔逊θ法。
cos 1 sin 1
2
2 2 0u 0 g 由u u u 2 u g (2 0 u 0 u u)
02
t
0
g (t )e i0 ( t ) {2 cos[ (t ) ] sin (t )}d u
地震工程学-傅里叶变换-反应谱计算
一、得到地震波数据 (1)二、地震波分析 (2)2.1时程曲线绘制 (2)2.2傅里叶谱的绘制 (3)2.3反应谱结果分析 (8)[参考文献] (17)附件说明 (17)一、得到地震波数据访问/网站，下载相应的附录说明文件，在这些文件中可以找到相关地震的信息，便于有效的规划检索词，可以看到题目需要的地震波RSN 编号为6,；根据所查到的RSN编号直接检索，如图1.2所示输入RSN号即可检索得到检索结果下载即可得到所需的地震波，如图1.3所示图1.1查询相应的编号图1.2检索条件具体地震波数据见附件1.二、地震波分析2.1时程曲线绘制根据得到的地震波数据，进行MATLAB程序编制，绘制的到相应的竖向地震、180度地震程曲线，270度地震的加速度时程曲线，分别为图2.1.1，2.1.2,2.1.3.图1.3检索结果图2.1.1竖向地震加速度时程曲线图2.1.2 180度地震加速度时程曲线图2.1.3 270度地震加速度时程曲线2.2傅里叶谱的绘制根据离散傅里叶的变换准则可以得到A k（k=0,1，…,N/2）、B k(k=1,….N/2-1)，进而计算得到相应的参数T∗√A k2+B k2；式2.12); 式2.2φ=acrtan(−B kA kf=k式2.3TS k=T×C k2式2.4式中N=1024，T=10.24s;在具体计算时采用了MATLAB中的FFT函数，并且对参考文献中[1]中的例题波进行了试算得到了理想的结果；故可证明程序中所使用的算法是没有问题的；得到的结果见下。
图2.2.1 竖向地震波的傅里叶幅值曲线图2.2.2 180度地震波的傅里叶幅值曲线图2.2.3 270度地震波的傅里叶幅值曲线图2.2.5 180度地震波的傅里叶相位曲线图2.2.6 270度地震波的傅里叶相位曲线图2.2.7竖向地震波的功率谱曲线2.3反应谱结果分析根据地震波是可以求得相应的不同周期的单自由度体系的各种反应谱的，这包括：位移反应谱，速度反应谱，准速度反应谱，加速度反应谱，绝对加速度反应谱，具体的原理是根据杜哈梅积分对单自由度进行积分获得，相应的公式如下[2]x(t)=|−1ωd ∫a(τ)e−εω(t−τ)sin (ωd(t−τ))tdτ|max式2.5ẋ(t)=|−ωωd ∫a(τ)e−εω(t−τ)cos (ωd(t−τ)+α)tdτ| max式2.6ẍ(t)=|−ω2ωd∫a(τ)e−εω(t−τ)sin(ωd(t−τ)+2α)tdτ| max式2.7其中，ωd=ω√(1−ε2),式中ε为阻尼比，在本题中分别取为0和0.05。
傅里叶红外光谱分析原理与方法
特征频率
• 而且，整个分子的红外吸收光谱与其分子结构有一一对应 关系
- 分子的身份证，指纹 • 特征吸收峰：代表基团存在并具有较强强度的吸收峰
这些特征频率是如何得到的？
• 理论计算 • 从大量的同类型物质进行对比总结 • 最有说服力的方法：
当两个相同原子和一个中心原子相连时
亚甲基
H CH
二氧化碳
O=C=O
其伸缩振动有两种方式：
如果两个相同（H）原子同时沿键轴；离开中心
（C）原子，则称为对称伸缩振动，用符号vs表示
如果一个（HⅠ ）原子移向中心（ C ）原子，
而另一个（HⅡ ）原子离开中心（C)原子，
则称为反对称伸缩振动，用符号vas表示
y
有效
有效
因此 线性分子的振动自由度为3N-5
非线性分子（H2O)有三种转动方式：转动 自由度=3
x
H
H
y
H
H
H zH
因此 非线性分子的振动自由度为3N-6
小结:
对于由N个原子组成的 分子,非线形者有3N-6个振 动自由度,线形者有3N-5个振动自由度。对于有M个自 由度的分子的振动状况，可认为是由M个互相独立的基
• 但二者各有优缺点。
一个光谱的例子
• 波长与透光度的关系 • 每一个吸收峰对应于一个振动模式
T1=0.9; T2=0.15
T1=0.99; T2=0.01
红外光谱原理
• 用一束红外光（连续波长）照射试样； • 若其频率相应的能量与某个分子的振动或转动
能量差相当时，就被分子吸收； ΔE=E’’-E’ =hv=hc/λ
反应谱和傅里叶谱+地震波选取
【拓展知识1-2】功率谱,反响谱和傅里叶谱,地震波选取,地震持续时间确定功率谱功率谱是功率谱密度函数的简称.对于一般情况的随机振动,具时间历程具有明显的非周期性,具有连续的多种频率成分,每种频率有对应的功率或能量, 用图像来表示这种关系,称为功率在频率域内的函数,简称功率谱密度.加速度功率谱是对地震动加速度时程进行快速傅里叶变换〔FFT〕得到的［1］.对于非平稳随机过程,功率谱密度的单位是G的平方/频率.G指的是随机过程.对于加速度功率谱,加速度的单位是m/s2,那么功率谱密度的单位是〔m/s2〕2/Hz, Hz的单位是1/s,故加速度功率谱密度的单位为m2/s3.加速度功率谱密度函数曲线下方的面积代表随机加速度的总方差,即加速度功率谱可以理解为“随机加速度方差的密度分度〞.参考文献［1］庄表中.随机振动入门.科学出版社,1981.反响谱和傅里叶谱反响谱〔earthquake response spectrum,是单自由度弹性系统对于某个实际地震加速度的最大反响〔可以是加速度、速度和位移〕和体系的自振特征〔自振周期或频率和阻尼比〕之间的函数关系.反响谱是地震工程中分析结构和设备在地震中的性能的非常有用的工具,因为许多主要表现为简单的振荡器〔也称为单自由度系统〕.因此,如果能找出结构的固有频率,那么建筑的峰值响应可以通过从地面响应谱中读取相应频率的值来估计.在地震区域的大多数建筑标准中,这个值构成了计算结构必须反抗的力的根底〔地震分析〕.如前所述,地面响应谱是在地球自由外表所做的响应图.如果建筑物的响应与地面运动〔共振〕的组成局部“协调〞,可能会发生重大的地震破坏,这些成分可以从响应谱中识别出来.傅里叶谱,全称为傅里叶振幅谱.地震波是在时间上连续的随机过程, 地震动记录仪是根据一定的采样频率得到该连续曲线上离散的点,想要复原这个曲线, 可以通过解N 元1次方程组,更简洁有效的方式是采用有限傅里叶级数来近似 原始的时间历程.对这个近似的函数进行物理意义的探讨, 傅里叶级数或者说傅里叶变换是将原始 的随机波分解成多个不同周期波的叠加.描述这些地震波分量的频率与振幅的关 系的直方图是一种“谱〞,对这种谱的纵坐标乘以T/2秒,得到傅里叶振幅谱, 简称为傅里叶谱［2］.傅里叶振幅谱可以使加速度谱,也可以使速度谱或位移谱.事实上,严格意义的傅里叶谱应该是直方图, 由于数据的个数是有限的,与各振 型相对应的各离散频率之间,存在的信息是位置的,因此,将举行的顶点用折线 连起来并没有实际意义.但一般大家仍采用折线作图,故成为一种默认.参考文献 ［2］大崎顺彦.地震动的谱分析入门地震波选取地震动具有强烈随机性,分析说明,结构的地震反响随输入地震波的不同而 差距很大,相差高达几倍甚至十几倍之多. 故要保证时程分析结果的合理性, 必 须合理选择输入地震波.归纳起来,选择输入地震波时应当考虑以下几方面的因素：峰值、频谱特性、地震动持时以及地震波数量,其中,前三个因素称为 地震动的三要素.1、峰值调整地震波的峰值一定程度上反映了地震波的强度,因此要求输入结构的地震波 峰值应与设防烈度要求的多遇地震或罕遇地震的峰值相当,否那么应按下式对该 地震波的峰值进行调整.其中,A'(t)和A'max 分别为地震波时程曲线与峰值,A'max 取设防烈度要求A'(t)A' max A max A(t) 式(1)的多遇或罕遇地震的地面运动峰值；A〔t〕和A max分别为原地震波时程曲线与峰值.2、频谱特性频谱即地面运动的频率成分及各频率的影响程度.它与地震传播距离、传播区域、传播介质及结构所在地的场地土性质有密切关系.地面运动的特性测定说明,不同性质的土层对地震波中各种频率成分的吸收和过滤的效果是不同的. 一般来说,同一地震,震中距近,那么振幅大,高频成分丰富,震中距远,那么振幅小,低频成分丰富.因此,在震中附近或岩石等坚硬场地土中,地震波中的短周期成分较多,在震中距很远或当冲积土层很厚而土质又较软时,由于地震波中的短周期成分被吸收而导致长周期成分为主.合理的地震波选择应从两个方面着手：1〕所输入地震波的卓越周期应尽可能与拟建场地的特征周期一致.2〕所输入地震波的震中距应尽可能与拟建场地的震中距一致.3、地震动持时地震动持时也是结构破坏、倒塌的重要因素.结构在开始受到地震波的作用时,只引起微小的裂缝,在后续的地震波作用下,破坏加大,变形积累,导致大的破坏甚至倒塌.有的结构在主震时已经破坏但没有倒塌,但在余震时倒塌, 就是由于震动时间长,破坏过程在屡次地震反复作用下完成,即所谓低周疲劳破坏.总之,地震动的持续时间不同,地震能量损耗不同,结构地震反响也不同.工程实践中确定地震动持续时间的原那么是：1〕地震记录最强烈局部应包含在所选持续时间内.2〕假设仅对结构进行弹性最大地震反响分析, 持续时间可取短些；假设对结构进行弹塑性最大地震反响分析或耗能过程分析,持续时间可取长些.3〕一般可考虑取持续时间为结构根本周期的5倍〜10倍.4、地震波数量输入地震波数量太少,缺乏以保证时程分析结果的合理性；输入地震波数量太多, 那么工作量较大.研究说明,在充分考虑以上三个因素的情况下,采用3条〜5条地震波可根本保证时程分析结果的合理性.。
傅里叶谱
k 0 k0 k N /2 k N /2
2 π( N / 2)m 2πkm N cos cos N N m0 0
(i)
式(i)实际上是式(4.9)第 1 式分别取 l =0 和 l = N /2 的情况。自身相同的项经相乘后相加 时，在式(4.9)中为 N /2。而对于稍有不同的式(i)，有趣的是，它不再为 N /2，而是 N 。 准备工作到此才大体完成。现返回本来的目的：求满足式(4.8)， 且由式(4.7)表示的 N 个 系数。
取它为某些特殊的偶数。以后就会知道，这样取是有好处的。 要决定一个函数，使它全部通过 N 个样本值 xm ，实际上有多种方法。这里采用一种三 角函数的方法来进行说明。 一般当
A0 , A1 , A2 , , Ak , B0 , B1 , B2 , , Bk ,
为常数时，把以下形式的级数称为三角级数，
N 1
综上，式(h)表示的三角函数乘积的和为
0 S N 2
k l k l
接着，来看当 k , l 为 l k ， l N / 2 1 的整数时， sin 同类乘积的和：
- 21 -
sin
m 0
N 1
2πlm 2πkm sin N N
其结果和式(h)对 cos 同类乘积的情况一样，当 l k 时为 0，仅当 l k 时为 N / 2 。 最后看 sin 与 cos 乘积的和：
t
2π t T
或根据式(4.1)进行如下代换：
t
2π t NΔt
则得到的还是三角级数：
- 18 -
A cos B sin NΔt NΔt
k 0 k k
地震工程学-傅里叶变换-反应谱计算
一、得到地震波数据 (1)二、地震波分析 (2)时程曲线绘制 (2)傅里叶谱的绘制 (3)反应谱结果分析 (8)[参考文献] (17)附件说明 (17)）一、得到地震波数据访问网站，下载相应的附录说明文件，在这些文件中可以找到相关地震的信息，便于有效的规划检索词，可以看到题目需要的地震波RSN 编号为6,；根据所查到的RSN编号直接检索，如图所示输入RSN号即可检索得到检索结果下载即可得到所需的地震波，如图所示图查询相应的编号图检索条件图检索结果具体地震波数据见附件1.二、地震波分析时程曲线绘制根据得到的地震波数据，进行MATLAB程序编制，绘制的到相应的竖向地震、180度地震程曲线，270度地震的加速度时程曲线，分别为图，图竖向地震加速度时程曲线图 180度地震加速度时程曲线傅里叶谱的绘制根据离散傅里叶的变换准则可以得到A k （k=0,1，…,N/2）、B k (k=1,….N/2-1)，进而计算得到相应的参数T2∗√T T T T ； 式φ=acrtan (−TT T T); 式T =TT 式T T =T ×T T 2 式式中N=1024，T=;在具体计算时采用了MATLAB 中的FFT 函数，并且对参考文献中[1]中的例题波进行了试算得到了理想的结果；故可证明程序中所使用的算法是没有问题的；得到的结果见下。
10图竖向地震波的傅里叶幅值曲线图 180度地震波的傅里叶幅值曲线图 270度地震波的傅里叶幅值曲线图竖向地震波的傅里叶相位曲线图 180度地震波的傅里叶相位曲线图竖向地震波的功率谱曲线图 180度地震波的功率谱曲线反应谱结果分析根据地震波是可以求得相应的不同周期的单自由度体系的各种反应谱的，这包括：位移反应谱，速度反应谱，准速度反应谱，加速度反应谱，绝对加速度反应谱，具体的原理是根据杜哈梅积分对单自由度进行积分获得，相应的公式如下[2]x(t)=|−1T T ∫T(T)T−TT(T−T)sin(T T(T−T))TTT|TTT式x(t)=|−TT T ∫T(T)T−TT(T−T)cos(T T(T−T)+T)TTT| TTT式x(t)=|−T2T T∫T(T)T−TT(T−T)sin(T T(T−T)+2T)TTT| TTT式2.7其中，T d=T√(1−T2),式中ε为阻尼比，在本题中分别取为0和。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
kx
k
c
cx
y(t)
此时结构加速度为 (xy)
x
m
m(xy)
kx
k
c
cx
y(t)
绝对位移：x+y，相对位移：x
m ( x y ) c x k 0 x
此地解震称反为应杜振哈动梅方积分程。
m x c x k x m y
x(t)0 tm f( )e0 (t )si n (t)d
埃森特罗地震记录南北分量时程（上）与傅立叶 振幅谱（左下）和相位谱（右下）
加速度/Gal
400 300 200 100
0 -100 -200 -300
0
10
20
30
40
50
时间/sec 160
傅立叶幅值谱/cm/sec
140
120
/2
Hale Waihona Puke /4 100相位/弧 度
0 80
/4
60 / 2
40
0
5
tk k t,a k a tk ,n n ,
可得
1 N1
i2N k n
F ae n
k
Nk0
N1 i2kn
ak Fne N n0
k0,1,2,,N1 k0,1,2,,N1
2
N t
N
2
N
2 t
解释：
1）t 称为采样间隔，时间步距；f 1 1 为频率间
隔，或频率的分辨率。
Nt T
反应谱表征了地震动加速度时间过程作用于单 自由度弹性体系的最大反应（加速度、速度和 位移）随体系的自振特性（周期、阻尼比）变 化的函数关系。
表示地震反应的物理量不同，反应谱的形 式也不同，如加速度、速度、位移等
位移反应谱； 速度反应谱； 加速度反应谱（绝对）
反应谱计算公式：
ut 1
t
ug e0t sin t d
3、频谱分析的结果可以方便地在结构抗震、基础 隔震等工程问题中应用。
二、反应谱（地震工程的灵魂）
回顾单自由度的地震反应分析，这是结构 地震反应分析最简单，最本质的分析模型。
m
k
c
x
m
k
c
y(t)
此时没有外力作用，结构随地面一起
运动，位移是（x+y），x，y分别为
相对位移和牵连位移。
x
m
m(xy)
2）频率的分辨率取决于时间步距，间隔越小，或时 间步距固定，则采样点（可以补零）越多，频率分 辨率越高。
T
f T
f T
时程 x(t)
傅立叶幅值谱 F(f)
3）时间分割和频率分割精度互相矛盾，是测不准 原理的体现。
4）由于 FnFNn（共轭），故 N 为对称点，称 2
为折叠频率，
可得
N
2
N
2 t
如何估计结构地震反应？
基于两个简化： 1）结构体系模型用单质点模型，计算简单 2）用一个不变的（静）力表示，便于设计 应用
方法：解单质点对地震动的反应 只选取加速度最大值 加速度乘以质量就是惯性力
与重力之比称为地震系数 F*W
反应谱就是以不同单质点自振频率（或周期） 为横坐标，反应最大值为纵坐标画出的曲线， 因为与频率（周期）对应，亦是谱的一种。
2）主要分量在0—10Hz之间
3）因为是离散计算，曲线毛刺多，可以平 滑化，如红线示意。
分解对研究地震动特性具有重要的意义
1、它定量揭示地震动的动力特性。以不同频率 分量的表现来研究地震动及其对结构的作用， 是动力分析的特点； 2、用严格的数理方法求解结构地震反应时，简 谐振动反应容易求解。对线性体系，先求解每个 简谐振动的反应，再叠加求得总体反应的方法， 比直接求解便利得多；
此处加标题
傅立叶谱与反应谱
眼镜小生制作
什么是频谱？
频谱是针对以时间为自变量的物理量变化 函数而言，例如将任意复杂的振动分解为不同 频率的简谐分量。
一般情况，变化非常复杂的物理量难以直接 定量分析，将复杂的现象分解为许多简单现象的 合成，化繁为简，化难为易。
如同白光能分解为各种频率的彩色光一样， 任何复杂的地震动也都可以分解为许多简谐振动， 这些简谐振动的振幅、初相位随频率变化。
fN/2 1/2t
若给定时间序列的采样间隔 t 则从采样后的离散信号中所能分辨出的最高频 率就是折叠频率：
fc 1(2t)
称为采样定理，折叠频率亦称奈奎斯特频率 （Nyquist frequency）
fR
例如加速度时程的时间步距为0.01sec，则离 散傅里叶谱的有效高频为50Hz。
加速度振幅谱特点： 1）开始（频率小)很小，迅速增大，然后 减小，高频和低频都小，意味某个频段幅 值大，这些振动分量强
10
15
20
25
30
频率/Hz
20
0
0
5
10
15
20
25
30
傅里叶谱是复数，由实部和虚部组成，它的模称 为幅值谱，幅角为相位谱
C () A () iB () C ()e i( )
C() A2()B2()
( ) a rc ta n [A ( )/B ( )]
实际上地震动时程无法用解析函数表示，因此要 用离散方法做数值计算
(-1)
将振动信号（或任意变化的函数）分解为 简谐振动（三角函数）的过程称为傅立叶 （Fourier）分解。得到振幅和相位随频率变化 的关系称为傅里叶谱，包括振幅谱和相位谱， 统称(大多指振幅谱)为频谱，完成分解的运算 称为傅立叶变换。
振动时程各不相同，傅里叶谱也互相各异， 有各自的特点，即频谱特性。振动信号的频谱 有重要的物理意义和应用。
一、傅里叶谱
在数学上可将任意变化的周期函数分解为正弦或 余弦变化的三角函数（称为简谐函数）之和，分解 是通过傅立叶展开或傅立叶变换完成的。
对随时间任意变化的物理量，分解具有物理意义： 例如对振动信号来说，时间变量的三角函数表示简 谐振动，这种分解就是将复杂的地震动化为无数简 谐振动之和。法国数学家傅立叶（J. Fourier）首 先研究将任意函数分解为三角级数之和的方法和条 件，并建立数学基础，后人发展并以其名字命名， 称为傅立叶级数。
0
ut 0
t
ug e0t cos t
0
d
ut
ug
t
02
t
ug e0t sin t 2d
0
对反应谱的注释：
反应谱是单自由度体系对地震动输入的地震 反应的最大值；
反应谱反映的是地震动的频谱特性； 反应谱还与体系的阻尼比有关。在反应谱计
时域和频域：时程和频谱
不同频率信号的时域图和频域图
复杂周期信号波形
从一般信号分析（数学化）的观点来看： 连续信号与离散信号 周期信号与非周期信号 周期信号：
ftTf(t) t
连续信号
f(t) 0
f(t)
f0
f1
t
t
0
f2
离散信号
f(tk)
(6)
(4.5)
(3) (1.5)
(2)
-1
t
01 2 3 4