VERHULST2018 - Cochlear transmission-line model (improved, incl. brainstem)

Usage

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

Input parameters

`insig`	the input signal to be processed. Each column is processed in parallel, so it is possible to run several simulations in parallel
`fs`	sampling rate (Hz)
`fc_flag`	list of frequencies specifying the probe positions along the basilar membrane, or `all` to probe all 1000 cochlear sections, or `abr` to probe 401 locations between 112 and 12000 Hz.
`numH`	number of high-SR fibers. Must be larger than zero.
`numM`	number of medium-SR fibers. Can be zero.
`numL`	number of low-SR fibers. Can be zero.
`ic`	calculate the IC responses

Output parameters

`cf`	Center frequencies (Hz) of the probed basiliar membrane sections.
`output`	Structure with the following fields:
`fs_an`	Sample rate (Hz) of the output.
`fs_abr`	Sample rate (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 channel) 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

This function computes the basilar membrane displacement.

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'`	velocity of the basilar membrane sections [time section channel] and input to the AN (modelled by `verhulst2018_ihctransduction`. Optional, provided only when called with flag `v`.
`'y'`	displacement of the basilar membrane sections [time section channel]. Can be disabled by the flag `no_y`.
`'oae'`	ottoacoustic emission as sound pressure at the middle ear. Can be disabled by the flag `oae`.

License

This model is licensed under the UGent Academic License. Further usage details are provided in the UGent Academic License which can be found in the AMT directory "licences" and at <https://raw.githubusercontent.com/HearingTechnology/Verhulstetal2018Model/master/license.txt>.

References:

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