0

我正在尝试用 Python 编写的树莓派构建一个 SPL 测量设备。我录制了一个声音，然后我做了 FFT，但我没有正确的对数值。

    fs_rate, signal = wavfile.read(filename)
    l_audio = len(signal.shape)
    if l_audio == 2:
        signal = signal.sum(axis=1) / 2
    N = signal.shape[0]
    secs = (N / float(fs_rate))
    Ts = 1.0/fs_rate # sampling interval in time
    t = scipy.arange(0, secs, Ts) # time vector as scipy arange field / numpy.ndarray
    FFT = abs(scipy.fft(signal))
    FFT_side = FFT[range(N//2)] # one side FFT range
    freqs = scipy.fftpack.fftfreq(signal.size, t[1]-t[0])
    fft_freqs = np.array(freqs)
    freqs_side = freqs[range(N//2)] # one side frequency range
    fft_freqs_side = np.array(freqs_side)
    fft_log = 94 + np.log10((FFT_side/(1.79*10**308)))       
    print("Type: %s" % (FFT_side.dtype))

上面你可以看到我正在使用的 FFT 代码。我有两个问题：

  1. FFT_side 是一个 float64 值，最大范围是多少？

我的处理器是 armv7l 32 位。我在互联网上发现这是最大范围 1.7976931348623157e+308，对吗？

  1. 如何在不知道麦克风两端电压的情况下获得 dB 图？我想做 value_read/max_float64 但这根本不起作用。
4

0 回答 0