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

justlikethisis

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

clc
4 `6 C$ ~+ Y+ I2 `clear;
8 Q! r$ L* _* ?0 F& w4 \* g
8 J+ a! U2 j3 f: m# s+ Ux = 0:0.0001:6*pi;
' Y5 S# D: P0 e( U4 \y1  = cos(100*pi*x+pi/4);

! z7 E: T# l: j  q5 w' sfigure(1)
. I* l# M/ _" L6 bplot(x,y1,'k');
%hold on;

%对上述信号进行采样
fs = 2000;
N=128;
dt = 1/fs;
$ U" G6 T, _  O' D. \TT = (0:N-1)*dt;
! L1 O$ |1 R9 V0 [7 YX = cos(100*pi*TT+pi/4);
/ q) b& F3 q  Sfigure(2)
stem(TT,X);
! A# w  |4 t% L% u2 V6 }- Y) m
2 b, a6 o; t3 `8 H%对采样的信号做fft
" I) v  O, y6 R9 b! D- LT=1/fs;
t=(0:N-1)*T;
  n) R6 N( j8 J! f* u- I' `p = X;
Y=fft(p);
! ^2 a6 Q  L6 d( ~* {8 PmagY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
/ ?- M2 [/ X( _9 d5 l0 l' `7 E# hangX = angle(Y(1:N/2));
' `' W5 l1 @( z7 d
2 `& o0 ?- {+ ]+ J; dfigure(3)
6 r# X- h+ l& I+ T% B) F2 Q! Yh=stem(f,magY,'fill','--');
/ q. k+ Z7 B( M+ h: q1 uset(h,'MarkerEdgeColor','red','Marker','*')
grid on
title('频谱图');
0 U0 S, r9 ]! ^$ n8 m% C7 c: Hxlabel('f (Hz)')
ylabel('幅值')
$ L, n4 |0 ~; U. O; W) E6 Q
! k+ U0 d, V/ M3 z2 R; o" hfigure(4)
; c# t1 u; s5 k) M6 p" th=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
, b5 a% ^6 |2 t6 kgrid on
title('相位谱');
xlabel('f (Hz)')
ylabel('相位')

相位谱

" h, A6 P1 I& u. I频谱
9 @1 G$ E6 d! H# d/ y1 \5 C) A1 x
) W, k2 J  I+ b) E1 B6 h' y: N, v7 m时域

shapeofyou888

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