How do I set thresholds and implement STA-LTA algorithm (short time average through long time average) on a seismic data (for instance "IndonesiaQ16.dat")?

14 visualizaciones (últimos 30 días)
D=load ('IndonesiaQ16.dat'); D1 = detrend(D); dn=zeros(1,360); l=0; for i=1:10:3600 m=0;t=1:10:3600; l=l+1; for j=0:9 m=D1(i+j)+m; end dn(l)=m/10; end plot(t,dn)
g=abs(D1);
plot(g)
dn=zeros(1,360);
l=0;
for i=1:10:3600
m=0;t=1:10:3600;
l=l+1;
for j=0:9
m=g(i+j)+m;
end
dn(l)=m/10;
end
plot(t,dn)

Respuestas (0)

Categorías

Más información sobre Seismology en Help Center y File Exchange.

Productos

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by