声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: a.gain

[FFT] 请问王老师关于apFFT问题

[复制链接]
发表于 2009-11-19 01:08 | 显示全部楼层

回复14楼kleeblatt 的贴子

你的问题
"我编写的一个函数为
x1(i)=2*cos(2.2*pi*2*(i-127)/128+10*pi/180),如果改成
x2(i)=2*cos(2.2*pi*2*i/128+10*pi/180)  相位结果就不对了,为什么呢?"

因这二个信号相位差 2.2*p*2*1/127*180/pi 度, 相位结果当然不对了

看初相位更清楚, 对笫二个信号 x2(i)=2*cos(2.2*pi*2*i/128+10*pi/180) 其初相位为10度

笫一个信号  x1(i)=2*cos(2.2*pi*2*(i-127)/128+10*pi/180),它的初相位不只是10度, 应是(即x1(i)式中i=0时的相位)mod(-127/128*2*pi*2.2*180/pi+10, 360)=304.1875度

所以x1(i)=2*cos(2.2*pi*2*(i-127)/128+10*pi/180)=2*cos(2.2*pi*2*i/128+304.1875*pi/180)

它和x2(i) 是同频2.2,同振幅2,不同初相(10度和304.1880度)的余弦信号.在任何时刻二者相位差304.1880-10=294.1875度


   你问: "每个包含x(0) 的组合移位 ——求和——DFT的结果 与 每个包含x(0)的组合移位——每个组合DFT——每组FFT结果对应求和得出的结果不同"

    每个包含x(0) 的组合移位相加的dft的 应等於每个包含x(0)的组合移位的DCT的复数之和, 这是DFT的线性性质, 你用小N 例如N=6为例,仔细分析, 不会错的. 从中还可看到泄漏互相抵消,相位互相补偿的细节.
  <
数字信号全相位谱分析与滤波技术>一书中第
3.4.1 数据实测分析(p37-38) N=4为例分析了整个过程.

[ 本帖最后由 zhwang554 于 2009-11-19 08:43 编辑 ]
回复 支持 反对
分享到:

使用道具 举报

发表于 2010-1-4 16:39 | 显示全部楼层

回复 7楼 xiaokang 的帖子

你好    我也想对实测的一些数据进行频率估计  不知道你用的什么方法?  还是就是我想用apfft时移相位差法   但是不知道最后程序应该怎么改    请多指教   谢谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 15:27 , Processed in 0.056497 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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