|
1 | 1 | function equateRmsWav(subjName) |
2 | 2 |
|
3 | | -if nargin<1 |
4 | | - subjName = input('Enter Subjects name: ','s'); |
5 | | -end |
| 3 | + if nargin < 1 |
| 4 | + subjName = input('Enter Subjects name: ', 's'); |
| 5 | + end |
6 | 6 |
|
7 | | -referenceDir = fullfile(pwd,'Static'); |
8 | | -subjectDir = fullfile(pwd,[ 'Motion/sub-',subjName ]); |
| 7 | + referenceDir = fullfile(pwd,'Static'); |
| 8 | + subjectDir = fullfile(pwd,['Motion/sub-', subjName]); |
9 | 9 |
|
10 | | -% reference folder (Static) |
11 | | -cd(referenceDir) |
12 | | -referenceWavFn = fullfile(pwd,'Static.wav'); |
| 10 | + % reference folder (Static) |
| 11 | + cd(referenceDir); |
| 12 | + referenceWavFn = fullfile(pwd, 'Static.wav'); |
13 | 13 |
|
14 | | -% Subject motion folder |
15 | | -cd(subjectDir) |
| 14 | + % Subject motion folder |
| 15 | + cd(subjectDir); |
16 | 16 |
|
17 | | -targetWavFn = ['sub-', subjName,'_LRL.wav']; |
18 | | -runFunction (referenceWavFn,targetWavFn) |
| 17 | + targetWavFn = ['sub-', subjName, '_LRL.wav']; |
| 18 | + runFunction (referenceWavFn, targetWavFn) |
19 | 19 |
|
20 | | -targetWavFn = [subjName,'_RLR.wav']; |
21 | | -runFunction (referenceWavFn,targetWavFn) |
| 20 | + targetWavFn = [subjName, '_RLR.wav']; |
| 21 | + runFunction (referenceWavFn, targetWavFn) |
22 | 22 |
|
23 | | -% Targets |
24 | | -cd(referenceDir) |
25 | | -referenceWavFn = fullfile(pwd,'Static_T.wav'); |
| 23 | + % Targets |
| 24 | + cd(referenceDir); |
| 25 | + referenceWavFn = fullfile(pwd, 'Static_T.wav'); |
26 | 26 |
|
27 | | -%target_wav_fn = ['M.wav']; |
28 | | -cd(subjectDir) |
| 27 | + %target_wav_fn = ['M.wav']; |
| 28 | + cd(subjectDir); |
29 | 29 |
|
30 | | -targetWavFn = [subjName,'_LRL_T.wav']; |
31 | | -runFunction (referenceWavFn,targetWavFn) |
| 30 | + targetWavFn = [subjName, '_LRL_T.wav']; |
| 31 | + runFunction (referenceWavFn,targetWavFn) |
32 | 32 |
|
33 | | -targetWavFn = [subjName,'_RLR_T.wav']; |
34 | | -runFunction (referenceWavFn,targetWavFn) |
| 33 | + targetWavFn = [subjName, '_RLR_T.wav']; |
| 34 | + runFunction (referenceWavFn, targetWavFn) |
35 | 35 |
|
36 | | -end |
| 36 | + end |
| 37 | + |
| 38 | + function runFunction (referenceWavFn, targetWavFn) |
| 39 | + % This Script takes a file (targetWavFn) and equates its rms with |
| 40 | + % another reference audio file (referenceWavFn) amd gives the equated |
| 41 | + % wav file as an output ('final_wave.wav') |
| 42 | + |
| 43 | + %reference_wav_fn = 'R_L.wav'; |
| 44 | + %target_wav_fn = 'L_R.wav'; |
| 45 | + |
| 46 | + % Get the rms of the original sound |
| 47 | + [referenceWav, referenceFs]= audioread(referenceWavFn); |
| 48 | + referenceRMS = rms(referenceWav); |
| 49 | + disp('rms of the reference wav file') |
| 50 | + disp(referenceRMS) |
| 51 | + |
| 52 | + % Get the rms for the edited combined sound (static) |
| 53 | + [targetWav, ~] = audioread(targetWavFn); |
| 54 | + targetRms = rms(targetWav); |
| 55 | + disp('rms of the target wav file') |
| 56 | + disp(targetRms) |
| 57 | + |
| 58 | + |
| 59 | + % correct for the rms differences in each channel |
| 60 | + finalWav = [targetWav(:,1)*(referenceRMS(1)/targetRms(1)) ... |
| 61 | + targetWav(:,2)*(referenceRMS(2)/targetRms(2))]; |
| 62 | + |
| 63 | + % check that the rms of the final is similar to the original |
| 64 | + finalRms = rms(finalWav); |
| 65 | + disp('rms of the final wav file'); |
| 66 | + disp(finalRms); |
| 67 | + |
| 68 | + audiowrite([targetWavFn(1:end-4), '_rms.wav'], finalWav, eferenceFs); |
| 69 | + |
| 70 | + %% plot the reference wav and final wav files |
| 71 | + figure() |
| 72 | + subplot(2, 1, 1) |
| 73 | + plot(referenceWav(:,1), 'r') |
| 74 | + hold on |
| 75 | + plot(referenceWav(:,2), 'b') |
| 76 | + title('Reference wav file') |
37 | 77 |
|
38 | | -function runFunction (referenceWavFn,targetWavFn) |
39 | | -% This Script takes a file (targetWavFn) and equates its rms with |
40 | | -% another reference audio file (referenceWavFn) amd gives the equated |
41 | | -% wav file as an output ('final_wave.wav') |
42 | | - |
43 | | -%reference_wav_fn = 'R_L.wav'; |
44 | | -%target_wav_fn = 'L_R.wav'; |
45 | | - |
46 | | -% Get the rms of the original sound |
47 | | -[referenceWav , referenceFs]= audioread(referenceWavFn); |
48 | | -referenceRMS = rms(referenceWav) ; |
49 | | -disp('rms of the reference wav file') |
50 | | -disp(referenceRMS) |
51 | | - |
52 | | -% Get the rms for the edited combined sound (static) |
53 | | -[targetWav, ~] = audioread(targetWavFn); |
54 | | -targetRms = rms(targetWav) ; |
55 | | -disp('rms of the target wav file') |
56 | | -disp(targetRms) |
57 | | - |
58 | | - |
59 | | -% correct for the rms differences in each channel |
60 | | -finalWav = [ targetWav(:,1)*(referenceRMS(1)/targetRms(1)) ... |
61 | | - targetWav(:,2)*(referenceRMS(2)/targetRms(2))] ; |
62 | | - |
63 | | -% check that the rms of the final is similar to the original |
64 | | -finalRms = rms(finalWav); |
65 | | -disp('rms of the final wav file') |
66 | | -disp(finalRms) |
67 | | - |
68 | | -audiowrite([targetWavFn(1:end-4), '_rms.wav'],finalWav,referenceFs) |
69 | | - |
70 | | -%% plot the reference wav and final wav files |
71 | | -figure() |
72 | | -subplot(2,1,1) |
73 | | -plot(referenceWav(:,1),'r') |
74 | | -hold on |
75 | | -plot(referenceWav(:,2),'b') |
76 | | -title('Reference wav file') |
77 | | - |
78 | | -subplot(2,1,2) |
79 | | -plot(finalWav(:,1),'r') |
80 | | -hold on |
81 | | -plot(finalWav(:,2),'b') |
82 | | -title('Final wav file') |
| 78 | + subplot(2, 1, 2) |
| 79 | + plot(finalWav(:,1), 'r') |
| 80 | + hold on |
| 81 | + plot(finalWav(:,2), 'b') |
| 82 | + title('Final wav file') |
83 | 83 |
|
84 | 84 | end |
0 commit comments