Las

Minor bugfixes (room acoustics, strength calculation) and additions (DAFFv15 examples)

applications/RoomAcoustics/AdditionalRoutines/ita_roomacoustics_sound_strength.m

@@ -81,12 +81,13 @@ if isa(reference, 'itaAudio')
     % upper inegration limit is intersection time
     [del1, snr, intersectionTime] = ita_roomacoustics_reverberation_time_lundeby(refGefiltert, 'freqRange', sArgs.freqRange, 'bandsPerOctave', sArgs.bandsPerOctave,  'shortRevTimeMode');
     idxEnd = round(refGefiltert.samplingRate * intersectionTime.freqData);
-    
     refGefiltertData = refGefiltert.timeData.^2 ;
+    idxEnd = min(idxEnd,length(refGefiltertData)-1);
+
     % TODO: checkw as besser ist max vs ita_start IR
     % %     [del idxStart] = max(refGefiltertData);
     idxStart = ita_start_IR(refGefiltert);
-    
+    idxStart = min(idxStart,length(refGefiltertData)-1);
     % integrate
     L_ref  = zeros(numel(freqVec),1);
     for iFreq = 1:numel(freqVec)
@@ -178,6 +179,10 @@ for iCh = 1:nChannels
         end
     
     % integrate
+    
+    idxEnd = min(idxEnd,rirBands.nSamples-1);
+    idxStart = min(idxStart,rirBands.nSamples-1);
+    
     rirBandsData = rirBands.timeData.^2 ;
     for iFreq = 1:numel(freqVec)
         if any(isnan([idxStart(iFreq) idxEnd(iFreq)]))