[...] = wierstorf2013(X,Y,phi,xs,src,L,resolution,method,...);
X  range of the xaxis [xmin,xmax] in m (for a single x point you have to specify [x x] 
Y  range of the yaxis [ymin ymax] in m (for a single y point you have to specify [y y] 
phi  orientation of the listener in rad (0 is in the direction of the xaxis 
xs  position of the point source in m / direction of the plane wave 
src 

L  length/diameter of the loudspeaker array 
resolution  resolution x resolution is the number of points the localization should be estimated in the listening area. The points are evenly distributed along the axes. 
method 

localization_error  
deviation from the desired direction, defined as perceived_direction  desired_direction / rad  
perceived_direction  
the direction of arrival the binaural model has estimated for our given setup / rad  
desired_direction  
the desired direction of arrival indicated by the source position xs / rad  
x  corresponding xaxis 
y  corresponding yaxis 
x0  position and directions of the loudspeakers in the form n x 6, where n is the number of loudspeakers 
wierstorf2013(X,Y,phi,xs,src,'wfs',L,nls,array) calculates the localization error for the defined wave field synthesis or stereophony setup. The localization error is defined here as the difference between the perceived direction as predicted by the dietz2011 binaural model and the desired direction given by xs. The loudspeaker setup for the desired reproduction method is simulated via HRTFs which are than convolved with white noise which is fed into the binaural model.
The following parameters may be passed at the end of the line of input arguments:
'resolution',resolution  
Resolution of the points in the listening area. Number of points is resoluation x resolution. If only one point in the listening area is given via single values for X and Y, the resolution is automatically set to 1.  
'nls',nls  Number of loudspeaker of your WFS setup. Default value is 2. 
'array',array  Array type to use, could be 'linear' or 'circle'. Default value is 'linear'. 
'hrtf',hrtf  HRTF database. This have to be in the TUBerlin matformat, see: https://dev.qu.tuberlin.de/projects/measurements/wiki/IRs_file_format Default HRTF set is the 3m one from TUBerlin measured with the KEMAR. 
'lookup',lookup  Lookup table to map ITD values to angles. This can be created by the itd2anglelookuptable function. Default value is the lookup table wierstorf2013itd2anglelookup.mat that comes with AMT. 
For the simulation of the wave field synthesis or stereophony setup this functions depends on the SoundFieldSynthesis Toolbox, which is available here: http://github.com/sfstoolbox/sfs. It runs under Matlab and Octave. The revision used to generate the figures in the corresponding paper is a8914700a4.
M. Dietz, S. D. Ewert, and V. Hohmann. Auditory model based direction estimation of concurrent speakers from binaural signals. Speech Communication, 53(5):592605, 2011. [ DOI  http ]
H. Wierstorf, M. Geier, A. Raake, and S. Spors. A free database of headrelated impulse response measurements in the horizontal plane with multiple distances. In Proceedings of the 130th Convention of the Audio Engineering Society, 2011.
H. Wierstorf, A. Raake, and S. Spors. Binaural assessment of multichannel reproduction. In J. Blauert, editor, The technology of binaural listening, chapter 10. Springer, BerlinHeidelbergNew York NY, 2013.