V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
V2EX 提问指南
reking11
V2EX  ›  问与答

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

  •  
  •   reking11 · 2015-10-20 12:06:18 +08:00 · 1258 次点击
    这是一个创建于 3173 天前的主题,其中的信息可能已经有所发展或是发生改变。
    最近在搞一个声波通信的项目,在检测采集到的声波频率的时候,有一段代码看了很久也看不懂,烦请各位大牛指点一下可否?
    关键代码如下:
    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 )的概率或者权重之类的东西???
    目前尚无回复
    关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   实用小工具   ·   5058 人在线   最高记录 6679   ·     Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 · 38ms · UTC 09:12 · PVG 17:12 · LAX 02:12 · JFK 05:12
    Developed with CodeLauncher
    ♥ Do have faith in what you're doing.