搜索
bottom↓
回复: 8

请问用 FFT 测量相位是怎麽做的?

[复制链接]

出0入0汤圆

发表于 2011-11-8 10:49:52 | 显示全部楼层 |阅读模式
请问像 AD5933 这种用傅利叶转换来测量两波形的相位差是怎麽做的?

我只知道会用到高速的 A/D 转换器做采样,但後面怎麽计算就不知道了,

有大大们可介绍一下哪里有作法讨论的网站吗?谢谢。

阿莫论坛20周年了!感谢大家的支持与爱护!!

知道什么是神吗?其实神本来也是人,只不过神做了人做不到的事情 所以才成了神。 (头文字D, 杜汶泽)

出0入0汤圆

 楼主| 发表于 2012-9-6 14:17:32 | 显示全部楼层
很久都没人回覆所以自己来顶一下,希望能够抛砖引玉。

出0入0汤圆

发表于 2012-9-6 18:24:57 来自手机 | 显示全部楼层
查下fft结果的意义,虚部除实部的反正切就是初相

出0入0汤圆

发表于 2012-9-6 18:27:19 来自手机 | 显示全部楼层
单频信号相位精度.001度不是问题

出0入0汤圆

发表于 2012-9-12 17:13:46 | 显示全部楼层
哥们 一个人孤单啊    都没人来。。。。。。。。。。。。。。。。。。。。我也做这个可以一起探讨下qq。925279861

出0入0汤圆

发表于 2012-9-12 19:30:59 | 显示全部楼层
我也小了解一下 呵呵

出0入0汤圆

发表于 2012-9-12 19:31:42 | 显示全部楼层
希望能详细点 谢谢again

出0入0汤圆

发表于 2012-9-12 21:17:56 | 显示全部楼层
理论这个东西还是去翻书吧!贴个实用matlab程序
clear all;clc   
N=1024;
fs=100000;
f=10000
a=1.2345;
ph=50.123;
dc=0.123;
t=0:N-1;
x=a*sin(2*pi*f*t/fs+ph*pi/180)+dc;
x=x+0.0005*randn(1,N);
xfft=fft(x,N);
s1=xfft(2);
S=abs(s1)*2/N;
y=abs(xfft);
[Y1, k]=max(y(2:N/2));
k=k+1;
z1=xfft(k)-(xfft(k-1)+xfft(k+1))/2;
z2=xfft(k+1)-(xfft(k)+xfft(k+2))/2;
ka=abs(z1)/abs(z2);
r=(2-ka)/(1+ka);
fo=(k-1+r)*fs/N
ao=2*pi*r*(1-r*r)*(abs(z1))/(N*sin(r*pi))
pho=(angle(z1)-pi*r)*180/pi+90
---------------------------------------------------
结果精度非常高
fo =

     9.999999905023071e+03


ao =

   1.234515159107427


pho =

  50.120900758689679

出0入0汤圆

发表于 2012-9-12 21:31:52 | 显示全部楼层
我也小了解一下
回帖提示: 反政府言论将被立即封锁ID 在按“提交”前,请自问一下:我这样表达会给举报吗,会给自己惹麻烦吗? 另外:尽量不要使用Mark、顶等没有意义的回复。不得大量使用大字体和彩色字。【本论坛不允许直接上传手机拍摄图片,浪费大家下载带宽和论坛服务器空间,请压缩后(图片小于1兆)才上传。压缩方法可以在微信里面发给自己(不要勾选“原图),然后下载,就能得到压缩后的图片。注意:要连续压缩2次才能满足要求!!】。另外,手机版只能上传图片,要上传附件需要切换到电脑版(不需要使用电脑,手机上切换到电脑版就行,页面底部)。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|Archiver|amobbs.com 阿莫电子技术论坛 ( 粤ICP备2022115958号, 版权所有:东莞阿莫电子贸易商行 创办于2004年 (公安交互式论坛备案:44190002001997 ) )

GMT+8, 2024-7-23 14:22

© Since 2004 www.amobbs.com, 原www.ourdev.cn, 原www.ouravr.com

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