function ex13()
xx = load("stoerung.dat");
t = xx(:,1);
x = xx(:,2);
dt = t(2) - t(1);     % sollten alle gleich sein
T = max(t);
N = length(x);

fN = 1/(2*dt)

F = fft(x);
power = abs(F(1:floor(N/2))).^2;
freq = [0:(N/2)-1]/T;
plot(freq,power);

xlabel("f [Hz]", "FontSize", 14);
ylabel("|F(\omega)|^2", "FontSize", 14);

% wo treten die Störfrequenzen im Aliasing auf?
f0 = 39.6;
fStoer = f0*(1:20);
index = find(fStoer > fN);
fAlias = 2*fN - fStoer(index);

plot(freq, power, "b-", fStoer, 0*fStoer, "g*", fAlias, 0*fAlias, "r*");
axis([200 500 0 1500])
xlabel("f [Hz]", "FontSize", 14);
legend("Spektrum", "Störfrequenzen", "Aliasfrequenzen")

F = getframe(gcf);
imwrite(F.cdata, "bild48.png");
