【地震】基于STALTA算法检测地震P波（含三维地震仪轨迹的可视化和估计、S波到达时间）附Matlab代码
✅作者简介：热爱科研的Matlab仿真开发者，擅长毕业设计辅导、数学建模、数据处理、程序设计科研仿真。
🍎完整代码获取 定制创新 论文复现点击：Matlab科研工作室
👇 关注我领取海量matlab电子书和数学建模资料
🍊个人信条：做科研，博学之、审问之、慎思之、明辨之、笃行之，是为：博学慎思，明辨笃行。
🔥 内容介绍
一、引言
地震波的准确检测和分析对于地震学研究、地震预警等方面至关重要。P 波作为地震发生时最先到达的地震波，其快速准确的检测能够为后续的地震应对措施争取宝贵时间。STALTA（Short - Term Average over Long - Term Average）算法是一种常用的地震波检测算法，它通过比较短时间和长时间内地震信号的平均振幅来识别地震波的到达。本文将详细介绍基于 STALTA 算法检测地震 P 波，并进一步探讨三维地震仪轨迹的可视化以及 S 波到达时间的估计方法。
二、STALTA 算法原理
基本概念
⛳️ 运行结果
📣 部分代码
%ip = importdata('F:\Earthquake data\Magnitude 7\All Mag 7 vt\TW_vt.txt'); % Data file path
Acc_EW = importdata('ADIB.HHE.dat');
Acc_NS = importdata('ADIB.HHN.dat');
Acc_ver = importdata('ADIB.HHZ.dat');
Fs = 200; %sampling frequency
%% Signal Pre-Processing
%Filter Design
digfilt = designfilt('lowpassiir', 'PassbandFrequency', 20, 'StopbandFrequency', 25, 'PassbandRipple', 1, 'StopbandAttenuation', 60, 'SampleRate', 200);
% Filtering Data
Acc_EW_filt = filter(digfilt,Acc_EW);
Acc_NS_filt = filter(digfilt,Acc_NS);
Acc_ver_filt = filter(digfilt,Acc_ver);
Fhp = 0.8; % high pass filter cutofff frequency
[b1,a1] = butter(3,Fhp/Fs,'high'); %3rd order high pass butterworth filter
fildat = filter(b1,a1,Acc_ver); % filtered acceleration data
vel = cumtrapz(fildat)./Fs; % Integrating acceleration data for velocity
[b2,a2] = butter(3,Fhp/Fs,'high'); %3rd order high pass butterworth filter
fildat1 = filter(b2,a2,vel); % filtered velocity data
dis = cumtrapz(fildat1)./Fs; % Integrating velocity data