Java .wav 文件频率分析 - 不正确的频率

Question

我正在做一个必须分析的项目 .wav文件，这对我来说基本上归结为能够显示所述文件中存在的频率。
我正在使用 WavFile 类来读取文件，然后我使用 JTransforms 对它们进行 FFT类(我实际上做了一个 realForward 因为我只给它实数)。

一切似乎都正常，直到我将数据输入 Excel:我在第一列中抛出一个行编号(我正在测试的文件为 1-8000)，然后将 FFT 的输出放入下一列。我提供的文件是一个简单的单频声音，频率为 440Hz，持续时间为 1 秒。

但是，在看到图表后，有一个问题:我有一个单频率峰值，这正是我所期望的，但峰值位于 880 位置，这是实际频率的两倍。有人可以向我解释为什么会这样吗？

奖金问题:为什么我会在 e-16 附近获得值?除了峰值之外的一切都应该是0，对吧？ (每次我得到的数据是 <= 1 时，我都会通过向文件写入 0 来解决这个问题 - 请参阅下面的代码)。也许是“噪音”？

代码:

有两个类。第一个，readWav习惯于阅读 .wav文件。第二个，wavFFT , 是实际对数据进行 FFT 的人。
readWav.java的代码:

import WavFile.*;
import java.io.*;

public class readWav {

    // getter methods
    public static long getWavFrames(File file)
    {
        // try loop to catch any exception
        try {
            // open the wav file
            WavFile wavFile = WavFile.openWavFile(file);

            // return the number of frames
            return wavFile.getNumFrames();

        } catch (Exception e) {
            System.err.println(e);

            // error value
            return -1;
        }

    }

    public static int getWavChannels(File file)
    {
        // try loop to catch any exception
        try {
            WavFile wavFile = WavFile.openWavFile(file);

             return wavFile.getNumChannels();

        } catch (Exception e) {
            System.err.println(e);

            // error value
            return -1;
        }
    }

    public static double[] getWavData(File file)
    {
        // try loop to catch any exception
        try {
            // open the file
            WavFile wavFile = WavFile.openWavFile(file);

            // use the getter method to get the channel number (should be mono)
            int numChannels = getWavChannels(file);

            // same, but with the frame getter method
            int numFrames = (int) getWavFrames(file); // possible data loss

            // create a buffer the size of the number of frames
            double[] buffer = new double[numFrames * numChannels];

            // Read frames into buffer
            wavFile.readFrames(buffer, numFrames);

            // Close the wavFile
            wavFile.close();

            return buffer;

        } catch (Exception e) {
            System.err.println(e);

            // throw an error, if this runs something went wrong in reading the .wav file
            throw new RuntimeException("[could not read wav file " + file + "]");
        }
    }


    // main method, solely for testing purposes
    public static void main(String[] args)
    {
        // test, everything seems to be working
        File fichier_son = new File("son/freq1.wav");
        double[] test = getWavData(fichier_son);
        for(int i = 0; i<test.length; i++){
            System.out.println(test[i]);
        }
    }

}

wavFFT.java 的代码:

import org.jtransforms.fft.DoubleFFT_1D;
import java.io.File;
import java.io.PrintWriter;
import java.io.IOException;

public class wavFFT {

    public static double[] realFFT(File file)
    {
        // Get the .wav data using the readWav class
        double[] data_to_fft = readWav.getWavData(file);

        /* Get the length of the array.
        Since we are feeding real numbers into the fft,
        the length of the array should be equal to the
        number of frames, which we get using the readWav class. */
        int n = (int) readWav.getWavFrames(file);

        // Make a new fft object
        DoubleFFT_1D fft = new DoubleFFT_1D(n);

        // Perform the realForward fft
        fft.realForward(data_to_fft);

        // Return the final data
        return data_to_fft;
    }


    public static void writeToFile(File in, File out) throws IOException
    {
        PrintWriter print_out = new PrintWriter(out);
        int i;
        double[] data_to_file = realFFT(in);

        for(i=0; i<data_to_file.length; i++){
            if(data_to_file[i] > 1){
                print_out.println(data_to_file[i]);
            } else {
                print_out.println(0);
            }

        }
        print_out.close();
    }

    // main method, solely for testing purposes
    public static void main(String[] args) {
        File fichier_son = new File("son/freq1.wav");
        double[] test = realFFT(fichier_son); 
        int i;

        for(i=0; i<test.length; i++){
            System.out.println(test[i]);
        }


        try{
            writeToFile(fichier_son, new File("datafiles/output.txt"));
        } catch (IOException e){
            System.out.println("error");
        }
    }

}

Answer 1

您没有说明如何解释使用上述代码生成的 Excel 中的结果。然而，一个可能的错误是误解了 FFT fft.realForward() 的输出。 - 这是一个复数数组，其实部和虚部占据连续的元素，如文档所示 here .如果您只是使用出现峰值的数组的索引，您的结果将相差两倍。请注意，此 FFT 实现仅计算 Nyqvist 速率(除此之外仅产生一个“别名”)。

其他注意事项:

您正在申请 rectangular window function到 sample 。坦率地说，我很惊讶其他垃圾箱中的出血小到 10e-16 .您的应用程序有更好的选择。

该窗口似乎是文件的全长，可能会导致巨大的 FFT。一般来说，2 次方 FFTS 的效率要高得多。您通常会看到在固定长度窗口上执行的频率分析。

FFT 实际上是一系列带通滤波器，它们是输入样本的加权和。对您看到的非零值的贡献只是浮点舍入误差。您将整个文件放在窗口中这一事实意味着有很多操作会对此有所贡献。