LI2020 - Sound externalization in reverberant spaces

Usage:

E = li2020(target_direct,template_direct,target_brir,template_brir)

`target`	binaural impulse response or head-related impulse response of the target sound. Matrix dimensions: time x receiver.
`template`	binaural impulse response or head-related impulse response of the template sound. Matrix dimensions: time x receiver. Note that the dimensions of target and template must be identical.
`target_brir`	target binaural room impulse response
`template_brir`	template binaural room impulse response

`E`	predicted degree of externalization

Examples:

Anechoic signal (rendered with HRTFs or direct sound part of BRIRs):

E = Li2020(tar_HRIR,tem_HRIR,[],[],'stim',stimulus,'fs',44100,'fsstim',44100,'flow',200, 'fhigh',16000, 'space',1);

Reverberant signal (rendered with BRIRs with unkonwn direct sound part):

E = Li2020([],[],tar_BRIR,tem_BRIR,'stim',stimulus,'fs',44100,'fsstim',44100,'flow',200, 'fhigh',16000, 'space',1);

Reverberant signal (rendered with BRIRs with konwn direct part (HRTFs)):

E = Li2020(tar_BRIR_dir,tem_BRIR_dir,tar_BRIR,tem_BRIR,'stim',stimulus,'fs',44100,'fsstim',44100,'flow',200, 'fhigh',16000, 'space',1);

The implementation is based on baumgartner2021.m

S. Li, R. Baumgartner, and J. Peissig. Modeling perceived externalization of a static, lateral sound image. Acta Acustica, 4(5), 2021.