Upload New File

freq_spectrum.m

+function [pxx,f,peak_f,max_f] = freq_spectrum(signal,fs)
+    [pxx,f] = periodogram(signal,[],[],fs); % Power spectrum
+    max_f = max(pxx); % Maximum power frequency
+    peak_f = f(pxx==max_f); %peak frequency
+end
\ No newline at end of file