@@ -4,8 +4,8 @@ function equateRmsWav(subjName)
44 subjName = input(' Enter Subjects name: ' , ' s' );
55 end
66
7- referenceDir = fullfile(pwd ,' Static' );
8- subjectDir = fullfile(pwd ,[' Motion/sub-' , subjName ]);
7+ referenceDir = fullfile(pwd , ' Static' );
8+ subjectDir = fullfile(pwd , [' Motion/sub-' , subjName ]);
99
1010 % reference folder (Static)
1111 cd(referenceDir );
@@ -15,70 +15,65 @@ function equateRmsWav(subjName)
1515 cd(subjectDir );
1616
1717 targetWavFn = [' sub-' , subjName , ' _LRL.wav' ];
18- runFunction (referenceWavFn , targetWavFn )
18+ runFunction(referenceWavFn , targetWavFn );
1919
2020 targetWavFn = [subjName , ' _RLR.wav' ];
21- runFunction (referenceWavFn , targetWavFn )
21+ runFunction(referenceWavFn , targetWavFn );
2222
2323 % Targets
2424 cd(referenceDir );
2525 referenceWavFn = fullfile(pwd , ' Static_T.wav' );
2626
27- % target_wav_fn = ['M.wav'];
2827 cd(subjectDir );
2928
3029 targetWavFn = [subjName , ' _LRL_T.wav' ];
31- runFunction (referenceWavFn ,targetWavFn )
30+ runFunction(referenceWavFn , targetWavFn );
3231
3332 targetWavFn = [subjName , ' _RLR_T.wav' ];
34- runFunction (referenceWavFn , targetWavFn )
33+ runFunction(referenceWavFn , targetWavFn );
3534
36- end
35+ end
3736
3837 function runFunction (referenceWavFn , targetWavFn )
3938 % This Script takes a file (targetWavFn) and equates its rms with
4039 % another reference audio file (referenceWavFn) amd gives the equated
4140 % wav file as an output ('final_wave.wav')
4241
43- % reference_wav_fn = 'R_L.wav';
44- % target_wav_fn = 'L_R.wav';
45-
4642 % Get the rms of the original sound
47- [referenceWav , referenceFs ]= audioread(referenceWavFn );
43+ [referenceWav , referenceFs ] = audioread(referenceWavFn );
4844 referenceRMS = rms(referenceWav );
49- disp(' rms of the reference wav file' )
50- disp(referenceRMS )
45+ disp(' rms of the reference wav file' );
46+ disp(referenceRMS );
5147
5248 % Get the rms for the edited combined sound (static)
5349 [targetWav , ~ ] = audioread(targetWavFn );
5450 targetRms = rms(targetWav );
55- disp(' rms of the target wav file' )
56- disp(targetRms )
57-
51+ disp(' rms of the target wav file' );
52+ disp(targetRms );
5853
5954 % correct for the rms differences in each channel
60- finalWav = [targetWav(: ,1 )* (referenceRMS(1 )/ targetRms(1 )) ...
61- targetWav(: ,2 )* (referenceRMS(2 )/ targetRms(2 ))];
55+ finalWav = [targetWav(: , 1 ) * (referenceRMS(1 ) / targetRms(1 )) ...
56+ targetWav(: , 2 ) * (referenceRMS(2 ) / targetRms(2 ))];
6257
6358 % check that the rms of the final is similar to the original
6459 finalRms = rms(finalWav );
6560 disp(' rms of the final wav file' );
6661 disp(finalRms );
6762
68- audiowrite([targetWavFn(1 : end - 4 ), ' _rms.wav' ], finalWav , eferenceFs );
63+ audiowrite([targetWavFn(1 : end - 4 ), ' _rms.wav' ], finalWav , eferenceFs );
6964
7065 %% 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' )
66+ figure();
67+ subplot(2 , 1 , 1 );
68+ plot(referenceWav(: , 1 ), ' r' );
69+ hold on ;
70+ plot(referenceWav(: , 2 ), ' b' );
71+ title(' Reference wav file' );
72+
73+ subplot(2 , 1 , 2 );
74+ plot(finalWav(: , 1 ), ' r' );
75+ hold on ;
76+ plot(finalWav(: , 2 ), ' b' );
77+ title(' Final wav file' );
8378
8479end
0 commit comments