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

justlikethisis

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

, ?* }5 ?3 C% a5 K" S. Pclc
6 [# g( ~5 J: v+ L* ~% nclear;

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

figure(1)
plot(x,y1,'k');
9 s8 T7 ?- [) p9 E5 ~, D1 p: a, v%hold on;

%对上述信号进行采样
fs = 2000;
/ ]& b* I- n, @- }8 f* uN=128;
) j  i( X) |$ L, r+ D! [1 ndt = 1/fs;
1 A2 v9 y" C; e6 p. G4 m+ pTT = (0:N-1)*dt;
X = cos(100*pi*TT+pi/4);
figure(2)
/ R. L' }* K5 dstem(TT,X);
8 o6 v# u. b. y& G% F  \; c
8 Q3 G7 F, v0 t4 X%对采样的信号做fft
T=1/fs;
t=(0:N-1)*T;
p = X;
- u$ }% }' R: }! Z3 |6 R0 EY=fft(p);
magY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
8 L! j4 j/ }, f* E+ DangX = angle(Y(1:N/2));

" s; f  W1 N+ n3 R' Sfigure(3)
" k& d8 G1 }; r* [% x3 h( o1 Th=stem(f,magY,'fill','--');
set(h,'MarkerEdgeColor','red','Marker','*')
grid on
title('频谱图');
2 U. w9 {0 C4 D' J& F! w, H' j) [9 dxlabel('f (Hz)')
, U8 J' k1 @0 B. u4 ^ylabel('幅值')
/ `1 L) Z4 |9 d% O' G
% h) Y( ^5 T! }8 Afigure(4)
h=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
grid on
title('相位谱');
xlabel('f (Hz)')
ylabel('相位')

相位谱

" V6 y- h0 b" G! ~, ?1 C5 W频谱

* `5 J# V% @5 v; m; x1 V( n时域
9 O: j3 F5 |4 E9 z4 d$ V3 b
  t6 z8 C! x, w; F5 {

shapeofyou888

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