声波通信 检测频率的一个问题

2015-10-20 12:06:18 +08:00
 reking11
最近在搞一个声波通信的项目,在检测采集到的声波频率的时候,有一段代码看了很久也看不懂,烦请各位大牛指点一下可否?
关键代码如下:
unsigned int freq = 16000; // 频率 16000hz
int fftIdx;
int NyquistMaxFreq = sampleRate/2.0; // 奈奎斯特 采样定律
fftIdx = freq / NyquistMaxFreq * fftLength; // 此处是计算 fft 后, 在幅度-频率坐标系当中, 频率 freq 对应的 x 轴坐标值???

double fftIdx_i, fftIdx_f;
fftIdx_f = modf(fftIdx, &fftIdx_i);
// 将 fftIdx 的整数部分放到 fftIdx_i, 小数部分放到 fftIdx_f;

signed char fft_l, fft_r;
float fft_l_fl, fft_r_fl;
float interpVal;

// fftData 是傅立叶变换后的频谱数据
fft_l = (fftData[(int)fftIdx_i] & 0xFF000000) >> 24;
fft_r = (fftData[(int)fftIdx_i + 1] & 0xFF000000) >> 24;
// 不知道这两句是想干嘛???

fft_l_fl = (float)(fft_l + 80) / 64.;
fft_r_fl = (float)(fft_r + 80) / 64.;
interpVal = fft_l_fl * (1. - fftIdx_f) + fft_r_fl * fftIdx_f;
// 还有这三句

#define CLAMP(min,x,max) (x < min ? min : (x > max ? max : x))
interpVal = sqrt(CLAMP(0., interpVal, 1.));
最后得到的参数 interpVal ,我猜测貌似是表示频率( freq = 16000 )的概率或者权重之类的东西???
1258 次点击
所在节点    问与答
0 条回复

这是一个专为移动设备优化的页面(即为了让你能够在 Google 搜索结果里秒开这个页面),如果你希望参与 V2EX 社区的讨论,你可以继续到 V2EX 上打开本讨论主题的完整版本。

https://www.v2ex.com/t/229476

V2EX 是创意工作者们的社区,是一个分享自己正在做的有趣事物、交流想法,可以遇见新朋友甚至新机会的地方。

V2EX is a community of developers, designers and creative people.

© 2021 V2EX