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

justlikethisis

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

/ z  P- M+ ]+ E& t5 [2 `- Mclc
clear;

7 Q- b; |8 |8 C: k: S8 A9 {" u5 U- lx = 0:0.0001:6*pi;
y1  = cos(100*pi*x+pi/4);
* d' I4 o# g0 w9 X
figure(1)
plot(x,y1,'k');
  s% u4 s5 z6 Q4 l%hold on;

%对上述信号进行采样
fs = 2000;
N=128;
dt = 1/fs;
: P5 }% m4 w# A  M9 l3 C# t" TTT = (0:N-1)*dt;
X = cos(100*pi*TT+pi/4);
figure(2)
stem(TT,X);
# q, s9 u4 k5 m% {" t) b
%对采样的信号做fft
. Y. N6 ^3 D  qT=1/fs;
3 |# r4 H! q# bt=(0:N-1)*T;
p = X;
Y=fft(p);
8 j7 Y# D8 W$ r! K& emagY=abs(Y(1:N/2))*2/N;%真实幅度
8 E- k* r% j' \- U8 {% mf=(1:N/2)'*fs/N;
, J: N) n/ y: l% X. _: p8 ^' XangX = angle(Y(1:N/2));
+ K& S) t3 F9 w+ {
figure(3)
h=stem(f,magY,'fill','--');
set(h,'MarkerEdgeColor','red','Marker','*')
grid on
title('频谱图');
xlabel('f (Hz)')
ylabel('幅值')
, ?, R/ k& l& B- o1 [* U
figure(4)
h=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
. k& I( ~, }# K4 mgrid on
4 q1 B1 x; B* m+ p2 I/ Htitle('相位谱');
xlabel('f (Hz)')
" v. R% C5 `9 ~7 ]* R1 bylabel('相位')

9 @7 s/ }. y- O8 b相位谱

频谱

2 k2 k3 B8 n8 J+ `时域
: F. ]6 _$ ]# G
( C0 v! t3 N( z- c) }/ J1 |; W9 g

shapeofyou888

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