matlab怎么画地震反应谱,地震工程学-反应谱和地震时程波的相互转化matlab编程...

一、地震波生成反应谱

1 所取的地震波为Elcentro地震波加速度曲线，如图1所示。

图1 Elcentro地震波加速度曲线

2 所调用的Matlab程序为：

% ***********读入地震记录*********** ElCentro;

Accelerate= ElCentro(:,1)*9.8067;%单位统一为m和s N=length(Accelerate);%N 读入的记录的量 time=0:0.005:(N-1)*0.005; %单位 s

%初始化各储存向量

Displace=zeros(1,N); %相对位移 Velocity=zeros(1,N); %相对速度 AbsAcce=zeros(1,N); %绝对加速度

% ***********A,B矩阵***********

Damp=0.02; %阻尼比0.02

TA=0.0:0.05:6; %TA=0.000001:0.02:6; %结构周期 Dt=0.005; %地震记录的步长

%记录计算得到的反应，MaxD为某阻尼时最大相对位移，MaxV为某阻尼最大相对速度，MaxA某阻尼时最大绝对加速度，用于画图

MaxD=zeros(3,length(TA)); MaxV=zeros(3,length(TA)); MaxA=zeros(3,length(TA)); t=1;

for T=0.0:0.05:6

NatualFrequency=2*pi/T ; %结构自振频率

DampFrequency=NatualFrequency*sqrt(1-Damp*Damp); %计算公式化简