annotate aux/loudnesstest.m @ 9:ee22af6610a3
Sanitised user ID, end message, bug fixes.
author |
Brecht De Man <b.deman@qmul.ac.uk> |
date |
Sun, 30 Nov 2014 21:44:03 +0000 |
parents |
ed0a8621b96a |
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)]); |