THE AUDITORY MODELING TOOLBOX

Applies to version: 1.6.0

View the code

Go to function

verhulst2018
Cochlear transmission-line model including a model of the brainstem

Usage:

output = verhulst2018(insig,fs,fc_flag)
[output, cf] = verhulst2018(insig,fs,fc_flag, ...)

Input parameters:

insig Acoustic signal (in Pa). Size: (time x ear).
fs Sampling rate (in Hz).
fc_flag Either a vector with frequencies specifying the probe positions along basilar-membrane sections, or 'all' to probe all 1000 cochlear sections, or 'abr' to probe 401 locations between 112 and 12000 Hz.

Output parameters:

cf Center frequencies (in Hz) of the probed basiliar membrane sections.
output

Structure with the following fields:

  • fs_an: Sample rate (in Hz) of the output.
  • fs_abr: Sample rate (in Hz) of the brainstem sections (IC, CN, W1, W3, and W5).
  • w1: Wave 1, output of the AN model.
  • w3: Wave 3, output of the CN model.
  • w5: Wave 5, output of the IC model.
  • an_summed: Sum of HSR, MSR and LSR responses (per ear) and the input to the CN (modelled by verhulst2015_cn). Provided by default. Can be disabled by the flag 'no_an'.
  • ihc: IHC receptor potential. Provided by default. Can be disabled by the flag 'no_ihc'.
  • cn: Detailed output of the CN. Provided by default. Can be disabled by the flag 'no_cn'.
  • ic: Detailed output of the IC. Provided by default. Can be disabled by the flag 'no_ic'.

Description:

verhulst2018(..) computes the basilar membrane displacement. The following additinal flags can be specified:

'an' Calculate the AN responses. Default.
'no_an' Disable the caclculation of the AN and following stages.
'cn' Calculate the CN responses. Default.
'no_cn' Disable the caclculation of the CN and following stages.
'ic' Calculate the IC responses. Default.
'no_ic' Disable the caclulation of the IC and following stages.
'no_mfb' Determines the output of W3 and W5: The output of only a single CN filter is returned, output.w3 and output.w5 are numeric. Default.
'mfb' Determines the output of W3 and W5: The output of the whole modulation filterbank is returned, output.w3 and output.w5 are cell arrays.
'progress' Display the calculation progress. Default.
'silent' No progress display. Opposite of 'progress'.
'no_debug' No display of debugging information. Default.
'debug' Display additional debugging information. Opposite of 'no_debug'.

verhulst2018(..) takes also many additional parameters, see arg_verhulst2018 for more details.

The output can optionally provide the following information:

anfH: Responses of the high-SR fibers. Optional, only when called with flag 'anfH'.

anfM: Responses of the medium-SR fibers. Optional, only when called with flag 'anfM'.

anfL: Responses of the low-SR fibers. Optional, only when called with flag 'anfL'.

v: Optional velocity of the basilar membrane sections. Size: (time x section x ear). This is also the input to the AN modelled by verhulst2018_ihctransduction. Provided only when called with flag 'v'.

y: Optional displacement of the basilar membrane sections. Size: (time x section x ear). Provided only when called with flag 'y'.

oae: Optional ottoacoustic emission as sound pressure at the middle ear. Provided only when called with flag 'oae'.

License

This model is licensed under the UGent Academic License. For non-commercial academic research, you can use this file and/or modify it under the terms of that license. Further usage details are provided in the in the AMT directory "licences".

References:

S. Verhulst, A. Altoè, and V. Vasilkov. Functional modeling of the human auditory brainstem response to broadband stimulation. Hearing Research, 360:55--75, 2018.