声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2050|回复: 5

[FFT] 求高手帮忙关于频谱质心实现的问题

[复制链接]
发表于 2008-12-7 16:49 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我想在matlab里面算出频谱质心,公式见附件1:
我在用matlab实现的时候不知道K取多大,我觉得FFT是关于中点对称的,所以,是否K该取FFT(N)的1/2N,还是取整个N????





[ 本帖最后由 weidaidi 于 2008-12-7 16:55 编辑 ]
公式.JPG
回复
分享到:

使用道具 举报

发表于 2008-12-7 21:20 | 显示全部楼层
对了...仔细看了一下你給的公式,仿佛 SC(n) 就恒= K ? 还是我看花眼了?
分母上的K可以提出来,然后 |dft(n,k)| 部分就约掉了?

分母应该是  K.|dft(n,k)|    or    k.|dft(n,k)| ?

fft在实数输入的时候,结果是对称的.但是能量是分布在N点上的,
从频谱质心的定义上来看,既然是能量分布的平均点,
应该K = N.
 楼主| 发表于 2008-12-8 21:32 | 显示全部楼层
非常感谢您的解答,公式写错了,应该是 k.|dft(n,k)|。
发表于 2008-12-9 08:42 | 显示全部楼层
我参看一下吉小军等 “AR谱结构特征提取及其对车辆目标的识别效果分析”(探测与控制学报 1999 21(2))中提到了谱重心(和楼主的谱质心是一个概念),计算时仅积分(累加)到fmax,所以楼主计算时把K取作N/2。
发表于 2008-12-9 09:28 | 显示全部楼层
之前没接触过谱质心,
刚才去网上找资料看了一下,
从语音处理的谱质心的物理意义上来看, songzy41说的是对的.
N/2后面的镜象频率不用计入.
K = N/2就可以了.
 楼主| 发表于 2008-12-9 15:30 | 显示全部楼层
谢谢大家了,这下知道怎么做了:victory:
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-19 01:01 , Processed in 0.087338 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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