Matlab: Spektrum aus Zeitsignal mit FFT

[ETStudent (Gast)]

Hallo,

ich schreibe ein Skript, mit dem man das Spektrum eines beliebigen 
Zeitsignals geplottet bekommt. Ich nutze dafür die FFT von Matlab. Ich 
habe aber das Problem, dass mein Skript ungenau ist. Wenn ich 
beispielsweise ein Cosinussignal mit 50Hz angeben, wird mir im Spektrum 
nicht genau bei 50Hz ein peak angezeigt. Was kann ich tun oder mache ich 
falsch?

Hier der Code:

%FFT Beispiel anders programmiert
clc;
clear all;

%zeitvektor definieren
Startzeit = 0;
Abtastfrequenz = 600000; %Abtastdauer
Schrittweite = 1/Abtastfrequenz;
Endzeit = 40*10^(-3);
t = Startzeit:Schrittweite:Endzeit;

%Anzahl der Abtastpunkte
Abtastpunkte = length(t);

%Definition des Zeitsignals
Frequenz = 50;
y = cos(2*pi*Frequenz*t)+cos(2*pi*2*Frequenz*t)+cos(2*pi*3*Frequenz*t);

%plot des Zeitsignals
subplot(2,1,1)
plot(t,y);
title('Signal')
xlabel('Zeit (Millisekunden)')
grid on;

%FFT
AnzahlMesswerte = length(t)
AnzahlStuetzpunkteFFT = 2^nextpow2(AnzahlMesswerte*10)
Y = fft(y,AnzahlStuetzpunkteFFT)/AnzahlMesswerte;

%Definition der Frequenzachse
f = Abtastfrequenz*linspace(0,1,AnzahlStuetzpunkteFFT);

%plot des Spektrums
subplot(2,1,2)
plot(f,abs(Y(1:AnzahlStuetzpunkteFFT)));
xlim([0 400]);
grid on;

[meiner (Gast)]

%Definition der Frequenzachse
f = Abtastfrequenz*linspace(0,1,AnzahlStuetzpunkteFFT)

frequenzachse geht nur bis zur halben abtastfrequenz, bzw von -f/2 bis 
+f/2

Wo ist denn dein peak?

[meiner (Gast)]

Mehr Messzeit schadet auch nicht

[Plopp.. (Gast)]

Ein kleine Ahnung zur Fouriertransformation und FFT im Speziellen sollte 
man bringen. Einfach eine Library aufrufen ist etwas an der unteren 
Grenze.

Als kleine Uebung :
-Was ist der Unterschied von einem Linienspektrum
   zu einem kontinuierlichen Spektrum ?
-Wodurch ist die Bandbreite und die Aufloesung bestimmt ?
-Was sind die Eigenfunktionen der Fouriertransformation ?