spektrum.m

function [Y, f] = spektrum(y, Fs)
% berechnet das Spektrum der Zeitreihe y (bei gegebener Samplerate fS)
% gibt die Spektralfunktion Y und die zugehörigen Frequenzwerte f zurück

N = length(y);
T = N/Fs;            % Zeitdauer von y
X = fft(y);
Y = abs(X(1:floor(N/2)))/N;
f = (0:(N/2)-1)/T;