THE AUDITORY MODELING TOOLBOX

Applies to version: 0.9.7

View the code

Go to function

ZIEGELWANGER2014 - Time of arrival estimates

Usage

[Obj,results]=ziegelwanger2014(data,estimation,outlierDetection,model,p0_onaxis)

Input parameters

Obj SOFA object
estimation (optional) select one of the estimation methods (1: Threshold-Detection, 2: Centroid of squared IR, 3: Mean Groupdelay, 4: Minimal-Phase Cross-Correlation (Max) (default), [TOAest]: pre-estimated TOAs)
outlierDetection (optional) detect outliers in estimated TOAs (0: off, 1: on (default values: [0.05;0.01]), [alpha r]: rejects outliers using the extreme Studentized deviance test with the significance level of ALPHA and upper bound of outlier rate R. )
model (optional) correct estimated toa, using geometrical TOA-Model (0: TOA estimated, 1: off-axis TOA modeled (default), 2: on-axis TOA modeled)
p0_onaxis (optional) startvalues for lsqcurvefit

Output parameters

Obj SOFA Object
results.toa
data matrix with time of arrival (TOA) for each impulse response (IR):
results.p_onaxis
estimated on-axis model-parameters
results.p_offaxis
estimated off-axis model-parameters

Description

Estimates the Time-of-Arrival for each measurement in Obj (SOFA) and corrects the results with a geometrical model of the head.

Requirements:

  1. SOFA API from http://sourceforge.net/projects/sofacoustics for Matlab (in e.g. thirdparty/SOFA)
  2. Optimization Toolbox for Matlab
  3. Data in hrtf/ziegelwanger2014

Examples:

To calculate the model parameters for the on-axis time-of-arrival model (p_onaxis) and for the off-axis time-of-arrival model (p_offaxis) for a given HRTF set (SOFA object, 'Obj') with the minimum-phase cross-correlation estimation, use:

[Obj,results]=ziegelwanger2014(Obj,4,1);

References:

H. Ziegelwanger and P. Majdak. Modeling the direction-continuous time-of-arrival in head-related transfer functions. J. Acoust. Soc. Am., 135:1278-1293, 2014.