view testPrograms/repeatTester.m @ 36:3ea506487b3b

changes to multiThreshold
author Ray Meddis <rmeddis@essex.ac.uk>
date Tue, 20 Sep 2011 14:04:47 +0100
parents 25d53244d5c8
children
line wrap: on
line source
restorePath=path;
addpath (['..' filesep 'MAP'],    ['..' filesep 'wavFileStore'], ...
    ['..' filesep 'utilities'])

global OMEParams DRNLParams IHC_cilia_RPParams IHCpreSynapseParams
global AN_IHCsynapseParams MacGregorParams MacGregorMultiParams
global dt ANdt  savedBFlist saveAN_spikesOrProbability saveMAPparamsName...
    savedInputSignal OMEextEarPressure TMoutput OMEoutput ARattenuation ...
    DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV...
    IHCoutput ANprobRateOutput ANoutput savePavailable ANtauCas  ...
    CNtauGk CNoutput  ICoutput ICmembraneOutput ICfiberTypeRates ...
    MOCattenuation 

signalCharacteristics.type='tones';
signalCharacteristicssignalCharacteristics.sampleRate=50000;
signalCharacteristics.duration= 1;
signalCharacteristics.rampDuration=0.05;     
signalCharacteristics.beginSilence=0.05;
signalCharacteristics.endSilence=0.05;                  
signalCharacteristics.toneFrequency=1000; 
signalCharacteristics.leveldBSPL=50;

showMapOptions.printModelParameters=0;   % prints all parameters
showMapOptions.showModelOutput=0;       % plot of all stages
showMapOptions.printFiringRates=0;      % prints stage activity levels
showMapOptions.showACF=0;               % shows SACF (probability only)
showMapOptions.showEfferent=0;          % tracks of AR and MOC
showMapOptions.surfProbability=0;       % 2D plot of HSR response 
showMapOptions.surfSpikes=0;            % 2D plot of spikes histogram
showMapOptions.ICrates=0;               % IC rates by CNtauGk

tic
fprintf('\n')
disp('Computing ...')

levels=80:10:100;
    figure(10), clf
    
for level=levels
signalCharacteristics.leveldBSPL=level;
% MAPrunner(MAPparamsName, AN_spikesOrProbability, ...
%     signalCharacteristics, paramChanges, showMapOptions)
MAPrunner('Normal', 'spikes', ...
    signalCharacteristics, {}, showMapOptions)
ARmin=min(ARattenuation);
disp([num2str(level) ':' num2str(ARmin)])
time=dt:dt:dt*length(ARattenuation);
hold on, plot(time,(1-ARattenuation)/(1-min(ARattenuation)))
pause(0.1)
end

path(restorePath)