从LomontFFT计算基频

本文关键字：计算 LomontFFT | 更新日期: 2023-09-27 17:53:16

我使用来自http://www.lomont.org/Software/Misc/FFT/LomontFFT.html的LomontFFT从信号的采样值中获得基频。为了测试基频是否正确确定，我使用了过去的一些样本(已知基频)。

下面是我编写的代码，用于调用LomontFFT算法并确定FFT:

private void buttonFFT_Click(object sender, EventArgs e)
    {
        //double fftavg = 0;
        double fftmax = 0;
        var fftData = new byte[512];
        double[] fftValues = Enumerable.Repeat(0.0, 512).ToArray();
                           Array.Copy(sapmledDoubleValuesADC1, fftValues, sapmledDoubleValuesADC1.Length);
var fftMethod = new Lomont.LomontFFT();
        fftMethod.RealFFT(fftValues, true);
        for (int i = 0; i < 512; i += 2)
        {
            double fftmag = Math.Sqrt((fftValues[i] * fftValues[i]) + (fftValues[i + 1] * fftValues[i + 1]));
            if (fftmag > fftmax)
                fftmax = fftmag;
            //fftavg += fftmag;
            //fftData[i] = (byte)fftmag;
            //fftData[i + 1] = fftData[i];
        }
        textBoxFundaFreq.Text = "Frey = " + fftmax.ToString();
        for (int x = 1; x < 512; x++)
        {
            this.chart2.Series[0].Points.AddXY(x, fftValues[x]);
        }
    }

但问题是频率的大小是错误的。FFT也不匹配，但这是可能的，因为有多个解决方案，但频率应该是相同的。这个算法已经被证明了很多年，所以它绝对没有错。我在调用代码时做错了什么吗?

(我在采样数据中只有实值)

从LomontFFT计算基频

基频提取/基音检测是不一种简单的算法。对于大多数输入信号(除了单个正弦/余弦波之外的任何信号)，FFT将显示几个峰值，通常最好估计这些峰值之间的距离。

进一步，您需要插值FFT箱以获得单个峰值的准确结果。

对于大多数应用程序，最好还是计算自相关函数(ACF)。

相关文章：