MATLAB 程序精确法求解反应谱 1. 反应谱的概念 反应谱是在 1932 年由 M.A.Biot 引入的，它是用来描述地面运动及其对结构的效应的一种实用工具。现在， 反应谱作为地震工程的核心概念，提供了一种方便的手段概括所有可能的线性单自由度体系对地面运动的某个 特定分量的峰值反应。它还提供了一种实用的方法，将结构动力学的知识应用于结构的设计以及建筑规范中侧 向力条文的制定。某个反应量的峰值作为体系的固有振动周期 Tn,（或者循环频率 fn）那样的相关参数的函数图形 ，称为该反应量的反应谱。每一个这样的图形针对的是有一个具有固定阻尼比的单自由度体系，多个具有不同阻尼比值范围。 2. 反应谱的计算 2.1 反应谱数值计算方法 计算反应谱的方法有很多，又卷积计算法，傅立叶变换法，线性加速度法，中点加速度法，精确法等。 2.2 精确法 本文中采用精确法做计算，该方法是 N.C.Nigam 和 P.C.Jennings 于 1969 年提出的，此法的出发点是把地面运动的加速度记录相邻点间的值用分段线性差值表示，从而获得地面运动的连续表达式。基于方程本身基础上进行，入误差，也不会产生任何的截断误差，所谓精确法就是指在这个意义上式精确的而然。正因为这种方法不会引起数值运算就可以达到采用其他方法需要较多次运算才能达到的精度。 “由于在sohu 博客上的文章发表后，陆续有问参考文献的邮件，因此将参考文献pdf 版放上来供大家学习、参考，地震动的谱分析入门 强震观测与分析原理 ” ResponSespectrumProgram（精确法求解） %%%%%%%%%%%%%%%%%%%%% 反应谱 精确法 程序 Begin With matlab6.5%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear % ***********读入地震记录*********** fid = fopen('CHI010.txt'); [Accelerate,count] = fscanf(fid,'%g'); %count 读入的记录的量 Accelerate=9.8*Accelerate'; %单位统一为 m 和 s time=0:0.005:(count-1)*0.005; %单位 s % ***********精确法计算各反应*********** %初始化各储存向量 Displace=zeros(1,count); %相对位移 Velocity=zeros(1,count); %相对速度 AbsAcce=zeros(1,count); %绝对加速度 % ***********A,B 矩阵*********** DampA=[0,0.05,0.1]; %三个阻尼比 TA=0.0:0.05:6; %TA=0.000001:0.02:6; %结构周期 Dt=0.005; %地震记录的步长 %记录计算得到的反应，MDis 为某阻尼时最大相对位移，MVel 为某阻尼 %时...
