Goertzel算法给出无限结果

Question

我有一个使用 Audacity 软件创建的 20hz - 1 振幅的正弦波。也只有500ms。

我正在使用以下算法来检测频率。

我只想检测音调幅度是否超过阈值并在 20 赫兹频率周期下给出肯定结果。

static float goertzel_mag(int numSamples,int TARGET_FREQUENCY,int SAMPLING_RATE, float* data)
{
    int     k,i;
    float   floatnumSamples;
    float   omega,sine,cosine,coeff,q0,q1,q2,magnitude,real,imag;

    float   scalingFactor = numSamples / 2.0;

    floatnumSamples = (float) numSamples;
    k = (int) (0.5 + ((floatnumSamples * TARGET_FREQUENCY) / SAMPLING_RATE));
    omega = (2.0 * M_PI * k) / floatnumSamples;
    sine = sin(omega);
    cosine = cos(omega);
    coeff = 2.0 * cosine;
    q0=0;
    q1=0;
    q2=0;

    for(i=0; i<numSamples; i++)
    {
        q0 = coeff * q1 - q2 + data[i];
        q2 = q1;
        q1 = q0;
    }

    // calculate the real and imaginary results
    // scaling appropriately
    real = (q1 - q2 * cosine) / scalingFactor;
    imag = (q2 * sine) / scalingFactor;

    magnitude = sqrtf(real*real + imag*imag);
    return magnitude;
}

调用函数

// If there's more packets, read them
        inCompleteAQBuffer->mAudioDataByteSize = numBytes;
        CheckError(AudioQueueEnqueueBuffer(inAQ,
                                           inCompleteAQBuffer,
                                           (sound->packetDescs?nPackets:0),
                                           sound->packetDescs),
                   "couldn't enqueue buffer");
        sound->packetPosition += nPackets;

        NSLog(@"number of packets %i",nPackets);

        float *data=(float*)inCompleteAQBuffer->mAudioData;
        int nn = sizeof(data)/sizeof(float);
        float gort = goertzel_mag(nn, 20, 44100, data);
        NSLog(@"gort:%f", gort);
        if (gort == INFINITY)
             NSLog(@"positive infinity");

函数内部的断点

输出

number of packets 8192
gort:36029896530591744.000000
number of packets 8192
gort:inf
positive infinity
number of packets 5666
gort:inf
positive infinity

为什么我得到的是 inf 结果？我不知道如何读取 return 值，我知道幅度必须始终为正值，但我创建的文件幅度为 1，难道我不应该得到 0 到 1 的结果吗？

编辑； Aduio 信息

afinfo 500ms.aiff
File:           500ms.aiff
File type ID:   AIFF
Num Tracks:     1
----
Data format:     1 ch,  44100 Hz, 'lpcm' (0x0000000E) 16-bit big-endian signed integer
                no channel layout.
estimated duration: 0.500000 sec
audio bytes: 44100
audio packets: 22050
bit rate: 705600 bits per second
packet size upper bound: 2
maximum packet size: 2
audio data file offset: 54
optimized
source bit depth: I16

Answer 1

我认为这些行有一个问题：

float *data=(float*)inCompleteAQBuffer->mAudioData;
int nn = sizeof(data)/sizeof(float);

我相信这是为了告诉你样本的数量。我没有信息或资源来重现您的代码，但可以通过以下方式重现错误：

#include <stdio.h>
#include <stdlib.h>

int main(void) {
    float *data=malloc(sizeof(float) * 10);
    printf ("Sizeof 'data' = %d\n", sizeof(data));
    return 0;
}

程序输出

Sizeof 'data' = 4

在我的 32 位编译中是数组指针的大小，不是数组。使用 sizeof(*data) 不会给你带来任何好处，因为它只会告诉你数据类型 float 的大小，而不是数组。

你无法从它的指针确定数组的大小或元素的数量，所以我的回答是，遗憾的是，你需要更多信息，也许 numBytes？或者 numBytes/sizeof(float)?

Goertzel算法给出无限结果

Goertzel algorithm giving infinite result

c

audio

signal-processing

dft

ios