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. 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. hearingresearch, 360:55--75, 2018.