How can i get frequency domain of an earthquake?
21 views (last 30 days)
I need to get the frequency domain of an earthquake. I have used "FFT" in Matlab, but I could not get correctly.
the earthquake is attached here called EQ.txt.
dt=0.01 sec time step;
Can you please correct that?
Thanks for your help.
Star Strider on 27 Jan 2021
Edited: Star Strider on 27 Jan 2021
acc = readmatrix('EQ.txt');
L = numel(acc);
Ts = 0.01;
Fs = 1/Ts;
Fn = Fs/2;
t = linspace(0, L, L)*Ts;
FTacc = fft(acc)/L;
Fv = linspace(0, 1, fix(L/2)+1)*Fn;
Iv = 1:numel(Fv);
ylabel('Acceleration Amplitude (Units)')
ylabel('Acceleration Power (Units^2)')
EDIT — (27 Jan 2021 at 18:50)
Corrected typographical errors.