声波通信 检测频率时这几句代码是什么意思？

Question

最近在搞一个声波通信的项目，在检测采集到的声波频率的时候，有一段代码看了很久也看不懂，烦请各位大牛指点一下可否？
关键代码如下：

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 是傅立叶变换后的频谱数据指针(int *型) 
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 ）的概率或者权重之类的东西？？？

Answer 1

看上去像是在利用插值计算freq=16000频点的信号功率

Answer 2

你好。我最近也在研究声波通信的代码。跟你的应该是一样的,这一块你弄懂了吗？相互学习啊