MATLAB: Spektrale Leistungsdichte

Das Leistungsdichtespektrum erhält man durch Fourier Transformation der 
AKF.
Also fft(...)

Da Matlab aber ausschließlich mit diskreten Stützstellen rechnen kann, 
wird man hier keine schöne Lösung als Formel erhalten, sondern höchstens 
eine numerische, und die Richtigkeit dieser hängt auch noch davon ab wie 
genau man das Eingangssignal beschrieben hat.

Die Funktion psd(...) dient zur SCHÄTZUNG der signalen Leistungsdichte 
bei gefensterten Zeitsignalen.

Wie müsste ich dann die fft Funktion an dem angegebenen Beispiel 
anwenden?

Da braucht man aber kein Matlab dazu. Bissel nachdenken reicht schon.

AKF: (sin(x)/x)^2

Spektrale Leistungsdichte ist Fourier-Transformierte der AKF.
Fourier-Transformierte von sin(x)/x ist ein Rechteck. Multiplikation im 
Zeitbereich (also bei der AKF) bedeutet Faltung im Frequenzbereich.
Fouriertransformierte von (sin(x)/x)^2 ist also Faltung von Rechteck mit 
Rechteck, was ein Dreieck ergibt.

Aber wie konnte man dies trotzdem mit Matlab grafisch darstellen?