Mercurial > hg > map
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)