以下是matlb源程序.docVIP

以下是matlab源程序 =========================================== %%%%%%%%%%%%%%人造地震波程序%%%%%%%%%%%%% % 本程序利用反应谱转换为功率谱方法模拟人造地震波 % 变量说明： % s0: 谱强度因子 dt: 地震波采样时间间隔 % kao: 地震波总时间长度 % Ts: 平稳地震动的持续时间 % t1,t2,C: 控制主震段首末时间和衰减快慢的参数 % keceg,omigag: 地表覆盖层的阻尼比和卓越频率 % omigah: 反映基岩特性的谱参数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% xg(1:1500,1)=zeros(1500,1); dt=0.02; kao=21.9; t1=1.5; t2=15; c=0.18; Domiga=2*pi*(25-1/6)/100; Tg=0.25; kece=0.05; alfa=0.9; for i=1:100 omiga=(2*pi/6)+Domiga*(i-0.5); T=2*pi/omiga; if Tlt;0.1 r=5.5*alfa*T+0.45*alfa; elseif (Tgt;=0.1Tlt;=Tg) r=alfa; else r=alfa*(Tg/T)^0.9; end gama=omiga/pi; M1=(t2-0.3*t1)/kao; M2=(exp(-2*kece*omiga*t2)-exp(-2*kece*omiga*t1))/(2*kece*omiga*kao); M3=(1-exp(-c*(kao+0.5*t1-t2)))/(c*kao); M=M1+M2+M3; Gomiga=4*kece*r^2/(pi*omiga*M*(sqrt(2*log(gama*kao))+0.5772/sqrt(2*log(gama*kao)))^2); ak=2*sqrt(Gomiga*Domiga); qrand=rand(1)*2*pi; for j=1:1500 xg(j)=xg(j)+ak*sin(omiga*j*dt+qrand); end end for i=1:1500 if i*dtlt;t1 fai=(i*dt/t1)^2; elseif (i*dtgt;=t1i*dtlt;t2) fai=1; else fai=exp(-c*(i*dt-t2)); end xg(i)=fai*xg(i); end subplot(2,1,1); plot(xg); %xg=xg./(max(abs(xg))); save ren xg -ascii; fans=fft(xg); fans=fans.*conj(fans); subplot(2,1,2); plot(fans);