b@5: % test loudness mod
b@5: 
b@6: tic; 
b@6: 
b@5: % target loudness and test file
b@5: target_loudness = -18;  % dBLU
b@6: filename = 'pro.wav';
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@6: audio = 10^(difference_loudness/20) .* audio;
b@6: disp(['Gain: ' num2str(10^(difference_loudness/20))]);
b@5: 
b@5: % measure and display loudness (should be equal to target loudness)
b@5: resulting_loudness = getloudness(audio, fs, 'ITU', 0);
b@6: disp(['Resulting loudness: ' num2str(resulting_loudness)]);
b@6: 
b@6: elapsed_time = toc;
b@6: disp(['Elapsed time: ' num2str(elapsed_time)]);