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

地震波位移、速度和加速度反应谱绘制问题：绘制El-Centro-NS地震波位移、速度和加速度的反应谱答：读取El-Centro-NS地震波的加速度数据，使用Matlab进行编程，利用Matlab 内部函数 lsim 求解，绘制出4 种不同阻尼比的情况下的反应谱图， 4 种 阻尼比分别为 0、 0.02、 0.05、 1Matlab 程序如下：%%%%%%画出加速度时程曲线%%%%%%clearformat compactt=xlsread('C:\Users\shaoqunl23\Desktop\el-centro 地震波南北向数据.xlsx','Sheetl','Al:A2675'); data=xlsread('C:\Users\shaoqunl23\Desktop\el-centro 地震波南北向数据.xlsx','Sheetl','Bl: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;-wA2 -2*kesi(i)*w];B=[0;1];C=[1 0; 0 1];D=[0 ;0];[y,Z]=lsim(A,B,C,D,data,t);d=Z(:,1);v=Z(:,2);a=-2*kesi(i)*w*v-wA2*d;Sd (i,j)=max(abs(d));Sv(i,j)=max(abs(v));Sa(i,j)=max(abs(a));endend%% 画图figure(1)plot(T,Sd(1,:),T,Sd(2,:),T,Sd(3,:),T,Sd(4,:)) title('displacement response spectrum') xlabel('T(s)');ylabel('Sd(m)');legend ('kesi=0','kesi=0.02','kesi=0.05','kesi=1'); pause;figure(2) plot(T,Sv(1,:),T,Sv(2,:),T,Sv(3,:),T,Sv(4,:)) title('velocity response spectrum') xlabel('T(s)');ylabel('Sv(m/s)');legend ('kesi=0','kesi=0.02','kesi=0.05','kesi=1'); pause;figure(3) plot(T,Sa(1,:),T,Sa(2,:),T,Sa(3,:),T,Sa(4,:)) title('acceleration response spectrum') xlabel('T(s)');ylabel('Sa(m/sA2)');legend ('kesi=0','kesi=0.02','kesi=0.05','kesi=1'); pause;close;070.40.30r20.524.55displacement response spectrumkesi=Q kesi=0.02 kesi=0.05 kesi=1251.61.4120.80.60,40.2velocity response spectrum0 Q.5 1 1P5 2 2.5 3 3,5 4 4P5 5252015105acceleration response 导pectrum0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5。