從Android獲取超聲波使用頻率

Question

我想要從任何Android設備獲取超聲波，例如頻率在18KHz和19KHz之間的超聲波。

我使用下面的代碼來計算頻率，但它似乎並沒有得到我正確的頻率。我得到的頻率停留在11 KHz和13KHz之間。

private void  calculateFrequency() 
{ 
    // 1 - Initialize audio 
    int channel_config = AudioFormat.CHANNEL_CONFIGURATION_MONO; 
    int format = AudioFormat.ENCODING_PCM_16BIT; 
    int sampleRate = 8000; 
    int bufferSize = 2048; 

    if (bufferSize < AudioRecord.getMinBufferSize(sampleRate, channel_config, format)) 
     bufferSize = AudioRecord.getMinBufferSize(sampleRate, channel_config, format); 
    AudioRecord audioInput = new AudioRecord(AudioSource.MIC, sampleRate, channel_config, format, bufferSize); 

    // 2 - Get sound 
    byte[] audioBuffer = new byte[bufferSize]; 
    audioInput.startRecording(); 
    int nbRead = audioInput.read(audioBuffer, 0, bufferSize); 
    audioInput.stop(); 
    audioInput.release(); 

    // 3 - Transform to double array 
    double[] micBufferData = new double[bufferSize]; 
    final int bytesPerSample = 2; // As it is 16bit PCM 
    final double amplification = 100.0; // choose a number as you like 
    for (int index = 0, floatIndex = 0; index < nbRead - bytesPerSample + 1; index += bytesPerSample, floatIndex++) { 
     double sample = 0; 
     for (int b = 0; b < bytesPerSample; b++) { 
      int v = audioBuffer[index + b]; 
      if (b < bytesPerSample - 1 || bytesPerSample == 1) { 
       v &= 0xFF; 
      } 
      sample += v << (b * 8); 
     } 
     double sample32 = amplification * (sample/32768.0); 
     micBufferData[floatIndex] = sample32; 
    } 

    // 4 - Create complex array 
    Complex[] fftTempArray = new Complex[bufferSize]; 
    for (int i=0; i<bufferSize; i++) 
    { 
     fftTempArray[i] = new Complex(micBufferData[i], 0); 
    } 

    // 5 - Calculate FFT 
    Complex[] fftArray = FFT.fft(fftTempArray); 

    // 6 - Calculate magnitude 
    double[] magnitude = new double[bufferSize/2]; 
    for (int i = 0; i < (bufferSize/2); i++) 
    { 
     magnitude[i] = Math.sqrt(fftArray[i*2].re() * fftArray[i*2].re() + fftArray[i*2].im() * fftArray[i*2].im()); 
    } 

    // 7 - Get maximum magnitude 
    double max_magnitude = -1; 
    for (int i = 0; i < bufferSize/2; i++) 
    { 
     if (magnitude[i] > max_magnitude) 
     { 
      max_magnitude = magnitude[i]; 
     } 
    } 

    // 8 - Calculate frequency 
    int freq = (int)(max_magnitude * sampleRate/bufferSize); 

    ((TextView) findViewById(R.id.textView1)).setText("FREQUENCY = " + freq + "Hz"); 
} 

我使用兩部手機：一個發送超聲波與this app，另一個得到這個超聲波。 我用這個question作爲起始點，我接受了FFT和複雜類。

我的代碼有什麼問題？

Answer 1

爲了得到正確的無鋸齒頻率估計值，必須使用比音頻輸入中的最高頻率高兩倍的採樣率（可能大10％至20％以避免濾波器滾降），因此是您想要找到的最高頻率的兩倍以上。

這是由於採樣定理所需的奈奎斯特率。

所以，如果你想找到一個19千赫的信號，就需要一個接近48000的採樣率。

Answer 2

步驟7和8是不完全正確 - 你需要爲了與最大的幅度爲使用FFT倉的指數確定頻率：

// 7 - Get maximum magnitude 
    double max_magnitude = -1; 
    int max_magnitude_index = -1; 
    for (int i = 0; i < bufferSize/2; i++) 
    { 
     if (magnitude[i] > max_magnitude) 
     { 
      max_magnitude = magnitude[i]; 
      max_magnitude_index = i; 
     } 
    } 

    // 8 - Calculate frequency 
    int freq = (int)(max_magnitude_imdex * sampleRate/bufferSize); 

並經@ hotpaw2，指出您的採樣率在8 kHz時太低 - 至少需要44.1 kHz，最好是48 kHz。