- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
三联反应谱绘制程序 同济大学 褚腾峰
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
您可能关注的文档
- 七年级下册语文生字及词语解释汇编.doc
- 七年级下第1课怎样确保骑车安全.doc
- 七年级下学期地理预习复习资料.doc
- 七年级下学期期中考地理试卷.doc
- 七年级下英语期末考试命题双向细目表.doc
- 七年级下语文第三单元导学案(教师版).doc
- 七年级下语文期中复习资料.doc
- 七年级新生入学考试语文试题.doc
- 七年级信息技术教学设计7-5.1 第一节 电子板报的设计.doc
- 七年级信息技术教学设计7-5.3 第三节 制作图文并茂的电子板报.doc
- 重庆新速达物业服务集团股份凯里公司招聘笔试题库2024.pdf
- 浙江杭州千岛湖泰众肉类食品有限公司招聘笔试题库2024.pdf
- 陕西新华出版传媒集团新华书店分公司招聘笔试题库2024.pdf
- 宁夏海原县肉牛产业发展集团有限公司招聘笔试题库2024.pdf
- 浙江温岭市化工轻工建筑材料有限公司招聘笔试题库2024.pdf
- 中航飞机汉中航空零组件制造有限公司招聘笔试题库2024.pdf
- 山东海阳市海昇矿产资源开发有限公司招聘笔试题库2024.pdf
- 山东潍坊寿光市公共交通运输有限公司招聘笔试题库2024.pdf
- 内蒙古察右前旗农业投资开发有限公司招聘笔试题库2024.pdf
- 线上清明祭英烈活动策划.pptx
文档评论(0)