diff multithreshold 1.46/testSynapse.m @ 0:f233164f4c86

first commit

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multithreshold 1.46/testSynapse.m	Fri May 27 13:19:21 2011 +0100
@@ -0,0 +1,88 @@
+function testSynapse
+% testSynapse tracks the quantity of available transmitter vesicles
+%  the computations are single channel using the first frequency
+%  in the targetFrequency box of the expGUI.
+% For, speed this function uses only probability and HSR fibers.
+% This cannot be changed because of the way AN_IHCsynapse is coded.This).
+
+global experiment  method IHCpreSynapseParams
+global  AN_IHCsynapseParams  stimulusParameters
+savePath=path;
+addpath (['..' filesep 'utilities'],['..' filesep 'MAP'])
+
+figure(6),clf
+plotColors='rbgkrbgkrbgkrbgkrbgkrbgk';
+set(gcf,'position',[5    32   264   243])
+
+sampleRate=5e4; dt=1/sampleRate;
+
+maskerLevels=-0:10:100;
+
+BFlist=stimulusParameters.targetFrequency(1);
+targetFrequency=stimulusParameters.targetFrequency;
+targetFrequency=targetFrequency(1); % only one frequency used
+
+showParams=0;
+useEfferent=0;  % no efferent
+
+silenceDuration=0.005;
+maskerDuration=0.1;
+recoveryDuration=0.15;
+rampDuration=0.004;
+
+maskerTime=dt:dt:maskerDuration;
+
+rampTime=dt:dt:rampDuration;
+ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ...
+    ones(1,length(maskerTime)-length(rampTime))];
+ramp=ramp.*fliplr(ramp);
+
+initialSilence=zeros(1,round(silenceDuration/dt));
+recoverySilence=zeros(1,round(recoveryDuration/dt));
+
+signal=sin(2*pi*targetFrequency'*maskerTime);
+signal= ramp.*signal;
+signal=[initialSilence signal  recoverySilence];
+
+levelCount=0;
+qtMatrix=[];
+for leveldB=maskerLevels
+    levelCount=levelCount+1;
+
+    amp=28e-6*10^(leveldB/20);
+    inputSignal=amp*signal;
+
+    AN_spikesOrProbability='probability';
+    MAPparamsName=experiment.name;
+    showPlotsAndDetails=0;
+
+    global savePavailable
+   
+        MAP1_14(inputSignal, 1/dt, targetFrequency, ...
+            MAPparamsName, AN_spikesOrProbability);
+
+    % ignore LSR channels (if any) at the top of the matrix
+    qt=savePavailable(end, :);
+
+    synapsedt=dt;
+    time=synapsedt:synapsedt:synapsedt*length(qt);
+
+    figure(6)
+    qtMatrix=[qtMatrix; qt];
+    plot(time,qt,  plotColors(levelCount))
+    hold on
+    xlim([0 max(time)])
+    ylim([0 AN_IHCsynapseParams.M])
+end
+
+set(gcf,'name','pre-synaptic available transmitter')
+title(['q - available vesicles:' num2str(BFlist) ' Hz'])
+legend(strvcat(num2str(maskerLevels')),'location','southeast')
+legend boxoff
+grid on
+
+figure(88), [c,H]=contour(time, maskerLevels,qtMatrix); clabel(c, H);
+set(gcf,'position',[ 276    31   264   243])
+xlabel('time'), ylabel('maskerLevels')
+
+path(savePath);