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)
Mostrar comentarios más antiguos
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)
1 comentario
Raj Bhattacharjee
el 21 de Sept. de 2020
is indonesiaQ16.dat available? just to understand the program better...
Respuestas (0)
Ver también
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!