4

第一次问虽然我已经访问了一段时间。

这是问题所在:

我目前正在尝试使用以下属性隔离 WAVE 数据文件中包含的信号的基频:

  • PCM 音频格式,即线性量化
  • 8000 赫兹采样率
  • 每个样本 16 位
  • 16000 字节速率

只有一个通道,没有交错。

获取字节值:

System.IO.FileStream WaveFile = System.IO.File.OpenRead(@"c:\tet\fft.wav");
        byte[] data = new byte[WaveFile.Length];
        WaveFile.Read(data,0,Convert.ToInt32(WaveFile.Length));

将其转换为双精度数组:

 for (int i = 0; i < 32768; i++)//this is only for a relatively small chunk of the file 
        {
           InReal[i] =BitConverter.ToDouble(data, (i + 1) * 8 + 44);
        }

并最终将其传递给变换函数。

       FFT FftObject = new FFT();
       FftObject.Transform(InReal, InImg, 0, 32768, out outReal, out outImg, false);

现在第一个问题,据我了解,wav 文件的 PCM 值应该在 -1 和 1 的边界内,但是当转换为 Double 时,我得到了这个值:

 2.65855908666825E-235
 2.84104982662944E-285
-1.58613492930337E+235
-1.25617351166869E+264
 1.58370933499389E-242
 6.19284549187335E-245
-2.92969500042228E+254
-5.90042665390976E+226
 3.11954507295188E-273
 3.06831908609091E-217
 NaN
 2.77113146323761E-302
 6.76597919848376E-306
-1.55843653898344E+291

这些是这些限制中的第一个数组,也是数组的其余部分。

我的结论是我有某种代码故障,但我似乎能够找到它。任何帮助,将不胜感激。

第二个问题,因为我只在响应向量中向 FFT 算法提供真实数据,我是否也应该只期望真实部分数据?

非常感谢。

4

1 回答 1

2

我终于能够找出问题所在了准备波形文件的函数中的代码。

public static Double[] prepare(String wavePath, out int SampleRate)

    {
        Double[] data;
        byte[] wave;
        byte[] sR= new byte[4];
        System.IO.FileStream WaveFile = System.IO.File.OpenRead(wavePath);
        wave = new byte[WaveFile.Length];
        data = new Double[(wave.Length - 44) / 4];//shifting the headers out of the PCM data;
        WaveFile.Read(wave,0,Convert.ToInt32(WaveFile.Length));//read the wave file into the wave variable
        /***********Converting and PCM accounting***************/
        for (int i = 0; i < data.Length - i * 4; i++)
        {
            data[i] = (BitConverter.ToInt32(wave, (1 + i) * 4)) / 65536.0;
            //65536.0.0=2^n,       n=bits per sample;
        }
        /**************assigning sample rate**********************/
        for (int i = 24; i < 28; i++)
        {
            sR[i-24]= wave[i];
        }
        SampleRate = BitConverter.ToInt32(sR,0);
        return data;
    }

您现在需要做的就是将采样率和返回的结果发送到您的 FFT 算法。代码未处理,因此请根据需要自行处理。我已经测试了电话录音、忙音、铃声和语音,它可以正常工作。

于 2012-12-01T17:12:21.330 回答