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

justlikethisis

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

$ O  Q% u/ y# {2 b& E& Yclc
1 H0 w0 Z" p8 j5 zclear;

x = 0:0.0001:6*pi;
( P3 w4 e! p! \8 k  Zy1  = cos(100*pi*x+pi/4);
  O1 p  A# ]4 J6 Z$ Y4 G' E
figure(1)
0 y/ O8 O, L0 Q9 g+ Aplot(x,y1,'k');
: j9 B* Z" x( H. H%hold on;

1 [- Q, M7 g: a%对上述信号进行采样
+ q/ M+ e* ~8 m  Q$ A0 gfs = 2000;
N=128;
dt = 1/fs;
TT = (0:N-1)*dt;
. L' m3 w7 l! e. BX = cos(100*pi*TT+pi/4);
: E: a4 E) N- f5 \, Yfigure(2)
, |& l) t9 p: {stem(TT,X);

; }, H: U9 I) _+ _" q! u%对采样的信号做fft
6 R' |2 g+ o( C6 p3 _2 ZT=1/fs;
9 S+ J* Q: W( o( W# Zt=(0:N-1)*T;
p = X;
Y=fft(p);
magY=abs(Y(1:N/2))*2/N;%真实幅度
3 U) `" ?. k; \0 R. F# ]" zf=(1:N/2)'*fs/N;
0 g7 [/ T6 T+ JangX = angle(Y(1:N/2));
5 d5 P; i! F3 Y$ ~- r$ X  S; J* L, J
' T1 q: a2 o& i" d- Kfigure(3)
" Z$ _# P' q$ @  x* X' ]) f' H! `h=stem(f,magY,'fill','--');
+ _( J8 J% f* I" u& [% [set(h,'MarkerEdgeColor','red','Marker','*')
grid on
; }8 q& L+ R* I9 M7 J1 B/ Mtitle('频谱图');
xlabel('f (Hz)')
4 `1 D  L; Q0 ?$ P8 y6 Mylabel('幅值')

$ j, E1 g, O6 G0 y" Z; U# Ffigure(4)
4 ~2 G3 t; f+ sh=stem(f,angX/pi*180,'fill','--');
! K  Q: ~0 s( x' x2 V, gset(h,'MarkerEdgeColor','green','Marker','*')
grid on
title('相位谱');
xlabel('f (Hz)')
9 |! ^' k* ^( C' |4 @ylabel('相位')

: ^7 ?% ]- P( E$ ?$ x  N; i相位谱

8 G* g/ h+ P2 B8 g' c; f; J频谱
2 @+ K4 _: G: f4 F, \. k
时域
1 E8 H6 _9 c0 B

shapeofyou888

如果只看相对相移，直接在相位谱的该频率处计算的二者之差。
" l: [9 m  q/ `+ r如果要计算各自初始的绝对相位，在有频谱泄露的情况下要还需做矫正。