Applies to version: 0.9.2

View the code

Go to function

WIERSTORF2013 - estimate the localization within a WFS or stereo setup


[...] = wierstorf2013(X,Y,phi,xs,src,L,resolution,method,...);

Input parameters

X range of the x-axis [xmin,xmax] in m (for a single x point you have to specify [x x]
Y range of the y-axis [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 x-axis
xs position of the point source in m / direction of the plane wave
source type
'ps' for a point source 'pw' for a plane wave
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.
reproduction setup
'wfs' for wave field synthesis 'setreo' for stereophony

Output parameters

 deviation from the desired direction, defined as perceived_direction - desired_direction / rad
 the direction of arrival the binaural model has estimated for our given setup / rad
 the desired direction of arrival indicated by the source position xs / rad
x corresponding x-axis
y corresponding y-axis
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 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 TU-Berlin mat-format, see: Default HRTF set is the 3m one from TU-Berlin 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 Sound-Field-Synthesis Toolbox, which is available here: 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):592-605, 2011. [ DOI | http ]

H. Wierstorf, M. Geier, A. Raake, and S. Spors. A free database of head-related 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 multi-channel reproduction. In J. Blauert, editor, The technology of binaural listening, chapter 10. Springer, Berlin-Heidelberg-New York NY, 2013.