% EXAMPLE for the usage of the combined audio quality model that is described
% in Fleßner et al. (2019). It is a reference-based audio quality model and based on
% monaural (intensity and amplitude modulation) and binaural (interaural level and
% phase differences, interaural vector strength) cues to predict subjective audio
% quality ratings. 

cur_dir=pwd;
addpath(genpath(cur_dir)); % include all folders from the current path
%% input signals

if verLessThan('matlab','8.0')
 % reference signal (clean)
[RefSig, fsRef] = wavread('Stimuli/guitar_ref.wav');
% test signal (processed)
[TestSig, fsTest] = wavread('Stimuli/guitar_midAnchor.wav');
else
 % reference signal (clean)
[RefSig, fsRef] = audioread('Stimuli/guitar_ref.wav');
% test signal (processed)
[TestSig, fsTest] = audioread('Stimuli/guitar_midAnchor.wav');
end

% compare sampling frequencies
if fsTest ~= fsRef,
    error('signals have different sampling frequencies')
else
    fs = fsTest;
end

% minlen = min([length(TestSig), length(RefSig)]);
% TestSig = TestSig(1:minlen,1);
% RefSig = RefSig(1:minlen,1);

% calculate objective perceptual measures
stOut = GPSMqBin(RefSig, TestSig, fs); % monaural model output
binQ = BAMQ(RefSig, TestSig, fs); % binaural model output

% **********************************************
% **********Combine outputs of BAMq and GPSMq***
% **********************************************
obj_meas_mon=stOut.out(:,1); % predicted monaural quality rating
obj_meas_bin=binQ(:,1); % predicted binaural quality rating
%*********f3************
mon_tmp=0.0528.*obj_meas_mon;
bin_tmp=0.0078*obj_meas_bin;
obj_meas=min(log10(mon_tmp),bin_tmp);

disp('***********************')
disp('****monaural measures**')
disp('***********************')
disp(stOut)
%        out: 5.7203
%     outFix: 1.7207
%     SNR_dc: 0.3909
% SNR_dc_fix: 1.0952
%     SNR_ac: 0.3909
% SNR_ac_fix: 1.0952
disp('***********************')
disp('***binaural measures***')
disp('***********************')
disp(binQ)
disp('***********************')
disp('****overall measures***')
disp('***********************')
disp(obj_meas)