Usage

[Obj,results]=ziegelwanger2013(Obj,method,model,p0_onaxis)

Input parameters

Output parameters

results.toa

data matrix with time of arrival (TOA) for each impulse response (IR):

System Message: WARNING/2 (<string>, line 42)

Block quote ends without a blank line; unexpected unindent.

dim 1: each toa in samples dim 2: each record channel

results.p_onaxis

estimated on-axis model-parameters

System Message: WARNING/2 (<string>, line 46)

Block quote ends without a blank line; unexpected unindent.

dim 1: [sphere-radius in m,

azimut of ear in radiants, elevation of ear in radiants, direction-independent delay in seconds]

System Message: WARNING/2 (<string>, line 50)

Definition list ends without a blank line; unexpected unindent.

dim 2: each record channel

results.p_offaxis

estimated off-axis model-parameters

System Message: WARNING/2 (<string>, line 53)

Block quote ends without a blank line; unexpected unindent.

dim 1: [sphere-radius in m,

xM in m, yM in m, zM in m, direction-independent delay in seconds, channel (starting at 1), azimut of ear in radiants, elevation of ear in radiants]

System Message: WARNING/2 (<string>, line 61)

Definition list ends without a blank line; unexpected unindent.

dim 2: each record channel

Description

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

The value of method is an integer choosing one of the following methods. XXX Explain for each method a little about how they work:

1. Threshold-Detection
2. Centroid of squared IR
3. Mean Groupdelay
4. Minimal-Phase Cross-Correlation (Max) (default)

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/ziegelwanger2013

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 method, use:

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