地震工程作业(完整资料）(可以直接使用，可编辑优秀版资料，欢迎下载）作业1绘制1９40ElCｅｎtro，N—S分量地震动的绝对加速度、相对速度和相对位移反应谱。地震动：在PEＥRGrouｎdMotｉoｎＤatabａsｅ自行下载经典的１9４0ＥlＣenｔrｏ，N—S分量.要求:在此模板内完成，A4纸打印。自编程序与软件(Biｓpｅｃ或Ｓeisｍosignａ等)计算反应谱进行对比。提交自编写程序。Ｍaｔｌab程序:clearfid=ｆopｅn('E：\Earｔhqｕake\Elcentro.ｔxｔ')；[Accｅlerａte，count］=fscaｎf（fid,’％g’）；%ｃount读入的记录的量Acceleｒaｔe=9．8＊Accelerａte’;%单位统一为m和ｓtｉmｅ=0:0。02：(ｃｏunｔ－1)＊0．02；%单位sＤispｌace=zeroｓ（１,ｃounｔ)；％相对位移Ｖelocｉtｙ＝zeｒos(1,coｕnt）;％相对速度AbｓAｃce=ｚeros（１，ｃount）；％绝对加速度DamｐA=[０。０0,0.02，０.０5]；%三个阻尼比ＴA=０．0：0。0２：４;%TＡ=0．000001：0．0２:4；结构周期Dt＝0。０2；％地震记录的步长％记录计算得到的反应，ＭDiｓ为某阻尼时最大相对位移%MVｅl为某阻尼时最大相对速度,MＡcｃ某阻尼时最大绝对加速度ＭDｉs＝zeroｓ（３，length(TA)）；MVel＝zeros（3,leｎgtｈ(TA）);MAcc=zｅros(3，lengｔh（ＴA））；ｊ=1;fｏrDaｍp＝[0．00，０.０5,0.1]t＝1;ｆorT＝０．0：0.０2：4Frcｙ=２＊pi/T；DaｍＦrｃｙ=Frcｙ＊ｓｑrｔ(1-Dａmp*Ｄaｍｐ)；e_ｔ=exp(－Daｍｐ＊Ｆｒｃｙ*Ｄt）；ｓ=sin（DaｍＦｒcy*Dt）；c=cos(DamFｒｃy＊Ｄt)；A=zｅrｏs（２,2）;A（1，1)=e_t＊(s*Damp/ｓｑrt(１-Damｐ*Ｄamp）＋ｃ）；Ａ（1，2)=e_t*s／ＤａmFrcy；A（2，1）=-Frcｙ＊e_ｔ＊ｓ/sｑｒt（1－Dａmｐ＊Ｄaｍｐ）;A（２，2)＝e_t*（－ｓ＊Ｄamp/sqｒt（１—Ｄaｍp＊Damp）+c）；d＿ｆ=(2＊Damｐ＾2—1)／（Ｆrcy^2＊Dt）;ｄ_３t=Damp/(Frcｙ^3＊Ｄｔ）；B=zeros(2，2）;Ｂ(1,１）=ｅ_ｔ*（（d_ｆ+Damp/Ｆrcy)＊s/DamFrcy＋（2＊d_3ｔ+1/Frcy＾2)*c）-2＊d_3ｔ；Ｂ(１，2）=-ｅ_t＊（ｄ＿f＊s/DａmFｒcy+２*ｄ_３t＊c）-1／Frｃy^2+2＊d_3t；B（２，１）=ｅ_t＊（（d＿f+Damp/Frcｙ）＊（c-Dａmp／sqrt（1—Damｐ^２）＊ｓ)—（2＊d_3t+１/Ｆｒcy＾2）＊（ＤａｍFｒcy*ｓ+Dａｍp＊Ｆｒcy*c））+１/（Frcy^2*Ｄt);B（２，2）＝ｅ_ｔ*(1／（Ｆrcy＾2＊Dt）＊ｃ+s*Ｄａmp/(Fｒcy*DamFrcｙ＊Ｄt））—１/（Ｆrｃｙ＾2＊Dt）;fori=1：（couｎt-1）Disｐｌace（i＋1）=A（1，1）＊Diｓplａce(i）+A(1,2)*Velｏciｔy（i)＋Ｂ（１，1）*Acｃeleｒate（i)＋B（1，2）*Acceｌeraｔe（i+1）；Ｖｅlociｔｙ(ｉ+1)=A（２,１)*Diｓplａce（ｉ）+A(２,2)＊Vｅlocity(i）+B（2，1）＊Aｃceｌｅrａｔe(i)+B(2,２）＊Acｃｅlerate(i+1）；AbsAcｃe(ｉ+1）=－2＊Daｍp*Fｒcy*Velocｉｔｙ（i+１)—Frcｙ＾2＊Diｓplaｃe(i＋1);enｄＭDis（j,t）=mａx(abs（Disｐlace））；MVel(j，t)=ｍaｘ（aｂｓ(Velｏcｉｔy)）；ifＴ＝=0．0MＡcc（ｊ,t)=max(abs(Aｃceｌeratｅ)）；elseＭAｃc（j，t)＝mａx(abs(AbｓAcｃe））；enｄＤiｓｐlaｃe＝zerｏｓ（1,cｏunt）；Veｌoｃity=zｅｒos（１，couｎt);AbsAcｃe=zｅroｓ(1，cｏunt)；ｔ=ｔ＋１；endｊ=j+1;ｅnｄｃｌoseallfigｕre％绘制位移反应谱plot(TA，ＭDis(1，:），'—b’，ＴA，ＭDis（２，：）,'—r’，TA，ＭDis(３，:）,':k’)tiｔle(’Displacem