反应谱.rar_冲击_冲击荷载_冲击谱_反应谱

%地震信号加速度反应谱 clear clc close all hidden format long %fni=input('地震信号加速度反应谱－输入数据文件名','s'); %fid=fopen(fni,'r') %sf=fscanf(fid,'%f',1); %读入采样频率值 %dp=fscanf(fid,'%d',1); %阻尼比 %hf=fscanf(fid,'%d',1); %计算最高截止频率 %df=fscanf(fid,'%d',1); %频率间隔 %fno=fscanf(fid,'%d',1); %读入输出数据文件名 %xg=fscanf(fid,'%f',[1 inf]); %读入地震波数据存成列向量 %status=fclose(fid); sf=500;dp=0.025;hf=50;df=0.09765625;fun=2;fno='out6_5.mat'; load y xg=y(1:1000); n=length(xg); m=round(hf/df)+1; f=0:df:(m-1)*df; t=0:1/sf:(n-1)/sf; for j=0:m-1 w=2*pi*df*j; wd=w*sqrt(1-dp*dp); e=exp(-t*w*dp); a=t.*wd; s=sin(a).*((1-2*dp*dp)/(1-dp*dp)); c=cos(a).*(2*dp/sqrt(1-dp*dp)); h=wd*e.*(s+c)/sf; r(j+1)=max(abs(conv(xg,h))); end subplot(2,1,1); plot(t,xg); xlabel('时间(s)'); ylabel('加速度(g)'); grid on; subplot(2,1,2); plot(f,r); xlabel('频率(Hz)'); ylabel('加速度(g)'); grid on; fid=fopen(fno,'w'); for k=1:m fprintf(fid,'%f%f\n',f(k),r(k)); end status=fclose(fid);