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)]); |