mftUntitled

从天然地震数据中提取瑞雷面波的群速度频散曲线

fid=fopen（‘C:\Documents and Settings\Administrator\桌面\SEGNAL.DAT‘‘r‘）;
x=fscanf（fid‘%f %f‘[2 inf]）;
y=x（2:）;
N=length（y）;
dt= 0.4153;
alpha=50.3;
beta=3.15;

n=0:N-1;t=n*dt;f=n/（N*dt）;


nfft=2^nextpow2（N）;
Y=fft（ynfft）; 


w=0.01:0.001:0.2;
i=size（w2）;
tt=zeros（1i）;

for j=1:i
    
    fmin=（1-（beta/alpha）^0.5）*w（j）;
    fmax=（1+（beta/alpha）^0.5）*w（j）;
      
    
    for（ m=0:N-1）
      
        if（m/（N*dt）>fmin&m/（N*dt）           yy（m+1）=Y（m+1）*exp（-alpha*（（（m+1）/（N*dt）-w（j））^2/w（j）^2））;
        
        elseif（m/（N*dt）>（1/dt-fmax）& m/（N*dt）<（1/dt-fmin））
           yy（m+1）=Y（m+1）*exp（-alpha*（（（m+1）/（N*dt）-1/dt+0.5*（fmin+fmax））^2/w（j）^2））;
       
         else
             yy（m+1）=0;
       
        end
      end 
  
    ZZ=abs（（ifft

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        994  2011-05-22 19:23  mftUntitled.m

     文件     100110  1996-11-08 11:01  SEGNAL.DAT

----------- ---------  ---------- -----  ----

               101104                    2