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")?
显示 更早的评论
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 个评论
Raj Bhattacharjee
2020-9-21
is indonesiaQ16.dat available? just to understand the program better...
回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Seismology 的更多信息
产品
2018-2-7
2020-9-21
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
