求助，傅里叶相位谱怎么看，如何从相位谱中得到相位信？

justlikethisis

如下，对一列没有初相的余弦信号和初相为pi/4的余弦信号，分别做傅里叶变换，怎么样从频谱图中看出相位信息，信号频率是50hz，采样频率是2000matlab代码如下：

clc
' U9 z6 A) U* D8 nclear;

x = 0:0.0001:6*pi;
y1  = cos(100*pi*x+pi/4);
/ o, O  {6 {( v6 r0 i% M0 U
7 y) N; C+ Q2 S) `8 Y+ e( Jfigure(1)
plot(x,y1,'k');
%hold on;
/ ^2 Y5 A- B7 x9 _/ ~3 O
%对上述信号进行采样
, \$ F7 s% O$ c# E( c( Q7 a3 Cfs = 2000;
N=128;
7 V- \1 S2 ?$ q+ H3 @dt = 1/fs;
3 s+ g4 P; I  }# aTT = (0:N-1)*dt;
1 D+ \" F9 c1 i% I. `' S% sX = cos(100*pi*TT+pi/4);
figure(2)
stem(TT,X);
: B9 S" a& O+ P0 m1 U
7 U0 v3 }9 v* Z+ U7 K  W) D& E%对采样的信号做fft
T=1/fs;
t=(0:N-1)*T;
p = X;
, Z6 M+ j: g$ M3 E) L, Q( Q, hY=fft(p);
magY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
angX = angle(Y(1:N/2));

7 K' d  Y! s4 o- n. vfigure(3)
h=stem(f,magY,'fill','--');
6 M% i* N) m" |2 Q( f; |2 Hset(h,'MarkerEdgeColor','red','Marker','*')
grid on
2 K- f# R9 p" g( G/ Y& otitle('频谱图');
+ r3 w2 \& C: d8 `0 h+ P  P0 {4 Pxlabel('f (Hz)')
ylabel('幅值')

figure(4)
9 `4 b) C/ R* ^+ I( hh=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
7 y$ z2 n' |8 j4 `grid on
title('相位谱');
2 p" w- R! Z9 H9 J- Bxlabel('f (Hz)')
  K  C/ e- }, C" mylabel('相位')
2 O; k- P  f7 I
相位谱

频谱

时域

shapeofyou888

如果只看相对相移，直接在相位谱的该频率处计算的二者之差。
如果要计算各自初始的绝对相位，在有频谱泄露的情况下要还需做矫正。