Как можно реализовать функцию MATLAB wavread() в Python?

Question

В Matlab:

[s, fs] = wavread(file); 

Поскольку значение каждого элемента в S находится в пределах от -1 до 1, поэтому мы импортировать scikits.audiolab используя свою wavread функцию.

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

Но когда я красный и тот же файл ввода Wav, значение х в Matlab и Python были совершенно разные. Как я могу получить тот же результат, что и в MATLAB?

p.s. WAV-файл - это простой 16-битный файл с одним каналом, отснятый на частоте 44100 Гц.

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');