Измерьте 532-1 ISO стационарную громкость свободного поля. Предположим, что уровень записи калиброван так, что тональный сигнал 1 кГц регистрируется как 100 дБ на SPL-метре.

[audioIn,fs] = audioread('WashingMachine-16-44p1-stereo-10secs.wav');

loudness = acousticLoudness(audioIn,fs)

loudness = 1×2

   28.2688   27.7643

Создайте два стационарных сигнала с эквивалентной степенью: розовый сигнал шума и белый сигнал шума.

fs = 48e3;
dur = 5;
pnoise = 2*pinknoise(dur*fs);
wnoise = rand(dur*fs,1) - 0.5;
wnoise = wnoise*sqrt(var(pnoise)/var(wnoise));

Функции acousticLoudness использование метода ISO 532-1 (Zwicker) по умолчанию и отсутствие выходных аргументов для построения графика громкости розового шума. Функции acousticLoudness снова, на этот раз с выходными аргументами, чтобы получить определенную громкость.

figure
acousticLoudness(pnoise,fs)

[~,pSpecificLoudness] = acousticLoudness(pnoise,fs);

Постройте график громкости для сигнала белого шума и затем получите определенные значения громкости.

figure
acousticLoudness(wnoise,fs)

[~,wSpecificLoudness] = acousticLoudness(wnoise,fs);

Вызовите acousticSharpness функция для сравнения резкости розового шума и белого шума.

pSharpness = acousticSharpness(pSpecificLoudness);
wSharpness = acousticSharpness(wSpecificLoudness);
fprintf('Sharpness of pink noise = %0.2f acum\n',pSharpness)

Sharpness of pink noise = 2.00 acum

fprintf('Sharpness of white noise = %0.2f acum\n',wSharpness)

Sharpness of white noise = 2.62 acum

Чтение в аудио файла.

[audioIn,fs] = audioread('JetAirplane-16-11p025-mono-16secs.wav');

Постройте график изменяющейся во времени акустической громкости в соответствии с ISO 532-1 и получите процентили. Слушайте аудиосигнал.

acousticLoudness(audioIn,fs,'SoundField','diffuse','TimeVarying',true)

sound(audioIn,fs)

Функции acousticLoudness снова с теми же входами и получить процентили. Напечатайте Nmax и N5 процентиля. Процентиль Nmax - это максимальная сообщаемая громкость. N5 процентиль - это громкость ниже, которая составляет 95% от сообщаемой громкости.

[~,~,perc] = acousticLoudness(audioIn,fs,'SoundField','diffuse','TimeVarying',true);
fprintf('Max loudness = %0.2f sones\n',perc(1))

Max loudness = 89.48 sones

fprintf('N5 loudness = %0.2f sones\n',perc(2))

N5 loudness = 81.77 sones

Чтение в аудио файла.

[audioIn,fs] = audioread('Turbine-16-44p1-mono-22secs.wav');

Функции acousticLoudness без выходных аргументов для построения графика определенной громкости. Примите коэффициент калибровки 0,15 и давление ссылки 21 микропаскали. Чтобы определить коэффициент калибровки, характерный для вашей аудиосистемы, используйте calibrateMicrophone функция.

calibrationFactor = 0.15;
refPressure = 21e-6;
acousticLoudness(audioIn,fs,calibrationFactor,'PressureReference',refPressure)

acousticLoudness позволяет вам задать промежуточное представление, уровни звукового давления, вместо входов во временной области. Это позволяет повторно использовать промежуточные расчеты SPL. Другое преимущество состоит в том, что если ваш физический SPL-счетчик не сообщает о громкости в соответствии с ISO 532-1 или ISO 531-2, можно использовать сообщенные 1/3-октавные SPL для вычисления совместимой со стандартом громкости.

Чтобы вычислить уровни звукового давления из аудиосигнала, сначала создайте splMeter объект. Вызовите splMeter объект со звуковым входом.

spl = splMeter("SampleRate",fs,"Bandwidth","1/3 octave", ...
    "CalibrationFactor",calibrationFactor,"PressureReference",refPressure, ...
    "FrequencyWeighting","Z-weighting","OctaveFilterOrder",6);

splMeasurement = spl(audioIn);

Вычислите средний уровень SPL, пропустив первые 0,2 секунды. Только держите полосы от 25 Гц до 12,5 кГц (первые 28 полос).

SPLIn = mean(splMeasurement(ceil(0.2*fs):end,1:28));

Используя вход SPL, вызовите acousticLoudness без выходных аргументов для построения графика определенной громкости.

acousticLoudness(SPLIn)

Настройте эксперимент, как указано схемой.

Создайте audioDeviceReader объект для чтения с микрофона и audioDeviceWriter объект для записи в динамик.

fs = 48e3;
deviceReader = audioDeviceReader(fs);
deviceWriter = audioDeviceWriter(fs);

Создайте audioOscillator объект для генерации синусоиды 1 кГц.

osc = audioOscillator("sine",1e3,"SampleRate",fs);

Создайте dsp.AsyncBuffer объект для буферизации данных, полученных из микрофона.

dur = 5;
buff = dsp.AsyncBuffer(dur*fs);

В течение пяти секунд проигрывайте синусоиду через динамик и записывайте с помощью микрофона. Пока аудиопотоки, обратите внимание на громкость, о которой сообщил ваш SPL-счетчик. После завершения считайте содержимое объекта буфера.

numFrames = dur*(fs/osc.SamplesPerFrame);
for ii = 1:numFrames
    audioOut = osc();
    deviceWriter(audioOut);
    
    audioIn = deviceReader();
    write(buff,audioIn);
end

SPLreading = 60.4;

micRecording = read(buff);

Чтобы вычислить коэффициент калибровки для микрофона, используйте calibrateMicrophone функция.

calibrationFactor = calibrateMicrophone(micRecording,deviceReader.SampleRate,SPLreading);

Функции acousticLoudness с записью микрофона, частотой дискретизации и коэффициентом калибровки. Громкость, о которой сообщают acousticLoudness - истинное измерение акустической громкости, заданное в 532-1.

loudness = acousticLoudness(micRecording,deviceReader.SampleRate,calibrationFactor)

loudness = 14.7902

Теперь можно использовать коэффициент калибровки, который вы определили, чтобы измерить громкость любого звука, который получен через ту же цепь записи микрофона.

Считывайте аудиосигнал.

[audioIn,fs] = audioread('TrainWhistle-16-44p1-mono-9secs.wav');

[~,specificLoudness] = acousticLoudness(audioIn,fs,'TimeVarying',true);

barkBins = 0.1:0.1:24;
hzBins = bark2hz(barkBins);

t = 0:2e-3:2e-3*(size(specificLoudness,1)-1);
surf(t,hzBins,sum(specificLoudness,3).','EdgeColor','interp')
set(gca,'YScale','log')
view([0 90])
axis tight
xlabel('Time (s)')
ylabel('Frequency (Hz)')
colorbar
title('Specific Loudness (sones/Bark)')

[~,specificLoudness] = acousticLoudness(audioIn,fs,'Method','ISO 532-2');

erbBins = 1.8:0.1:38.9;
hzBins = erb2hz(erbBins);

semilogx(hzBins,specificLoudness)
xlabel('Frequency (Hz)')
ylabel('Loudness (sones)')
title('Specific Loudness')
grid on

Чтение в аудио файла.

[x,fs] = audioread('WashingMachine-16-44p1-stereo-10secs.wav');

ISO 532-2 позволяет вам задать пользовательский ответ наушника при вычислении громкости. Создайте матрицу 30 на 2, где первый столбец является частотой, а второй - отклонением наушника от плавного ответа.

tdh = [   0,    80,   100,  200,  500,  574,  660,  758,  871, 1000, 1149,  1320,  1516,  1741,  2000, ...
       2297,  2639,  3031, 3482, 4000, 4500, 5000, 5743, 6598, 7579, 8706, 10000, 12000, 16000, 20000; ...
        -50, -15.3, -13.8, -8.1, -0.5,  0.4,  0.8,  0.9,  0.5,  0.1, -0.8,  -1.5,  -2.3,  -3.2,  -3.9, ...
       -4.2,  -4.3,  -4.3, -3.9, -3.2, -2.3, -1.1, -0.3,   -2, -5.4,   -9, -12.1, -15.2,   -30,   -50 ].';

Вычислите громкость по ISO 532-2. Задайте SoundField как earphones и ответ наушника как матрица, которую вы только что создали.

acousticLoudness(x,fs,'Method','ISO 532-2','SoundField','earphones','EarphoneResponse',tdh)

Создайте dsp.AudioFileReader объект для чтения в кадр за кадром аудиосигнала. Задайте продолжительность системы координат 50 мс. Это будет системой координат продолжительность, в течение которой вы вычисляете стационарную громкость.

fileReader = dsp.AudioFileReader('Engine-16-44p1-stereo-20sec.wav');

frameDur = 0.05;
fileReader.SamplesPerFrame = round(fileReader.SampleRate*frameDur);

Создайте audioDeviceWriter объект для записи аудио на устройство выхода по умолчанию.

deviceWriter = audioDeviceWriter('SampleRate',fileReader.SampleRate);

Создайте timescope объект для отображения стационарной громкости с течением времени.

scope = timescope( ...
    'SampleRate',1/frameDur, ...
    'YLabel','Loudness (sones)', ...
    'ShowGrid',true, ...
    'PlotType','Stairs', ...
    'TimeSpanSource','property', ...
    'TimeSpan',20, ...
    'AxesScaling','Auto', ...
    'ShowLegend',true);

В цикле:

while ~isDone(fileReader)
    audioIn = fileReader();
    loudness = acousticLoudness(audioIn,fileReader.SampleRate);
    deviceWriter(audioIn);
    scope(loudness)
end
release(fileReader)
release(deviceWriter)
release(scope)