Bandpass filter of set of Data
34 views (last 30 days)
I want to filter some seismic data I have like an array (27*4352) elements and I want to filter this data with Pass-Band filter so only Frequency between 100-2000 Hz will access I was reading a lot of codes and functions like Butterworth or filter(filter) that could help but I'm not able to perform the filter untill now? any suggestions to do that!
I use trying something but not sure!
% read the data set in mtx files
filename = ('shot_1.mtx')
seismogram = readSeismogram(filename);
Fs = 1/DT;
T = 0.135;
% design filter
[b_pass, a_pass] = butter(2,[100,2000]/(Fs/2), 'bandpass');
plot(filtfilt(b_pass, a_pass, seismogram));