comparison Code/Descriptors/Matlab/Common/AudioSpectralBandPower.m @ 4:92ca03a8fa99 tip

Update to ICASSP 2013 benchmark

+function [ASBP] = AudioSpectralBandPower(raw, fs, noOfFrames, frameLength)
+% this file normalises the frequency spectrum and calculates the spectral
+% power in frequency band 100 - 1000 Hz and in the top 
+x = buffer( raw, frameLength );
+myColor = 'rx';
+
+hiEdgeSamples = getMaxFreq; % the highest fundamental
+loEdgeSamples = getMinFreq; % the lowest fundamental
+
+hiFreqBin = ceil(((frameLength/2)/(fs/2))*hiEdgeSamples);
+loFreqBin = floor(((frameLength/2)/(fs/2))*loEdgeSamples);
+
+bandWidth = hiFreqBin - loFreqBin;
+
+for n = 1:noOfFrames
+    spectrum = abs( fft( x(:,n) ));
+    spectrum = spectrum( 1 : frameLength/2 );
+    % normalise
+%     spectrum = spectrum/max(spectrum);
+    % fundamental frequency band
+    lBand = ( loFreqBin : hiFreqBin );
+    hBand = ( frameLength/2 - ( bandWidth ) : frameLength/2 );
+    lBandPower = sum( spectrum( lBand ));
+    hBandPower = sum( spectrum( hBand ));
+%     gm(n) = nthroot( prod( spectrum ), frameLength );
+%     am(n) = sum( spectrum )/frameLength;
+    ASBP(n) = lBandPower/hBandPower;
+%     subplot(421); plot( x(:,n) );
+%     subplot(422); plot( spectrum ); hold on; 
+%     plot( loFreqBin : hiFreqBin, spectrum( lBand ), 'r' );
+%     plot( hBand, spectrum( hBand ), 'r' ); hold off
+%     
+% %     look at the frame with reference to the surrounding 10 frames
+%     if( n>5 )
+%         startPlot = (n-5)*frameLength+1;
+%         endPlot = (n+5)*frameLength;
+%     else
+%         startPlot = 1;
+%         endPlot = 10*frameLength;
+%     end
+%     subplot(412); plot( raw( startPlot : endPlot ) ); hold on; 
+%     plot( ((n-1)*frameLength+1) - startPlot : (n*frameLength) - startPlot, x(:,n), 'r' ); hold off;
+%     
+%     subplot(425); plot( n, lBandPower, myColor ); hold on;%plot( n, gm(n), myColor ); hold on;
+%     subplot(426); plot( n, hBandPower, myColor ); hold on;%plot( n, am(n), myColor ); hold on;
+%     subplot(414); plot( n, lBandPower / hBandPower, myColor ); hold on;%plot( n, ASF(n), myColor ); hold on;
+% end
+end
+
+
+
+
+