三联反应谱的绘制程序.docx

三联反应谱绘制程序 同济大学 褚腾峰clear;clc;% 三联反应谱的绘制程序% step1 坐标轴的设置figureaxis(0.02,50,0.2,100);%横纵坐标界限 set(gca, XTickLabel,0.02,0.05,0.1,0.2,. 0.5,1,2,5,10,20,50,. XTick,0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.2 0.3 0.4 0.5 0.6. 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 ,XScale,log); set(gca,YTickLabel,0.2,0.5,1,2,5,. ,10,20,50,100,. YTick,0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 . 9 10 20 30 40 50 60 70 80 90 100 ,YScale,log); xlabel(Period T_n(s); ylabel(Pseudo-volocity V(in/s); hold ongrid onset(gca,gridlinestyle,-,minorgridlinestyle,-);% step2 垂直位移刻度的斜率为-1的斜线xinitial=0.02;yinitial=0.2;c=0.001:0.001:0.01 0.02:0.01:0.1 0.2:0.1:1 2:1:10 20:10:100 200:100:1000;n=length(c);a=zeros(1,2);b=zeros(1,2);a(1,1)=xinitial;b(1,2)=yinitial;for i=1:n z=2*pi*c(i);b(1,1)=z/xinitial;a(1,2)=z/yinitial;if c(i)=0.2 plot(a,b,b-,linewidth,0.8);else plot(a,b,k-); endhold onend% step3 垂直伪加速度刻度的斜率为1的斜线xlast=50;yinitial=0.2;d=0.001:0.001:0.01 0.02:0.01:0.1 0.2:0.1:1 2:1:10 20:10:100 200:100:1000;a=zeros(1,2);b=zeros(1,2);m=length(d);for j=1:mz=2*pi*d(j);b(1,2)=xlast*z;a(1,1)=yinitial/z;b(1,1)=yinitial;a(1,2)=xlast; if d(j)=0.2 plot(a,b,b-,linewidth,0.8); else plot(a,b,k-); endendbox on% step 4 添加斜向刻度值cc=0.2;dd=0.01 0.1 1 5 10 50 100;nn=length(dd);for k=1:nnx=sqrt(cc/dd(k);x=double(x);y=2*pi*sqrt(cc*dd(k);y=double(y);text(x,y,num2str(dd(k)/10),rotation,45,BackgroundColor,1 1 1);endee=0.2;ff=0.01 0.1 1 5 10 50 100;mm=length(ff)for k=1:mmx=sqrt(ff(k)/ee);x=double(x);y=2*pi*sqrt(cc*dd(k);y=double(y);text(x,y,num2str(dd(k),rotation,-45,BackgroundColor,1 1 1);end% step 5 输入数据画图load pseudo_sudu.txt; %加载伪速度反应谱数据psd=pseudo_sudu/2.54;%单位转化，由厘米转为英寸t=0:0.05:50-0.05;loglog(t,psd,linewidth,1.5); %双对数坐标text(2,1.8,Displacement D/in,rotation,45,Backgr