b@5: % test loudness mod
b@5: 
b@5: % target loudness and test file
b@5: target_loudness = -18;  % dBLU
b@5: filename = '';
b@5: 
b@5: % import sound file
b@5: [audio, fs] = audioread(filename); 
b@5: 
b@5: % measure and display loudness
b@5: initial_loudness = getloudness(audio, fs, 'ITU', 0);
b@5: disp(['Initial loudness: ' num2str(initial_loudness)]);
b@5: 
b@5: % apply gain
b@5: difference_loudness = target_loudness - initial_loudness;
b@5: disp(['Difference in loudness: ' num2str(difference_loudness)]);
b@5: audio = 10^(difference_loudness/10) .* audio;
b@5: disp(['Gain: ' num2str(10^(difference_loudness/10))]);
b@5: 
b@5: % measure and display loudness (should be equal to target loudness)
b@5: resulting_loudness = getloudness(audio, fs, 'ITU', 0);
b@5: disp(['Resulting loudness: ' num2str(resulting_loudness)]);