Hi,
data = readtable("Data.csv");
t = table2array(data(:,1)); % Time Vector
s = table2array(data(:,2)); % Signal Vector
Ts = mean(diff(t)); % Sampling Interval
Fs = 1/Ts; % Sampling Frequency
Fn = Fs/2; % Nyquist Frequency
L = numel(t); % Signal Length (Number Of Samples)
FTs = fft(s - mean(s))/L; % Fourier Transform
Fv = linspace(0, 1, fix(L/2)+1)*Fn; % Frequency Vector
Iv = 1:numel(Fv); % Index Vector
figure
plot(Fv, abs(FTs(Iv))*2)
grid
xlabel('Frequency')
ylabel('Amplitude')
centroid = spectralCentroid(abs(FTs(Iv))*2 , Fv,'SpectrunType' , 'magnitude');
Hope this helps.
