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

三联反应谱绘制程序 同济大学 褚腾峰 clear;clc; %% 三联反应谱的绘制程序 %% step1 坐标轴的设置 figure axis([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 on grid on set(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-); end hold on end %% 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:m z=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-); end end box on %% step 4 添加斜向刻度值 cc=0.2; dd=[0.01 0.1 1 5 10 50 100]; nn=length(dd); for k=1:nn x=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]); end ee=0.2; ff=[0.01 0.1 1 5 10 50 100]; mm=length(ff) for k=1:mm x=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,BackgroundColor,[1 1 1]); text(0.33,2.9,Pseudo-acceler