找回密码
 注册
关于网站域名变更的通知
查看: 558|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-8-6 13:36 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
如下,对一列没有初相的余弦信号和初相为pi/4的余弦信号,分别做傅里叶变换,怎么样从频谱图中看出相位信息,信号频率是50hz,采样频率是2000matlab代码如下:
) Y  j& j* M# W& B; C- n% G6 {; [) Q+ P  k2 L( Q) N
clc
8 X$ M6 I2 x3 H$ N! F. l% G) V  ?clear;8 m# Q, o6 n0 D" E

0 n* Y) C' E  v6 D3 a+ ?0 w* fx = 0:0.0001:6*pi;9 {  N9 O  A% ~, k; c9 _
y1  = cos(100*pi*x+pi/4);
0 h' ~5 {% e; Y4 i* F! I; K2 l, O; [3 P: \! z7 L: E  J1 o3 a; n
figure(1), X$ |% U8 n4 G6 _- r5 O& a# v
plot(x,y1,'k');
2 ~2 K& a! c- ~1 n* r, @) p%hold on;
' m+ ?3 ~  z+ M( {! W( b0 i- Q
, e' P( b' X# w7 m/ s  l3 v2 H%对上述信号进行采样  d4 L2 u4 ^' y0 B. k8 k; p) I
fs = 2000;! i  E( H" g3 X: y
N=128;7 c2 q* x1 U* @
dt = 1/fs;
9 W2 u% A) i1 STT = (0:N-1)*dt;" H8 _# f% f5 c6 d; C6 T/ ?) X
X = cos(100*pi*TT+pi/4);8 H; e/ @4 }1 m% ?7 q+ ~6 q
figure(2)3 Z/ T* l8 K) ^7 @7 y0 ?
stem(TT,X);
; Y4 b' {& f) A' ~3 c3 e8 I- K* Z' P' g& I+ B: \
%对采样的信号做fft
0 N: K7 [8 G. W' v1 t. z: c& o% hT=1/fs;
, U; R+ ~, o& h, L. |t=(0:N-1)*T;  `( d& h6 }. a9 M  H* Z
p = X;6 I! y( Z" Q! n0 L$ J
Y=fft(p);
  a! B" O" h. E3 K+ `3 s- _- E: nmagY=abs(Y(1:N/2))*2/N;%真实幅度  `8 `+ M) F7 |' [0 H" f) Y; ^
f=(1:N/2)'*fs/N;
6 o+ v; [; l) P8 H1 N; vangX = angle(Y(1:N/2));
" I! p8 s$ l7 a, m$ h
, |  V, Y  c5 v" H. H4 bfigure(3)
& O# |) C4 x& c, m5 k8 y8 q! u& w& s  kh=stem(f,magY,'fill','--');
6 `' U! Z" {( zset(h,'MarkerEdgeColor','red','Marker','*')
! p* B0 @5 e. ^4 J1 o6 [grid on' b; f* U/ a4 A& p' l! H
title('频谱图');1 J  i+ ?, _5 A$ f: ~& f
xlabel('f (Hz)')
1 q( D$ o6 _# {+ @4 X2 ~7 G) \/ Dylabel('幅值')
0 C- U% M1 B/ h( C. K7 ]! y) O/ H* R
figure(4)# B! e6 ]+ J1 J& u# D2 M, m* k
h=stem(f,angX/pi*180,'fill','--');
' L3 x3 U5 G' }  n8 sset(h,'MarkerEdgeColor','green','Marker','*')
1 j5 _7 ~$ H  Q) egrid on
* }' |4 ^- a; I9 h2 y! [title('相位谱');0 ~& `/ W4 r; Y2 B$ c! |% S3 k) c
xlabel('f (Hz)')
9 l. ?  i" a& Q7 ?, {ylabel('相位')

, H( L2 b9 }0 v* R: q! ~. @; U
, E7 X4 S% j0 @7 }
相位谱
0 J5 ~4 ]4 s0 _. r8 V. \ ( z' Z6 v/ K# }4 H7 d0 M5 J0 @  ]; ?
频谱
* f( [" E: ~0 s  F6 r
# Y8 w2 c% \. G$ m7 `- X+ F时域
! f( _2 X- _+ K2 A# @3 e
  C3 e2 G' F  z: M6 V8 k3 X
  • TA的每日心情
    开心
    2022-12-5 15:37
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2020-8-6 14:18 | 只看该作者
    如果只看相对相移,直接在相位谱的该频率处计算的二者之差。
    ) w; S/ X7 K7 @* S/ a如果要计算各自初始的绝对相位,在有频谱泄露的情况下要还需做矫正。
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-24 23:13 , Processed in 0.187500 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表