Фазовый спектр синусоидальной волны с Matlab

Question

Я пытаюсь построить фазовый спектр синусоидальной волны с 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]'); 

Answer 1

Фаза частотного компонента плохо определена (и бессмысленна), когда величина настолько мала. Попробуйте случайным образом изменить фазу для каждого спектрального компонента, преобразовать обратно во временную область (с ifft), и вы все равно сможете восстановить синусоиду.