dawn@0: function [ vuv ] = calculate_VoicedUnvoicedDecision( sampleName, x, fs, frameLength, noOfFrames, silentFrames )
dawn@0: 
Dawn@4: DEBUG = 1; % do we want to plot the output?
dawn@0: 
dawn@0: y = buffer( x, frameLength );
dawn@0: 
dawn@0: % % open original annotation file for PhD speech
dawn@0: % fileName = [sampleName '.txt'];
dawn@0: % % read metrics from file
dawn@0: % fileID = fopen( fileName );
dawn@0: % data = fscanf( fileID, '%d', inf );
dawn@0: % annotatedPitch = data(3:4:end);
dawn@0: % fclose( fileID );
dawn@0: 
dawn@0: % open HNR file
dawn@0: fileName = [sampleName '_HNR.txt'];
dawn@0: % read metrics from file
dawn@0: fileID = fopen( fileName );
dawn@0: data = fscanf( fileID, '%f', inf );
dawn@0: harmonic2noise = data(2:2:end);
dawn@0: fclose( fileID );
dawn@0: 
dawn@0: largestH2NValue = max( [max(harmonic2noise) abs(min(harmonic2noise))] );
dawn@0: %normalise
dawn@0: harmonic2noise = harmonic2noise/largestH2NValue;
dawn@0: %round to 1dp
dawn@0: harmonic2noise= (round(harmonic2noise*10))/10;
dawn@0: 
dawn@0: % open audio power file
dawn@0: fileName = [sampleName '_AP.txt'];
dawn@0: % read metrics from file
dawn@0: fileID = fopen( fileName );
dawn@0: data = fscanf( fileID, '%f', inf );
dawn@0: audioPower = data(2:2:end);
dawn@0: fclose( fileID );
dawn@0: maxPower = max(audioPower);
dawn@0: 
dawn@0: 
dawn@0: minFreq = ceil(fs/getVariables('getMinFreq')); %default minimum fundatmental frequency
dawn@0: maxFreq = ceil(fs/getVariables('getMaxFreq')); %default maximum fundamental frequency
dawn@0:     
dawn@0: % open pitch file
dawn@0: fileName = [sampleName '_YIN_pitch.txt'];
dawn@0: % read metrics from file
dawn@0: fileID = fopen( fileName );
dawn@0: data = fscanf( fileID, '%f', inf );
dawn@0: pitch = data(2:2:end);
dawn@0: fclose( fileID );
dawn@0: 
dawn@0: % open band power file
dawn@0: % fileName = [sampleName '_ASBP_norm.txt'];
dawn@0: % % read metrics from file
dawn@0: % fileID = fopen( fileName );
dawn@0: % data = fscanf( fileID, '%f', inf );
dawn@0: % bandPowerRatio = data(2:2:end);
dawn@0: % fclose( fileID );
dawn@0: % bandPowerRatioRMS =   sqrt( mean(( bandPowerRatio/max(bandPowerRatio) ).^2 ));
dawn@0: % 
dawn@0: 
dawn@0: 
dawn@0: HNRThresh = 0;
dawn@0: if(DEBUG) hold off; end
dawn@0: minLength = min([length(harmonic2noise) length(audioPower) length(pitch)]);
dawn@0: for( i=1:minLength )
dawn@0: 
dawn@0: 
dawn@0:      if( silentFrames(i) == 0 ) % harmonic2noise(i) == 0 )
dawn@0:          %silent
dawn@0:         if(DEBUG) plot( ((i-1)*frameLength)+1:(i*frameLength), y(:,i), 'y' );hold on; end;
dawn@0:         vuv(i) = 0;
dawn@0:      elseif( ((harmonic2noise(i) <= HNRThresh))... %|| (bandPowerRatio(i)/max(bandPowerRatio) < 0.001)) ...
dawn@0:             || (pitch(i) > minFreq) ...
dawn@0:             || (pitch(i) < maxFreq)) %...
dawn@0: %             || (audioPower(i)/maxPower < 0.0005))
dawn@0:         % unvoiced
dawn@0:         if(DEBUG) plot( ((i-1)*frameLength)+1:(i*frameLength), y(:,i), 'r' );hold on; end
dawn@0:         vuv(i) = 2;
dawn@0:      else
dawn@0:          %voiced
dawn@0:         if(DEBUG) plot( ((i-1)*frameLength)+1:(i*frameLength), y(:,i), 'g' );hold on; end
dawn@0:         vuv(i) = 1;
dawn@0:     end
dawn@0: end
dawn@0: 
dawn@0: if(DEBUG)
dawn@0:     
dawn@0:     plot( frameLength/2:frameLength:(frameLength*length(harmonic2noise)), harmonic2noise, 'b' );hold on;
dawn@0: %     plot( frameLength/2:frameLength:(frameLength*length(audioPower)), audioPower/maxPower, 'm' );hold on;
dawn@0:     plot( frameLength/2:frameLength:(frameLength*length(pitch)), pitch/max(pitch), 'w' );hold on;
dawn@0: %     plot( frameLength/2:frameLength:(frameLength*length(bandPowerRatio)), bandPowerRatio/max(bandPowerRatio), 'c' );hold on;
dawn@0: 
dawn@0:     L=line([0 length(x)], [maxFreq/max(pitch) maxFreq/max(pitch)]);
dawn@0:     set(L,'color',[1 0 0]);
dawn@0:     L=line([0 length(x)], [minFreq/max(pitch) minFreq/max(pitch)]);
dawn@0:     set(L,'color',[1 0 0]);
dawn@0: %     L=line([0 length(x)], [0.0001 0.0005]);
dawn@0: %     set(L, 'color', [0 1 1] );
dawn@0:     %HNR threshold
dawn@0:     L=line([0 length(x)], [HNRThresh HNRThresh]);
dawn@0:     set(L, 'color', [0 0 1] );
dawn@0: end
dawn@0: 
dawn@0: