STA-LTA 长短视窗处理地震信号算法

STA-LTA 长短视窗处理地震信号算法

clc;
clear;
Ap=0.3;fp=300;Np=300;h=-150;dt=0.0005;
As=1;fs=200;Ns=450;
for i=300:1000
 A1（i）=Ap*exp（h*（i*dt-Np*dt）） *sin（2*pi*fp*（i*dt-Np*dt））;
end

for i=450:1000
     A2（i）=-As*exp（h*（i*dt-Ns*dt）） *sin（2*pi*fs*（i*dt-Ns*dt））;
end

for i=1:1000
s（i）=A1（i）+A2（i）;
end



 

NOISE=randn（size（s））;

Y=s+0.03*NOISE;


A=abs（Y）;



for i=180:1000
    LTA1（i）=0;
    for n=i-60+1:i
    LTA1（i）=LTA1（i）+A（n）*A（n）;
    end 
end


for i=180:1000
    LTA3（i）=0;
    for n=i-180+1:i
    LTA3（i）=LTA3（i）+A（n）*A（n）;
    end 
end





for i=180:1000
    STA（i）=0;
    for n=i-20+1:i
   STA（i）=STA（i）+A（n）*A（n）;
    end 
end
sta=STA./20;
r1=3.*STA./LTA1;
lta1=LTA1./60;
r3=9.*STA./LTA3;
lta3=LTA3./180;
THR（1:1000）=3;
figure（1）
subplot 511 ;plot（Y）;xlim（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1436  2016-04-25 21:08  duanshic.m
     文件         933  2016-04-25 21:21  jianquanmer.m
     文件         588  2016-04-25 21:20  jiaquanxishu（1）.m
     文件        1228  2016-04-25 17:28  lta.m
     文件        1224  2016-04-25 21:34  CHANG.m