地震波位移、速度和加速度反应谱绘制.docVIP

地震波位移、速度和加速度反应谱绘制 问题：绘制EI-Centro-NS地震波位移、速度和加速度的反应谱 答：读取EI-Centro-NS地震波的加速度数据，使用Matlab进行编程，利用Matlab内部函数lsim求解，绘制出4种不同阻尼比的情况下的反应谱图，4种阻尼比分别为0、0.02、0.05、1。 Matlab程序如下： %%%%%%画出加速度时程曲线%%%%%% clear format compact t=xlsread(C:\Users\shaoqun123\Desktop\el-centro地震波南北向数据.xlsx,Sheet1,A1:A2675); data=xlsread(C:\Users\shaoqun123\Desktop\el-centro地震波南北向数据.xlsx,Sheet1,B1:B2675); %% 计算反应谱 kesi=[ 0 0.02 0.05 0.1]; T=0.1:0.1:5; Sd=zeros(4,length(T)); Sv=zeros(4,length(T)); Sa=zeros(4,length(T)); for i=1:length(kesi) for j=1:length(T) w=2*pi./T(j); A=[0 1;-w^2 -2*kesi(i)*w]; B=[0;1