|
rmeddis@35
|
1 restorePath=path;
|
|
rmeddis@35
|
2 addpath (['..' filesep 'MAP'], ['..' filesep 'wavFileStore'], ...
|
|
rmeddis@35
|
3 ['..' filesep 'utilities'])
|
|
rmeddis@35
|
4
|
|
rmeddis@35
|
5 global OMEParams DRNLParams IHC_cilia_RPParams IHCpreSynapseParams
|
|
rmeddis@35
|
6 global AN_IHCsynapseParams MacGregorParams MacGregorMultiParams
|
|
rmeddis@35
|
7 global dt ANdt savedBFlist saveAN_spikesOrProbability saveMAPparamsName...
|
|
rmeddis@35
|
8 savedInputSignal OMEextEarPressure TMoutput OMEoutput ARattenuation ...
|
|
rmeddis@35
|
9 DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV...
|
|
rmeddis@35
|
10 IHCoutput ANprobRateOutput ANoutput savePavailable ANtauCas ...
|
|
rmeddis@35
|
11 CNtauGk CNoutput ICoutput ICmembraneOutput ICfiberTypeRates ...
|
|
rmeddis@35
|
12 MOCattenuation
|
|
rmeddis@35
|
13
|
|
rmeddis@35
|
14 signalCharacteristics.type='tones';
|
|
rmeddis@35
|
15 signalCharacteristicssignalCharacteristics.sampleRate=50000;
|
|
rmeddis@35
|
16 signalCharacteristics.duration= 1;
|
|
rmeddis@35
|
17 signalCharacteristics.rampDuration=0.05;
|
|
rmeddis@35
|
18 signalCharacteristics.beginSilence=0.05;
|
|
rmeddis@35
|
19 signalCharacteristics.endSilence=0.05;
|
|
rmeddis@35
|
20 signalCharacteristics.toneFrequency=1000;
|
|
rmeddis@35
|
21 signalCharacteristics.leveldBSPL=50;
|
|
rmeddis@35
|
22
|
|
rmeddis@35
|
23 showMapOptions.printModelParameters=0; % prints all parameters
|
|
rmeddis@35
|
24 showMapOptions.showModelOutput=0; % plot of all stages
|
|
rmeddis@35
|
25 showMapOptions.printFiringRates=0; % prints stage activity levels
|
|
rmeddis@35
|
26 showMapOptions.showACF=0; % shows SACF (probability only)
|
|
rmeddis@35
|
27 showMapOptions.showEfferent=0; % tracks of AR and MOC
|
|
rmeddis@35
|
28 showMapOptions.surfProbability=0; % 2D plot of HSR response
|
|
rmeddis@35
|
29 showMapOptions.surfSpikes=0; % 2D plot of spikes histogram
|
|
rmeddis@35
|
30 showMapOptions.ICrates=0; % IC rates by CNtauGk
|
|
rmeddis@35
|
31
|
|
rmeddis@35
|
32 tic
|
|
rmeddis@35
|
33 fprintf('\n')
|
|
rmeddis@35
|
34 disp('Computing ...')
|
|
rmeddis@35
|
35
|
|
rmeddis@35
|
36 levels=80:10:100;
|
|
rmeddis@35
|
37 figure(10), clf
|
|
rmeddis@35
|
38
|
|
rmeddis@35
|
39 for level=levels
|
|
rmeddis@35
|
40 signalCharacteristics.leveldBSPL=level;
|
|
rmeddis@35
|
41 % MAPrunner(MAPparamsName, AN_spikesOrProbability, ...
|
|
rmeddis@35
|
42 % signalCharacteristics, paramChanges, showMapOptions)
|
|
rmeddis@35
|
43 MAPrunner('Normal', 'spikes', ...
|
|
rmeddis@35
|
44 signalCharacteristics, {}, showMapOptions)
|
|
rmeddis@35
|
45 ARmin=min(ARattenuation);
|
|
rmeddis@35
|
46 disp([num2str(level) ':' num2str(ARmin)])
|
|
rmeddis@35
|
47 time=dt:dt:dt*length(ARattenuation);
|
|
rmeddis@35
|
48 hold on, plot(time,(1-ARattenuation)/(1-min(ARattenuation)))
|
|
rmeddis@35
|
49 pause(0.1)
|
|
rmeddis@35
|
50 end
|
|
rmeddis@35
|
51
|
|
rmeddis@35
|
52 path(restorePath)
|
