Mercurial > hg > ape
comparison aux/loudnesstest.m @ 5:85bff3d1b6fe
Equalise loudness in prepare2listen.m
author | Brecht De Man <b.deman@qmul.ac.uk> |
---|---|
date | Sun, 23 Nov 2014 23:24:35 +0000 |
parents | |
children | ed0a8621b96a |
comparison
equal
deleted
inserted
replaced
4:b28ffd29e6e1 | 5:85bff3d1b6fe |
---|---|
1 % test loudness mod | |
2 | |
3 % target loudness and test file | |
4 target_loudness = -18; % dBLU | |
5 filename = ''; | |
6 | |
7 % import sound file | |
8 [audio, fs] = audioread(filename); | |
9 | |
10 % measure and display loudness | |
11 initial_loudness = getloudness(audio, fs, 'ITU', 0); | |
12 disp(['Initial loudness: ' num2str(initial_loudness)]); | |
13 | |
14 % apply gain | |
15 difference_loudness = target_loudness - initial_loudness; | |
16 disp(['Difference in loudness: ' num2str(difference_loudness)]); | |
17 audio = 10^(difference_loudness/10) .* audio; | |
18 disp(['Gain: ' num2str(10^(difference_loudness/10))]); | |
19 | |
20 % measure and display loudness (should be equal to target loudness) | |
21 resulting_loudness = getloudness(audio, fs, 'ITU', 0); | |
22 disp(['Resulting loudness: ' num2str(resulting_loudness)]); |