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

justlikethisis

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

" M! H4 G$ `# V" _% i% Z4 Zclc
' O3 r! |( T. {) `+ zclear;

x = 0:0.0001:6*pi;
y1  = cos(100*pi*x+pi/4);

figure(1)
plot(x,y1,'k');
3 c" x, R# g! G, L%hold on;
9 N5 m/ @2 m! |* y; o7 ~& K
8 a$ l: q% {% a7 c%对上述信号进行采样
: T* ^/ K9 K3 L; d& i, kfs = 2000;
N=128;
dt = 1/fs;
TT = (0:N-1)*dt;
X = cos(100*pi*TT+pi/4);
& P* V2 F* W" bfigure(2)
9 k' P' }" g/ s, r* }4 n+ J/ Lstem(TT,X);

%对采样的信号做fft
T=1/fs;
" W8 v) r8 E6 P' W7 Ot=(0:N-1)*T;
9 e& {# D) y2 v1 J+ ]- xp = X;
Y=fft(p);
magY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
angX = angle(Y(1:N/2));

3 a5 S; c+ X3 i: ^( l) z0 u4 Rfigure(3)
h=stem(f,magY,'fill','--');
set(h,'MarkerEdgeColor','red','Marker','*')
) k' D8 R0 X" `3 P5 P) Kgrid on
title('频谱图');
) l  [# o7 B3 _0 j: Y$ z) T6 D! Nxlabel('f (Hz)')
2 S/ _+ G- l# Z" C$ q+ Kylabel('幅值')
( _$ E. O! K  U4 p/ u) T. h- p
+ T4 f5 I# O2 }7 l( T9 K; y9 hfigure(4)
h=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
6 F. h: A) n( h, k& Y5 vgrid on
" V2 G( E2 Y8 k9 G; ]2 etitle('相位谱');
; O/ C$ k8 H- G% J; V3 K; W. b* {, zxlabel('f (Hz)')
ylabel('相位')
9 |& K7 K! ~9 s; i+ @- X
相位谱
  O+ |% h! p  k( a6 T
3 \/ B* h, f8 a; N& B5 s频谱
, o/ N+ ]( O( U2 e, y4 w
时域
" `; h7 I' t) y/ E1 z" V
4 Z6 v5 a: X1 _" g/ n

shapeofyou888

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