EDA365电子论坛网

标题: 求助，傅里叶相位谱怎么看，如何从相位谱中得到相位信？ [打印本页]
作者: justlikethisis    时间: 2020-8-6 13:36
标题: 求助，傅里叶相位谱怎么看，如何从相位谱中得到相位信？
如下，对一列没有初相的余弦信号和初相为pi/4的余弦信号，分别做傅里叶变换，怎么样从频谱图中看出相位信息，信号频率是50hz，采样频率是2000matlab代码如下：4 v: O; `+ q8 x- Q6 {4 h

% k' l! V' Q- ]/ ^- n! vclc4 R. T9 |+ T% K0 D* o- ?  O1 y$ l9 ^7 r
clear;1 d( u8 \5 s! A9 x% I* N
6 u6 w( ^* \  [$ N* @
x = 0:0.0001:6*pi;
% s9 N& `$ e( }* k7 I) My1  = cos(100*pi*x+pi/4);' N( q, W7 G. u4 y9 z8 a8 A
8 _: [% c' k7 e0 y; w6 k
figure(1)1 i7 x) A1 g0 W& M
plot(x,y1,'k');
* y! Y6 C* u- P2 v" u' ?1 X%hold on;3 ?  W4 Y0 n4 R

- T2 E: r; y! W6 _% e, h2 v%对上述信号进行采样  v5 Q3 d8 f6 e  @. t5 X7 [2 {
fs = 2000;
& _) ]; I4 Q6 EN=128;& M' [; M- X; ]# K6 V' ?0 R( i
dt = 1/fs;0 f, J2 `& o6 W  |9 j$ b7 v
TT = (0:N-1)*dt;; G8 Q# I5 C  X' e1 ~
X = cos(100*pi*TT+pi/4);
0 j/ c4 x5 y' D  c0 H) O' afigure(2)
5 f9 x# Y( z' @2 D4 _; n! p& j) sstem(TT,X);
7 S8 `4 B# l4 ]( f* {. Z% O8 d* s1 B$ [7 X4 W
%对采样的信号做fft! u. k8 h! g3 d# h# {
T=1/fs;$ e1 N. Y" K4 U1 @$ `. [# D3 O* W% A
t=(0:N-1)*T;
: \, {) t. a( U# o( Wp = X;
  z# Q+ t- G& oY=fft(p);
* E  P8 {& K- BmagY=abs(Y(1:N/2))*2/N;%真实幅度9 [' U% s0 N& R$ o
f=(1:N/2)'*fs/N;7 C: r  ?/ a* e7 Q* u/ ]1 B/ O
angX = angle(Y(1:N/2));6 f1 a( Z* y( [) E6 N
4 N% r/ N! t4 I0 U, \
figure(3)# `) f% U$ ?! ]
h=stem(f,magY,'fill','--');( T. G. o9 a, p6 N, b1 z' N
set(h,'MarkerEdgeColor','red','Marker','*')3 m" b+ ]2 C+ X- F' ]$ x
grid on- ~; \) u6 u) ]- L  K. K( q
title('频谱图');
+ V' G! A: O+ N  [& Q, C/ exlabel('f (Hz)')
. a- y5 i% [* Zylabel('幅值')" u" M9 d/ T' k0 Q
' x: R- F( ^3 p
figure(4)# o+ q2 X& ^) e/ s7 B7 p% X
h=stem(f,angX/pi*180,'fill','--');/ Y* Q3 T0 C9 e' D* q/ g/ @' c
set(h,'MarkerEdgeColor','green','Marker','*')# J5 f6 n7 {* ]- E; Z5 a
grid on/ D$ g5 ~! O" t7 f2 T0 c, M8 b
title('相位谱');/ ^# Y  _7 p/ Y/ a6 g$ |
xlabel('f (Hz)')
' v/ m4 q" i# B5 T! w) e9 ^ylabel('相位')
5 [7 Q# }! H9 j' j- n8 h  C5 H
" s. h/ a! k7 e6 ]/ V2 B9 h相位谱! k3 t$ [; w0 D5 W% Z' W

! y* `: o$ J- |频谱
- Y3 d7 R* A+ k; {# O# F " q$ T  W# w7 j( ^' f6 F
时域
% `$ j  ~" N- N3 B% Q 9 f( C' n8 Q$ K+ P1 u& F
作者: shapeofyou888    时间: 2020-8-6 14:18
如果只看相对相移，直接在相位谱的该频率处计算的二者之差。6 v( H: a7 F  W5 p/ \& I) G
如果要计算各自初始的绝对相位，在有频谱泄露的情况下要还需做矫正。



欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2