% 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 (ILDs, IPDs), interaural vector strength (IVS)) 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


%% monaural GPSMq
stOut = GPSMqBin(RefSig, TestSig, fs); % monaural model output
%% binaural BAM-Q
[binQ, ILDdiff, ITDdiff, IVSdiff] = BAMQpc(RefSig, TestSig, fs);
%% Combine outputs of BAMq and GPSMq
[obj_meas]=combine_binQ_OPM(stOut.out(:,1), binQ(:,1));


disp('***********************')
disp('****monaural measures**')
disp('***********************')
disp(stOut)
%     SNR_dc: 0.0179
%     SNR_ac: 0.0225
% SNR_dc_fix: 0.0179
% SNR_ac_fix: 0.0225
%        out: 132.6339
%     outFix: 132.7494
disp('***********************')
disp('***binaural measures***')
disp('***********************')
disp(['binQ: ',num2str(binQ)])
%     binQ: 100

disp('***********************')
disp('****overall measures***')
disp('***********************')
disp(['overall_measure: ',num2str(obj_meas)])
%     overall measure: 0.7800