FFT - Phase

Ja, da weiß man schon eher was Du meinst.
An den Diagrammen ist alles in Ordnung, alles ist normal. Trotzdem Du einen reinen Sinus von 10Hz eingegeben hast, besteht das Amplitudenspekrum nicht nur aus dieser einen Spektrallinie. Das sieht man in der logarithmischen Darstellung:
   
Ich nehme an, das hängt mit Digitalisierungsfehlern zusammen, gegen die man nichts machen kann. Die Phase berechnet sich zu
Phase = -arctan(Sinusamplitude/Cosinusamplitude).
Bei 10Hz kommt der Werte -Pi/2 richtig heraus, bei den anderen Frequenzen entsteht irgendein zufälliger Wert durch der Division von 2 Werten der Größenordnung 10^-18.

Dagegen kannst Du nur mit Nachbearbeitung der Arrays etwas machen. Z.B so, daß, wenn ein Amplitudenwert unterhalb einer Schranke liegt, der Phasenwert, weil uninteressant und weil er sich auf etwas quasi nicht Existierendes bezieht, nicht mehr dargestellt wird.

Tip(hat nichts mit dem Problem zu tun):
Als Samplezahl immer 2^n Werte verwenden, sonst kann die schnelle FFT nicht angewendet werden (z.B in Deinem Beipiel fs=1024 Hz und #s=1024).

Gruß Ludwig