annotate aux/loudnesstest.m @ 6:ed0a8621b96a

prepare2listen small fixes, manual, standard target loudness (-23 dBLUFS)
author Brecht De Man <b.deman@qmul.ac.uk>
date Mon, 24 Nov 2014 16:55:55 +0000
parents 85bff3d1b6fe
children 0014c50188da
rev   line source
b@5 1 % test loudness mod
b@5 2
b@6 3 tic;
b@6 4
b@5 5 % target loudness and test file
b@5 6 target_loudness = -18; % dBLU
b@6 7 filename = 'pro.wav';
b@5 8
b@5 9 % import sound file
b@5 10 [audio, fs] = audioread(filename);
b@5 11
b@5 12 % measure and display loudness
b@5 13 initial_loudness = getloudness(audio, fs, 'ITU', 0);
b@5 14 disp(['Initial loudness: ' num2str(initial_loudness)]);
b@5 15
b@5 16 % apply gain
b@5 17 difference_loudness = target_loudness - initial_loudness;
b@5 18 disp(['Difference in loudness: ' num2str(difference_loudness)]);
b@6 19 audio = 10^(difference_loudness/20) .* audio;
b@6 20 disp(['Gain: ' num2str(10^(difference_loudness/20))]);
b@5 21
b@5 22 % measure and display loudness (should be equal to target loudness)
b@5 23 resulting_loudness = getloudness(audio, fs, 'ITU', 0);
b@6 24 disp(['Resulting loudness: ' num2str(resulting_loudness)]);
b@6 25
b@6 26 elapsed_time = toc;
b@6 27 disp(['Elapsed time: ' num2str(elapsed_time)]);