diff --git a/semestralka/1.m b/semestralka/1.m index 7c61de6..1a3ca05 100644 --- a/semestralka/1.m +++ b/semestralka/1.m @@ -9,13 +9,24 @@ data = int16(str2double(nums)); X = DFT_priama(double(data)) n = length(X) -freq = (0:n-1)/n +Fs = 200; +freq = (0:n-1)*Fs/n figure; plot(double(data)); figure; plot(freq, X) figure; +x = double(data); +x = x - mean(x); +plot(x); % centrovane data +figure; plot(freq, abs(X)) figure; plot(freq, angle(X)) + +magX = abs(X); +[pks, locs] = findpeaks(magX, 'NPeaks', 6, 'SortStr', 'descend'); + +f0 = freq(locs); +disp(f0);