seimology - EW=fft(y3,nt) TD=1/fc+0.05*R Twin=2*TD fmax=50

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
Sheet1 Page 1 nt= y=randn(1,nt) ymean=mean(y) y1=y-ymean ystd=std(y1) y2=y1./ystd e=0.2 et=0.05 a=exp(1)/e b=-e*log(et)/(1+e*(log(e)-1)) t=0:nt Wt=a*t^b.*exp(-c*t) y3=y2.*Wt R=28.2 %km Bs=3.4 %km/s Dxi=1e7 %dyn/cm2 ro=2.67 % Mw=7.74 Mo=10^((Mw+10.7)*3/2) fc=4900000*Bs*(Dxi/Mo)^(1/3)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: EW=fft(y3,nt) TD=1/fc+0.05*R Twin=2*TD fmax=50 nto=round(log(2*fmax*Twin)/log(2)) nw=nt/2 dt=Twin/nt df=1/Twin Am=abs(EW)*dt f=df*(1:nw) S1=mean(Am^2) SF=sqrt(S1) SA=Am./SF SObs=SA*f Ts=ifft(SObs,nw) Tse=Ts(1:nw/2)*df...
View Full Document

This note was uploaded on 10/26/2008 for the course CIVIL 602 taught by Professor Duta during the Spring '08 term at Air Force Institute of Technology, Ohio.

Ask a homework question - tutors are online