MAP

%

% When the demonstration is satisfactory, freeze it by renaming it <demoxx>

restorePath=path;

addpath (['..' filesep 'MAP'],    ['..' filesep 'wavFileStore'], ...

    ['..' filesep 'utilities'])

% or

% NB probabilities are not corrected for refractory effects

%% #3 pure tone, harmonic sequence or speech file input

signalType= 'tones';

sampleRate= 100000;

% toneFrequency= 250:250:8000;    % harmonic sequence (Hz)

rampDuration=.005;              % seconds

% or

%% #4 rms level

% signal details

%% #5 number of channels in the model

BFlist=round(logspace(log10(lowestBF), log10(highestBF), numChannels));

%   or specify your own channel BFs

% BFlist=toneFrequency;

%  *after* the MAPparams file has been read

% This example declares only one fiber type with a calcium clearance time

% constant of 80e-6 s (HSR fiber) when the probability option is selected.

% paramChanges={'AN_IHCsynapseParams.ANspeedUpFactor=5;', ...

%     'IHCpreSynapseParams.tauCa=86e-6;'};

% paramChanges={'DRNLParams.rateToAttenuationFactorProb = 0;'};

%% delare 'showMap' options to control graphical output

showMapOptions.printModelParameters=1;   % prints all parameters

showMapOptions.showModelOutput=1;       % plot of all stages

showMapOptions.printFiringRates=1;      % prints stage activity levels

showMapOptions.showACF=0;               % shows SACF (probability only)

showMapOptions.showEfferent=1;          % tracks of AR and MOC

showMapOptions.surfProbability=1;       % 2D plot of HSR response 

if strcmp(AN_spikesOrProbability, 'spikes')

    % avoid nonsensical options

    showMapOptions.surfProbability=0;

    showMapOptions.showACF=0;

end

if strcmp(signalType, 'file')

    % needed for labeling plot

%% Generate stimuli

switch signalType

    case 'tones'

        inputSignal=createMultiTone(sampleRate, toneFrequency, ...

inputSignal=inputSignal.*ramp;

% add 10 ms silence

inputSignal= [silence inputSignal silence];