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

justlikethisis

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

, Z/ q3 [8 `: ]8 Z) Z5 lclc
1 y! D8 ~' X% T' b" ?0 s, |clear;

0 M' j- t' [+ _: {7 o& }- ?x = 0:0.0001:6*pi;
y1  = cos(100*pi*x+pi/4);

; M5 y5 t' Q' gfigure(1)
) `2 f" ~/ a; \; t/ R; |  r# nplot(x,y1,'k');
%hold on;
* N7 G7 g/ G" U- _) u
%对上述信号进行采样
% n/ e1 l2 N& r+ y" k# _( K  sfs = 2000;
N=128;
6 f0 i; H! i' V1 Gdt = 1/fs;
* M& M& s/ n: [& N/ R* a/ Y7 hTT = (0:N-1)*dt;
2 l7 M- t# U4 k+ TX = cos(100*pi*TT+pi/4);
figure(2)
stem(TT,X);

%对采样的信号做fft
T=1/fs;
t=(0:N-1)*T;
+ h: Z/ z  S+ Cp = X;
Y=fft(p);
. e5 O/ s( t/ GmagY=abs(Y(1:N/2))*2/N;%真实幅度
f=(1:N/2)'*fs/N;
" p  {' m! t& mangX = angle(Y(1:N/2));

8 O7 w4 `! ?4 hfigure(3)
! N7 D4 ?1 S8 f& ih=stem(f,magY,'fill','--');
set(h,'MarkerEdgeColor','red','Marker','*')
/ D7 |* B4 ^* B6 f0 M6 zgrid on
+ O% n4 V! p4 A6 C8 a+ {  Y8 Ytitle('频谱图');
6 L  b: z: q4 G% A: kxlabel('f (Hz)')
0 k0 M2 [% Z& Vylabel('幅值')
4 H) r1 c! I" r
figure(4)
h=stem(f,angX/pi*180,'fill','--');
set(h,'MarkerEdgeColor','green','Marker','*')
' ?# n0 Z8 w# [6 p7 Ggrid on
title('相位谱');
8 y$ Z+ m, b4 N( r! Txlabel('f (Hz)')
ylabel('相位')

* m6 ?3 z1 G' l- }$ I( T( \相位谱

# n2 x# S: X+ f0 _! k频谱
- f* h, d: g5 e
时域
; m/ G2 z8 d/ K' C
! r! k# [" l- l2 h. I  |

shapeofyou888

如果只看相对相移，直接在相位谱的该频率处计算的二者之差。
- h& K$ Z* {( Q3 X4 \' C, z% h如果要计算各自初始的绝对相位，在有频谱泄露的情况下要还需做矫正。