C语言使用fft库函数求wav文件的时域信号功率谱

使用C语言中的fft库函数可以方便地求出WAV文件的时域信号功率谱。以下是一个简单的示例代码：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "fft.h" // 包含FFT库函数头文件

#define N 1024 // FFT窗口大小

int main()
{
    FILE *fp;
    short buf[N];
    double fft_real[N], fft_imag[N], power[N/2+1];
    int i;

    // 读取WAV文件
    fp = fopen("test.wav", "rb");
    if (fp == NULL) {
        printf("Cannot open file!\n");
        exit(1);
    }

    // 跳过WAV文件头
    fseek(fp, 44, SEEK_SET);

    // 循环读取WAV文件中的数据，并进行FFT计算
    while (fread(buf, sizeof(short), N, fp) == N) {
        // 将时域信号复制到FFT缓冲区
        for (i = 0; i < N; i++) {
            fft_real[i] = (double)buf[i];
            fft_imag[i] = 0.0;
        }

        // 进行FFT计算
        fft(fft_real, fft_imag, N);

        // 计算功率谱
        for (i = 0; i <= N/2; i++) {
            power[i] = fft_real[i]*fft_real[i] + fft_imag[i]*fft_imag[i];
        }

        // 输出功率谱
        for (i = 0; i <= N/2; i++) {
            printf("%d %lf\n", i, power[i]);
        }
    }

    fclose(fp);

    return 0;
}

需要注意的是，上述代码中的fft.h头文件中包含了FFT库函数的定义，需要根据实际情况进行修改。