请教：VC下实现apFFT，但幅度校正不准确，高手帮忙看看

fudeyun

我在VC下实现了apFFT算法，但是幅度校正不准确，情况如下：
#define  DATA_PER_T    (512)   //采样点数
i 取 - DATA_PER_T+1 到 2* DATA_PER_T – 1的范围。
Fsig是输入的信号频率。
测试信号：2.356+8.365 * cos(2 * M_PI * Fsig * (i-DATA_PER_T)/DATA_PER_T - M_PI * 1/6 )
所有数据采用float型，采用apFFT/apFFT算法，对相位、频率、幅度进行校正。
测试结果发现，相位和频率校正非常准确，但是幅度校正值较差。实际幅度为8.365，当频点在.9到.4时（比如36.9、36.2、36.4等等），幅度校正结果在8.36左右，当从.4到.8变化时，幅度由8.358105递减到8.335485，最大误差有0.38%。这是什么原因？是否有什么改进的方法让精度高一些？还是我的算法有什么问题？

幅度的校正算法如下：
当p大于0.5时，也就是频点大于.5时，a2取最大点前面的一个，否则取最大点的幅值。
p为相位差，在0-1之间。
a = 2*M_PI*p*(1-p*p)/sin(M_PI*p);
a2 = (a*a)*a2/2;
a2即为校正的幅值。

详细测试结果见附件。初始相角值不正确，实际计算的不是按照下面的相角，另作了处理的。

[ 本帖最后由 fudeyun 于 2009-11-20 17:05 编辑 ]

zhwang554

你的程序对,但p大于0.5时,不在峰值处取值, 取最大点前面的一个,有时影响精度

    可参见whoru 发表于 2009-6-23 16:33 的贴子"求教：apFFT频谱校正问题！'中的程序

http://forum.vibunion.com/thread-83661-1-1.html
     这个程序始终在峰值处取值

      另外apfft程序中的hanning和hanning卷积窗 可改为 hann和hanning卷积窗 一试

[ 本帖最后由 zhwang554 于 2009-11-21 01:07 编辑 ]

fudeyun

多谢赐教，非常感谢！

zhwang554

可参见 振动论坛上  下日志
<apfft谱分析中的振幅校正精度>
zhwang554     2009-12-01 01:40
http://home.vibunion.com/blog-62061-17986.html

[ 本帖最后由 zhwang554 于 2009-12-2 22:23 编辑 ]

wzy75101

王老师就是厉害啊！中国的学术强人！

wzy75101

感谢楼主的提问？为中国的学术做出了贡献！

wzy75101

APFFT是一个很好的算法，向王老师的科研团队致敬！