%% Acoustic Sharpness of Audio Signal
% Compute the acoustic sharpness of turbine noise. Assume it is stationary and was recorded in a diffuse sound field.

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

sharpness = acousticSharpness(audioIn,fs,'SoundField','diffuse');

fprintf('Acoustic sharpness = %0.2f acum\n',sharpness)

%
% Copyright 2019 The MathWorks, Inc.