MATLAB 函数 wavread() 如何在 Python 中实现？

Question

在 Matlab 中：

[s, fs] = wavread(file);  

由于 s 中每个元素的值都在 -1 和 1 之间，所以我们使用其 wavread 函数导入 scikits.audiolab。

from scikits.audiolab import wavread
s, fs, enc = wavread(filename)

但是当我对同一个输入wav文件进行reduce时，Matlab和Python中s的值是完全不同的。如何获得与 MATLAB 中相同的 s 输出？

附言wav 文件是简单的 16 位单声道文件，以 44100Hz 采样。

Answer 1

Matlab wavread() 函数默认返回标准化数据，即将所有数据缩放到 -1 到 +1 的范围内。如果您的音频文件是 16 位格式，那么原始数据值将在 -32768 到 +32767 的范围内，这应该与 scikits.audiolab.wavread() 返回的范围相匹配。

要将这些数据标准化到 -1 到 +1 的范围内，您只需将数据数组除以具有最大幅度的值（在这种情况下使用 numpy 模块）：

from scikits.audiolab import wavread
import numpy

s, fs, enc = wavread(filename)  # s in range -32768 to +32767
s = numpy.array(s)
s = s / max(abs(s))             # s in range -1 to +1

请注意，您还可以将 'native' 选项与 Matlab 函数一起使用以返回未标准化的数据值，正如 horchler 在他的评论中所建议的那样。

[s, fs] = wavread(file, 'native');