I have a sine series and calculated the amplitude, power and phase spectra for that. I am not sure about my code and don't know how to read the angle from phase spectra, My program is like this::
close all;clear all;clc; t=0:0.5:100; f=10; fi= 30*pi/180; x=sin(2*pi*f*t+fi); figure plot(t,x) Y=fft(x,128); Amp = (abs(Y)); % Amplitude Pyy= (Y.* conj(Y))/128; % Power f=(0:100)/128; figure; plot(f,Amp(1:101)) figure; plot (f,Pyy(1:101)) Phase= (imag(Y)./real (Y)); Pa=angle(Y)*180/pi; %P2=(P1)*180/pi; plot (Pa')