understanding the given algorithm
5 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
In the given code ,why there are two filtering process.This is a part of STA/LTA algorithm for earthquake detection which was posted in matlab
clc
clear all
%% Data Inputs
ip = importdata("G:\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\Nepal\ADIB.HHN.dat"); % Data file path
Acc_EW = importdata("G:\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\Nepal\ADIB.HHN.dat");
Acc_NS = importdata("G:\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\Nepal\ADIB.HHN.dat");
Acc_ver = importdata("G:\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\STA-LTA-Algorithm-and-Seismometer-Trajectory-visualization-in-3D-master\Nepal\ADIB.HHN.dat");
Fs = 100; %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 for displacement
peakToPeakRange = max(fildat) - min(fildat);
dt = 1/Fs; %sampling time
nt = length(fildat); % length of the input signal
time = (1:nt).*dt; % time duration of the input signal
1 commentaire
Réponses (0)
Voir également
Catégories
En savoir plus sur Filter Analysis dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!