Matlabで正弦波の位相スペクトルをプロットしようとしています。以下にコードを添付します。評判がないため、図を添付できません。マグニチュードスペクトルは良さそうですが、位相スペクトルは正しくないようで、まるでノイズのようです。理由はわかりますか？
clear all;
fs=8000;
l=1000;
t=1/fs*(1:l);
x1=sin(2*pi()*1000*t);
spec_x1=fft(x1,1000);
magnitude=2*abs(spec_x1)/l;
phase=angle(spec_x1)*180/pi;
figure
plot(fs/2*linspace(0,1,500),magnitude(1:500));
title('Magnitude spectrum');
xlabel('F[Hz]');
ylabel('Magnitude');
figure
plot(fs/2*linspace(0,1,500),phase(1:500));
title('Phase spectrum');
xlabel('F[Hz]');
ylabel('Phase [degrees]');