如何从FFT幅值数组获取对应频率及解决索引始终为0问题

一、幅值索引转对应频率 #

FFT结果的每个索引对应的频率可以用公式直接计算：

频率 = (索引值 × 采样率) / FFT长度

其中：

• 采样率：你采集麦克风音频时使用的采样频率（比如常见的44100Hz、48000Hz）
• FFT长度：你的real和imaginary数组的原始长度（因为你的toMagnitudes只取了前一半数据，所以索引范围是0到FFT长度/2 - 1，对应频率范围是0到采样率/2，也就是奈奎斯特频率）

举个代码示例，假设采样率是44100Hz，FFT长度是real.length，计算索引i对应的频率：

int sampleRate = 44100;
int fftLength = real.length;
double frequency = (i * sampleRate) / (double)fftLength;

二、最大幅值索引始终为0的原因及解决 #

索引0对应的是直流分量（0Hz），出现这种情况大概率是音频信号带有直流偏移，或者采集到的有效音频信号太弱。解决方法如下：

1. 去除音频信号的直流偏移 #

在执行FFT之前，先计算原始音频样本的平均值，然后每个样本减去这个平均值，去掉直流分量：

// 假设rawAudio是你采集到的原始音频float数组
float[] rawAudio = ...;
float dcOffset = 0;
for (float sample : rawAudio) {
    dcOffset += sample;
}
dcOffset /= rawAudio.length;

// 去除直流偏移
for (int i = 0; i < rawAudio.length; i++) {
    rawAudio[i] -= dcOffset;
}

// 再对处理后的rawAudio执行FFT得到real和imaginary数组

2. 检查麦克风采集状态 #

• 确认麦克风权限已开启，能正常采集声音
• 测试时靠近麦克风输入明显音频（比如说话、播放固定频率声音），避免采集静态噪声

3. 可选：添加FFT窗口函数 #

如果去除直流偏移后仍有问题，可以给音频数据加汉明窗，减少频谱泄漏的影响：

// 给rawAudio加汉明窗
for (int i = 0; i < rawAudio.length; i++) {
    float hamming = (float)(0.54 - 0.46 * Math.cos(2 * Math.PI * i / (rawAudio.length - 1)));
    rawAudio[i] *= hamming;
}

你的现有代码整理 #

主代码片段 #

final double[] magnitudes = toMagnitudes(real, imaginary);
System.out.println("Largest index: " + maxIndex(magnitudes));

辅助函数 #

public static int maxIndex(double[] array) {
    int maxIndex = 0;
    for (int i = 0; i < array.length; i++) {
        if (array[i] > array[maxIndex]) {
            maxIndex = i;
        }
    }
    return maxIndex;
}

private static double[] toMagnitudes(final float[] realPart, final float[] imaginaryPart) {
    final double[] powers = new double[realPart.length / 2];
    for (int i = 0; i < powers.length; i++) {
            powers[i] = Math.sqrt(realPart[i] * realPart[i] + imaginaryPart[i] * imaginaryPart[i]);
    }
    return powers;
}

内容的提问来源于stack exchange，提问作者SprintKeyz