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

justlikethisis

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

+ X' l  O# s" @8 iclc
& M3 q! g& O( [2 Iclear;

2 ^4 ]8 p2 s6 }/ P" mx = 0:0.0001:6*pi;
y1  = cos(100*pi*x+pi/4);

figure(1)
plot(x,y1,'k');
1 ~! O3 [6 I- Q3 n0 C%hold on;

+ E4 u5 Q% h3 j* t1 F, H: E%对上述信号进行采样
1 Y  Q" l" i. i# ~fs = 2000;
N=128;
) X  F1 V4 x' J2 ~1 odt = 1/fs;
TT = (0:N-1)*dt;
, t2 K. B% D: I$ s% @/ l+ w$ nX = cos(100*pi*TT+pi/4);
figure(2)
stem(TT,X);
/ t+ ~- {( |% Y6 J1 A
" K2 K. `# o; m$ c; V  m%对采样的信号做fft
T=1/fs;
t=(0:N-1)*T;
6 ^/ U& e8 s. w3 b, bp = X;
, }  z7 w, n6 Y1 L& W* S; MY=fft(p);
magY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
' d5 P! |! _# A% P# hangX = angle(Y(1:N/2));

* |1 }5 Y. }+ `! _' G" P: W, t& Yfigure(3)
4 M2 D" N; W2 m6 c* X2 lh=stem(f,magY,'fill','--');
set(h,'MarkerEdgeColor','red','Marker','*')
0 s+ h: n2 d4 ]* n4 r5 X# E  Z2 o  B- ?grid on
title('频谱图');
2 k+ w) I" P9 F* o$ ]5 y: g: zxlabel('f (Hz)')
7 d, r3 G  Q6 f2 ?3 S* Fylabel('幅值')
0 w3 |  i) ?& I; A) g) i) E. Z; Y
* p5 N2 r! ?3 m1 E) M4 ofigure(4)
/ h" r: ~2 Y4 Q; W$ g: {7 j8 ^h=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
7 U9 |$ Y6 K. a* ~9 Ugrid on
8 E, |! o. ~$ T) w3 ?6 [title('相位谱');
xlabel('f (Hz)')
+ h/ w: @2 B: {1 ?+ A$ C( yylabel('相位')
1 Z# j3 ?# T6 T2 H
7 i( j. m) n+ h0 w3 T& _% X相位谱

频谱
7 p$ ^* _; [4 k' I: p4 m
5 r' g5 ^, F0 d0 e# |& `时域

shapeofyou888

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